Re: Re: about color key and plot a rectangular

["Rongzhen Gao" <gaorz@xxxxxxxxxxx>]

Hi, Jaison Kurian,

 

You did me a big favor. Thanks.

Why i couldn't find out it in the Ferret user guide.

 

Question 1: sp rm -f map.gif    >>>>i wonder why you add this command befor "frame/file=map.gif".

Question 2: pls see the attached jnl file. when i plot the upper view it is ok, but when i plot the upper and lower view together there is a question. As if the etopo5 data file closed after the first plot by itself, and when plot the second figure it does not exist. It is so strange!  Have you ever encountered it?

--------------------------------error massage----------------------------------  

yes? go d:\wind
Re-defining viewport V11
 **ERROR: unknown data set: D=1
CONTOUR/FILL/nolabel/nokey/lev=(-10000,0,1000)/pal=bluescale/set rose[d=1]
Command file, command group, or REPEAT execution aborted
----------------------------------               

 

Gao Rongzhen

---

Hi Gao Rongzhen,

                  I made few modifications to your script and is 

working fine for me (on my pc, Fer5.8). Have a trial with the modified

script (please see the attachment color_key1.jnl).

 

   You can get a box overlayed on your plot in many ways. Two methods are

illustrated in color_key1.jnl. 

 

 1. use box.jnl 

 

        usage : GO box xlo xhi ylo yhi [pen_number]

    

        for details see this file /usr/local/fer5.8/go/box.jnl

    

        this is the most easy method    

 

 2. use plot command

 

        define the corner points and plot it using the qualifiers

 

           /ov/vs/line

 

 Please let me know if you need more explanation.

 

With Regards 

 

Jaison

 

 

On Tue, 11 Jan 2005, Rongzhen Gao wrote:

 

>  Hi, All

>  

>  I'm trying to plot topography by fill command, why the color of the shading key is not continuous? 

>  Pls look at the attached figure. My problem is why the color from -9000 to -10000 is gray not blue? Is there a way to solve it?

>  

>  My second question is how to plot a rectangular over the Marianas trench?  Thanks in advance!

>  

>  Ferret version is 5.5.

>  Gao Rongzhen

>  

>  ---------------------------map.jnl  i used-------------------------------------

>  \ cancel mode verify

>    cancel region;  cancel data/all;cancel var/all;cancel sym/all

>  

>    def view/xlim=0.05,0.95/ylim=0.2,0.8 v11

>    def sym lon1=105;def sym lon2=165

>    def sym lat1=0;  def sym lat2=40

>    def sym dx=10;  def sym dy=5

>    def sym pall=bluescale

>    def sym levv=(-10000,0,1000)

>    set window/size=1/aspect=1.294

>  

>    set data etopo5

>    set view v11

>    set region/x=($lon1):($lon2)/y=($lat1):($lat2)

>    fill/nolabel/nokey/lev=($levv)/pal=($pall)/set rose

>    ppl xaxis ($lon1),($lon2),($dx); ppl yaxis ($lat1),($lat2),($dy)

>    ppl labset,,0.16,0.16;ppl axlsze 0.16,0.16;ppl tics .1,.12,.1,.12;ppl axlint 1,1  

>    let ylo = ($ppl$yorg);let yhi = ylo + ($ppl$ylen);let xlo = ($ppl$xorg) + ($ppl$xlen) + .2; let xhi = xlo + .2

>    ppl shakey 1, 1, 0.12, 0, 5, 6, `xlo`, `xhi`, `ylo`, `yhi`; ppl fill;ppl shaset reset 

>    contour/nolabel/line=2/lev=(-10000,-9000,1000)/over/set rose;ppl contour/over

>    ppl xaxis ($lon1),($lon2),($dx); ppl yaxis ($lat1),($lat2),($dy);  ppl labset,,0.16,0.16; ppl axlsze 0.16,0.16;ppl tics .1,.12,.1,.12;ppl axlint 1,1

>    ppl ylab LATITUDE;ppl xlab LONGITUDE;

>    ppl contour

>    go fland 5  

>  

>   

>    frame/file=map.gif

>  

 

-- 

___________________________________________________

 

    Jaison Kurian                           

    Centre for Atmospheric and Oceanic Sciences

    Indian Institute of Science

    B A N G A L O R E   560 012

    Ph: +91-80-3942505

        +91-80-3600450

    Fax:+91-80-3600865

___________________________________________________

 

---

Attachment: wind.jnl
Description: Binary data

Attachment: face-2.GIF
Description: GIF image