[MITgcm-support] Initial hydrog/bathy-Files and CS-grid
Martin.Losch at awi.de
Mon Mar 5 03:13:55 EST 2012
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:
> 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 Brunetti
> Institute for Environmental Sciences (ISE)
> University of Geneva
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
More information about the MITgcm-support