! -------------------------------------------
! Purpose: read in L_r and C from Chelton et al. 1998
!          from ascii data file produced with F90 program "raw2grid.f90".
!          Output to netCDF file
! Usage:   start ferret (see ferret.wrc.noaa.gov/Ferret)
!	   type: go map2netCDF.jnl
!
! Coded by Rob Scott
! do with it whatever you want.
! -------------------------------------------
DEFINE SYMBOL outfile = rossby_radius.cdf

DEFINE AXIS/x=0.5E:359.5E:1/UNIT=longitude/MODULO x_axis
DEFINE AXIS/y=89.5S:89.5N:1/UNIT=latitude y_axis

DEFINE GRID/x=x_axis/y=y_axis grid_2d

FILE/GRID=grid_2d/ORDER=xy/COLUMN=6/VAR=jj,ii,lat,lon,c,lr rossby_radius_map.dat

SET VARIABLE/BAD=999/UNITS=degrees/TITLE="latitude" lat
SET VARIABLE/BAD=999/UNITS=degrees/TITLE="longitude" lon
SET VARIABLE/BAD=-1/UNITS="m/s"/TITLE="Rossby phase speed" c
SET VARIABLE/BAD=-1/UNITS=km/TITLE="Rossby radius" lr

SAVE/FILE=($outfile) lr, c