[gradsusr] GrADS help, plotting shaded contour precip/ptype products over a basemap/background
Jeff Chabot
jsc219 at gmail.com
Wed Nov 10 23:11:28 EST 2010
Hello again GrADS Users,
I solved the background/basemap under a shaded contour. And here is how I
did it:
Please note: A Google base map did not work, either due to a format issue or
a dimension issue.
1.) Create a basemap in GrADS. There are several ways to do this, I found
two methods.
a.) Use the Library Script: basemap.gs tool to mask out land/oceans
regions to create a basemap and use the printim command export as an image.
b.) The other method was an idea from Stephen (thanks Stephen) to use
shapefiles to produce the map.
I have attached sample code (.gs files) for each option.
2.) When you create your product, use the following printim command to use
the background basemap under your shaded contour:
'printim ptype_NE_'.png png -t 0 x570 y435 -b NE_ptype.png white'
Be sure to use the exact same size file, same lat/lon, scale, parea, etc.
I used the following lines for this example, both the basemap and the final
product:
'set mpdset hires'
'set mproj scaled'
'set parea 0.5 10.5 1. 7.5'
'set lat 36.5 51.5'
'set lon '-90.0 -63'
I have also attached the blank basemap and an example of the final product.
Thanks again to all,
Jeff
On Wed, Nov 10, 2010 at 1:45 PM, <gradsusr-request at gradsusr.org> wrote:
> Send gradsusr mailing list submissions to
> gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
> gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
> gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
> 1. vertical cross sections (Jeffrey Duda)
> 2. Re: gradsusr Digest, Vol 9, Issue 23
> (Juan Jose Amides Figueroa Urbano)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 10 Nov 2010 16:02:14 -0600
> From: Jeffrey Duda <jdduda at iastate.edu>
> Subject: [gradsusr] vertical cross sections
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <AANLkTimafbMUp=OGARs6rRKknMyfWdL7q29T5QGs_W8u at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Grads users,
> Is there a way to plot a vertical cross section that isn't aligned
> meridionally or zonally (i.e., posessing some sort of NW-SE or SW-NE etc
> orientation)?
>
> Jeff Duda
>
> --
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101110/8938ced1/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Wed, 10 Nov 2010 16:50:14 -0600
> From: Juan Jose Amides Figueroa Urbano <juanjosefigueroa at gmail.com>
> Subject: Re: [gradsusr] gradsusr Digest, Vol 9, Issue 23
> To: gradsusr at gradsusr.org
> Message-ID:
> <AANLkTimTMFJCUvmwhNkhjkWp0XoSd7zf825VTrRvVf-9 at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Thanks for the help!
>
> On Tue, Nov 9, 2010 at 5:28 AM, <gradsusr-request at gradsusr.org> wrote:
>
> > Send gradsusr mailing list submissions to
> > gradsusr at gradsusr.org
> >
> > To subscribe or unsubscribe via the World Wide Web, visit
> > http://gradsusr.org/mailman/listinfo/gradsusr
> > or, via email, send a message with subject or body 'help' to
> > gradsusr-request at gradsusr.org
> >
> > You can reach the person managing the list at
> > gradsusr-owner at gradsusr.org
> >
> > When replying, please edit your Subject line so it is more specific
> > than "Re: Contents of gradsusr digest..."
> >
> >
> > Today's Topics:
> >
> > 1. Re: netCDF control file (Heiner K?rnich)
> > 2. Cuestion (Juan Jose Amides Figueroa Urbano)
> > 3. Re: Cuestion (RafaelCB)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Tue, 9 Nov 2010 14:33:08 +0100
> > From: Heiner K?rnich <heiner at misu.su.se>
> > Subject: Re: [gradsusr] netCDF control file
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <AANLkTi=THDPRVWHVfG0Kq0nRHnC4mBigtHQqPH=uqC7W at mail.gmail.com>
> > Content-Type: text/plain; charset="windows-1252"
> >
> > Hi Deniz,
> >
> > ok. When you tried to open the file with sdfopen and it failed, this
> means
> > that the file is not consistent with the coards-standard. So we need a
> > ctl-file and xdfopen.
> >
> > You need the name of all coordinates of the file. ncdump -h would give
> > that.
> > When you wrote "dimension1", this is just an example for the first
> > dimension. Probably, this should be the name of your longitude-variable
> in
> > matlab, eg. "long" or "x" or so.
> >
> > "fwrite" in matlab follows the logic of the programming language c. You
> > have
> > to open the file first with fopen. It is explained in the matlab-manual.
> > Here is an example from me:
> > % open the file 'ures.bin' in binary write-mode.
> > fid=fopen('ures.bin','wb');
> > % write binary data from the field winds to file in real*4 %precision
> > fwrite(fid,winds,'real*4');
> >
> > Regards,
> > Heiner
> >
> > On Tue, Nov 9, 2010 at 1:37 PM, Deniz Demirhan Bari <bari at iap-kborn.de
> > >wrote:
> >
> > > Dear Heiner,
> > >
> > > Grads 2.0 in our server is now changed and it is compiled with netCDF
> > > support.
> > >
> > > When I tried:
> > >
> > > >sdfopen ures.nc
> > >
> > > it says:
> > >
> > > >Scanning self-describing file: ures.nc
> > >
> > > gadsdf: SDF file has no discernable X coordinate.
> > >
> > > When I try:
> > >
> > > >xdfopen ures.ctl
> > >
> > > >Scanning Descriptor File: ures.ctl
> > >
> > > gadsdf: Lon dimension dimension1 is not an SDF dimension.
> > >
> > > SDF Descriptor file ures.ctl was not successfully opened & parsed.
> > >
> > >
> > >
> > > ---Matlab fwrite, I tried it but couldn?t manage.
> > >
> > > Thank you
> > >
> > > Deniz
> > >
> > > *Von:* gradsusr-bounces at gradsusr.org [mailto:
> > gradsusr-bounces at gradsusr.org]
> > > *Im Auftrag von *Heiner K?rnich
> > > *Gesendet:* Montag, 8. November 2010 17:07
> > >
> > > *An:* GrADS Users Forum
> > > *Betreff:* Re: [gradsusr] netCDF control file
> > >
> > >
> > >
> > > Hi Deniz,
> > >
> > > when you write recognized "sdfopen" or "xdfopen", do you mean that
> grads
> > > does not know the command? Is your grads-version compiled with
> > > netcdf-support?
> > >
> > > It is certainly a good thing to create netcdf-files, as they are
> > > self-descriptive. But still you could write binary data from matlab
> with
> > the
> > > command fwrite:
> > > http://www.mathworks.com/help/techdoc/ref/fwrite.html
> > >
> > > Then you have to follow the logic of grads binary files, ie. first
> > > longitude, then latitude, then height, then time, as described in the
> > > manual.
> > >
> > > So far,
> > > Heiner
> > >
> > > On Mon, Nov 8, 2010 at 4:25 PM, Deniz Demirhan Bari <bari at iap-kborn.de
> >
> > > wrote:
> > >
> > > Hello Heiner,
> > >
> > > Thank you for your mail. I tried 'sdfopen ures.nc' but neither grads
> > 2.0
> > > nor grads 1.8 recognized sdfopen. Then I prepared a control file as
> > below,
> > > by looking at the ?CDO sinfo ures.nc?. Then tried ?xdfopen ures.ctl?
> > but
> > > none of the grads versions recognized xdfopen again.
> > >
> > > *
> > >
> > > DSET ^ures.nc
> > >
> > > DTYPE netcdf
> > >
> > > UNDEF -9e+33
> > >
> > > XDEF dimension1 96 LINEAR 0.000000 3.750000
> > >
> > > YDEF dimension2 48 LEVELS -87.159 -83.479 -79.777 -76.070 -72.362
> > -68.652
> > >
> > > -64.942 -61.232 -57.521 -53.810 -50.099 -46.389
> > >
> > > -42.678 -38.967 -35.256 -31.545 -27.833 -24.122
> > >
> > > -20.411 -16.700 -12.989 -9.278 -5.567 -1.856
> > >
> > > 1.856 5.567 9.278 12.989 16.700 20.411
> > >
> > > 24.122 27.833 31.545 35.256 38.967 42.678
> > >
> > > 46.389 50.099 53.810 57.521 61.232 64.942
> > >
> > > 68.652 72.362 76.070 79.777 83.479 87.159
> > >
> > > ZDEF dimension3 114 LINEAR 1 1
> > >
> > > TDEF dimension4 1 LINEAR 00:00Z00jan0000 1mn
> > >
> > > TITLE ures.grb T31 grid
> > >
> > > VARS 1
> > >
> > > ures 114 1,100 ures
> > >
> > > ENDVARS
> > >
> > > *
> > >
> > >
> > >
> > > My data is three dimensionaI and if I try to save my data in an ASCII
> > file
> > > in MATLAB it gives a two dimensional output. I also couldn?t solve this
> > > problem
> > >
> > > Thank you
> > >
> > > Deniz
> > >
> > >
> > >
> > >
> > >
> > > *Von:* gradsusr-bounces at gradsusr.org [mailto:
> > gradsusr-bounces at gradsusr.org]
> > > *Im Auftrag von *Heiner K?rnich
> > > *Gesendet:* Montag, 8. November 2010 14:06
> > > *An:* GrADS Users Forum
> > > *Betreff:* Re: [gradsusr] netCDF control file
> > >
> > >
> > >
> > > Hi Deniz,
> > >
> > > you don't nead a ctl-file to open a netcdf-file in grads. If your file
> > > follows the COARDS-standard, then you can open it with 'sdfopen
> ures.nc
> > '.
> > > If grads won't recognise it, you might want to try an older
> > grads-version.
> > > We had some problems here with grads 2.0 and netcdf.
> > >
> > > If grads still won't open your file, you can write a ctl-file and try
> to
> > > open that with the command 'xdfopen'. The ctl-file looks a bit
> different
> > > from standard ctl-files, look here:
> > > http://www.iges.org/grads/gadoc/gradcomdxdfopen.html
> > >
> > > Then, you would need the output from ncdump in order to know what is
> > inside
> > > the netcdf-file. Alternatively you can use the climate data operators
> > (CDO)
> > > with the command "cdo sinfo ures.nc".
> > > Or in matlab, if you have the mexnc-package installed, you read the
> > > information with
> > > info=nc_info('ures.nc')
> > > Then you can access all details in info.
> > >
> > > Finally, if it is only you who will use the data, it might be simpler
> to
> > > write binary output in matlab and use a standard ctl.
> > >
> > > Kind regards from Stockholm,
> > > Heiner
> > >
> > > On Mon, Nov 8, 2010 at 10:11 AM, Deniz Demirhan Bari <
> bari at iap-kborn.de>
> > > wrote:
> > >
> > > Dear Grads users,
> > >
> > > I have calculated residual wind components in MATLAB and saved them in
> > > netCDF file. Now I want to read them in grads but I cannot create a
> > netcdf
> > > control file.
> > >
> > > I tried to get the output of ncdump but my computer doesn?t recognize
> > ncdump although I installed it.
> > >
> > > I am using the latest version of Grads on a Suse Linux computer.
> > >
> > > My attempt to create a ctl file is below:
> > >
> > > *
> > >
> > > DSET ^ures.nc
> > >
> > > DTYPE netCDF
> > >
> > > XDEF 96 LINEAR 0.000000 3.750000
> > >
> > > YDEF 48 LEVELS -87.159 -83.479 -79.777 -76.070 -72.362 -68.652
> > >
> > > -64.942 -61.232 -57.521 -53.810 -50.099 -46.389
> > >
> > > -42.678 -38.967 -35.256 -31.545 -27.833 -24.122
> > >
> > > -20.411 -16.700 -12.989 -9.278 -5.567 -1.856
> > >
> > > 1.856 5.567 9.278 12.989 16.700 20.411
> > >
> > > 24.122 27.833 31.545 35.256 38.967 42.678
> > >
> > > 46.389 50.099 53.810 57.521 61.232 64.942
> > >
> > > 68.652 72.362 76.070 79.777 83.479 87.159
> > >
> > > ZDEF 114 LINEAR 1 1
> > >
> > > TDEF 1 LINEAR 00:00Z00jan0000 1mn
> > >
> > > OPTIONS yrev
> > >
> > > UNDEF -9e+33
> > >
> > > VARS 1
> > >
> > > ures 114 1,100 ures
> > >
> > > ENDVARS
> > >
> > > *
> > >
> > > Thanks in advance
> > >
> > > Deniz DEMIRHAN BARI
> > >
> > > Leibniz-Institut f?r Atmosph?renphysik e.V.
> > > an der Universit?t Rostock
> > > Schlossstra?e 6
> > > 18225 K?hlungsborn
> > >
> > >
> > >
> > >
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > >
> > > --
> > > Heiner K?rnich
> > > Dept. of Meteorology Tel: +46 8 164333
> > > Stockholms University, SE-106 91 Stockholm, Sweden
> > > Email: heiner at misu.su.se www.misu.su.se/~heiner/<http://www.misu.su.se/%7Eheiner/>
> <http://www.misu.su.se/%7Eheiner/>
> > <http://www.misu.su.se/%7Eheiner/>
> > >
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> > >
> > >
> > > --
> > > Heiner K?rnich
> > > Dept. of Meteorology Tel: +46 8 164333
> > > Stockholms University, SE-106 91 Stockholm, Sweden
> > > Email: heiner at misu.su.se www.misu.su.se/~heiner/<http://www.misu.su.se/%7Eheiner/>
> <http://www.misu.su.se/%7Eheiner/>
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> >
> >
> > --
> > Heiner K?rnich
> > Dept. of Meteorology Tel: +46 8 164333
> > Stockholms University, SE-106 91 Stockholm, Sweden
> > Email: heiner at misu.su.se www.misu.su.se/~heiner/<http://www.misu.su.se/%7Eheiner/>
> <http://www.misu.su.se/%7Eheiner/>
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> >
> http://gradsusr.org/pipermail/gradsusr/attachments/20101109/ca3bf19e/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 2
> > Date: Tue, 9 Nov 2010 09:16:11 -0600
> > From: Juan Jose Amides Figueroa Urbano <juanjosefigueroa at gmail.com>
> > Subject: [gradsusr] Cuestion
> > To: gradsusr at gradsusr.org
> > Message-ID:
> > <AANLkTinzWXC9M-NGiuTXZ=qg9cu2=ayukRurzhjkA_ga at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Hello, i new in GrADS, i need a list of available diagnostics or
> productos
> > to see in GrADS from WRF model, i try whit rainc for rain and pressure
> ("d
> > rainc" and "d pressure") can you help me?
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> >
> http://gradsusr.org/pipermail/gradsusr/attachments/20101109/bf5c6f70/attachment-0001.html
> >
> > ------------------------------
> >
> > Message: 3
> > Date: Tue, 9 Nov 2010 13:33:45 -0200
> > From: RafaelCB <rafael.ciram at gmail.com>
> > Subject: Re: [gradsusr] Cuestion
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <AANLkTi=nRm3XEFHb_kqMe+w-KoLtYSEns5sb3b5y8v_y at mail.gmail.com<nRm3XEFHb_kqMe%2Bw-KoLtYSEns5sb3b5y8v_y at mail.gmail.com>
> <nRm3XEFHb_kqMe%2Bw-KoLtYSEns5sb3b5y8v_y at mail.gmail.com<nRm3XEFHb_kqMe%252Bw-KoLtYSEns5sb3b5y8v_y at mail.gmail.com>
> >
> > >
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Hello
> > In WRF has two types of precipitation
> > rainc 0 0 TOTAL CUMULUS Accumulated Precipitation (cm)
> > rain 0 0 TOTAL GRID SCALE Accumulated Precipitation (cm)
> > You need to add the two variables to have the total precipitation that
> > time.
> >
> > * Precipitation
> > set gxout shaded
> > d rain+rainc
> >
> > * Pressure
> > contour set gxout
> > d slp
> >
> > hug.
> >
> > 2010/11/9 Juan Jose Amides Figueroa Urbano <juanjosefigueroa at gmail.com>
> >
> > > Hello, i new in GrADS, i need a list of available diagnostics or
> > productos
> > > to see in GrADS from WRF model, i try whit rainc for rain and pressure
> > ("d
> > > rainc" and "d pressure") can you help me?
> > >
> > > _______________________________________________
> > > gradsusr mailing list
> > > gradsusr at gradsusr.org
> > > http://gradsusr.org/mailman/listinfo/gradsusr
> > >
> > >
> >
> >
> > --
> > Rafael Censi Borges
> > T?c.Meteorologia - Epagri/Ciram
> > +55 (48)3239-8064/8066
> > http://ciram.epagri.sc.gov.br
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL:
> >
> http://gradsusr.org/pipermail/gradsusr/attachments/20101109/f5053ea2/attachment.html
> >
> > ------------------------------
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> > End of gradsusr Digest, Vol 9, Issue 23
> > ***************************************
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101110/9d09d91b/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 9, Issue 29
> ***************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101110/96dfd457/attachment-0003.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ne.png
Type: image/png
Size: 7774 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20101110/96dfd457/attachment-0006.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: precip_NE_2.png
Type: image/png
Size: 15104 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20101110/96dfd457/attachment-0007.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nebasemap.gs
Type: application/octet-stream
Size: 419 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20101110/96dfd457/attachment-0006.obj
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nebackmap.gs
Type: application/octet-stream
Size: 1456 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20101110/96dfd457/attachment-0007.obj
More information about the gradsusr
mailing list