[gradsusr] generate seasonal climatology, means and anomalies from monthly data
Mubashar Dogar
mubashardogar at gmail.com
Wed Dec 8 14:50:03 EST 2010
Hi saiguran,
First you need to convert your data in to required season and then just find
the climatology by averaging that season throughout the period. Then just
subtract this climatology from each season and plot them in form of time
series.
Regards,
Mubashar
On Tue, Dec 7, 2010 at 2:22 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. Re: PDEF problem, WRF (Matei Georgescu)
> 2. Re: PDEF problem, WRF (Eric Altshuler)
> 3. hello (namrata rajput)
> 4. script for moist static energy (Rashed Mahmood)
> 5. generate seasonal climatology, means and anomalies from
> monthly data (saiguran loisulie)
> 6. Re: script for moist static energy (Jeffrey Duda)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 6 Dec 2010 14:43:23 -0700
> From: Matei Georgescu <Matei.Georgescu at asu.edu>
> Subject: Re: [gradsusr] PDEF problem, WRF
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <AANLkTi=cTDDiG9CK3aCh5X3KudHpHAV-_-z_UY22x=2V at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> I've always used ARWPost to do this (creates Grads-readable file and
> associated control file) ...
>
> Please check here:
>
> http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/users_guide_chap9.htm#_ARWpost_2
>
> hope this helps,
> matei
>
>
> On Mon, Dec 6, 2010 at 2:37 PM, Song Feng <sfeng2 at unlnotes.unl.edu> wrote:
>
> > Hi All:
> >
> > I'd like to display the WRF output (in NetCDF format) using GrADS. The
> > projection used is Lambert Conformal. I am trying to create a 'CTL' file
> for
> > the output. The following are the PDEF options for the CTL file:
> >
> > *PDEF isize jsize LCC latref lonref iref jref Struelat Ntruelat slon** dx
> > dy*
> >
> > I can find the information for isize, jsize, latref, lonref, iref, jref,
> dx
> > and dy. *The problem is how can I find the information for "Struelat
> > Ntruelat slon"?* I am bugging this for a while. Can anybody give me some
> > clues on how to solve this problems?
> >
> >
> > The following are the WRF ncdump output:
> >
> > dimensions:
> > time = UNLIMITED ; // (1095 currently)
> > yc = 109 ;
> > xc = 134 ;
> > bnds = 2 ;
> > variables:
> > char Lambert_Conformal ;
> > Lambert_Conformal:grid_mapping_name =
> > "lambert_conformal_conic" ;
> > Lambert_Conformal:standard_parallel = 30., 60. ;
> > Lambert_Conformal:longitude_of_central_meridian = -97. ;
> > Lambert_Conformal:latitude_of_projection_origin = 47.5 ;
> > Lambert_Conformal:false_easting = 3325000. ;
> > Lambert_Conformal:false_northing = 2700000. ;
> > double yc(yc) ;
> > yc:long_name = "y-coordinate in Cartesian system" ;
> > yc:standard_name = "projection_y_coordinate" ;
> > yc:axis = "Y" ;
> > yc:units = "m" ;
> > double xc(xc) ;
> > xc:long_name = "x-coordinate in Cartesian system" ;
> > xc:standard_name = "projection_x_coordinate" ;
> > xc:axis = "X" ;
> > xc:units = "m" ;
> > double lon(yc, xc) ;
> > lon:long_name = "longitude" ;
> > lon:standard_name = "longitude" ;
> > lon:units = "degrees_east" ;
> > lon:axis = "X" ;
> > double lat(yc, xc) ;
> > lat:long_name = "latitude" ;
> > lat:standard_name = "latitude" ;
> > lat:units = "degrees_north" ;
> > lat:axis = "Y" ;
> > double time(time) ;
> > time:long_name = "time" ;
> > time:standard_name = "time" ;
> > time:units = "days since 1968-01-01 00:00:0.0" ;
> > time:axis = "T" ;
> > time:calendar = "365_day" ;
> > time:bounds = "time_bnds" ;
> > double time_bnds(time, bnds) ;
> > float tasmax(time, yc, xc) ;
> > tasmax:long_name = "Maximum Daily Surface Air
> Temperature"
> > ;
> > tasmax:standard_name = "air_temperature" ;
> > tasmax:units = "K" ;
> > tasmax:_FillValue = 1.e+20f ;
> > tasmax:missing_value = 1.e+20f ;
> > tasmax:coordinates = "lon lat level" ;
> > tasmax:grid_mapping = "Lambert_Conformal" ;
> > tasmax:cell_methods = "time: maximum (interval: 24
> hours)"
> > ;
> > tasmax:original_name = "T2" ;
> > tasmax:original_units = "K" ;
> > tasmax:level_desc = "2 m" ;
> > double level ;
> > level:standard_name = "height" ;
> > level:units = "m" ;
> > level:positive = "up" ;
> > level:axis = "Z" ;
> >
> > // global attributes:
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
> *======================================
> Matei Georgescu, Ph.D.*
> *School of Mathematical and Statistical Sciences and
> Center for Environmental Fluid Dynamics*
> *Arizona State University*
> *Phone*: (480) 965-3770
> *Fax*: (480) 965-8119
> *E-mail*: Matei.Georgescu at asu.edu
> *==================================**====*
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101206/8fe0a17e/attachment-0001.html
>
> ------------------------------
>
> Message: 2
> Date: Mon, 6 Dec 2010 17:03:11 -0500 (EST)
> From: Eric Altshuler <ela at cola.iges.org>
> Subject: Re: [gradsusr] PDEF problem, WRF
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <1040925276.11305.1291672991875.JavaMail.root at mail.iges.org>
> Content-Type: text/plain; charset=utf-8
>
> Hello Song Feng,
>
> >From the ncdump output, Struelat=30, Ntruelat=60, and slon=-97. These
> should be the same as the corresponding values you specified in your
> namelist.wps when you ran WPS.
>
> Best regards,
>
> Eric L. Altshuler
> Assistant Research Scientist
> Center for Ocean-Land-Atmosphere Studies
> 4041 Powder Mill Road, Suite 302
> Calverton, MD 20705-3106
> USA
>
> E-mail: ela at cola.iges.org
> Phone: (301) 902-1257
> Fax: (301) 595-9793
>
> ----- Original Message -----
> From: "Song Feng" <sfeng2 at unlnotes.unl.edu>
> To: gradsusr at gradsusr.org
> Sent: Monday, December 6, 2010 4:37:43 PM
> Subject: [gradsusr] PDEF problem, WRF
>
>
> Hi All:
>
> I'd like to display the WRF output (in NetCDF format) using GrADS. The
> projection used is Lambert Conformal. I am trying to create a 'CTL' file for
> the output. The following are the PDEF options for the CTL file:
>
> PDEF isize jsize LCC latref lonref iref jref Struelat Ntruelat slon dx dy
>
> I can find the information for isize, jsize, latref, lonref, iref, jref, dx
> and dy. The problem is how can I find the information for "Struelat Ntruelat
> slon"? I am bugging this for a while. Can anybody give me some clues on how
> to solve this problems?
>
>
> The following are the WRF ncdump output:
>
> dimensions:
> time = UNLIMITED ; // (1095 currently)
> yc = 109 ;
> xc = 134 ;
> bnds = 2 ;
> variables:
> char Lambert_Conformal ;
> Lambert_Conformal:grid_mapping_name = "lambert_conformal_conic" ;
> Lambert_Conformal:standard_parallel = 30., 60. ;
> Lambert_Conformal:longitude_of_central_meridian = -97. ;
> Lambert_Conformal:latitude_of_projection_origin = 47.5 ;
> Lambert_Conformal:false_easting = 3325000. ;
> Lambert_Conformal:false_northing = 2700000. ;
> double yc(yc) ;
> yc:long_name = "y-coordinate in Cartesian system" ;
> yc:standard_name = "projection_y_coordinate" ;
> yc:axis = "Y" ;
> yc:units = "m" ;
> double xc(xc) ;
> xc:long_name = "x-coordinate in Cartesian system" ;
> xc:standard_name = "projection_x_coordinate" ;
> xc:axis = "X" ;
> xc:units = "m" ;
> double lon(yc, xc) ;
> lon:long_name = "longitude" ;
> lon:standard_name = "longitude" ;
> lon:units = "degrees_east" ;
> lon:axis = "X" ;
> double lat(yc, xc) ;
> lat:long_name = "latitude" ;
> lat:standard_name = "latitude" ;
> lat:units = "degrees_north" ;
> lat:axis = "Y" ;
> double time(time) ;
> time:long_name = "time" ;
> time:standard_name = "time" ;
> time:units = "days since 1968-01-01 00:00:0.0" ;
> time:axis = "T" ;
> time:calendar = "365_day" ;
> time:bounds = "time_bnds" ;
> double time_bnds(time, bnds) ;
> float tasmax(time, yc, xc) ;
> tasmax:long_name = "Maximum Daily Surface Air Temperature" ;
> tasmax:standard_name = "air_temperature" ;
> tasmax:units = "K" ;
> tasmax:_FillValue = 1.e+20f ;
> tasmax:missing_value = 1.e+20f ;
> tasmax:coordinates = "lon lat level" ;
> tasmax:grid_mapping = "Lambert_Conformal" ;
> tasmax:cell_methods = "time: maximum (interval: 24 hours)" ;
> tasmax:original_name = "T2" ;
> tasmax:original_units = "K" ;
> tasmax:level_desc = "2 m" ;
> double level ;
> level:standard_name = "height" ;
> level:units = "m" ;
> level:positive = "up" ;
> level:axis = "Z" ;
>
> // global attributes:
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> ------------------------------
>
> Message: 3
> Date: Tue, 7 Dec 2010 00:51:41 -0800 (PST)
> From: namrata rajput <namrata_rajput47 at yahoo.com>
> Subject: [gradsusr] hello
> To: gradsusr at gradsusr.org
> Message-ID: <742465.55854.qm at web55408.mail.re4.yahoo.com>
> Content-Type: text/plain; charset="utf-8"
>
> hey there, m using AIRS data for the analysis related to my area of
> research, as this data are in hdf form how should i proceed further?
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101207/3a8a05bd/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Tue, 7 Dec 2010 18:12:57 +0800
> From: Rashed Mahmood <rashidcomsis at gmail.com>
> Subject: [gradsusr] script for moist static energy
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <AANLkTim0MdwjVp3L97UCkYmoioK6H6byns6ZBFGOXixq at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi all
>
> Did anyone calculate moist static energy and/or dry static energy.
> If yes, can you please share the script(s) related to it.
>
> Thanks
>
> Rashed
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101207/ca7051a1/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Tue, 7 Dec 2010 06:40:56 -0800 (PST)
> From: saiguran loisulie <saiguranloi at yahoo.com>
> Subject: [gradsusr] generate seasonal climatology, means and anomalies
> from monthly data
> To: gradsusr at gradsusr.org
> Message-ID: <764978.92085.qm at web111710.mail.gq1.yahoo.com>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi,
> I want to generate seasonal climatology, means and anomalies from monthly
> data.
> May I have a script or explanations on how to do it in GrADS!
> Thank you!
>
> Saiguran Loisulie
>
> Climate Systems Analysis Group
>
> Environmental and Geographical Science Department
>
> University of Cape Town
>
> Private Bag X3
>
> Rondebosch
>
> 7701
>
> SOUTH AFRICA
>
>
>
> CELL 0748 068 230
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101207/c4e7ddc5/attachment-0001.html
>
> ------------------------------
>
> Message: 6
> Date: Tue, 7 Dec 2010 09:28:59 -0600
> From: Jeffrey Duda <jdduda at iastate.edu>
> Subject: Re: [gradsusr] script for moist static energy
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <AANLkTin0gTHnaU69SEsKPpBR5zeRB0-dJp15ARehtNu+ at mail.gmail.com<AANLkTin0gTHnaU69SEsKPpBR5zeRB0-dJp15ARehtNu%2B at mail.gmail.com>
> >
> Content-Type: text/plain; charset="iso-8859-1"
>
> Rashed,
> The equations for moist and dry static energy are quite simple and all you
> need to do to compute them is run this line:
>
> 'define dse = gravity*height + Cp*T' (for dry static energy)
> 'define mse = dse + L*q' (for moist static energy)
>
> where "gravity" and "height" represent variables for the gravitational
> constant (9.8 m/s2) and geopotential height, Cp = 1005 J/Kg*K is the heat
> constant for dry air at constant pressure, T is temperature (in K), L is
> heat of phase change (generally vaporization), and q is moisture quantity
> (mixing ratio).
>
> Hope this helps.
>
> Jeff Duda
>
> On Tue, Dec 7, 2010 at 4:12 AM, Rashed Mahmood <rashidcomsis at gmail.com
> >wrote:
>
> > Hi all
> >
> > Did anyone calculate moist static energy and/or dry static energy.
> > If yes, can you please share the script(s) related to it.
> >
> > Thanks
> >
> > Rashed
> >
> > _______________________________________________
> > gradsusr mailing list
> > 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>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20101207/336e0f38/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 10, Issue 20
> ****************************************
>
--
Muhammad Mubashar Ahmad Dogar
Scientific Officer (Climatology Section),
Global Change Impact Studies Centre (GCISC),
NCP complex, Quaid-e-Azam University Campus,
Shahdra road, Islamabad, Pakistan.
Tel: +92519230226
Mob:+923315144169
email: mubashardogar at yahoo.com
mubashar.ahmad at gcisc.org.pk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101208/fe08dd57/attachment-0003.html
More information about the gradsusr
mailing list