[MITgcm-devel] suggestion llc/cubed sphere
Patrick Heimbach
heimbach at MIT.EDU
Thu Nov 23 07:29:24 EST 2006
Hi Martin,
these angles are actually in the grid file:
> ncdump grid.t001.nc
... (snip) ...
double AngleCS(Y, X) ;
AngleCS:description = "Cos of grid orientation angle
at cell cen
ter" ;
double AngleSN(Y, X) ;
AngleSN:description = "Sin of grid orientation angle
at cell cen
ter" ;
I think they're also in the .mitgrid file,
at least from looking at code ini_curvilinear_grid.F
Finally, Chris has started the heroic effort to write matlab
script for plotting vectors. Maybe he can send u latest version.
Cheers
-Patrick
PS:
We're making progress on LLC grid runs, see
http://mitgcm.org/~heimbach/llc_setup/
for 15-level runs (still problems with 50-level runs)
On Nov 23, 2006, at 7:09 AM, <Martin.Losch at awi.de> wrote:
> Hi there,
>
> I have been trying desperately to write a matlab scripts that plots
> velocities in the same way that llc_pcol for cs_pcol plot
> skalars for the lat-lon-cap grid and the cubed sphere grid
> respectively. Now I realize that this is only possible, if I have the
> proper grid information: the angles between the grid and the lat-
> lon directions. I believe they are supposed to be in angleCS
> and angleSN (angleCosC and angleSinC in the code). Unfortunately
> they are not in the tile???.mitgrid by default. Maybe it
> is a good idea to alway have these angle in the code and the
> grid.t???.nc files. Or is there a simple way to compute them?
> Or is my entire effort to plot velocities on the cubed-sphere
> obsolete because there is matlab function that does this
> already?
>
> Martin
>
>
> Martin Losch
> Alfred Wegener Institute
> Postfach 120161, 27515 Bremerhaven, Germany;
> Tel./Fax: ++49(0471)4831-1872/1797
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20061123/4569ff1d/attachment.htm>
More information about the MITgcm-devel
mailing list