Re: Curvilinear data, mp_fland, and mp_graticule...

[Jonathan Callahan <callahan@xxxxxxxxxxxxx>]

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