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

[Maura BRUNETTI]

Dear Martin, 

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 
global_ocean.cs32x15 experiment).
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). 

Cheers, 

Maura


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.
> 
> 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
> 
> 
> 
> 
> ------------------------------
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 
> 
> End of MITgcm-support Digest, Vol 105, Issue 2
> **********************************************
> 
> 
--
Maura Brunetti
Institute for Environmental Sciences (ISE)
University of Geneva
Switzerland
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120306/49629e1c/attachment.htm>