[MITgcm-support] Re : Initial hydrog/bathy-Files and CS-grid
Maura.Brunetti at unige.ch
Tue Mar 6 05:05:35 EST 2012
thank you, I agree that initial conditions are our 'responsibility' and this is why I
am trying to reconstruct and reproduce those in the tutorials (in particular, the
The problem was on the cubed sphere configuration and the fact that I am not familiar with it!
I think that now I have solved most of the problems, starting from the matlab scripts in
MITgcm_contrib (for example mk_input.m, where the matlab function 'interp2' is used also
for cs configuration).
On Mar 5, 2012, Martin Losch wrote:
> Hi Maura,
> did you get an answer to your question?
> The input fields for bathymetry and initial conditions (and for surface forcing) are your "responsibility". The interpolation by the exf-package could be used, but it's really bi-linear or bi-cubic interpolation, which you can more easily do offline with e.g. matlab (or octave). Use the grid information for a script like this (I am using matlab syntax):
> xc = rdmds('XC');
> yc = rdmds('YC');
> % 2d field fld2d on lon-lat grid:
> fld = interp2(lon,lat,fld2d,xc(:),yc(:),method); % method = 'linear','cubic', etc see help interp2
> fld = reshape(fld,size(xc));
> there are more complicated and robust methods, especially bathymetry probably needs to be tweaked a little.
> On Feb 23, 2012, at 4:29 PM, Maura BRUNETTI wrote:
> > Hi,
> > I would like to use the cubed sphere configuration for paleoclimate simulations. For this I need to create bathyFile, hydrogThetaFile and hydrogSaltFile on the cubed sphere for different past bathymetries.
> > In other words, I would like to know how to create files analogous to bathy_cs32.bin, lev_S_cs_15k.bin, lev_T_cs_15k.bin in MITgcm/verification/global_ocean.cs32x15/input for different bathymetries.
> > I know that the exf package allows us to interpolate on the curvilinear grid, and indeed I use this package for interpolating all the surface boundary conditions (uvwind, swdown, precip, etc) on the cubed-sphere grid. Is it possible to use this package for interpolating bathyFile and hydrogTheta/SaltFile on the cubed sphere? How do I need to set parameters in data and data.exf? Or what do you suggest to use?
> > Many thanks,
> > Maura
> > --
> > Maura Brunetti
> > Institute for Environmental Sciences (ISE)
> > University of Geneva
> > Switzerland
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-support
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> End of MITgcm-support Digest, Vol 105, Issue 2
Institute for Environmental Sciences (ISE)
University of Geneva
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the MITgcm-support