[Thread Prev][Thread Next][Index]

Re: polar stuff



----------
X-Sun-Data-Type: text
X-Sun-Data-Description: text
X-Sun-Data-Name: text
X-Sun-Charset: us-ascii
X-Sun-Content-Lines: 106


Thanks for you help Ansley,  Now I don't have a problem with the 
actual plot, because I now define my grid correctly, but I run into problems
with go mp_land when using sterographic_south (but not with north)
Here is my script and error...
set data "Have_con_stero.cdf"
set data "Havex1_stero.cdf"
let diff=H_avexp[d=Havex1_stero]-H_ave[d=Have_con_stero]
set region/y=-90:0/z=100
load diff
show grid
 Default grid for DEFINE VARIABLE is ABSTRACT
 Last successful data access was on grid GTJ21
    GRID GTJ1
name		axis		#pts	start			end
LONGITUDE	LONGITUDE	96mr 	0E			3.75w		
LATITUDE	LATITUDE	73 r 	90S			90N
P 		HEIGHT (mbar)	20 i-	0.08			1000
TANNAVE	  	TIME		12 r 	16-JAN 		00.00	16-DEC	00.00
show variable diff
 DIFF=H_avexp[d=Havex1_stero]-H_ave[d=Have_con_stero]
set grid diff
go stereographic_south
go margins 0.2,0.2,1,0.6
fill/noaxis/set_up diff[z=100,l=1:12@ave], x_page, y_page
ppl fill
go mp_land
**ERROR:illegal limits: "MASKED_X_PAGE" is not in the range Y=-90:0
     Axis extremes are Y=0.5;10000000
plot/set/vs/nolab/line=1/i=1:22585/overlay masked_x_page, masked_y_page
Command file, command group, or REPEAT execution aborted

So, is there something else I should be using for the land mask when using
southern projections?
Thanks
Alice

----- Begin Included Message -----


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement