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