[Thread Prev][Thread Next][Index]

Re: vector and shade in curvilinear coordinates



Dear Eric,

thanks for the script. I already thought that there would be such a solution, but did not know how to find out the details.

Funny that the vector command worked anyway.

Ciao,
Arne


Eric Maisonnave wrote:
Dear Arne,

please find in attachment the ferret script we use in Cerfacs to visualize
ORCA 2D fields with ferret. nav_lon,nav_lat variables must contain lon,lat
fields.

Thanks to Olivier Marti (LSCE).


Eric





On Tue, 15 Jul 2003, Arne Biastoch wrote:


Date: Tue, 15 Jul 2003 09:49:05 +0200
From: Arne Biastoch <abiastoch@ifm.uni-kiel.de>
To: Ferret mailing list <ferret_users@ferret.wrc.noaa.gov>
Subject: vector and shade in curvilinear coordinates


Dear Ferreters,

plotting output of the ORCA2 model (global standard configuration of the
OPA model) I tried to make use of the 3-argument versions of the shade
and vector commands. The data are on a tripolar grid, latitude and
longitude are therefore 2D fields (see attached grid.gif).

Plotting vectors using
	yes? vector u,v,lon,lat
works fine (land masks can easily and correctly be identified, see
vector.gif).

But doing a shade (or fill) using
	yes? shade temp,lon,lat
gives distorted fields, only a small central regions seems to be ok (see
shade.gif), outer regions are "fringed". I played around with the axes
(modulo etc.), without any enhancement.

The actual construction of the figures is different for both commands:
vectors are correctly drawn point by point, starting with indices 1,1 of
lat,lon. In contrast the shade command is build up from south to north.

Any ideas from the experts and/or other ORCA/ferret users?

Cheers,
	Arne


--

Dr. Arne Biastoch

Institute for Marine Research         phone: ++49 431 600-4013
FB1 Ocean Circulation and Climate     fax  : ++49 431 600-1515
Duesternbrooker Weg 20                email: abiastoch@ifm.uni-kiel.de
24105 Kiel, Germany

http://www.ifm.uni-kiel.de/fb/fb1/tm/data/pers/abiastoch.html


Eric MAISONNAVE

C.E.R.F.A.C.S.
42, Av. G. Coriolis
F-31057 Toulouse Cedex

tel: (+33) 5 61 19 30 12
fax: (+33) 5 61 19 30 00



------------------------------------------------------------------------

! NOAA/PMEL TMAP
! Program FERRET ! Version 5.22 - 08/23/00
! 23-Jul-01 13:56 LET ii = I[I=1:182]+J[J=1:149]*0 ; LET jj = I[I=1:182]*0+J[J=1:149]
LET/QUIET/TITLE="Longitude"/UNIT"=Degrees_east" lon1 = IF ii LE 52 AND nav_lon GE 0 THEN nav_lon-360.0 ELSE nav_lon
LET/QUIET/TITLE="Longitude"/UNIT"=Degrees_east" lon2 = IF ii LE 19 AND jj GE 98 AND jj LE 120 THEN lon1+360 ELSE lon1
SHADE/HLIM=-285:85 $1, lon2, nav_lat
--

Dr. Arne Biastoch

Institute for Marine Research         phone: ++49 431 600-4013
FB1 Ocean Circulation and Climate     fax  : ++49 431 600-1515
Duesternbrooker Weg 20                email: abiastoch@ifm.uni-kiel.de
24105 Kiel, Germany

http://www.ifm.uni-kiel.de/fb/fb1/tm/data/pers/abiastoch.html



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement