General

• doall.m is a master script that does call all interpolation scripts shown below
• archivedirectory.m contains a readable path to the data tree. All raw data are hold in seperate subdirectories (e.g. etopo5, hellerman/raw etc.)

--> please adopt variable ARCHIVEDIR
• input data are taken from the IRI/LDEO Climate Data Library. They are also help locally on ice.ucsd.edu in /export/ice3/biastoch/DATA
• setfmt.m specifies format of the output files (e.g. real*4, real*8, IEEE). These parameter have to match the specifications of the machine that runs the model. For options refer to routines rdda.m or wrda.m.
• setgrid.m contains grid parameter that have to be set:

lon_lo (westernmost longitude), lon_L (length of x-axis), dlon (increment)
lat_lo, lat_L, dlat similar
dz  (vertical grid spacing)

Topography

• topo.m interpolates ETOPO5 onto the model grid and allows interactive manipulations (uses extract_etopo5.m, grid_sphere.m, finegrid.m, gen_hxhy.m, xyrecur_fv.m, edtopo.m)
• To edit an existing topography:

topo=rdda('bathymetry.bin',[nxc nyc],1);  % nxc and nyc are x/y dimensions topo=edtopo(topo);

• gen_pmask.m generates mask for tracer points (uses hfac.m)
• Output files:

bathymetry.bin	bathymetry [x,y]
pmask.bin	mask at pressure points [x,y,z]

• tiling.m shows a splitting of the domain on multiprocessors. tnx and tny specify thereby the dimensions of the subdomains.

Initial Data (Levitus et al. 1994)

• levit.m calls exintlevitclim.m (interpolation of climatology, uses extract_levitus_clim.m and xyexpand.m) and exintlevitmonthly.m (iinterpolation of monhtly surface data, uses  extract_levitus_monthly.m and xyexpand.m) for interpolation onto the model grid. Arguments "temp" or "salt" are given to the routines.
• please note that input data for temperatures are potential temperatures. There is no conversion made and the model requires potential temperatures!
• both scripts interpolate gaps and make tests that all ocean points are filled with data. If check2dmask.m and check3dmask.m run through the files are ready and safe to put into the model.
• Output files:

lev_clim_temp.bin	annual mean temperature [x,y,z]	deg C
lev_clim_salt.bin	annual mean salinity [x,y,z]	ppm
lev_monthly_sstemp.bin	monthly sea surface temperature at Levitus z=0 [x,y,t]	deg C
lev_monthly_sssalt.bin	monthly sea surface salinity at Levitus z=0 [x,y,t]	ppm
lev_monthly_temp.bin	monthly sst, volume interpolated to first model level [x,y,t]	deg C
lev_monthly_salt.bin	monthly sss, volume interpolated to first model level [x,y,t]	ppm

• Data can be read by FERRET using the following script:

define axis/x=-98:45:1/units=deg xax               ! create x axis define axis/y=-18:81:1/units=deg yax               ! create y axis define axis/z=1:31:1 zaxtmp                        ! number of vertical dimensions define grid/z=zaxtmp zgridtmp                      ! create temporary z axis file/var=dz/grid=zgridtmp/columns=31  deltaz.dat   ! dz is read from a file define axis/z/from_data/name=zax/units=m/edges/depth dz[k=@iin] ! create z axis define grid/x=xax/y=yax/z=zax xyzgrid              ! create grid file/var=tmp/grid=xyzgrid/ez/format=stream/type=r8 lev_clim_temp.bin ! read data let/TITLE="TEMPERATURE"/UNIT="deg C" temp=ignore0(tmp) ! ignore zeros

Otherwise the conversion can be made together with all other input data by the script convert.jnl(see below).

Surface Forcing (COADS, DaSilva et al. 1994)

dasilva.m interpolates COADS data onto the model grid.

Output files:

dasilva_taux.bin	monthly zonal wind stress [x,y,t]	N/m2
dasilva_tauy.bin	monthly meridional wind stress [x,y,t]	N/m2
dasilva_qnet.bin	monthly net surface heat flux [x,y,t]	W/m2
dasilva_emp.bin	monthly evaporation - precipitation [x,y,t]	m/s

The following additional data sets can be provided: Hellerman and Rosenstein wind data, data from the global state estimation (Iteration 60).

Conversion to NetCDF

• conversion of binary data into NetCDF format can simply be done by using the FERRET script convert.jnl:

> source /net/ice/ecco/ferret/ferret_paths                    <-- if not already done at startup > ferret         NOAA/PMEL TMAP         Program FERRET (V500beta1.1)         Version 5.00 - 05/26/99         11-Nov-99 16:03  yes? go convert      .... yes? exit >

• after successfull conversion the following NetCDF files are created:

bathymetry.cdf	bathymetry and pmask
lev_clim.cdf	annual data from Levitus
lev_surf.cdf	monthly data from Levitus
dasilva.cdf	monthly data from COADS

• please note that the NetCDF contain full grid information. So FERRET takes care of nonuniform grid spacings when plotting or calculating.