[MITgcm-support] exf_interp

[Dimitris Menemenlis]

Georges, as currently written, exf_interp works with a latitude/longitude coordinate systems only.

What you need to do is to interoplate/extrapolate your forcing field to a regular
latitude/longitude grid and then describe the coordinates of this grid in the
data.exf file.  An example data.exf file for exf_interp is here:
MITgcm/verification/global_with_exf/input/data.exf

Specifically you need to provide:
C     lon_0, lat_0       :: lon and lat of sw corner of global input grid
C     lon_inc            :: scalar x-grid increment
C     lat_inc            :: vector y-grid increments
C  nlon, nlat       :: size in x & y direction of input file to read

Note that lon_inc is a scalar, i.e., zonal increments are constant,
while lat_inc is a vector, i.e., meridional incrments can vary with latitude.

Hope this helps, Dimitris

On Jan 4, 2010, at 1:39 PM, Georges Djoumna wrote:

> 
> Happy new year to all!
> 
> Hi Dimitris,
> 
> I don't understand quite well your sentence " for ex_interp, you need to provide input forcing fields on a regular latitude/longitude". Does exf_interp work only in the latitude- longitude framework? I have observation at three points in a Cartesian grid (I also have their latitude/longitude coordinate) I'm using a regular Cartesian grid in the data file. My understanding is that I have to convert the Cartesian grid into the latitude/longitude, am I right?
> 
> How to set a value for  forcing_lon0 (vstress_lon0)? are there some rules to follow?
>  
> Thanks.
> 
> Georges
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100107/68b90030/attachment.htm>