[MITgcm-support] initializing a grid

Jean-Michel Campin jmc at mit.edu
Fri Feb 24 14:11:34 EST 2017


Hi Matt,

I would suggest to set the longitude origin (xgOrigin = ... ) to -20. something 
(i.e., 20.W); Normally, it should just work fine. 
For instance, exf_interp.F does care of a 360 shift (line 288 - 297:
> C--   put xG in interval [ lon_0 , lon_0+360 [
and just to be on the safe side, you could try a short run with exf_debugLev=3
(and pkg/debug compiled) to check that everything is right.

Cheers,
Jean-Michel

On Fri, Feb 24, 2017 at 10:03:23AM -0800, Matthew Mazloff wrote:
> Hi Jean-Michel et al
> 
> I am setting up a grid that is from 20W to 20E. Since MITgcm deals in degrees E, what is recommended?
> 
> From looking at the code it isn???t clear to me what would happen if I told the MITgcm to have a grid 340 to 380???..
> 
> One option is to just shift this by 20 degrees east, and pretend the grid is from 0 to 40 instead of 340 to 360 to 20. This complicates any on the fly interpolation but I can deal with that. Is this the easiest?
> 
> Thanks!
> Matt
> 
> 
> 
> 
> 
> 
> 
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support



More information about the MITgcm-support mailing list