to test GSW

nicojourdain · Mar 6, 2017 · 6262a8f · 6262a8f
1 parent fdfa2c4
commit 6262a8f
Showing 1 changed file with 26 additions and 0 deletions.
diff --git a/test_GSW.f90 b/test_GSW.f90
@@ -0,0 +1,26 @@
+program test
+
+use netcdf
+
+use gsw_mod_kinds
+use gsw_mod_netcdf
+use gsw_mod_toolbox
+
+use gsw_mod_error_functions, only : gsw_error_code, gsw_error_limit
+
+IMPLICIT NONE
+
+REAL*8,DIMENSION(2,2) :: salout, sal, dep, lat, lon
+
+call gsw_saar_init (.true.)
+
+sal(:,:)=35.0
+dep(:,:)=10.0
+lat(:,:)=45.0
+lon(:,:)=175.0
+
+salout = gsw_sa_from_sp(sal,dep,lon,lat)
+
+write(*,*) salout
+
+end program test