[MITgcm-devel] notes from today: 2003-10-22

[Dimitris Menemenlis]

Hi Ed, going over this morning's todo list:

>> - read in data sets and interpolate (EXF has the ability to
>>   interpolate on the fly with bi-cubics???)

Yes, exf_interp.F will interpolate using either bilinear or bicubic.
input grid must be spherical-polar and global since it assumes
wrap-around (but it will also work for limited domains, if you make sure
there is sufficient padding (1 for bilinear, 2 for bicubic) on all
sides).  Zonal spacing must be even but meridional spacing can be
arbitrary.  There is a verification experiment (global_with_exf, see
data.exf and data.exf_clim for examples of how to specify input grids). 
It only works for scalar fields at this time though, i.e., no vector
rotation.  But hope to implement that while at MIT next week.

>> add a netcdf package

For MITgcm_contrib/ocean_inversion_project I am using netcdf output
files, see for example, code/ptracers_init.F,
write_netCDF/write_nc_MaskAreaBathy.F, write_netCDF/netcdf.inc, and
write_netCDF/linux_ia32_g77_netcdf

D.

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Insitute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720