[MITgcm-support] Velocity cubed sphere + pickup seaice

[Elja Huibregtse]

Thanks, Martin and Dimitris, for your useful replies.

My study object is a moon of Jupiter, thus the radius and other parameters differ from those used for the Earth. I changed all these parameters in set_defaults.f. Trying to calculate the values of angleCosC and angleSinC, I found out that the radius stored in RAC (calculated from: r_earth=sqrt(sum(sum(RAC))./(4.*pi)) in cubeCalcAngle.m, line 23), is still the radius of the Earth. 

Now, I think that all the 'coordinate' files, like RAC, RAS, RAZ etc. still contain the configuration of the Earth, which will, of course, influence my results. Is it possible to adapt those files to the correct configuration?

Thanks!
Elja  


-----Oorspronkelijk bericht-----
Van: mitgcm-support-bounces at mitgcm.org namens Martin Losch
Verzonden: wo 6/4/2008 8:54
Aan: mitgcm-support at mitgcm.org
Onderwerp: Re: [MITgcm-support] Velocity cubed sphere + pickup seaice
 
Hi Elja,

in addition, the orientation of u and v-velocity is along grid  
coordinates, which will give you funny results, if you try to plot  
them on some sort of map-projection. After interpolating, you need to  
rotate the velocities according to the values of angleCosC and  
angleSinC. Unfortunately, these fields are often not set, because the  
model does not need them for hydrostatic integrations (so that they  
are 1 and 0). There is a matlab script to compute the angles in  
<utils/matlab/cs_grid/cubeCalcAngle.m>

I am also attaching a matlab function that does the interpolation and  
takes care of the edges (for me, so I am not sure that this will work  
for you or anyone else, too).

Martin

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20080604/6c288f31/attachment.htm>