test_GSW.f90

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