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