<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>Dear Andrew,<BR>&nbsp;<BR>Thank you for your reply. Finally I could get the result and the script is as follows. I just changed r1= subwrd(result, 4) to r1= subwrd(result,1 4).<BR>&nbsp;<BR>&nbsp;<BR><div>'reinit'<br>'sdfopen H:/prate.sfc.mon.mean.nc'<br>'set lat 6.66'<br>'set lon 99.375'<br>'d ave(prate,t=517,t=636)'<br>r1= subwrd(result,1 4)<br>say r1<br>'close 1'</div><div>&nbsp;</div><div>Regards,<BR>Thareh<BR><br><br>&gt; Message: 1<br>&gt; Date: Sat, 4 Oct 2014 23:21:07 +0200<br>&gt; From: Andrew Friedman &lt;andfried@gmail.com&gt;<br>&gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; Message-ID: &lt;630FFD43-E9AF-487C-8833-3DDA857BB752@gmail.com&gt;<br>&gt; Content-Type: text/plain; charset=windows-1256<br>&gt; <br>&gt; Hi Thareh,<br>&gt; <br>&gt; To clarify:  What is the result of this line (cutting out the ?subwrd?? line) for different lat and lon values?<br>&gt; &gt; 'd ave(prate,t=517,t=636)?<br>&gt; <br>&gt; Does it give the correct answer then?<br>&gt; <br>&gt; I?m not very familiar with the ?subwrd? command, but perhaps that could be causing your issue.<br>&gt; <br>&gt; Andrew<br>&gt; <br>&gt; On Oct 4, 2014, at 7:37 PM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; <br>&gt; &gt; Dear Andrew,<br>&gt; &gt; Thanks for your help. Unfortunately I couldn't do that. <br>&gt; &gt; I wrote this script to find an average of 10-year sst data from ncep data in every special longitude and latitude.<br>&gt; &gt;  <br>&gt; &gt; 'reinit'<br>&gt; &gt; 'sdfopen H:/prate.sfc.mon.mean.nc'<br>&gt; &gt; 'set lat 6.66'<br>&gt; &gt; 'set lon 99.375'<br>&gt; &gt; 'd ave(prate,t=517,t=636)'<br>&gt; &gt; r1= subwrd(result,4)<br>&gt; &gt; say r1<br>&gt; &gt; 'close 1'<br>&gt; &gt;  <br>&gt; &gt; But when I change the lon and lat values I always get the same value "3," which is a wrong value.<br>&gt; &gt;  <br>&gt; &gt; Could you help me to write a correct script. <br>&gt; &gt; regards,<br>&gt; &gt; Thareh<br>&gt; &gt; <br>&gt; &gt;  <br>&gt; &gt; &gt; From: gradsusr-request@gradsusr.org<br>&gt; &gt; &gt; Subject: gradsusr Digest, Vol 56, Issue 6<br>&gt; &gt; &gt; To: gradsusr@gradsusr.org<br>&gt; &gt; &gt; Date: Thu, 2 Oct 2014 12:00:02 -0400<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; Send gradsusr mailing list submissions to<br>&gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; To subscribe or unsubscribe via the World Wide Web, visit<br>&gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; or, via email, send a message with subject or body 'help' to<br>&gt; &gt; &gt; gradsusr-request@gradsusr.org<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; You can reach the person managing the list at<br>&gt; &gt; &gt; gradsusr-owner@gradsusr.org<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; When replying, please edit your Subject line so it is more specific<br>&gt; &gt; &gt; than "Re: Contents of gradsusr digest..."<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; <br>&gt; &gt; &gt; Today's Topics:<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; 1. Re: Average of a parameter (Andrew Friedman)<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; <br>&gt; &gt; &gt; ----------------------------------------------------------------------<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; Message: 1<br>&gt; &gt; &gt; Date: Thu, 2 Oct 2014 17:53:50 +0200<br>&gt; &gt; &gt; From: Andrew Friedman &lt;andfried@gmail.com&gt;<br>&gt; &gt; &gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; Message-ID: &lt;DF6903BD-691F-4969-A2B9-EB71CAAC8608@gmail.com&gt;<br>&gt; &gt; &gt; Content-Type: text/plain; charset=windows-1256<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; Hi Thareh,<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; Check out the documentation on ASCII data in GrADS to save the averages in a table as you?re looking for.<br>&gt; &gt; &gt; http://www.iges.org/grads/gadoc/ascii.html<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; A related option is to use the OpenGrADS recipe fprintf.gs:<br>&gt; &gt; &gt; http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; Andrew<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; On Oct 2, 2014, at 4:58 PM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; Hi,<br>&gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; Unfortunately, It didn't work for my aim. In fact, I would like to have the average over the time but in every grid point (according to different lon and lat) separately. To show what kind of average I need is like a matrix as follows. each number is the averaged data during 10 years (t= 517 to t=643), based on Ncep data set. Can anyone help?<br>&gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; 99.375        101.25        103.125        105        106.875        108.75        110.625        112.5        114.375        116.25        118.125<br>&gt; &gt; &gt; &gt; 6.66        6.38E-05        1.02E-04        5.36E-05        5.58E-05        6.27E-05        7.12E-05        8.22E-05        8.61E-05        6.89E-05        4.79E-05        4.94E-05<br>&gt; &gt; &gt; &gt; 4.76        5.74E-05        1.11E-04        1.03E-04        6.50E-05        6.75E-05        7.09E-05        8.32E-05        7.41E-05        3.92E-05        1.06E-04        1.07E-04<br>&gt; &gt; &gt; &gt; 2.85        1.30E-04        5.39E-05        1.22E-04        6.92E-05        6.71E-05        6.29E-05        6.39E-05        1.13E-04        9.78E-05        1.06E-04        4.05E-05<br>&gt; &gt; &gt; &gt; 0.95        1.30E-04        1.31E-04        1.36E-04        6.95E-05        6.53E-05        5.40E-05        1.23E-04        1.07E-04        9.03E-05        1.12E-04        1.19E-04<br>&gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; I have 10*12=120 of tables like this for 10 years with 12 month. I would like to write a script that read the data from a file in each grid point in every month during 10 years and then calculate the average, and show the average in a text file.<br>&gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; Regards, <br>&gt; &gt; &gt; &gt; Thareh<br>&gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; From: gradsusr-request@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; Subject: gradsusr Digest, Vol 56, Issue 4<br>&gt; &gt; &gt; &gt; &gt; To: gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 22:45:15 -0400<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Send gradsusr mailing list submissions to<br>&gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; To subscribe or unsubscribe via the World Wide Web, visit<br>&gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; or, via email, send a message with subject or body 'help' to<br>&gt; &gt; &gt; &gt; &gt; gradsusr-request@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; You can reach the person managing the list at<br>&gt; &gt; &gt; &gt; &gt; gradsusr-owner@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; When replying, please edit your Subject line so it is more specific<br>&gt; &gt; &gt; &gt; &gt; than "Re: Contents of gradsusr digest..."<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Today's Topics:<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; 1. Re: Average of a parameter (Andrew Friedman)<br>&gt; &gt; &gt; &gt; &gt; 2. Re: Average of a parameter (Mingxuan Chen)<br>&gt; &gt; &gt; &gt; &gt; 3. Re: Average of a parameter (Joao Nobre)<br>&gt; &gt; &gt; &gt; &gt; 4. How to plot X-Y axis in log scale (Jian-Jian Wang)<br>&gt; &gt; &gt; &gt; &gt; 5. Correlation (Jian-Jian Wang)<br>&gt; &gt; &gt; &gt; &gt; 6. Correlation (Jian-Jian Wang)<br>&gt; &gt; &gt; &gt; &gt; 7. Moist Potential Vorticity (PHILBERT LUHUNGA)<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ----------------------------------------------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 1<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 18:21:42 +0200<br>&gt; &gt; &gt; &gt; &gt; From: Andrew Friedman &lt;andfried@gmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID: &lt;E643D2EE-87B0-4CE5-B290-7C9002C06F4F@gmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset=windows-1256<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Hi Thareh,<br>&gt; &gt; &gt; &gt; &gt; As written, your script averages over space, not time.<br>&gt; &gt; &gt; &gt; &gt; If you are trying to average over time, you should use an average of the form<br>&gt; &gt; &gt; &gt; &gt; ?ave(prate,t=517,t=643)'<br>&gt; &gt; &gt; &gt; &gt; Also, if you?re trying to calculate a spatial average, I suggest using ?aave?.<br>&gt; &gt; &gt; &gt; &gt; Andrew<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; On Oct 1, 2014, at 5:08 PM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; Dear Andrew Friedman,<br>&gt; &gt; &gt; &gt; &gt; &gt; Thank you for your help.<br>&gt; &gt; &gt; &gt; &gt; &gt; In fact I would like to calculate average of sst during 10 years in special area for example lon=100 120 and lat= 0 7. Then I want to put the result in a text file.<br>&gt; &gt; &gt; &gt; &gt; &gt; I wrote this script:<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; 'reinit'<br>&gt; &gt; &gt; &gt; &gt; &gt; 'sdfopen H:/prate.sfc.mon.mean.nc'<br>&gt; &gt; &gt; &gt; &gt; &gt; 'set t 517 643'<br>&gt; &gt; &gt; &gt; &gt; &gt; 'd ave(ave(prate,lon=100,lon=120),lat=0,lat=7)'<br>&gt; &gt; &gt; &gt; &gt; &gt; r1= subwrd(result,4)<br>&gt; &gt; &gt; &gt; &gt; &gt; say r1<br>&gt; &gt; &gt; &gt; &gt; &gt; 'close 1'<br>&gt; &gt; &gt; &gt; &gt; &gt; vars= r1<br>&gt; &gt; &gt; &gt; &gt; &gt; res=write('filename3.txt', vars ,append)<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; but unfortunately I just got one number but I must see too many data. Is it possible to have average of SST during 10 years in special longitude and latitude. I mean I would like to have sst1 ( lon1&amp;lat1), sst2 (lon2&amp;lat2), ..... which averaged during 10 years.<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; Regards,<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; Thareh<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Message: 2<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 11:10:44 +0200<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; From: Andrew Friedman &lt;andfried@gmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Message-ID: &lt;85BC0AE9-7894-4048-9A70-9126D5993A61@gmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset=windows-1256<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Hi Thareh,<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Check out the ?ave? command: http://www.iges.org/grads/gadoc/gradfuncave.html<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Best, Andrew<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; On Oct 1, 2014, at 9:16 AM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; Dear All,<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; Could anyone tell me how I can take the average of a parameter during 10 years? For example I would like to calculate 10-year average of SST in Grads.<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; Best regards,<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; Thareh <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 2<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 12:07:45 -0600<br>&gt; &gt; &gt; &gt; &gt; From: Mingxuan Chen &lt;chen@atmos.colostate.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID: &lt;7F53D3A4-970D-4C0E-ACB0-4B1919866950@atmos.colostate.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset=windows-1256<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; See you at Parking lot.<br>&gt; &gt; &gt; &gt; &gt; On Oct 1, 2014, at 10:21 AM, Andrew Friedman &lt;andfried@gmail.com&gt; wrote:<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; Hi Thareh,<br>&gt; &gt; &gt; &gt; &gt; &gt; As written, your script averages over space, not time.<br>&gt; &gt; &gt; &gt; &gt; &gt; If you are trying to average over time, you should use an average of the form<br>&gt; &gt; &gt; &gt; &gt; &gt; ?ave(prate,t=517,t=643)'<br>&gt; &gt; &gt; &gt; &gt; &gt; Also, if you?re trying to calculate a spatial average, I suggest using ?aave?.<br>&gt; &gt; &gt; &gt; &gt; &gt; Andrew<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; On Oct 1, 2014, at 5:08 PM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; Dear Andrew Friedman,<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; Thank you for your help.<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; In fact I would like to calculate average of sst during 10 years in special area for example lon=100 120 and lat= 0 7. Then I want to put the result in a text file.<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; I wrote this script:<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; 'reinit'<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; 'sdfopen H:/prate.sfc.mon.mean.nc'<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; 'set t 517 643'<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; 'd ave(ave(prate,lon=100,lon=120),lat=0,lat=7)'<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; r1= subwrd(result,4)<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; say r1<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; 'close 1'<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; vars= r1<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; res=write('filename3.txt', vars ,append)<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; but unfortunately I just got one number but I must see too many data. Is it possible to have average of SST during 10 years in special longitude and latitude. I mean I would like to have sst1 ( lon1&amp;lat1), sst2 (lon2&amp;lat2), ..... which averaged during 10 years.<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; Regards,<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; Thareh<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Message: 2<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Date: Wed, 1 Oct 2014 11:10:44 +0200<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; From: Andrew Friedman &lt;andfried@gmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Message-ID: &lt;85BC0AE9-7894-4048-9A70-9126D5993A61@gmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Content-Type: text/plain; charset=windows-1256<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Hi Thareh,<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Check out the ?ave? command: http://www.iges.org/grads/gadoc/gradfuncave.html<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; Best, Andrew<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; On Oct 1, 2014, at 9:16 AM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; Dear All,<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; Could anyone tell me how I can take the average of a parameter during 10 years? For example I would like to calculate 10-year average of SST in Grads.<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; Best regards,<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; Thareh <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt;&gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; <br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt;&gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 3<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 23:04:29 +0300<br>&gt; &gt; &gt; &gt; &gt; From: Joao Nobre &lt;j.pnobre@hotmail.com&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID: &lt;SNT148-W402759BEBC5A3BEBD441158FB80@phx.gbl&gt;<br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset="windows-1252"<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; First, you should go to this site and next you can make a plot on GrADS, using the normal Scripts or Codes, if you prefer!Name of site http://www.esrl.noaa.gov/psd/cgi-bin/db_search/SearchMenus.pl<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; From: andfried@gmail.com<br>&gt; &gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 11:10:44 +0200<br>&gt; &gt; &gt; &gt; &gt; &gt; To: gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt; Subject: Re: [gradsusr] Average of a parameter<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; Hi Thareh,<br>&gt; &gt; &gt; &gt; &gt; &gt; Check out the ?ave? command: http://www.iges.org/grads/gadoc/gradfuncave.html<br>&gt; &gt; &gt; &gt; &gt; &gt; Best, Andrew<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; On Oct 1, 2014, at 9:16 AM, T. haghroosta &lt;haghroosta@hotmail.com&gt; wrote:<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Dear All,<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Could anyone tell me how I can take the average of a parameter during 10 years? For example I would like to calculate 10-year average of SST in Grads.<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Best regards,<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; Thareh <br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; -------------- next part --------------<br>&gt; &gt; &gt; &gt; &gt; An HTML attachment was scrubbed...<br>&gt; &gt; &gt; &gt; &gt; URL: http://gradsusr.org/pipermail/gradsusr/attachments/20141001/c000739a/attachment-0001.html <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 4<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 21:28:19 +0000<br>&gt; &gt; &gt; &gt; &gt; From: Jian-Jian Wang &lt;jjwang@umd.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: [gradsusr] How to plot X-Y axis in log scale<br>&gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID: &lt;C270AC74-A859-4A1D-9653-F0114FD66246@umd.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset="us-ascii"<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Hi,<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; I would like to make a scatterplot in log scale. But seems Grads can only do 1-D at Z dimension?<br>&gt; &gt; &gt; &gt; &gt; Any suggestions?<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Thanks!<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Jay<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 5<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 21:32:21 +0000<br>&gt; &gt; &gt; &gt; &gt; From: Jian-Jian Wang &lt;jjwang@umd.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: [gradsusr] Correlation<br>&gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID: &lt;90F31DC7-BE47-4DF8-936C-70F475C8BBDA@umd.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset="us-ascii"<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Hi,<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; When make a scatter plot like:<br>&gt; &gt; &gt; &gt; &gt; 'set gxout scatter'<br>&gt; &gt; &gt; &gt; &gt; 'd var1 var2'<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; It has no problem. Then when I use:<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; d scorr(var1, var2, lon=-125, lon=-65, lat=20, lat=50)<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; It doesn't show anything. what did I do wrong? <br>&gt; &gt; &gt; &gt; &gt; Thanks.<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Jay<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 6<br>&gt; &gt; &gt; &gt; &gt; Date: Wed, 1 Oct 2014 21:46:05 +0000<br>&gt; &gt; &gt; &gt; &gt; From: Jian-Jian Wang &lt;jjwang@umd.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: [gradsusr] Correlation<br>&gt; &gt; &gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID: &lt;68443CE6-7209-4BF1-BC68-F972828F2030@umd.edu&gt;<br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset="us-ascii"<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; 'reinit'<br>&gt; &gt; &gt; &gt; &gt; 'open npp_swath_test.ctl'<br>&gt; &gt; &gt; &gt; &gt; 'open /user/parkin/validation/STAR/SWATH/NMQ/nmq_test.ctl'<br>&gt; &gt; &gt; &gt; &gt; 'd scorr(nmq.2, atm.1, lon=-125, lon=-65, lat=20, lat=50)'<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; When I run the script in Grads, it won't give me any number.<br>&gt; &gt; &gt; &gt; &gt; But if I do it line by line like:<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; open npp_swath_test.ctl<br>&gt; &gt; &gt; &gt; &gt; open /user/parkin/validation/STAR/SWATH/NMQ/nmq_test.ctl<br>&gt; &gt; &gt; &gt; &gt; d scorr(nmq.2, atm.1, lon=-125, lon=-65, lat=20, lat=50)<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; it shows: <br>&gt; &gt; &gt; &gt; &gt; Result value = 0.426107<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; why didn't it show the Result value when running the script?<br>&gt; &gt; &gt; &gt; &gt; Thank you in advance.<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Jay<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Message: 7<br>&gt; &gt; &gt; &gt; &gt; Date: Thu, 2 Oct 2014 02:42:33 +0000 (UTC)<br>&gt; &gt; &gt; &gt; &gt; From: PHILBERT LUHUNGA &lt;philuhunga@yahoo.com&gt;<br>&gt; &gt; &gt; &gt; &gt; Subject: [gradsusr] Moist Potential Vorticity<br>&gt; &gt; &gt; &gt; &gt; To: "gradsusr@gradsusr.org" &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; &gt; &gt; Message-ID:<br>&gt; &gt; &gt; &gt; &gt; &lt;1916954210.38931.1412217753822.JavaMail.yahoo@jws10646.mail.bf1.yahoo.com&gt;<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Content-Type: text/plain; charset="utf-8"<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; Hi I am writing a script to compute moist potential vorticity using entropy temperature Theta_s, and virtual potential temperature Theta_v. Please help me to check the script if it is coded right. When I drive it, it gives me results but am not sure the way I have coded if it is right. could anyone cross check it?*our script to compute MPV700hPa.<br>&gt; &gt; &gt; &gt; &gt; *###################Day1################################<br>&gt; &gt; &gt; &gt; &gt; *'open $DATA_DIR/WRF_9KM_Nonhy_Topo${RunDay}.ctl'<br>&gt; &gt; &gt; &gt; &gt; 'open 201112200000_arw_wrfout_d01.ctl'<br>&gt; &gt; &gt; &gt; &gt; ?<br>&gt; &gt; &gt; &gt; &gt; 'q file''q dims''pi=3.14159'<br>&gt; &gt; &gt; &gt; &gt; 'dtr=pi/180'<br>&gt; &gt; &gt; &gt; &gt; 'a=6.37122*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; 'omega=7.2921*pow(10,-5)'<br>&gt; &gt; &gt; &gt; &gt; 'g=9.81'<br>&gt; &gt; &gt; &gt; &gt; 'R=287'<br>&gt; &gt; &gt; &gt; &gt; 'Cp=1004.5'<br>&gt; &gt; &gt; &gt; &gt; 'L=2.453*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; *'L=2.453e6'<br>&gt; &gt; &gt; &gt; &gt; 'Rv=461'<br>&gt; &gt; &gt; &gt; &gt; * tk??? = Temp in Kelvin<br>&gt; &gt; &gt; &gt; &gt; * mixr? = water vapour Mixing ratio<br>&gt; &gt; &gt; &gt; &gt; * u???? = U-wind in m/s<br>&gt; &gt; &gt; &gt; &gt; * v???? = V-wind in m/s<br>&gt; &gt; &gt; &gt; &gt; * define ug=UGRDprs?? U-wind in m/s in pressure level<br>&gt; &gt; &gt; &gt; &gt; * define vg=VGRDprs? V-wind in m/s in pressure level<br>&gt; &gt; &gt; &gt; &gt; *'define ua=U-ug'<br>&gt; &gt; &gt; &gt; &gt; *'define va=V-vg'*PV @700hPa<br>&gt; &gt; &gt; &gt; &gt; 'set lev 800'<br>&gt; &gt; &gt; &gt; &gt; 'define f=2*omega*sin(lat*dtr)'<br>&gt; &gt; &gt; &gt; &gt; 'define p1=lev*100'<br>&gt; &gt; &gt; &gt; &gt; 'define tk1=TMPprs'<br>&gt; &gt; &gt; &gt; &gt; 'define u1=UGRDprs'<br>&gt; &gt; &gt; &gt; &gt; 'define v1=VGRDprs'<br>&gt; &gt; &gt; &gt; &gt; 'define rh1=rhprs'<br>&gt; &gt; &gt; &gt; &gt; 'define thetap1=tk1*pow((100000/p1),0.286)'<br>&gt; &gt; &gt; &gt; &gt; 'define es=6.1173*100*exp((L/Rv)*(1/273.16 - 1/tk1))'<br>&gt; &gt; &gt; &gt; &gt; 'define ep1=rh1*es/100'<br>&gt; &gt; &gt; &gt; &gt; 'WVmixr1=0.622*ep1/(p1)'<br>&gt; &gt; &gt; &gt; &gt; 'define T=tk1'<br>&gt; &gt; &gt; &gt; &gt; 'define VirT1=T*(1+0.61*WVmixr1)'<br>&gt; &gt; &gt; &gt; &gt; 'thetav1=VirT1*pow((100000/p1),0.2857))'<br>&gt; &gt; &gt; &gt; &gt; 'thetas1=thetap1*exp(5.87*WVmixr1)''set lev 600'<br>&gt; &gt; &gt; &gt; &gt; 'define p2=lev*100'<br>&gt; &gt; &gt; &gt; &gt; 'define tk2=TMPprs'<br>&gt; &gt; &gt; &gt; &gt; 'define u2=UGRDprs'<br>&gt; &gt; &gt; &gt; &gt; 'define v2=VGRDprs'<br>&gt; &gt; &gt; &gt; &gt; 'define rh2=rhprs'<br>&gt; &gt; &gt; &gt; &gt; 'define thetap2=tk2*pow((100000/p2),0.286)'<br>&gt; &gt; &gt; &gt; &gt; 'define es=6.1173*100*exp((L/Rv)*(1/273.16 - 1/tk2))'<br>&gt; &gt; &gt; &gt; &gt; 'define ep2=rh2*es/100'<br>&gt; &gt; &gt; &gt; &gt; 'WVmixr2=0.622*ep2/(p2)'<br>&gt; &gt; &gt; &gt; &gt; 'define T=tk2'<br>&gt; &gt; &gt; &gt; &gt; 'define VirT2=T*(1+0.61*WVmixr2)'<br>&gt; &gt; &gt; &gt; &gt; 'thetav2=VirT2*pow((100000/p2),0.2857))'<br>&gt; &gt; &gt; &gt; &gt; 'thetas2=thetap2*exp(5.87*WVmixr2)'<br>&gt; &gt; &gt; &gt; &gt; ********************PV at 700hPa level<br>&gt; &gt; &gt; &gt; &gt; 'set lev 700'<br>&gt; &gt; &gt; &gt; &gt; 'dp=p2-p1'<br>&gt; &gt; &gt; &gt; &gt; 'dudp=(u2-u1)/dp'<br>&gt; &gt; &gt; &gt; &gt; 'dvdp=(v2-v1)/dp'<br>&gt; &gt; &gt; &gt; &gt; 'define p0=lev*100'<br>&gt; &gt; &gt; &gt; &gt; 'define tk0=TMPprs'<br>&gt; &gt; &gt; &gt; &gt; 'define u0=UGRDprs'<br>&gt; &gt; &gt; &gt; &gt; 'define v0=VGRDprs'<br>&gt; &gt; &gt; &gt; &gt; 'vort0=hcurl(u0,v0)'<br>&gt; &gt; &gt; &gt; &gt; 'define Totvort=vort0+f''define alpha0=R*tk0/p0'<br>&gt; &gt; &gt; &gt; &gt; 'define rh0=rhprs'<br>&gt; &gt; &gt; &gt; &gt; 'theta0=tk0*pow((100000/p0),0.2857))'<br>&gt; &gt; &gt; &gt; &gt; 'define es=6.1173*100*exp((L/Rv)*(1/273.16 - 1/tk0))'<br>&gt; &gt; &gt; &gt; &gt; 'de<br>&gt; &gt; &gt; &gt; &gt; fine ep0=rh0*es/100'<br>&gt; &gt; &gt; &gt; &gt; 'WVmixr0=0.622*ep0/(p0)'<br>&gt; &gt; &gt; &gt; &gt; 'define T=tk0'<br>&gt; &gt; &gt; &gt; &gt; 'define VirT0=T*(1+0.61*WVmixr0)'<br>&gt; &gt; &gt; &gt; &gt; 'thetav0=VirT0*pow((100000/p0),0.2857))'<br>&gt; &gt; &gt; &gt; &gt; 'thetas0=theta0*exp(5.87*WVmixr0)'<br>&gt; &gt; &gt; &gt; &gt; 'dy=cdiff(lat,y)*dtr*a'<br>&gt; &gt; &gt; &gt; &gt; 'dx=cdiff(lon,x)*dtr*a*cos(lat*dtr)'<br>&gt; &gt; &gt; &gt; &gt; 'define dthetas0dx=cdiff(thetas0,x)/dx'<br>&gt; &gt; &gt; &gt; &gt; 'define dthetas0dy=cdiff(thetas0,y)/dy'<br>&gt; &gt; &gt; &gt; &gt; 'define dthetas0dp=(thetas2-thetas1)/dp'<br>&gt; &gt; &gt; &gt; &gt; 'define Totdthedtas=dthetas0dx+dthetas0dy+dthetas0dp''define dthetav0dx=cdiff(thetav0,x)/dx'<br>&gt; &gt; &gt; &gt; &gt; 'define dthetav0dy=cdiff(thetav0,y)/dy'<br>&gt; &gt; &gt; &gt; &gt; 'define dthetav0dp=(thetav2-thetav1)/dp'<br>&gt; &gt; &gt; &gt; &gt; 'define Totdthedtav=dthetav0dx+dthetav0dy+dthetav0dp''term11= -g*Totvort*dthetav0dp*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; 'term22=-g*dudp*dthetav0dy*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; 'term33=g*dVdp*dthetav0dx*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; 'PV7002=term1+term2+term3''term1= -g*Totvort*dthetas0dp*pow(10,6)''term2=-g*dudp*dthetas0dy*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; 'term3=g*dVdp*dthetas0dx*pow(10,6)'<br>&gt; &gt; &gt; &gt; &gt; 'PV700=term1+term2+term3'*'term1=Totvort*dthetas0dp*e6)'<br>&gt; &gt; &gt; &gt; &gt; *'term2=dudp*dthetas0dy*e6'<br>&gt; &gt; &gt; &gt; &gt; *'term3=-dVdp*dthetas0dx*e6'<br>&gt; &gt; &gt; &gt; &gt; 'PV700=term1+term2+term3''set t 7'<br>&gt; &gt; &gt; &gt; &gt; 'set gxout shaded'<br>&gt; &gt; &gt; &gt; &gt; 'set mpdset hires'<br>&gt; &gt; &gt; &gt; &gt; 'set display color white'<br>&gt; &gt; &gt; &gt; &gt; 'set vpage 1.0 8.5 2.0 8.5'<br>&gt; &gt; &gt; &gt; &gt; 'c'<br>&gt; &gt; &gt; &gt; &gt; *'d mp'd<br>&gt; &gt; &gt; &gt; &gt; 'd theta_s'<br>&gt; &gt; &gt; &gt; &gt; *'exec colors'<br>&gt; &gt; &gt; &gt; &gt; *'run colorset.gs'<br>&gt; &gt; &gt; &gt; &gt; *'set ylint 12'<br>&gt; &gt; &gt; &gt; &gt; *'set xlint 12'<br>&gt; &gt; &gt; &gt; &gt; 'cbarhov2.gs'<br>&gt; &gt; &gt; &gt; &gt; 'draw title 700hPa-Moist Potential Vorticity ''printim 700_mv.png'<br>&gt; &gt; &gt; &gt; &gt; 'set grads off'<br>&gt; &gt; &gt; &gt; &gt; *'quit'*'define MV=TotVort*ept(tmpprs,rhprs.prs)'<br>&gt; &gt; &gt; &gt; &gt; ??AS EVER<br>&gt; &gt; &gt; &gt; &gt; LUHUNGA<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; --------------------<br>&gt; &gt; &gt; &gt; &gt; Philbert Modest Luhunga<br>&gt; &gt; &gt; &gt; &gt; University of Pretoria?<br>&gt; &gt; &gt; &gt; &gt; Department of Geography,Geoinformatics and Meteorology<br>&gt; &gt; &gt; &gt; &gt; Private Bag X20 Hatfield 0028 South Africa<br>&gt; &gt; &gt; &gt; &gt; Tel +27 (0) 12 420 5164<br>&gt; &gt; &gt; &gt; &gt; Fax +27 (0) 12 420 6385<br>&gt; &gt; &gt; &gt; &gt; Mobile:+ 27826228060<br>&gt; &gt; &gt; &gt; &gt; Email address: philuhunga@yahoo.com<br>&gt; &gt; &gt; &gt; &gt; -------------- next part --------------<br>&gt; &gt; &gt; &gt; &gt; An HTML attachment was scrubbed...<br>&gt; &gt; &gt; &gt; &gt; URL: http://gradsusr.org/pipermail/gradsusr/attachments/20141002/51d6cc01/attachment.html <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; <br>&gt; &gt; &gt; &gt; &gt; End of gradsusr Digest, Vol 56, Issue 4<br>&gt; &gt; &gt; &gt; &gt; ***************************************<br>&gt; &gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; <br>&gt; &gt; &gt; <br>&gt; &gt; &gt; <br>&gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; _______________________________________________<br>&gt; &gt; &gt; gradsusr mailing list<br>&gt; &gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; &gt; &gt; <br>&gt; &gt; &gt; <br>&gt; &gt; &gt; End of gradsusr Digest, Vol 56, Issue 6<br>&gt; &gt; &gt; ***************************************<br>&gt; &gt; _______________________________________________<br>&gt; &gt; gradsusr mailing list<br>&gt; &gt; gradsusr@gradsusr.org<br>&gt; &gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; ------------------------------<br>&gt; <br>&gt; _______________________________________________<br>&gt; gradsusr mailing list<br>&gt; gradsusr@gradsusr.org<br>&gt; http://gradsusr.org/mailman/listinfo/gradsusr<br>&gt; <br>&gt; <br>&gt; End of gradsusr Digest, Vol 56, Issue 10<br>&gt; ****************************************<br><BR></div>                                               </div></body>
</html>