[Thread Prev][Thread Next][Index]

Re: How to plot the topography in Y-Z plate?



On Thu, 2 Mar 2000, Yang Kun wrote:

> There is one more question. 
> Now, for example,I use 
>    go fland_levitus z=1:5000 y=90s:90n
> It gives me a topography. But, what's the longitude it uses as
> the value of x? 
> I read the script but can not understand it well, but I guess
> it used some default lon. Since I want the lon. average 
> topography, how to modify your script to do this?

The script works if the longitude (or a range of longitudes) is defined.
It tried the following and it worked.

yes? go fland_levitus z=0:1000 gray overlay x=30:30 y=-90:90

The script will not work for

yes? go fland_levitus z=1:5000 y=90s:90n

FERRET reports an error "*** NOTE: 3rd argument of GO fland_levitus must
be "overlay", "basemap", or "transparent" ."

I am not sure how to do this for an average over a range of longitudes.
Since the file levitus_mask.nc has a MASK defined, I guess it is necessary
to define a new variable, say, MASK_AVE, to get the new mask over the
required region.  Then, you could use MASK_AVE instead.  Only, an average
won't work here because you need to define land at a given y-z point
even if there is land at just one longitude.

I have never done this, but think it should work. The FERRET gurus will
correct me if I am wrong :).

Regards, Shankar.

> > ------------------------------------------------------------------------
> > ! NOTE: for a contour plot in the XY plane BELOW the surface, you must supply
> > ! the first argument to this GO script.
> > ! usage: GO fland_levitus k=N  (where N is the required level)
> > ! or     GO fland_levitus z=Z  (where Z is the required depth)
> > !        GO fland_levitus      (in other cases)
> > \cancel mode verify
> > ! fland_levitus.jnl - a land mask for the LEVITUS dataset (overlay or "basemap")
> > ! 1/97
> > 
> > ! Description: Plot/overlay an an area-filled continental map using LEVITUS mask
> > 
> > ! usage:               arg1      arg2        arg3         arg4      arg5
> > !  GO fland_levitus [z=lo:hi] [palette] [overlay_style] [x=lo:hi] [y=lo:hi]
> > 
> > ! arg1  - depth limits - specify as for any FERRET depth range
> > !                        use either k=*:* or z=*:*
> > ! arg2  - palette - may be any palette on your system.  Normally solid colors
> > !	           like black, gray,red, etc. (see Fpalette '*') default=gray
> > ! arg3  - overlay_style - "basemap", "transparent", or "overlay" (default)
> > ! arg4  - longitude limits - specify as for any FERRET longitude range
> > ! arg5  - latitude  limits - specify as for any FERRET latitude  range
> >  
> > ! examples:
> > !	  GO fland_levitus 	          ! default - color gray, overlay
> > !       GO fland_levitus k=1            ! MASK for surface
> > !	  GO fland_levitus k=1 red        ! palette red.spk, overlay
> > !	  GO fland_levitus z=0:1000 red basemap ! color red, basemap
> > !	  GO fland_levitus z=1000 gray transparent ! color gray, transparent mask
> > !       GO fland k=1 gray overlay x=30e:100w y=20s:20n
> > 
> > ! Note : The transparent option creates a see-through overlay land mask
> > 
> > ! check the overlay/basemap/transparent option
> > query/ignore $3%overlay|overlay|basemap|transparent<3rd argument of GO fland_levitus must be "overlay", "basemap", or "transparent"%
> > 
> > ! save data set
> > set data/save
> > 
> > ! save region
> > define region/default levitus
> > 
> > ! set data for MASK
> > use levitus_mask.nc
> > 
> > set region/l=1
> > 
> > shade/nolab/nokey/levels="(-0.5,0.5,1.0)"$1" |*>/*"/pal=$2"gray"$3"/overlay|basemap> |*>/overlay"$4" |*>/*"$5" |*>/*" $3"MASK|transparent>if mod(i+j,2) eq 0 then MASK|*>MASK"
> > 
> > ! cancel MASK data set
> > cancel data levitus_mask.nc
> > 
> > ! restore region
> > set region levitus
> > 
> > ! restore original data set
> > set data/restore
> > 
> > set mode/last verify
> > 
> > ! END of fland_levitus.jnl
> > ----------------------------------------------------------------------------
> > 
> > On Wed, 1 Mar 2000, Yang Kun wrote:
> > 
> > > Hi, Dear ferret_users,
> > > 
> > > I am a new ferret_beginner.
> > > I contoured the temperature in Y-Z plate and want to
> > > fill the land in depth, how can I do it? I use a
> > > variable in my netcdf file as a mask which use missing
> > > value(-999.0) as land and 1 as sea water.
> > > Thank you for the help.
> > > 
> > > yours,
> > > Kun
> > > 
> > > =====================================================
> > > Kun Yang,   Researching Assistant
> > > Dept. of Oceanography, Texas A&M Uni. College Station
> > > Mailing Addr:     806F, Oran Circle, Bryan, TX, 77801
> > > Tel:  (409)-268-0963 (h), (409)-845-3238 (o)
> > > Email addr:    yangk@ocean.tamu.edu; yangk@netease.com
> > > =====================================================
> > > 
> > 
> > D. Shankar                           Ph. (O): [91](832) 22 6253 * 4400/4224
> > Physical Oceanography Division,                         22 1322 * 4400/4224
> > National Institute of Oceanography,      (R): [91](832) 23 8208      
> > Dona Paula,                          Fax    : [91](832) 22 3340/1360/9102
> > Goa 403 004,                         email  : shankar@csnio.ren.nic.in
> > India.                                        shankar@darya.nio.org
> > 
> > 
> 
> =====================================================
> Kun Yang,   Researching Assistant
> Dept. of Oceanography, Texas A&M Uni. College Station
> Mailing Addr:     806F, Oran Circle, Bryan, TX, 77801
> Tel:  (409)-268-0963 (h), (409)-845-3238 (o)
> Email addr:    yangk@ocean.tamu.edu; yangk@netease.com
> =====================================================
> 

D. Shankar                           Ph. (O): [91](832) 22 6253 * 4400/4224
Physical Oceanography Division,                         22 1322 * 4400/4224
National Institute of Oceanography,      (R): [91](832) 23 8208      
Dona Paula,                          Fax    : [91](832) 22 3340/1360/9102
Goa 403 004,                         email  : shankar@csnio.ren.nic.in
India.                                        shankar@darya.nio.org



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement