[Thread Prev][Thread Next][Index]

Curvilinear data, mp_fland, and mp_graticule...




I've been hammering away at trying to get a picture of some curvilinear data 
properly overlaid with mp_graticule, and it just ain't working. The graticule 
lines quit at 0E, as does the mp_fland.

Also, the image longitudinal extent is that of the data itself, regardless of
how I tweak the arguments to mp_lambert_cyl. The jnl file:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
use "/data/pcm/B06.33/ocn/B06.33.ocn.2010-01.nc"
use "/fs/cgd/data0/strandwg/ferret/data/pcm1_grid.nc"
!
let plto = temp[d=1,l=1,k=1]
!
set mem/siz=64
!
let xmin =   20; let xmax = 380 ; let ymin = -90; let ymax =  90
let xdel =   30; let ydel =  15
let mid_lon = `xmin+(xmax-xmin)/2` ; let mid_lat = `ymin+(ymax-ymin)/2`
!
go mp_lambert_cyl `mid_lon` `mid_lat`
let mp_lambda = lon[d=pcm1_grid.nc] * deg2rad
let mp_phi    = lat[d=pcm1_grid.nc] * deg2rad
set grid `ignore0(plto),return=grid`
!
can region; set reg/j=1:288
shade/set/noax/nolabel/lev=(-4,34,2)/title="" plto*mp_mask, x_page, y_page
ppl shade
ppl axset 1 1 1 1
!
! background
go mp_lambert_cyl `mid_lon` `mid_lat`
set region/x=`xmin`:`xmax`/y=`ymin`:`ymax`
go mp_graticule `xmin` `xmax` `xdel` `ymin` `ymax` `ydel` 
go mp_fland 60 black
go mp_land 1
!
frame/file=POP_CE.gif
!
quit
-------------------------------------------------------------------------------

And POP_CE.gif:

Attachment: POP_CE.gif
Description: POP_CE.gif

--
Gary Strand
strandwg@ucar.edu


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement