[Thread Prev][Thread Next][Index]
Re: Curvilinear data, mp_fland, and mp_graticule...
Gary,
Part of this is the same kind of problem described in the FAQ:
http://ferret.wrc.noaa.gov/Ferret/FAQ/graphics/working_with_go_fland.html
The mp_fland.jnl script has the same arguments as fland.jnl so you
should be able to fix the land mask with
> go mp_fland 60 black " " " " x=0:360
I'd also try playing round with setting the region specifically to
20:380 when you plot your data because those are the settings you pass
to the map projection script.
If you don't get it working you can send me a single XY slice of your
data files and I'll have a go at it.
-- Jon Callahan
Gary Strand wrote:
>
> 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:
>
> ------------------------------------------------------------------------
> Name: POP_CE.gif
> POP_CE.gif Type: GIF Image (image/gif)
> Encoding: base64
> 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