[Thread Prev][Thread Next][Index]
Re: about color key and plot a rectangular
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
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement