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
Our Free Plans just got better! | Auth0
You asked, we delivered! Auth0 is excited to expand our Free and Paid plans to include more options so you can focus on building, deploying, and scaling applications without having to worry about your security. Auth0 now, thank yourself later.
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 .................................................