<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Dear Jeff,<BR>&nbsp;<BR>I did a little change, but it didn't work again. In fact, I didn't get what you mean exactely. The errors I got&nbsp;are : <BR>&nbsp;<BR>Define error: invalid expression<BR>close error: only the last file may be closed<BR>error from pow: too many or two few arguments<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; two arguments expected<BR>operation error: error from pow function<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; error from sqrt function<BR>&nbsp;<BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'></span>&nbsp;<BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'>t=1<br>pow=0<br>'reinit'<br>while(t&lt;=4)<br>'sdfopen file1.nc'<br>lon1= 100<br>lat1= 21<br>TK.1= TK<br>'close 1'<br>'open file2.ctl'<br>TK.2= TK<br>'pow1 = pow((TK.1-TK.2))'<br>'pow = pow + pow1'<br>t=t+1<br>endwhile<br>'RMSE= sqrt(pow/t)'<br>
</span><BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'></span>&nbsp;<BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'>How can I also show the RMSE value on the screen?</span><BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'></span>&nbsp;<BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'>Regards,<br style="mso-special-character: line-break;">
<br style="mso-special-character: line-break;">
T.</span><BR><span style='color: rgb(42, 42, 42); line-height: 115%; font-family: "Calibri","sans-serif"; font-size: 11pt; mso-ascii-theme-font: minor-latin; mso-fareast-font-family: Calibri; mso-fareast-theme-font: minor-latin; mso-hansi-theme-font: minor-latin; mso-bidi-font-family: Arial; mso-bidi-theme-font: minor-bidi; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;'></span><font color="#2a2a2a" face="Calibri"></font><br>&nbsp;<BR><div>&nbsp;<br>&gt; <br>&gt; Today's Topics:<br>&gt; <br>&gt;    1. Re: How to define this formula? (Jeff Duda)<br>&gt; <br>&gt; <br>&gt; ----------------------------------------------------------------------<br>&gt; <br>&gt; Message: 1<br>&gt; Date: Sun, 3 Feb 2013 12:29:00 -0600<br>&gt; From: Jeff Duda &lt;jeffduda319@gmail.com&gt;<br>&gt; Subject: Re: [gradsusr] How to define this formula?<br>&gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; Message-ID:<br>&gt;         &lt;CAAig09DMieY67G5iX33RA-uy_he_ZYVXwwcB3qqqOTktvt4ZrQ@mail.gmail.com&gt;<br>&gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt; <br>&gt; T.,<br>&gt; Whenever you have multiple files open simultaneously, you need to use a .n<br>&gt; appended to the variable name to reference the file.  It seems in your<br>&gt; case, the variable in the second file is also TK, but you would need to<br>&gt; refer to it as TK.2.  Also, you have open statements in your while loop,<br>&gt; but no close statements.  Thus you will continue to open new files that<br>&gt; will need successively higher reference numbers (.3, .4, etc).  Either use<br>&gt; an index variable outside of quotes as the reference OR close each file at<br>&gt; the end of each iteration of the loop.  Finally, remember to enclose<br>&gt; statements in quotes if you would enter them at the Grads command prompt<br>&gt; (like setting lat and lon values).<br>&gt; <br>&gt; Jeff<br>&gt; <br>&gt; On Sun, Feb 3, 2013 at 12:14 PM, T. haghroosta &lt;haghroosta@hotmail.com&gt;wrote:<br>&gt; <br>&gt; &gt;  Dear Jeff Duda,<br>&gt
; &gt;<br>&gt; &gt; I did your suggestion as follows. As an extra explanation, I would like to<br>&gt; &gt; calculate RMSE for a special parameter like TK during four days in a<br>&gt; &gt; special Lon and Lat. I also have two files, observed and simulated. So, I<br>&gt; &gt; wrote this script, but it didn't work. Could you please help me? Another<br>&gt; &gt; question how can I show the last output (RMSE)?<br>&gt; &gt;<br>&gt; &gt;  t=1<br>&gt; &gt; pow=0<br>&gt; &gt; 'reinit'<br>&gt; &gt; while(t&lt;=4)<br>&gt; &gt; 'sdfopen observed.nc'<br>&gt; &gt; lon= 100<br>&gt; &gt; lat= 21<br>&gt; &gt; TK1=TK<br>&gt; &gt; 'close 1'<br>&gt; &gt; 'open simulated.ctl'<br>&gt; &gt; TK2=TK<br>&gt; &gt; 'pow1=pow(TK1(t)-TK2(t))<br>&gt; &gt; 'pow = pow + pow1<br>&gt; &gt; t=t+1<br>&gt; &gt; endwhile<br>&gt; &gt; 'RMSE= sqrt(pow/t)'<br>&gt; &gt;<br>&gt; &gt;  Regards,<br>&gt; &gt;<br>&gt; &gt; T.<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Today's Topics:<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; 1. How to define this formula? (T. haghroosta)<br>&gt; &gt; &gt; 2. Re: How to define this formula? (Jeff Duda)<br>&gt; &gt; &gt; 3. Re: lon lat to km (wendi harjupa)<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: Sun, 3 Feb 2013 05:24:48 +0000<br>&gt; &gt; &gt; From: "T. haghroosta" &lt;haghroosta@hotmail.com&gt;<br>&gt; &gt; &gt; Subject: [gradsusr] How to define this formula?<br>&gt; &gt; &gt; To: &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; Message-ID: &lt;DUB121-W214D8BDEBA3D102D60AD37DD020@phx.gbl&gt;<br>&gt; &gt; &gt; Content-Type: text/plain; charset="windows-1256"<br>&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; Dear All, Could you please tell me how can I define this equation which<br>&gt; &gt; I have attached to the email. I got confused completely, I did somethings<br>&gt; &gt;
 that don't work.<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Regards, T.<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;<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;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; -------------- next part --------------<br>&gt; &gt; &gt; An HTML attachment was scrubbed...<br>&gt; &gt; &gt; URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130203/f3a4e!<br>&gt; &gt; 59a/attachment-0001.html<br>&gt; &gt; &gt; -------------- next part ----------- ---<br>&gt; &gt; &gt; A non-text attachment was scrubbed...<br>&gt; &gt; &gt; Name: D.jpg<br>&gt; &gt; &gt; Type: image/jpeg<br>&gt; &gt; &gt; Size: 6713 bytes<br>&gt; &gt; &gt; Desc: not available<br>&gt; &gt; &gt; Url :<br>&gt; &gt; http://gradsusr.org/pipermail/gradsusr/attachments/20130203/f3a4e59a/attachment-0001.jpg<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Message: 2<br>&gt; &gt; &gt; Date: Sun, 3 Feb 2013 01:53:22 -0600<br>&gt; &gt; &gt; From: Jeff Duda &lt;jeffduda319@gmail.com&gt;<br>&gt; &gt; &gt; Subject: Re: [gradsusr] How to define this formula?<br>&gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; Message-ID:<br>&gt; &gt; &gt; &lt;CAAig09DhORYVv_MLG58R5SEAv9pypU41igt_7XwYt97v9Ljocw@mail.gmail.com&gt;<br>&gt; &gt; &gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; That looks like standard deviation or RMSE. Either way, assume your array<br>&gt; &gt; &gt; is VAR. Then it's<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; C = sqrt(amean(asum(pow(VAR -<br>&gt; &gt; &gt; amean(VAR,x=x1,x=x2,y=y1,y=y2)),2),x=x1,x=x2,y=y1,y=y2))<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; OR<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; C = sqrt(asum(V! AR - amean(VAR,x=x1,x=x2,y=y1,y=y2))/N), where N is the<br>&gt; &gt;<br>&gt; &gt; &gt;
 number of grid points over which you're summing.<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; You'll have to insert values for x1, x2, y1, and y2 yourself. You can use<br>&gt; &gt; &gt; lat and lon instead of y and x, though.<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; If that fails, try breaking it into parts:<br>&gt; &gt; &gt; 'average = amean(VAR,x=x1,x=x2,y=y1,y=y2)'<br>&gt; &gt; &gt; 'diffsquaredsum = asum(pow(VAR-average,2),x=x1,x=x2,y=y1,y=y2)'<br>&gt; &gt; &gt; 'SD (or) RMSE = sqrt(diffsquaredsum/N)'<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Jeff Duda<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; On Sat, Feb 2, 2013 at 11:24 PM, T. haghroosta &lt;haghroosta@hotmail.com<br>&gt; &gt; &gt;wrote:<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Dear All,<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Could you please tell me how can I define this equation which I have<br>&gt; &gt; &gt; &gt; attached to the email. I got confused completely, I did somethings that<br>&gt; &gt; &gt; &gt; don't work.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Regards,<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; T.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; **** ** ** ** ** ** ** ** ** ** ** ** ** ** ** *!<br>&gt; &gt; *************************<br>&gt; &gt; &gt; &gt; ** **** ********<br>&gt; &gt; &gt; &gt;*&gt; &gt; _______________________________________________<br>&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; &gt;<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; --<br>&gt; &gt; &gt; Jeff Duda<br>&gt; &gt; &gt; Graduate research assistant<br>&gt; &gt; &gt; University of Oklahoma School of Meteorology<br>&gt; &gt; &gt; Center for Analysis and Prediction of Storms<br>&gt; &gt; &gt; -------------- next part --------------<br>&gt; &gt; &gt; An HTML attachment was scrubbed...<br>&gt; &gt; &gt; URL:<br>&gt; &gt; http://gradsusr.org/pipermail/gradsusr/attachments/20130203/0cd5f463/attachment-0001.html<b
r>&gt; &gt; &gt;<br>&gt; &gt; &gt; ------------------------------<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Message: 3<br>&gt; &gt; &gt; Date: Sun, 3 Feb 2013 21:55:53 +0900<br>&gt; &gt; &gt; From: wendi harjupa &lt;wendiharjupa@gmail.com&gt;<br>&gt; &gt; &gt; Subject: Re: [gradsusr] lon lat to km<br>&gt; &gt; &gt; To: GrADS Users Forum &lt;gradsusr@gradsusr.org&gt;<br>&gt; &gt; &gt; Message-ID:<br>&gt; &gt; &gt; &lt;CA+u6P9METKDdi6QQkH-FZT7mMzcFSAR17MeuvY-hiOkwYQn6Ow@mail.gmail.com&gt;!<br>&gt; &gt; &gt; Content-Type: text/plain; charset="iso-8859-1"<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Dear Mr. Jeff<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Thank you very much for your guidance as far..<br>&gt; &gt; &gt; But, I am very sorry for asking you many times..<br>&gt; &gt; &gt; since it is very difficult for me..<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; I tried to change some parts of script, but still this did not run..<br>&gt; &gt; &gt; I got error message as below :<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Error from CDIFF: Specified dimension non varying<br>&gt; &gt; &gt; Operation Error: Error from cdiff function<br>&gt; &gt; &gt; Error ocurred at column 31<br>&gt; &gt; &gt; DEFINE error: Invalid expression.<br>&gt; &gt; &gt; *** glibc detected *** grads: corrupted double-linked list: 0x089e2918<br>&gt; &gt; &gt; ***<br>&gt; &gt; &gt; Would you mind to explain more??<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Thank you very much<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; Best regards,<br>&gt; &gt; &gt; wendi<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; On Sun, Feb 3, 2013 at 1:46 AM, Jeff Duda &lt;jeffduda319@gmail.com&gt; wrote:<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Use this code:<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; 'reinit'<br>&gt; &gt; &gt; &gt; 'open try2.ctl'!<br>&gt; &gt; &gt; &gt; 'set grads off'<br>&gt; &gt; &gt; &gt; 'set mproj latlon'<br>&gt; &gt; &gt; &gt; 'set mpdset hires'<br>&gt; &gt; &gt; &gt; 'set gxout shaded'<br>&gt; &gt; &gt; &gt; 'set grid on'<br>&gt; &gt; &gt; &gt; 'set lon 100 101'<br>&gt; &gt; &gt; &gt; 'set lat -1 0'<br>&
gt; &gt; &gt; &gt; 'define dx =<br>&gt; &gt; *6.371e3**cos(-1*(3.14159/180))*cdiff(100,0)*(3.14159/180)'<br>&gt; &gt; &gt; &gt; 'define dy = *6.371e3**cdiff(-1,0)*(3.14159/180))'<br>&gt; &gt; &gt; &gt; 'define xdist = ((lon-100)/0.004492369)*(dx/2)'<br>&gt; &gt; &gt; &gt; 'define ydist = ((lat-0)/0.004522556)*(dy/2)'<br>&gt; &gt; &gt; &gt; 'set z 5'<br>&gt; &gt; &gt; &gt; 'set t 1'<br>&gt; &gt; &gt; &gt; 'd ref'<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; The square brackets were just to show different groupings. When it<br>&gt; &gt; comes<br>&gt; &gt; &gt; &gt; time to code it, you must use only parentheses.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; The above should work.<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; Jeff<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; On Sat, Feb 2, 2013 at 1:41 AM, wendi harjupa &lt;wendiharjupa@gmail.com<br>&gt; &gt; &gt;wrote:<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt;&gt; Dear Mr. Jeff<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Thank you very much for your explanation,<br>&gt; &gt; &gt; &gt;&gt; but I am very sorry, ! I still do not get it,<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; I made new script<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; 'reinit'<br>&gt; &gt; &gt; &gt;&gt; 'open try2.ctl'<br>&gt; &gt; &gt; &gt;&gt; 'set grads off'<br>&gt; &gt; &gt; &gt;&gt; 'set mproj latlon'<br>&gt; &gt; &gt; &gt;&gt; 'set mpdset hires'<br>&gt; &gt; &gt; &gt;&gt; 'set gxout shaded'<br>&gt; &gt; &gt; &gt;&gt; 'set grid on'<br>&gt; &gt; &gt; &gt;&gt; 'set lon 100 101'<br>&gt; &gt; &gt; &gt;&gt; 'set lat -1 0'<br>&gt; &gt; &gt; &gt;&gt; 'define dx =<br>&gt; &gt; *6.371e3**cos(-1*(3.14159/180))*cdiff(100,0)*(3.14159/180)'<br>&gt; &gt; &gt; &gt;&gt; 'define dy = *6.371e3**cdiff(-1,0)*(3.14159/180))'<br>&gt; &gt; &gt; &gt;&gt; 'define xdist = ((lon-100)/0.004492369)*(dx/2]'<br>&gt; &gt; &gt; &gt;&gt; 'define ydist = ((lat-0)/0.004522556)*dy/2]'<br>&gt; &gt; &gt; &gt;&gt; 'set z 5'<br>&gt; &gt; &gt; &gt;&gt; 'set t 1'<br>&gt; &gt; &gt; &gt;&gt; 'd ref
'<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; I checked grads document and I found that 6.371e3 is radius of the<br>&gt; &gt; earth<br>&gt; &gt; &gt; &gt;&gt; (6371 to 6374 km)*.*<br>&gt; &gt; &gt; &gt;&gt; *<br>&gt; &gt; &gt; &gt;&gt; *<br>&gt; &gt; &gt; &gt;&gt; *I tried to change "*6.371e3" to! "0e10" which is mean 0 to 10 km)<br>&gt; &gt; &gt; &gt;&gt; But still the progr am did not run normally.<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; I changed the value of lon and lat in script with 100 and -1.<br>&gt; &gt; &gt; &gt;&gt; but still I got error.<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Please help me to solve it.<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; Thank you very much<br>&gt; &gt; &gt; &gt;&gt; Best regards,<br>&gt; &gt; &gt; &gt;&gt; wendi<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; On Sat, Feb 2, 2013 at 2:21 AM, Jeff Duda &lt;jeffduda319@gmail.com&gt;<br>&gt; &gt; wrote:<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt; Wendi,<br>&gt; &gt; &gt; &gt;&gt;&gt; You need to run those two define commands at the top of my previous<br>&gt; &gt; &gt; &gt;&gt;&gt; email before attempting to define the xdist and ydist fields. Also,<br>&gt; &gt; since<br>&gt; &gt; &gt; &gt;&gt;&gt; I can see your control file information, I can see that dlon is<br>&gt; &gt; 0.004492369<br>&gt; &gt; &gt; &gt;&gt;&gt; and dlat is 0.004522556. So substitute those two values into your<br>&gt; &gt; &gt; &gt;&gt;&gt; expressions for xdist and ydist also.<br>&gt; &gt; &gt; &gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt; Jeff<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; &gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt; On Fri, Feb 1, 2013 at 4:36 AM, wendi harjupa &lt;<br>&gt; &gt; wendiharjupa@gmail.com&gt;wrote:<br>&gt; &gt; &gt; &gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Dear Mr. Jeff<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; I have tried to use the scipt which you suggested,<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; But, still I could not do it. Would you help
 me more...<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; here I send my gs file :<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; 'open try2.ctl'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set grads off'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set mproj latlon'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set mpdset hires'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set gxout shaded'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set grid on'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set lon 100 101'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set lat 0 1'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'define xdist = [(lon 99.7)/dlon]*[dx/2]'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'define ydist = [(lat -1.02)/dlat]*[dy/2]'<br>&gt; &gt; &gt; &gt;&gt;! &gt;&gt; 'set z 5'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; 'set t 1'<br>&gt; &gt; &gt; &gt;&gt; &gt;&gt; 'd ref'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Below is my ctl file :<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; DSET 1013.dat<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; TITLE RADAR<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; OPTIONS TEMPLATE LITTLE_ENDIAN<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; UNDEF -1000<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; XDEF 321 LINEAR 99.688120996 0.004492369<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; YDEF 321 LINEAR -1.083708891 0.004522556<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; ZDEF 39 LINEAR 1.0 0.5<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; TDEF 15 LINEAR 13:00Z10apr2004 4mn<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; VARS 2<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; ref 39 99 Reflectivity(dBZ)<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; vel 39 99 Doppler Velocity(m/s)<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; ENDVARS<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Would you check my gs file, is there something that I have to<br>&gt; &gt; revise??<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Thank you very much.<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; 
&gt; &gt; &amp;g! t;&gt;&gt;&gt; Best regards,<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; wendi<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;&gt; On Fri, Feb 1, 2013 at 10:34 AM, Jeff Duda &lt;jeffduda319@gmail.com<br>&gt; &gt; &gt;wrote:<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; These commands will help.<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; 'define dx =<br>&gt; &gt; 6.371e3*cos(lat*(3.14159/180))*cdiff(lon,x)*(3.14159/180)'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; 'define dy = 6.371e3*cdiff(lat,y)*(3.14159/180))'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; They define grid increments in terms of spatial distance rather<br>&gt; &gt; than<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; geographic distance. If you're trying to change the axes and tick<br>&gt; &gt; marks,<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; you'll have to define an origin and define your x and y distances<br>&gt; &gt; based on<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; that. For instance, if you set your lower left domain corner<br>&gt; &gt; (latitude 0<br>&gt; &gt; &gt; &gt;&gt;&gt;&amp;gt! ;&gt; and longitude 100), then you can define your coordinate<br>&gt; &gt; distance s using<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; 'define xdist = [(lon-100)/dlon]*[dx/2]'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; 'define ydist = [(lat-0)/dlat]*[dy/2]'<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; I think this should work, but I haven't tested it myself, so I<br>&gt; &gt; don't<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; know. dlon and dlat are the step values from the XDEF and YDEF<br>&gt; &gt; lines of<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; the control file, respectively.<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; Jeff Duda<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&
gt;&gt;&gt; On Thu, Jan 31, 2013 at 7:02 AM, wendi harjupa &lt;<br>&gt; &gt; wendiharjupa@gmail.com<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; &gt; wrote:<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Dear Sir/Madam<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; I need help<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; I have problem with grads command,<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; I am not well understand wi! th command gr2XY,<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Actually I want to convert the lon 100 to 101, and lat 0 to 1 to<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; distance (km)<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; let think that lon 100 is 0 km, and lat 0 is 0 km,<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; what command I can use to convert lon and lat to km?<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Thank you very much for your help<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; --<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; -------------<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Wendi Harjupa. ST<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Shimane University Graduate School of Engineering Sciences<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Department of Electronic Control Systems<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; Remote Sensing Laboratory<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt; s119447<br>&gt; &gt; &gt; ! &gt;&gt;&gt;&gt;&gt;&gt; Cellphone : +81-080-4268-6676<br>&gt; &gt; &gt; &gt;&amp;gt ;&gt;&gt;&gt;&gt; email : wendi@rslab.riko.shimane-u.ac.jp<br>&gt; &gt;<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt;&gt;<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;&gt;<br>&gt; &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;&gt; --<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; Jeff Duda<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; Graduate research assistant<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; University of Oklahoma School of Meteorology<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; Center for Analysis and Prediction of Storms<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; gradsusr mailing list<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt;&gt; gradsusr@gradsusr.org<br>&gt; &gt; &gt; &gt;&gt;&amp;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;<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;&gt; Wendi Harjupa. ST<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Shimane University Graduate School of Engineering Sciences<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Department of Electronic Control Systems<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Remote Sensing Laboratory<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; s119447<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; Cellphone : +81-080-4268-6676<br>&gt; &gt; &gt; &gt;&gt;&gt;&gt; email : wendi@rslab.riko.shimane-u.ac.jp<br>&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;&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; &amp;! gt; &gt;&gt;&gt; --<br>&gt; &gt; &gt; &gt;&gt;&gt; Jeff Duda<br>&gt; &gt; &gt; &gt;&gt;&amp;gt ; Graduate research assistant<br>&gt; &gt;<br>&gt; &gt; &gt; &gt;&gt;&gt; University of Oklahoma School of Meteorology<br>&gt; &gt; &gt; &gt;&gt;&gt; Center for Analysis and Prediction of Storms<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;&gt;<br>&gt; &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; Wendi Harjupa. ST<br>&gt; &gt; &gt; &gt;&gt; Shimane University Graduate School of Engineering Sciences<br>&gt; &gt; &gt; &gt;&gt; Department of Electronic Control Systems<br>&gt; &gt; &gt; &gt;&gt; Remote Sensing Laboratory<br>&gt; &gt; &gt; &gt;&gt; s119447<br>&gt; &gt; &gt; &gt;&gt; Cellphone : +81-080-4268-6676<br>&gt; &gt; &gt; &gt;&gt; email : wendi@rslab.riko.shimane-u.ac.jp<br>&gt; &gt; &gt; &gt;&gt;<br>&gt; &gt; &gt; &gt;&gt; _______________________________________________<br>&gt; &gt; &gt; &gt;&gt; gradsusr mailing! list<br>&gt; &gt;<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;<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt; &gt; --<br>&gt; &gt; &gt; &gt; Jeff Duda<br>&gt; &gt; &gt; &gt; Graduate research assistant<br>&gt; &gt; &gt; &gt; University of Oklahoma School of Meteorology<br>&gt; &gt; &gt; &gt; Center for Analysis and Prediction of Storms<br>&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; &gt;<br>&gt; &gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; --<br>&gt; &gt; &gt; -------------<br>&gt; &gt; &gt; Wendi Harjupa. ST<br>&gt; &gt; &gt; Shimane University Graduate School of Engineering Sciences<br>&gt; &gt; &gt; Department of Electronic Control Systems<br>&gt; &gt; &gt; Remote Sensing Laboratory<br>&gt; &gt; &gt; s119447<br>&gt; &gt; &gt; Cellphone : +81-080-4268-6676<br>&gt; &gt; &gt; email : wendi@rslab.riko.shimane-u.ac.jp<br>&gt; &gt; &gt; -------------- next part --------------<br>&gt; &gt; &gt; An HTML ! attachment was scrubbed...<br>&gt; &gt; &gt; URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130203/3be6f25a/attachment.html<br>&gt; &gt; &gt;<br>&gt; &gt; &gt; ------------------------------<br>&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 36, Issue 8<br>&gt; &gt; &gt; ***************************************<br>&gt; &gt; *<br>&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; &gt;<br>&gt; &gt;<br>&gt; <br>&gt; <br>&gt; -- <br>&gt; Jeff Duda<br>&gt; Graduate research assistant<br>&gt; University of Oklahoma School of Meteorology<br>&gt; Center for Analysis and Prediction of Storms<br>&gt; -------------- next part --------------<br>&gt; An HTML attachment was scrubbed...<br>&gt; URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130203/d8f96332/attachment.html <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 36, Issue 11<br>&gt; ****************************************<br></div>                                               </div></body>
</html>