Browse Code
SOSIE is Only a Surface Interpolation Environment, a versatile tool that allows fast and high quality interpolation of geophysical 2D fields from a gridded domain to another.
Project Samples
Categories
Earth SciencesLicense
GNU General Public License version 2.0 (GPLv2)Follow SOSIE
Other Useful Business Software
MongoDB Atlas runs apps anywhere
MongoDB Atlas gives you the freedom to build and run modern applications anywhere—across AWS, Azure, and Google Cloud. With global availability in over 115 regions, Atlas lets you deploy close to your users, meet compliance needs, and scale with confidence across any geography.
Rate This Project
Login To Rate This Project
User Reviews
-
As it seems there is no bug tracker I write here. Does not compile with gfortran 5.4.0 ubuntu 16.04 (fortran format problem). The patch to correct it is: diff -Naur sosie-2.6.4.orig/src/mask_drown_field.f90 sosie-2.6.4/src/mask_drown_field.f90 --- sosie-2.6.4.orig/src/mask_drown_field.f90 2016-04-26 10:39:11.000000000 +0200 +++ sosie-2.6.4/src/mask_drown_field.f90 2017-03-10 16:59:58.926326000 +0100 @@ -170,8 +170,8 @@ - READ(cv_how_far, '(i)') i_how_far - IF ( l_drwn_f ) WRITE (*, '("How far (points) from sea to extrapolate into land ",i)') i_how_far + READ(cv_how_far, '(i8)') i_how_far + IF ( l_drwn_f ) WRITE (*, '("How far (points) from sea to extrapolate into land ",i8)') i_how_far READ(cv_ewp, '(i2)') iewper
-
A great interpolation library. May I ask a question. Does the drown method only works for global with with east-west periodicity? Would you please have a look? In the MODULE MOD_DROWN, if k_ew=-1, there would be an error. As the following code will not run, so the vim_per and vip_per will only have 0 values. IF (k_ew >= 0) THEN vim_per = (/ ni-k_ew , ni-1 /) vip_per = (/ 2 , 1+k_ew /) END IF Then, jip =0, array of maskv will overflow DO jci = 1, 2 jim = vim_per(jci) ! ji-1 ji = ivi(jci) ! first ji = 1, then ji = ni jip = vip_per(jci) ! ji+1 .................................................
