<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Hi Jeff, <br>thank you. It is working for the dp= constant (50) in my case.<br><br>Can you tell me how to go about doing it if I have the pressure levels where dp is not constant. <br>Thanks again<br>geeta<br><br><br><br>> From: gradsusr-request@gradsusr.org<br>> Subject: gradsusr Digest, Vol 8, Issue 42<br>> To: gradsusr@gradsusr.org<br>> Date: Mon, 25 Oct 2010 10:02:52 -0400<br>> <br>> Send gradsusr mailing list submissions to<br>> gradsusr@gradsusr.org<br>> <br>> To subscribe or unsubscribe via the World Wide Web, visit<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> or, via email, send a message with subject or body 'help' to<br>> gradsusr-request@gradsusr.org<br>> <br>> You can reach the person managing the list at<br>> gradsusr-owner@gradsusr.org<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. Query Reg Extracting the data in Vertical Levels (Geeta Geeta)<br>> 2. correlation at each grid poing (dhaval.prajapati@iccsir.org)<br>> 3. Re: Query Reg Extracting the data in Vertical Levels<br>> (Jeffrey Duda)<br>> 4. Creating "Lines" Shapefiles (Stephen McMillan)<br>> 5. (no subject) (SRAVANI A)<br>> 6. Re: Max Value Coordinates (Charles Seman)<br>> 7. (no subject) (Javier G. Corripio)<br>> <br>> <br>> ----------------------------------------------------------------------<br>> <br>> Message: 1<br>> Date: Mon, 25 Oct 2010 08:53:09 +0000<br>> From: Geeta Geeta <geeta124@hotmail.com><br>> Subject: [gradsusr] Query Reg Extracting the data in Vertical Levels<br>> To: grads grads <gradsusr@gradsusr.org><br>> Message-ID: <BLU152-w18DEC934AD8300BC78E01D8A410@phx.gbl><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> <br>> Dear All, <br>> I am trying to extract the SPECIFIC HUMIDITY for all the levels in vertical (19 Levels from 1000 to 100hPa) from the Model. <br>> So I am Using tvrh2q function of Grads. Following is the ctl file. <br>> <br>> 'open test.ctl'<br>> *'set lev 1000 100'<br>> 'set t 1'<br>> * while (lev< 100)<br>> *zlev=1000<br>> 'set gxout fwrite'<br>> 'set fwrite test.dat'<br>> zlev =1000<br>> 'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'<br>> 'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))'<br>> ************<br>> zlev =950<br>> 'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'<br>> 'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))' <br>> *'disable fwrite'<br>> *****************<br>> zlev =900<br>> 'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'<br>> 'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))' <br>> 'disable fwrite'<br>> ***********************<br>> <br>> * endwhile<br>> <br>> I am Unable to Extract the data Using the While Statement in Grads, so I am defining the Zlev separately . <br>> <br>> Kindly Suggest. <br>> geeta<br>> <br>> <br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/86067578/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 2<br>> Date: Mon, 25 Oct 2010 10:00:47 +0000<br>> From: dhaval.prajapati@iccsir.org<br>> Subject: [gradsusr] correlation at each grid poing<br>> To: gradsusr@gradsusr.org<br>> Message-ID: <W58022227159441288000847@webmail25><br>> Content-Type: text/plain; charset="utf-8"<br>> <br>> <br>> Hi,<br>> <br>> I have IMD rainfall data for India ( 1 deg. by 1 deg. Resolution) for 1951 to 2007 in binary format . I have Global Climate Index for the same peiode in text format. Now I want to do correlation analysis at each grid point for rainfall and Climate index and want to plot it in grads.<br>> How to do it ?<br>> <br>> Thanks in Advance,<br>> <br>> DhavalP. Prajapati<br>> Indian Centre for Climate and Societal Impacts Research (ICCSIR)<br>> Ahmedabad Education Society compound,<br>> Opp. H.L. Commerce College,<br>> Navrangpura,<br>> Ahmedabad 380 009,<br>> Gujarat,<br>> India<br>> Contact No. : +91(079)40045462<br>> <br>> <br>> <br>> <br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/2132346d/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 3<br>> Date: Mon, 25 Oct 2010 09:03:19 -0500<br>> From: Jeffrey Duda <jdduda@iastate.edu><br>> Subject: Re: [gradsusr] Query Reg Extracting the data in Vertical<br>> Levels<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <AANLkTik3A7wrCTsc0HhZJT72xtdu4T+j+3FHfaDTCcFh@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> Geeta,<br>> The documentation on the tvrh2q function says that it will assume it is at<br>> the pressure level of the current dimension environment. I don't see any<br>> "set lev" commands that are uncommented, so one of the reasons you are not<br>> getting the data is probably because your pressure level is not set<br>> correctly. Nonetheless, I should think the while loop should work just<br>> fine, but you need to be careful how you set your limits. If you want to go<br>> from low levels to upper levels, your while statement should look something<br>> like<br>> <br>> i = 1 ; zlev = 1000<br>> while (zlev > 100)<br>> ...<br>> ...<br>> ...<br>> zlev = zlev - dp*i (where dp is your pressure change between<br>> levels...hopefully it is constant. If it isn't there's another way to do<br>> this)<br>> i = i + 1<br>> endwhile<br>> <br>> Try that and see if it works.<br>> <br>> Jeff Duda<br>> <br>> On Mon, Oct 25, 2010 at 3:53 AM, Geeta Geeta <geeta124@hotmail.com> wrote:<br>> <br>> > Dear All,<br>> > I am trying to extract the SPECIFIC HUMIDITY for all the levels in vertical<br>> > (19 Levels from 1000 to 100hPa) from the Model.<br>> > So I am Using tvrh2q function of Grads. Following is the ctl file.<br>> ><br>> > 'open test.ctl'<br>> > *'set lev 1000 100'<br>> > 'set t 1'<br>> > * while (lev< 100)<br>> > *zlev=1000<br>> > 'set gxout fwrite'<br>> > 'set fwrite test.dat'<br>> > zlev =1000<br>> > 'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'<br>> > 'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))'<br>> > ************<br>> > zlev =950<br>> > 'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'<br>> > 'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))'<br>> > *'disable fwrite'<br>> > *****************<br>> > zlev =900<br>> > 'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'<br>> > 'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))'<br>> > 'disable fwrite'<br>> > ***********************<br>> ><br>> > * endwhile<br>> ><br>> > I am Unable to Extract the data Using the While Statement in Grads, so I am<br>> > defining the Zlev separately .<br>> ><br>> > Kindly Suggest.<br>> > geeta<br>> ><br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org<br>> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> ><br>> ><br>> <br>> <br>> -- <br>> Jeff Duda<br>> Iowa State University<br>> Meteorology Graduate Student<br>> 3134 Agronomy Hall<br>> www.meteor.iastate.edu/~jdduda<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/d800f45c/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 4<br>> Date: Mon, 25 Oct 2010 11:04:04 -0400<br>> From: Stephen McMillan <smcmillan@planalytics.com><br>> Subject: [gradsusr] Creating "Lines" Shapefiles<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID:<br>> <AANLkTinzsXXMzdhj3XaFmpEEtjeRiJ9-JLETbRYV_zQq@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> I think it's a great addition to have both the draw- and create-shapefiles<br>> capabilities in GrADS now, and have successfully done both from the latest<br>> release.<br>> <br>> As I understand the online documentation, shapefiles can be created using<br>> only the "display" command, and that "the shape types currently supported<br>> are points and lines...." with lines being contour lines. Are there any<br>> plans to support shapefile creation using the "draw" command, such as for<br>> creating straight line segments or polygons, as it does for drawing<br>> shapefiles?<br>> <br>> Perhaps the developers or another user knows of a way to write shapefiles<br>> derived from drawn (not displayed) map features. The only way I can think<br>> of is to somehow incorporate the desired line/polygon features within a<br>> gridded file, perhaps similar to the way "lat" and "lon" are done, except<br>> crossing lat/lon lines.<br>> <br>> I am using the OpenGrADS version 2.0.a9.oga.1 on a Windows XP.<br>> <br>> Thanks<br>> Stephen McMillan<br>> <br>> ***************************************************<br>> The information contained in this e-mail message <br>> is intended only for the use of the recipient(s) <br>> named above and may contain information that is <br>> privileged, confidential, and/or proprietary. <br>> If you are not the intended recipient, you may not<br>> review, copy or distribute this message. If you have<br>> received this communication in error, please notify <br>> the sender immediately by e-mail, and delete the original message.<br>> ***************************************************<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/847cf9c4/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 5<br>> Date: Mon, 25 Oct 2010 21:30:10 +0530<br>> From: SRAVANI A <srisravani.a@gmail.com><br>> Subject: [gradsusr] (no subject)<br>> To: gradsusr@gradsusr.org<br>> Message-ID:<br>> <AANLkTimzkFgLdkb=BZK_V+x197YzMXdF4tv_7YpD2mt2@mail.gmail.com><br>> Content-Type: text/plain; charset="iso-8859-1"<br>> <br>> HI<br>> grads user i had one problem i m using .nc wind files to find the psi and<br>> khi . now i want to partial diffrention wth respective to time so plese<br>> kindly help me<br>> <br>> <br>> -- <br>> A.Sravani<br>> Research Scholar<br>> Department of Meteorology & Oceanography<br>> Andhra university<br>> Visakhapatnam<br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/3d6e3b8c/attachment-0001.html <br>> <br>> ------------------------------<br>> <br>> Message: 6<br>> Date: Mon, 25 Oct 2010 12:43:51 -0400<br>> From: Charles Seman <Charles.Seman@noaa.gov><br>> Subject: Re: [gradsusr] Max Value Coordinates<br>> To: GrADS Users Forum <gradsusr@gradsusr.org><br>> Message-ID: <4CC5B3C7.9080106@noaa.gov><br>> Content-Type: text/plain; charset="windows-1252"<br>> <br>> Andrew and Jeff and GrADS Users,<br>> <br>> Kun-Hsuan Chou provided the reference code for this technique (which can <br>> be found in the attached file hurricane_tracking.txt, taken from an <br>> email exchange posted to the GrADS Listserv by Kun-Hsuan Chou and Arturo <br>> Caracas Uribe in Oct, 2004).<br>> <br>> Chuck<br>> <br>> Jeffrey Duda wrote:<br>> > Andrew,<br>> > I had posted a request on this forum about this exact question in <br>> > August. I believe it was Charles Seaman who provided the answer. The <br>> > meat of the code is<br>> ><br>> > 'd maxloc(max('ans',x='x1',x='x2'),y='y1',y='y2')'<br>> > yline = sublin(result,2)<br>> > ymax = subwrd(yline,4)<br>> > 'd maxloc(max('ans',y='y1',y='y2'),x='x1',x='x2')'<br>> > xline = sublin(result,2)<br>> > xmax = subwrd(xline,4)<br>> ><br>> ><br>> > where "ans" refers to the field for which you want to find the <br>> > max/min, and x1, x2 refer to x-grid values or longitude points and y1, <br>> > y2 refer to y-grid values or latitude points. xmax, ymax are the <br>> > coordinates of the maximum value. To find the minimum, just replace <br>> > maxloc and max with minloc and min.<br>> ><br>> > Jeff Duda<br>> ><br>> > On Fri, Oct 22, 2010 at 7:43 PM, Andrew Revering <andy@f5data.com <br>> > <mailto:andy@f5data.com>> wrote:<br>> ><br>> > I'd like to have the coordinates of our 'max value'? I know about<br>> > maxloc(), and tried various examples I found in forums,<br>> > documentation and email conversations, but I can't seem to get it<br>> > to perform as I expect.<br>> ><br>> > <br>> ><br>> > When it supposedly works, I just see some white lines on the<br>> > bottom and right of my image (not where the max value is)? is<br>> > there a way to actually get lat/lon coordinates of the max value?<br>> ><br>> > <br>> ><br>> > I should note that in my situation there will be the teeniest dot<br>> > of plotted value, the rest of the grid is blank. I don't know if<br>> > that matters in these functions.<br>> ><br>> > <br>> ><br>> > Andrew Revering<br>> ><br>> > Convective Development, Inc.<br>> ><br>> > http://www.convectivedevelopment.com/<br>> ><br>> > <br>> ><br>> > <br>> ><br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org <mailto:gradsusr@gradsusr.org><br>> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> ><br>> ><br>> ><br>> ><br>> > -- <br>> > Jeff Duda<br>> > Iowa State University<br>> > Meteorology Graduate Student<br>> > 3134 Agronomy Hall<br>> > www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda><br>> > ------------------------------------------------------------------------<br>> ><br>> > _______________________________________________<br>> > gradsusr mailing list<br>> > gradsusr@gradsusr.org<br>> > http://gradsusr.org/mailman/listinfo/gradsusr<br>> > <br>> <br>> -- <br>> <br>> Please note that Charles.Seman@noaa.gov should be considered my NOAA<br>> email address, not cjs@gfdl.noaa.gov.<br>> <br>> ********************************************************************<br>> Charles Seman Charles.Seman@noaa.gov<br>> U.S. Department of Commerce / NOAA / OAR<br>> Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547<br>> 201 Forrestal Road fax: (609) 987-5063<br>> Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/<br>> ********************************************************************<br>> <br>> "The contents of this message are mine personally and do not reflect any<br>> official or unofficial position of the United States Federal Government,<br>> the United States Department of Commerce, or NOAA."<br>> <br>> -------------- next part --------------<br>> An embedded and charset-unspecified text was scrubbed...<br>> Name: hurricane_tracking.txt<br>> Url: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/b5abbf46/attachment-0001.txt <br>> <br>> ------------------------------<br>> <br>> Message: 7<br>> Date: Mon, 25 Oct 2010 20:06:58 +0200<br>> From: "Javier G. Corripio" <jgcorripio@mac.com><br>> Subject: [gradsusr] (no subject)<br>> To: gradsusr@gradsusr.org<br>> Message-ID: <FA9D704E-063E-4948-94B7-CC2A550DB315@mac.com><br>> Content-Type: text/plain; charset="utf-8"<br>> <br>> Dear Grads Users,<br>> <br>> I wonder if there is a simple way to add labels to KML contour lines, <br>> the output of:<br>> <br>> set kml -ln fnameroot<br>> <br>> Thanks a lot.<br>> <br>> Javier<br>> <br>> P.S. Sorry if it has already been asked, I couldn't find an list <br>> search function, so I have downloaded the monthly text files and grep <br>> all "kml"s but I couldn't find an answer<br>> <br>> -- <br>> **************************************<br>> Javier G. Corripio<br>> email: jgcorripio@mac.com<br>> url: http://www.meteoexploration.com<br>> <br>> <br>> ?<br>> <br>> -------------- next part --------------<br>> An HTML attachment was scrubbed...<br>> URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/6678ffa7/attachment.html <br>> -------------- next part --------------<br>> A non-text attachment was scrubbed...<br>> Name: tmk.png<br>> Type: image/png<br>> Size: 3534 bytes<br>> Desc: not available<br>> Url : http://gradsusr.org/pipermail/gradsusr/attachments/20101025/6678ffa7/attachment.png <br>> <br>> ------------------------------<br>> <br>> _______________________________________________<br>> gradsusr mailing list<br>> gradsusr@gradsusr.org<br>> http://gradsusr.org/mailman/listinfo/gradsusr<br>> <br>> <br>> End of gradsusr Digest, Vol 8, Issue 42<br>> ***************************************<br> </body>
</html>