[gradsusr] Map projection difficulty

[Jacob Carley]

Hello GrADS users,

I'm attempting to plot some data (already projected on a Lambert grid), and
am having problems where my plot cuts off the southern portion of of my
domain (which I have a feeling is related to the SW corner of my domain).
I've tried setting mpdraw off as well as trying various types of projections
but I end up with the same/similar problem each time.  Any help would be
greatly appreciated.  My script and CTL file header are below and I've
provided a link to the example plotting problem.  I'm also a new GrADS user,
so if this is a newbie mistake, I apologize in advance!

Link to problem plot: http://web.ics.purdue.edu/~jrcarley/LCC_0100_03.png


GrADS plotting routine:

'open LCC_0100_03.ctl'
'set grads off'
'set frame on'
'set clab off'
'set mpvals 266 272 29.0 34.0'
'set mproj latlon'
'set mpdset hires'
'set gxout shaded'
'run radar.gs'
'set rbcols 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40'
'set clevs 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75'
'd REFD1000m'
'run cbarnew.gs'
'set gxout vector'
'set arrscl 0.5 15'
'set arrlab on'
'd skip(UGRD10m,10,10);VGRD10m'
'draw title 1km AGL Reflectivity and 10m wind'
'printim LCC_0100_03.png white'


Control file header:

dset LCC_0100_03.grb
index LCC_0100_03.grb.idx
undef 9.999E+20
title LCC_0100_03.grb
*  produced by grib2ctl v0.9.12.5p39a
dtype grib 255
pdef 501 401 lcc 31.459000 -93.762000 1 1 50.000000 50.000000 -111.000000
1050 900
xdef 501 linear -93.762000 0.0111902692769373
ydef 401 linear 31.459000 0.00818181818181818
tdef 1 linear 13Z24apr2010 1mo
*  z has 46 levels, for prs
zdef 46 levels
1000 975 950 925 900 875 850 825 800 775 750 725 700 675 650 625 600 575 550
525 500 475 450 425 400 375 350 325 300 275 250 225 200 175 150 125 100 75
70 50 30 20 10 7 5 2
vars 203



Thanks!
Jacob
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100920/a46daaac/attachment-0003.html