[MITgcm-support] Initial hydrog/bathy-Files and CS-grid

[Martin Losch]

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.

Martin

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