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

To: Steve Hankin <hankin@xxxxxxxxxxxxx>
Subject: Re: How to plot the topography in Y-Z plate?
From: Yang Kun <yangk@xxxxxxxxxxxxxx>
Date: Sun, 5 Mar 2000 16:05:03 -0600 (CST)
Cc: D Shankar <shankar@xxxxxxxxxxxxx>, ferret_users <ferret_users@xxxxxxxxxxxxxxxxxxx>
In-reply-to: <38BFDD6B.D1052444@pmel.noaa.gov>
Sender: owner-ferret_users@xxxxxxxxxxxxxxxxxxx
Dear Shankar and Steve

Thanks a lot!

I finally solved the problem by using the second way as Steve
suggested. It works well. 

As for fland_levitus.jnl, the built-in script to plot is like:

shade/nolab/nokey/levels="(-0.5,0.5,1.0)"/z=1:5000/pal=gray/overlay/x=180e:180e/y=90s:90n LEVI_MASK

this can be modified as
shade/nolab/nokey/levels="(-0.5,0.5,1.0)"/z=1:5000/pal=gray/overlay/y=90s:90n
LEVI_MASK[x=120e:50w@ave]

It can work as this way.

I used what looks like in two_dee_plot_of_1d to do this but using 
my_depth as Steve suggested, and it works. As for the first solution,
because the model resolution is not high, so the model topography
looks not so nice. 

Thanks again!!

yours,
Kun Yang

 
On Fri, 3 Mar
2000, Steve Hankin wrote:

> Hi Yang and Shankar,
> 
> Here are a  couple of additional techniques to think about.  I am assuming that your variable ("my_data") currently has missing value flags (-999.0) over land points
> 
> Simplest:
>     yes? LET one_zero_mask = MISSING(0*my_data, 1)
>     yes? SHADE/OVERLAY/PALETTE=black/NOLABEL/LEV=(0.5,1.5,1) one_zero_mask
>     or use "FILL" instead of "SHADE" for a smoother appearance
> 
> Another nice technique for overlaying a bathymetry mask is illustrated in
> 
>     yes? GO two_dee_plot_of_1d
> 
>     A GIF image using this technique is attached to this email. For this approach you need to
>     have the depths as a 1d field. An example would be
> 
>     yes? use etopo60
>     yes? LET my_depths = rose[x=160W,d=etopo60]
> 
>     If you want your depth mask to represent average depth from 160E to 160W simply use
> 
>     yes? LET my_depths = rose[x=160E:160W@AVE,d=etopo60]
> 
>     - steve
> 
> =========================================================
> 
> Yang Kun wrote:
> 
> > Dear Shankar
> >
> > Thanks for your scripts!
> > I generated a levitus_mask,and it works well.
> > I also generated a mask file with my model data, it also
> > works but the resolution is somewhat too coarse.
> >
> > 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?
> >
> > Thanks a lot!
> >
> > yours,
> > Kun
> >
> > On Thu, 2 Mar
> > 2000, D Shankar
> > wrote:
> >
> > > Hi,
> > > The following journal file is what we use to overlay a land-sea mask for
> > > Levitus data.  It works for vertical sections too.
> > > 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
> > =====================================================
> 
> --
> Steve Hankin
> NOAA/PMEL, 7600 Sand Point Way NE, Seattle, WA 98115-0070
> ph. (206) 526-6080 -- FAX (206) 526-6744
> 
> 
> 

=====================================================
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
=====================================================
References:
- Re: How to plot the topography in Y-Z plate?
  - From: Steve Hankin
Previous by thread: Re: How to plot the topography in Y-Z plate?
Next by thread: Vector into array...
[Thread Prev][Thread Next][Index]
Dept of Commerce / NOAA / OAR / PMEL / TMAP