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

[Steve Hankin <hankin@xxxxxxxxxxxxx>]

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