[gradsusr] Map projection difficulty
Jacob Carley
jrcarley at gmail.com
Mon Sep 20 16:15:30 EDT 2010
Hi Wesley and Jeff,
Thanks for your help. Tweaking the xdef/ydef fixed it!
-Jacob
On Mon, Sep 20, 2010 at 3:59 PM, Wesley Ebisuzaki <Wesley.Ebisuzaki at noaa.gov
> wrote:
> Hi Jacob,
>
> The program grib2ctl.pl uses a very simple method for determine the
> display grid. To see the corners,
> you'll have to tweek the definition of the display grid. Look for the
> "xdef" and "ydef" lines in the
> ctl file and change it to something that you like better. There is no
> "magic" in the numbers;
> you can change the nx/ny, x0/y0 and dx/dy.
>
> xdef nx linear x0 dx x0, dx are in degrees
>
> Wesley Ebisuzaki
>
>
> Jacob Carley wrote:
> > Hi Jeff,
> >
> > Thanks for your reply. The northwest corner is fine, but the GrADS
> > won't plot the part of the grid that is any farther south (like the SE
> > corner). I suppose a better example would be:
> >
> > http://web.ics.purdue.edu/~jrcarley/LCC_0100_03_ex2.png<http://web.ics.purdue.edu/%7Ejrcarley/LCC_0100_03_ex2.png>
> > <http://web.ics.purdue.edu/%7Ejrcarley/LCC_0100_03_ex2.png>
> >
> > Thanks,
> > Jacob
> >
> > On Mon, Sep 20, 2010 at 3:27 PM, Jeffrey Duda <jdduda at iastate.edu
> > <mailto:jdduda at iastate.edu>> wrote:
> >
> > Jacob,
> > I also use WRF data that is on a Lamber Conformal projection, but
> > I don't have those problems when I view the data in Grads. The
> > picture you attached looks like the northwest corner of your
> > domain, and by the look of it alone, it doesn't look like anything
> > is wrong with what you or Grads is doing. Did you mean to send a
> > different image?
> >
> > Jeff Duda
> >
> > On Mon, Sep 20, 2010 at 2:07 PM, Jacob Carley <jrcarley at gmail.com
> > <mailto:jrcarley at gmail.com>> wrote:
> >
> > 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<http://web.ics.purdue.edu/%7Ejrcarley/LCC_0100_03.png>
> > <http://web.ics.purdue.edu/%7Ejrcarley/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 <http://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 <http://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
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> >
> >
> > --
> > Jeff Duda
> > Iowa State University
> > Meteorology Graduate Student
> > 3134 Agronomy Hall
> > www.meteor.iastate.edu/~jdduda<http://www.meteor.iastate.edu/%7Ejdduda>
> > <http://www.meteor.iastate.edu/%7Ejdduda>
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org <mailto:gradsusr at gradsusr.org>
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100920/e860e5d6/attachment-0003.html
More information about the gradsusr
mailing list