Re: about color key and plot a rectangular

[Jaison Kurian <jaison@xxxxxxxxxxxxxxxxxx>]

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
___________________________________________________
\ 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/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
 
 contour/nolabel/line=2/lev=(-10000,-9000,1000)/over/set rose
   ppl ylab LATITUDE;ppl xlab LONGITUDE;
 ppl contour/over
 go fland 5

! mark mariana trench by a box in red color using box.jnl

  go box 141.6,143.2,10.8,11.9 2

! OR mark mariana trench by a box in red color using "plot" command

  ! let mx = {141.6, 143.2, 143.2, 141.6, 141.6}
  ! let my = { 10.8,  10.8,  11.9,  11.9,  10.8}
  ! plot/ov/vs/nolab/line/color=2 mx, my

 sp rm -f map.gif
 frame/file=map.gif