[gradsusr] Map projection difficulty
Wesley Ebisuzaki
Wesley.Ebisuzaki at noaa.gov
Mon Sep 20 15:59:42 EDT 2010
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>
>
> 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>
>
>
> 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>
>
> _______________________________________________
> 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
>
More information about the gradsusr
mailing list