Simba,<br>To add to Wesley&#39;s, if you have a grib file or a binary file (eg written by fortran) you can just edit an existing ctl file (without using grib2ctl) and put in the relevant info, xdef, ydef, zdef, tdef vars, levels  etc.<br>
<br>Fir installing grads on a linux box just do &quot;yum install grads&quot; then you do not have to worry about libraries and so on.<br>Zilore<br><br><div class="gmail_quote">On Fri, Apr 29, 2011 at 10:21 AM,  <span dir="ltr">&lt;<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>&gt;</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 &#39;help&#39; 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 &quot;Re: Contents of gradsusr digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: How to use grib2ctl option for one chosen variable<br>
      (Wesley Ebisuzaki)<br>
   2. Installing opengrads in x86_64 GNU/Linux (Roberto Mera)<br>
   3. Re: Installing opengrads in x86_64 GNU/Linux (James T. Potemra)<br>
   4. Re: Installing opengrads in x86_64 GNU/Linux (saeed bayat)<br>
   5. Re: how to find typhoon center via minloc function<br>
      (sushant puranik)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 28 Apr 2011 09:55:20 -0400<br>
From: Wesley Ebisuzaki &lt;<a href="mailto:Wesley.Ebisuzaki@noaa.gov">Wesley.Ebisuzaki@noaa.gov</a>&gt;<br>
Subject: Re: [gradsusr] How to use grib2ctl option for one chosen<br>
        variable<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:4DB971C8.9080904@noaa.gov">4DB971C8.9080904@noaa.gov</a>&gt;<br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Simba,<br>
<br>
    No option for grib2ctl/g2ctl.   The solution depends on why you want<br>
to make<br>
the ctl file for a subset.<br>
<br>
(1)  grib2ctl is choosing the wrong vertical coordinate<br>
<br>
    You have to make a new grib file with only the variables that you want.<br>
    Simple task for wgrib.<br>
<br>
(2) Vertical coordinate is ok.<br>
<br>
     Edit the ctl file by removing all the unwanted fields.  Then change<br>
the &quot;vars N&quot; line<br>
     so that N is the new number of fields in the control file.<br>
<br>
   Wesley Ebisuzaki<br>
<br>
Chidzambwa Simba wrote:<br>
&gt; Dear All,<br>
&gt; Is there an option that you can use with &#39;grib2ctl&#39; to create a data<br>
&gt; descriptor file for just one or two variables from a group of<br>
&gt; variables? For example for just the u and v wind components?<br>
&gt;<br>
&gt; Regards<br>
&gt; Simba<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; +----------------------------------------------------------------+<br>
&gt;<br>
&gt; Think B4U Print<br>
&gt;<br>
&gt; 1 ream of paper = 6% of a tree and 5.4kg CO2 in the atmosphere<br>
&gt;<br>
&gt; 3 sheets of A4 paper = 1 litre of water<br>
&gt;<br>
&gt; +----------------------------------------------------------------+<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; ------------------------------------------------------------------------<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 28 Apr 2011 16:09:30 -0700<br>
From: Roberto Mera &lt;<a href="mailto:rjmera@coas.oregonstate.edu">rjmera@coas.oregonstate.edu</a>&gt;<br>
Subject: [gradsusr] Installing opengrads in x86_64 GNU/Linux<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:4DB9F3AA.1060809@coas.oregonstate.edu">4DB9F3AA.1060809@coas.oregonstate.edu</a>&gt;<br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
  All,<br>
<br>
I was trying to install opengrads in a x86_64 GNU/Linux machine and I<br>
got the message below. Where should I put/get this libraries? How do I<br>
compile grads so that it works?<br>
<br>
Thanks,<br>
<br>
Robert<br>
<br>
opengrads<br>
<br>
               Welcome to the OpenGrADS Bundle Distribution<br>
               --------------------------------------------<br>
<br>
For additional information enter &quot;opengrads -h&quot;.<br>
<br>
Starting<br>
&quot;/home/server/pi/homes/rjmera/grads-2.0.a9.oga.1/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/grads<br>
-HC 1  &quot; ...<br>
<br>
/home/server/pi/homes/rjmera/grads-2.0.a9.oga.1/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/grads:<br>
error while loading shared libraries: libXaw.so.8: cannot open shared<br>
object file: No such file or directory<br>
<a href="http://maui.OCE.ORST.EDU" target="_blank">maui.OCE.ORST.EDU</a> 69:<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 28 Apr 2011 13:38:38 -1000<br>
From: &quot;James T. Potemra&quot; &lt;<a href="mailto:jimp@hawaii.edu">jimp@hawaii.edu</a>&gt;<br>
Subject: Re: [gradsusr] Installing opengrads in x86_64 GNU/Linux<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:4DB9FA7E.9060209@hawaii.edu">4DB9FA7E.9060209@hawaii.edu</a>&gt;<br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Roberto:<br>
<br>
Try this, from the INSTALL notes:<br>
<br>
&gt; grads: error while loading shared libraries: libXaw.so.7: cannot open<br>
&gt; shared obj<br>
&gt; ect file: No such file or directory<br>
&gt;<br>
&gt; for libXaw (or any other library) is because your system lacks this<br>
&gt; standard shared library. We have include many of these shared<br>
&gt; libraries used by GrADS under<br>
&gt;<br>
&gt;    Contents/Linux/Versions/2.0.a5.oga.3/x86_64/libs<br>
&gt;<br>
&gt; If it complains about some missing shared, copy each missing shared<br>
&gt; library (and only these) from this directory to<br>
&gt;<br>
&gt;    Contents/Linux/Versions/2.0.a5.oga.3/x86_64/gex<br>
&gt;<br>
&gt; and try again from Contents/.<br>
<br>
Jim<br>
<br>
<br>
On 4/28/11 1:09 PM, Roberto Mera wrote:<br>
&gt;    All,<br>
&gt;<br>
&gt; I was trying to install opengrads in a x86_64 GNU/Linux machine and I<br>
&gt; got the message below. Where should I put/get this libraries? How do I<br>
&gt; compile grads so that it works?<br>
&gt;<br>
&gt; Thanks,<br>
&gt;<br>
&gt; Robert<br>
&gt;<br>
&gt; opengrads<br>
&gt;<br>
&gt;                 Welcome to the OpenGrADS Bundle Distribution<br>
&gt;                 --------------------------------------------<br>
&gt;<br>
&gt; For additional information enter &quot;opengrads -h&quot;.<br>
&gt;<br>
&gt; Starting<br>
&gt; &quot;/home/server/pi/homes/rjmera/grads-2.0.a9.oga.1/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/grads<br>
&gt; -HC 1  &quot; ...<br>
&gt;<br>
&gt; /home/server/pi/homes/rjmera/grads-2.0.a9.oga.1/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/grads:<br>
&gt; error while loading shared libraries: libXaw.so.8: cannot open shared<br>
&gt; object file: No such file or directory<br>
&gt; <a href="http://maui.OCE.ORST.EDU" target="_blank">maui.OCE.ORST.EDU</a> 69:<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Fri, 29 Apr 2011 09:08:34 +0430<br>
From: saeed bayat &lt;<a href="mailto:saeedbayat7276@gmail.com">saeedbayat7276@gmail.com</a>&gt;<br>
Subject: Re: [gradsusr] Installing opengrads in x86_64 GNU/Linux<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;<a href="mailto:BANLkTik42mdM_zzcR6cCsXT5GAaAzdmtmw@mail.gmail.com">BANLkTik42mdM_zzcR6cCsXT5GAaAzdmtmw@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Hi Robert<br>
<br>
just try this command for open &quot;opengrads&quot;<br>
<br>
first: copy the &quot;Content&quot; folder to outside of &quot;opengrads&quot; folder, then<br>
<br>
/home/ which path it is/Content/opengrads-------------------&gt; Enter<br>
<br>
I hope this will hwlp you<br>
<br>
Abdolhalim bayat<br>
M.Sc student of Climatology<br>
Ferdowsi University of mashhad(Iran)<br>
<br>
On Fri, Apr 29, 2011 at 3:39 AM, Roberto Mera<br>
&lt;<a href="mailto:rjmera@coas.oregonstate.edu">rjmera@coas.oregonstate.edu</a>&gt;wrote:<br>
<br>
&gt;  All,<br>
&gt;<br>
&gt; I was trying to install opengrads in a x86_64 GNU/Linux machine and I<br>
&gt; got the message below. Where should I put/get this libraries? How do I<br>
&gt; compile grads so that it works?<br>
&gt;<br>
&gt; Thanks,<br>
&gt;<br>
&gt; Robert<br>
&gt;<br>
&gt; opengrads<br>
&gt;<br>
&gt;               Welcome to the OpenGrADS Bundle Distribution<br>
&gt;               --------------------------------------------<br>
&gt;<br>
&gt; For additional information enter &quot;opengrads -h&quot;.<br>
&gt;<br>
&gt; Starting<br>
&gt;<br>
&gt; &quot;/home/server/pi/homes/rjmera/grads-2.0.a9.oga.1/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/grads<br>
&gt; -HC 1  &quot; ...<br>
&gt;<br>
&gt;<br>
&gt; /home/server/pi/homes/rjmera/grads-2.0.a9.oga.1/Contents/Linux/Versions/2.0.a9.oga.1/x86_64/grads:<br>
&gt; error while loading shared libraries: libXaw.so.8: cannot open shared<br>
&gt; object file: No such file or directory<br>
&gt; <a href="http://maui.OCE.ORST.EDU" target="_blank">maui.OCE.ORST.EDU</a> 69:<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20110429/65cefb32/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20110429/65cefb32/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Fri, 29 Apr 2011 18:02:14 +0530<br>
From: sushant puranik &lt;<a href="mailto:sushantpuranik@gmail.com">sushantpuranik@gmail.com</a>&gt;<br>
Subject: Re: [gradsusr] how to find typhoon center via minloc function<br>
To: GrADS Users Forum &lt;<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
Message-ID: &lt;BANLkTinDFzmd=<a href="mailto:ck1Zc_O7KdK0O1JJPikEQ@mail.gmail.com">ck1Zc_O7KdK0O1JJPikEQ@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<br>
Hello Charles<br>
I want to make a jet stream plot. Can i make it with the help of hurricane<br>
center script by changing minloc and min to maxloc and max.<br>
any suggestion?<br>
<br>
thanks<br>
<br>
Sushant<br>
<br>
On Wed, Jan 5, 2011 at 1:01 AM, Charles Seman &lt;<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>&gt;wrote:<br>
<br>
&gt; Jie,<br>
&gt;<br>
&gt; Attached is a file File_Error_1_Grads_1.9-Plot_tc_shi.gs.txt containing the<br>
&gt; text from an email exchange to the GrADS Listserv between Joe Covert and<br>
&gt; Diane Stokes on 9/13/2004 -- on 4/27/06, found via a Google search for &quot;<br>
&gt; <a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a>&quot;: <a href="http://caos.iisc.ernet.in/gslib/plot_tc_shi.gs" target="_blank">http://caos.iisc.ernet.in/gslib/plot_tc_shi.gs</a><br>
&gt;<br>
&gt; The attached GrADS script <a href="http://plot_hurricane_center3.gs" target="_blank">plot_hurricane_center3.gs</a> (along with pause.gsf,<br>
&gt; aGrADS script function used by the plot script) is a locally modified<br>
&gt; version derived from Joe Covert&#39;s script &quot;<a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a>&quot; (<br>
&gt; <a href="mailto:Joe.Covert@noaa.gov">Joe.Covert@noaa.gov</a>) which was posted to the GrADS Listserv on 9/13/2004<br>
&gt; (his script was named <a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a>) and code from hurricane_tracking.txt<br>
&gt; (taken from an email exchange posted to the GrADS Listserv by Kun-Hsuan Chou<br>
&gt; and Arturo Caracas Uribe in Oct, 2004)...<br>
&gt;<br>
&gt; Please find below code from hurricane_tracking.txt (taken from an email<br>
&gt; exchange posted to the GrADS Listserv by Kun-Hsuan Chou and Arturo Caracas<br>
&gt; Uribe in Oct, 2004) illustrating a technique for finding max wind in a 2D<br>
&gt; field... and below that some code from attached script<br>
&gt; <a href="http://plot_hurricane_center3.gs" target="_blank">plot_hurricane_center3.gs</a> to find the min sea-level pressure in a 2D<br>
&gt; field...<br>
&gt;<br>
&gt; code from Kun-Hsuan Chou to find maximum wind speed from<br>
&gt; hurricane_tracking.txt:<br>
&gt; ---<br>
&gt; &#39;d maxloc(max(mag(u,v),lon=120,lon=130),lat=15,lat=25)&#39;<br>
&gt;  line=sublin(result,2)<br>
&gt;  ygrd=subwrd(line,4)<br>
&gt; &#39;d maxloc(max(mag(u,v),lat=15,lat=25),lon=120,lon=130)&#39;<br>
&gt;  line=sublin(result,2)<br>
&gt;  xgrd=subwrd(line,4)<br>
&gt; &#39;set x &#39;xgrd<br>
&gt; lonval = subwrd(result,4)<br>
&gt; &#39;set y &#39;ygrd<br>
&gt; latval = subwrd(result,4)<br>
&gt; &#39;q w2xy &#39;lonval&#39; &#39;latval<br>
&gt; xpos = subwrd(result,3)<br>
&gt; ypos = subwrd(result,6)<br>
&gt; &#39;draw mark 1 &#39;xpos&#39; &#39;ypos&#39; .2&#39;<br>
&gt; ---<br>
&gt;<br>
&gt; sample code from <a href="http://plot_hurricane_center3.gs" target="_blank">plot_hurricane_center3.gs</a> to find hurricane center (see<br>
&gt; script for supporting code):<br>
&gt; ---<br>
&gt; *<br>
&gt; * find minimum &quot;pressure&quot; within box area (x1,x2), (y1,y2)...<br>
&gt; *<br>
&gt;     &#39;set x &#39;x1<br>
&gt;     &#39;set y &#39;y1<br>
&gt;     &#39;set z 1&#39;<br>
&gt;     &#39;set t &#39;tt<br>
&gt;     &#39;set gxout print&#39;<br>
&gt;     nxp = x2-x1+1<br>
&gt; *      say &#39;    nxp = &#39;nxp<br>
&gt;     nyp = y2-y1+1<br>
&gt; *      say &#39;    nyp = &#39;nyp<br>
&gt;     &#39;d minloc(min(psl,y=&#39;y1&#39;,y=&#39;y2&#39;),x=&#39;x1&#39;,x=&#39;x2&#39;)&#39;<br>
&gt;       rec=sublin(result,nxp+3)<br>
&gt;       xc=subwrd(rec,1)<br>
&gt;     &#39;d minloc(min(psl,x=&#39;x1&#39;,x=&#39;x2&#39;),y=&#39;y1&#39;,y=&#39;y2&#39;)&#39;<br>
&gt;       rec=sublin(result,nyp+3)<br>
&gt;       yc=subwrd(rec,1)<br>
&gt;     say<br>
&gt;     say &#39;location of minimum &quot;pressure&quot;...&#39;<br>
&gt;     say<br>
&gt;     say &#39; (xc,yc) = (&#39;xc&#39;,&#39;yc&#39;)&#39;<br>
&gt; *<br>
&gt; * find &quot;world&quot; coordinates of (xc,yc) and convert &quot;world&quot; coordinates<br>
&gt; * to &quot;xy&quot; coordinates for plotting track of hurricane center...<br>
&gt; *<br>
&gt;     &#39;set x &#39;xc<br>
&gt;     lonval = subwrd(result,4)<br>
&gt;     &#39;set y &#39;yc<br>
&gt;     latval = subwrd(result,4)<br>
&gt;     say &#39; (lonval,latval) = (&#39;lonval&#39;,&#39;latval&#39;)&#39;<br>
&gt;     &#39;q w2xy &#39;lonval&#39; &#39;latval<br>
&gt;     xpos = subwrd(result,3)<br>
&gt;     ypos = subwrd(result,6)<br>
&gt;     say &#39; (xpos,ypos) = (&#39;xpos&#39;,&#39;ypos&#39;)&#39;<br>
&gt; *<br>
&gt; * write (xpos,ypos) to output file...<br>
&gt; *<br>
&gt;     res = write (&#39;&#39;<a href="http://hurricane.nf" target="_blank">hurricane.nf</a>&#39;&#39;,&#39;t = &#39;tt&#39; xpos = &#39;xpos&#39; ypos = &#39;ypos&#39;&#39;)<br>
&gt;<br>
&gt;     pause()<br>
&gt; ---<br>
&gt;<br>
&gt; Please let me know if you have any questions.<br>
&gt;<br>
&gt; Hope this helps,<br>
&gt; Chuck<br>
&gt;<br>
&gt; Jie TANG wrote:<br>
&gt;<br>
&gt;&gt;<br>
&gt;&gt; hi,grads folks ,<br>
&gt;&gt;  I am using grads v1.9, trying to find the center of typhoon via minloc<br>
&gt;&gt; function.<br>
&gt;&gt; the key script is shown as below:<br>
&gt;&gt; slplat=min(slp,lat=19,lat=26)<br>
&gt;&gt; tclon=minloc(slplat,lon=119,lon=126)<br>
&gt;&gt; slplon=min(slp,lon=119,lon=126)<br>
&gt;&gt;  tclat=minloc(slplon,lat=19,lat=26)<br>
&gt;&gt;<br>
&gt;&gt; but grads tell me that &quot;function not found min &quot; and when i changed my<br>
&gt;&gt; script to be :<br>
&gt;&gt; tclat=min(minloc(slpt,lon=119,lon=126), lat=19,lat=26)<br>
&gt;&gt; tclon=min(minloc(slpt, lat=19,lat=26 )lon=119,lon=126)<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; how can i finx my scrpit ? thank you .:)<br>
&gt;&gt; --<br>
&gt;&gt;<br>
&gt;&gt; TANG Jie<br>
&gt;&gt; Email: <a href="mailto:totangjie@gmail.com">totangjie@gmail.com</a> &lt;mailto:<a href="mailto:totangjie@gmail.com">totangjie@gmail.com</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Tel: 0086-2154896104<br>
&gt;&gt; Shanghai Typhoon Institute,China<br>
&gt;&gt; ------------------------------------------------------------------------<br>
&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; gradsusr mailing list<br>
&gt;&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt;&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;<br>
&gt; --<br>
&gt;<br>
&gt; Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> should be considered my NOAA<br>
&gt; email address, not <a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>.<br>
&gt;<br>
&gt; ********************************************************************<br>
&gt; Charles Seman                                <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
&gt; U.S. Department of Commerce / NOAA / OAR<br>
&gt; Geophysical Fluid Dynamics Laboratory         voice: <a href="tel:%28609%29%20452-6547" value="+16094526547">(609) 452-6547</a><br>
&gt; 201 Forrestal Road                              fax: <a href="tel:%28609%29%20987-5063" value="+16099875063">(609) 987-5063</a><br>
&gt; Princeton, NJ  08540-6649            <a href="http://www.gfdl.noaa.gov/%7Ecjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>
&gt; ********************************************************************<br>
&gt;<br>
&gt; &quot;The contents of this message are mine personally and do not reflect any<br>
&gt; official or unofficial position of the United States Federal Government,<br>
&gt; the United States Department of Commerce, or NOAA.&quot;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; worked great ... thanks.  joe<br>
&gt;<br>
&gt; Diane Stokes wrote:<br>
&gt;<br>
&gt; &gt; Joe,<br>
&gt; &gt;<br>
&gt; &gt; Read works fine for me in 1.9.<br>
&gt; &gt;<br>
&gt; &gt; I think this is the issue where &#39;pull&#39; in 1.9 throws in a carriage<br>
&gt; &gt; return.  The read statement is not getting the exact filename.<br>
&gt; &gt;<br>
&gt; &gt; Try adding:<br>
&gt; &gt;   file=sublin(file,1)<br>
&gt; &gt; after each:<br>
&gt; &gt;   pull file<br>
&gt; &gt;<br>
&gt; &gt; Diane<br>
&gt; &gt;<br>
&gt; &gt; Joe Covert wrote:<br>
&gt; &gt;<br>
&gt; &gt;&gt; There appears to be a problem with version 1.9 of GrADS &quot;read&quot; function.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; The <a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a> script (see below) uses the read function to open<br>
&gt; the<br>
&gt; &gt;&gt; track data file and the program gives the subject error message &quot;File<br>
&gt; &gt;&gt; Error 1&quot;.<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; jc<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; *  Script to draw an hurricane-track plot.<br>
&gt; &gt;&gt; *  Does little error checking on the input file.<br>
&gt; &gt;&gt; *  Assumes the input file is set up as follows:<br>
&gt; &gt;&gt; *<br>
&gt; &gt;&gt; *    Line 1:  Title<br>
&gt; &gt;&gt; *    Line 2:  Drawing primitives for marks: marktype size<br>
&gt; &gt;&gt; *    Line 3:  Drawing primitives for lines: color style thickness<br>
&gt; &gt;&gt; *    Line 4:  Starting hour and the interval of plotting points<br>
&gt; &gt;&gt; *             e.g., 0 6 means that track starts at 0 hour and mark<br>
&gt; &gt;&gt; *                   will be plotted every 6 hours.<br>
&gt; &gt;&gt; *    Rest of lines:  hour  long.  lat.<br>
&gt; &gt;&gt; *             e.g.,   0    -70.5  25.0<br>
&gt; &gt;&gt; *                     6    -71.8  25.2<br>
&gt; &gt;&gt; *                            :<br>
&gt; &gt;&gt; *                            :<br>
&gt; &gt;&gt; *<br>
&gt; &gt;&gt; *  Also assumes that a file has been opened (any file, doesn&#39;t<br>
&gt; &gt;&gt; *  matter -- the set command doesn&#39;t work until a file has been<br>
&gt; &gt;&gt; *  opened).<br>
&gt; &gt;&gt; *<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; function main()<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; *  &#39;clear&#39;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   &#39;open dummy.ctl&#39;<br>
&gt; &gt;&gt;   &#39;set lat 20 50&#39;<br>
&gt; &gt;&gt;   &#39;set lon -90 -30&#39;<br>
&gt; &gt;&gt;   &#39;set mpdset hires&#39;<br>
&gt; &gt;&gt;   &#39;set poli on&#39;<br>
&gt; &gt;&gt;   &#39;draw map&#39;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; say &#39;Enter File Name: (type q to stop)&#39;<br>
&gt; &gt;&gt; pull fname<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; while (fname != &#39;q&#39;)<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; *  Read the 1st record: Title<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   ret = read(fname)<br>
&gt; &gt;&gt;   rc = sublin(ret,1)<br>
&gt; &gt;&gt;   if (rc&gt;0)<br>
&gt; &gt;&gt;       say &#39;File Error 1&#39;<br>
&gt; &gt;&gt;       return<br>
&gt; &gt;&gt;   endif<br>
&gt; &gt;&gt;   title = sublin(ret,2)<br>
&gt; &gt;&gt;   say title<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; *  Read the drawing primitives<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   ret = read(fname)<br>
&gt; &gt;&gt;   rc = sublin(ret,1)<br>
&gt; &gt;&gt;   if (rc&gt;0)<br>
&gt; &gt;&gt;      say &#39;File Error 2&#39;<br>
&gt; &gt;&gt;      return<br>
&gt; &gt;&gt;   endif<br>
&gt; &gt;&gt;   dpline = sublin(ret,2)<br>
&gt; &gt;&gt;   marktype = subwrd(dpline,1)<br>
&gt; &gt;&gt;   marksize = subwrd(dpline,2)<br>
&gt; &gt;&gt;   ret = read(fname)<br>
&gt; &gt;&gt;   rc = sublin(ret,1)<br>
&gt; &gt;&gt;   if (rc&gt;0)<br>
&gt; &gt;&gt;      say &#39;File Error 3&#39;<br>
&gt; &gt;&gt;      return<br>
&gt; &gt;&gt;   endif<br>
&gt; &gt;&gt;   dpline = sublin(ret,2)<br>
&gt; &gt;&gt;   lcolor = subwrd(dpline,1)<br>
&gt; &gt;&gt;   lstyle = subwrd(dpline,2)<br>
&gt; &gt;&gt;   lthick = subwrd(dpline,3)<br>
&gt; &gt;&gt;   say &#39; marktype, marksize, lcolor, lstyle and lthick:&#39;<br>
&gt; &gt;&gt;   say &#39; &#39;marktype &#39; &#39; marksize &#39; &#39; lcolor &#39; &#39; lstyle &#39; &#39; lthick<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; * Read starting hour and the interval hours of plotting points<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   ret = read(fname)<br>
&gt; &gt;&gt;   rc = sublin(ret,1)<br>
&gt; &gt;&gt;   if (rc&gt;0)<br>
&gt; &gt;&gt;      say &#39;File Error 4&#39;<br>
&gt; &gt;&gt;      return<br>
&gt; &gt;&gt;   endif<br>
&gt; &gt;&gt;   dhour = sublin(ret,2)<br>
&gt; &gt;&gt;   start = subwrd(dhour,1)<br>
&gt; &gt;&gt;   jump = subwrd(dhour,2)<br>
&gt; &gt;&gt;   say &#39; starting hour and the interval hours of plotting points:&#39;<br>
&gt; &gt;&gt;   say &#39;  &#39;start&#39; &#39;jump<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; *  Read all data points<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   ret = read(fname)<br>
&gt; &gt;&gt;   rc = sublin(ret,1)<br>
&gt; &gt;&gt;   while (rc = 0)<br>
&gt; &gt;&gt;       loc = sublin(ret,2)<br>
&gt; &gt;&gt;       hour = subwrd(loc,1)<br>
&gt; &gt;&gt;       dlong.hour = subwrd(loc,2)<br>
&gt; &gt;&gt;       dlat.hour = subwrd(loc,3)<br>
&gt; &gt;&gt; *      prompt &#39; hour &#39; hour &#39; are read,&#39;<br>
&gt; &gt;&gt; *      say &#39; long=&#39; dlong.hour &#39;    lat=&#39; dlat.hour<br>
&gt; &gt;&gt;       ret = read(fname)<br>
&gt; &gt;&gt;       rc = sublin(ret,1)<br>
&gt; &gt;&gt;   endwhile<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   if (rc!=2 &amp; rc!=0)<br>
&gt; &gt;&gt;          say &#39;File Error 5, rc=&#39; rc<br>
&gt; &gt;&gt;          return<br>
&gt; &gt;&gt;   endif<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   endhour = hour<br>
&gt; &gt;&gt;   say &#39; endhour=&#39; endhour<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; * Plotting<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   &#39;set line &#39;lcolor&#39; &#39;lstyle&#39; &#39;lthick<br>
&gt; &gt;&gt;   &#39;query w2xy &#39;dlong.start&#39; &#39;dlat.start<br>
&gt; &gt;&gt;   xprev = subwrd(result,3)<br>
&gt; &gt;&gt;   yprev = subwrd(result,6)<br>
&gt; &gt;&gt;   &#39;draw mark &#39;marktype&#39; &#39;xprev&#39; &#39;yprev&#39; &#39;marksize<br>
&gt; &gt;&gt;   next = start+jump<br>
&gt; &gt;&gt;   while (next &lt;= endhour)<br>
&gt; &gt;&gt; *      say &#39; &#39;dlong.start&#39; &#39;dlat.start<br>
&gt; &gt;&gt;       &#39;query w2xy &#39;dlong.next&#39; &#39;dlat.next<br>
&gt; &gt;&gt;       xnext = subwrd(result,3)<br>
&gt; &gt;&gt;       ynext = subwrd(result,6)<br>
&gt; &gt;&gt;       &#39;draw line &#39;xprev&#39; &#39;yprev&#39; &#39;xnext&#39; &#39;ynext<br>
&gt; &gt;&gt; *      say &#39; &#39;xprev&#39; &#39;yprev&#39; &#39;xnext&#39; &#39;ynext<br>
&gt; &gt;&gt;       &#39;draw mark &#39;marktype&#39; &#39;xnext&#39; &#39;ynext&#39; &#39;marksize<br>
&gt; &gt;&gt;       next = next+jump<br>
&gt; &gt;&gt;       xprev = xnext<br>
&gt; &gt;&gt;       yprev = ynext<br>
&gt; &gt;&gt;   endwhile<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   say fname &#39; is working fine.&#39;<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; * read in the filename to be plotted next<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt;   say &#39;Enter File Name: (type q to stop)&#39;<br>
&gt; &gt;&gt;   pull fname<br>
&gt; &gt;&gt;<br>
&gt; &gt;&gt; endwhile<br>
&gt; &gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt; *----------------------------------------------------------------------<br>
&gt; *<br>
&gt; *  Adapted from <a href="http://plot_hurricane_center2c.gs" target="_blank">plot_hurricane_center2c.gs</a>, which was adapted from:<br>
&gt; *<br>
&gt; *  (1) &quot;plot_hurricane_center2[a][b].gs&quot; (from &quot;<a href="http://plot_hurricane_center.gs" target="_blank">plot_hurricane_center.gs</a><br>
&gt; &quot;),<br>
&gt; *      locally modified versions derived from Joe Covert&#39;s script<br>
&gt; *      &quot;<a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a>&quot; (<a href="mailto:Joe.Covert@noaa.gov">Joe.Covert@noaa.gov</a>), which was posted to the<br>
&gt; *      GrADS Listserv on 9/13/2004 (his script was named <a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a>)<br>
&gt; *      on 4/27/06, found via a Google search for &quot;<a href="http://plot_tc_shi.gs" target="_blank">plot_tc_shi.gs</a>&quot;:<br>
&gt; *      <a href="http://caos.iisc.ernet.in/gslib/plot_tc_shi.gs" target="_blank">http://caos.iisc.ernet.in/gslib/plot_tc_shi.gs</a><br>
&gt; *<br>
&gt; *  (2) code in &quot;hurricane_tracking.txt&quot; (taken from an email exchange<br>
&gt; *      posted to the GrADS Listserv by Kun-Hsuan Chou and<br>
&gt; *      Arturo Caracas Uribe in Oct, 2004)<br>
&gt; *<br>
&gt; *<br>
&gt; *  To use: grads -l<br>
&gt; *<br>
&gt; *----------------------------------------------------------------------<br>
&gt;<br>
&gt; function main()<br>
&gt;<br>
&gt;  &#39;reinit&#39;<br>
&gt;<br>
&gt; ************************************************************************<br>
&gt; *  Allow external GrADS functions<br>
&gt; ************************************************************************<br>
&gt;<br>
&gt;  rc = gsfallow(&quot;on&quot;)<br>
&gt;<br>
&gt; ************************************************************************<br>
&gt; *  Define dataset file information and plotting limits for base map...<br>
&gt; ************************************************************************<br>
&gt;<br>
&gt; **  nc_file1 = &#39;&#39;<br>
&gt; **  nc_file2 = &#39;&#39;<br>
&gt; **  nc_file3 = &#39;&#39;<br>
&gt; **  nc_file = &#39;&#39;nc_file1&#39; &#39;nc_file2&#39; &#39;nc_file3&#39;&#39;<br>
&gt;  nc_file = &#39;your_netCDF_file&#39;<br>
&gt;<br>
&gt;  lon1 = 298 ; lon2 = 308<br>
&gt;  lat1 =  20 ; lat2 =  25<br>
&gt;<br>
&gt; ************************************************************************<br>
&gt; *  Legend captions, colors, and marks for the plots<br>
&gt; ************************************************************************<br>
&gt;<br>
&gt; **  nmodels = 1<br>
&gt; **  model.1 = &#39;isabel_ras&#39;<br>
&gt; **  model.2 = &#39;isabel_isotke&#39;<br>
&gt; **  model.3 = &#39;isabel_my25&#39;<br>
&gt;<br>
&gt; **  titles.1 = model.1 ; colors.1 = 2 ; marks.1 = 2<br>
&gt; **  titles.2 = model.2 ; colors.2 = 3 ; marks.2 = 2<br>
&gt; **  titles.3 = model.3 ; colors.3 = 4 ; marks.3 = 2<br>
&gt; ***<br>
&gt; ***  Define the legend plotter input...<br>
&gt; ***<br>
&gt; **  titles = &#39;-t&#39; ; colors = &#39;-c&#39; ; marks = &#39;-m&#39; ; lines = &#39;-l&#39;<br>
&gt; **  n=1<br>
&gt; **  while ( n &lt;= nmodels )<br>
&gt; **    titles = &#39;&#39;titles&#39; &quot;&#39;titles.n&#39;&quot;&#39;<br>
&gt; **    colors = &#39;&#39;colors&#39; &#39;colors.n<br>
&gt; **    marks = &#39;&#39;marks&#39; &quot;&#39;marks.n&#39;&quot;&#39;<br>
&gt; **    lines = &#39;&#39;lines&#39; 1&#39;<br>
&gt; **    n=n+1<br>
&gt; **  endwhile<br>
&gt;<br>
&gt; **  legend_info = &#39;&#39;colors&#39; &#39;lines&#39; &#39;marks&#39; &#39;titles&#39; -p&#39;<br>
&gt;<br>
&gt;  nmodels = 1<br>
&gt;  model.1 = &#39;title_name_for_your_netCDF_dataset&#39;<br>
&gt;<br>
&gt;  titles.1 = model.1 ; colors.1 = 2 ; marks.1 = 2<br>
&gt; *<br>
&gt; *  Define the legend plotter input...<br>
&gt; *<br>
&gt;  titles = &#39;-t&#39; ; colors = &#39;-c&#39; ; marks = &#39;-m&#39; ; lines = &#39;-l&#39;<br>
&gt;  n=1<br>
&gt;  while ( n &lt;= nmodels )<br>
&gt;    titles = &#39;&#39;titles&#39; &quot;&#39;titles.n&#39;&quot;&#39;<br>
&gt;    colors = &#39;&#39;colors&#39; &#39;colors.n<br>
&gt;    marks = &#39;&#39;marks&#39; &quot;&#39;marks.n&#39;&quot;&#39;<br>
&gt;    lines = &#39;&#39;lines&#39; 1&#39;<br>
&gt;    n=n+1<br>
&gt;  endwhile<br>
&gt;<br>
&gt;  legend_info = &#39;&#39;colors&#39; &#39;lines&#39; &#39;marks&#39; &#39;titles&#39; -p&#39;<br>
&gt;<br>
&gt; *  ...other drawing primitives<br>
&gt;<br>
&gt;  marksize = 0.1<br>
&gt;<br>
&gt;  lstyle = 1 ; lthick = 5<br>
&gt;<br>
&gt; ************************************************************************<br>
&gt; *  Use &quot;cbar_line&quot; or &quot;cbar_line_box&quot;?<br>
&gt; ************************************************************************<br>
&gt;<br>
&gt;  legend_plotter = &#39;cbar_line&#39;<br>
&gt; *  legend_plotter = &#39;cbar_line_box&#39;<br>
&gt;<br>
&gt; ************************************************************************<br>
&gt; *  Plot page limits<br>
&gt; ************************************************************************<br>
&gt;<br>
&gt;  plot_vpage = &#39;0 11 0 8.5&#39;<br>
&gt;  plot_area  = &#39;1 10 1 7.5&#39;<br>
&gt;<br>
&gt;  xl = subwrd(plot_area,1) ; xr = subwrd(plot_area,2)<br>
&gt;  yb = subwrd(plot_area,3) ; yt = subwrd(plot_area,4)<br>
&gt;<br>
&gt;  xc = xl + (xr-xl)/2<br>
&gt;<br>
&gt; ************************************************************************<br>
&gt; *  GrADS metafile output...<br>
&gt; ************************************************************************<br>
&gt;<br>
&gt;  &#39;enable print hurricane_track.gx&#39;<br>
&gt;<br>
&gt; *-----------------------------------------------------------------------<br>
&gt; *  Define and plot track of minimum surface pressure for each dataset...<br>
&gt; *-----------------------------------------------------------------------<br>
&gt;<br>
&gt;  nf=1<br>
&gt;  while ( nf &lt;= nmodels )<br>
&gt;<br>
&gt; * Open dataset file and define parameters...<br>
&gt;<br>
&gt;    fname = subwrd(nc_file,nf)<br>
&gt;    &#39;sdfopen &#39;fname<br>
&gt;    &#39;set dfile &#39;nf<br>
&gt;    &#39;q file&#39;<br>
&gt;    rec3 = sublin(result,3) ; binary_file = subwrd(rec3,2)<br>
&gt;    say<br>
&gt;    say<br>
&gt;    say &#39;*** &#39;binary_file&#39; ***&#39;<br>
&gt;    say<br>
&gt;    rec5 = sublin(result,5)<br>
&gt;    nx = subwrd(rec5,3) ; ny = subwrd(rec5,6) ; nz = subwrd(rec5,9) ; nt =<br>
&gt; subwrd(rec5,12)<br>
&gt;<br>
&gt;    if( nf = 1 )<br>
&gt;<br>
&gt; * Plot base map...<br>
&gt;<br>
&gt;      plot_base_map(plot_vpage,plot_area,lon1,lon2,lat1,lat2)<br>
&gt;<br>
&gt;    endif<br>
&gt;<br>
&gt; * Define locations of &quot;hurricane center&quot; for plotting points, and write<br>
&gt; * locations out to ASCII text file for read-in for plotting track later...<br>
&gt;<br>
&gt;    <a href="http://hurricane.nf" target="_blank">hurricane.nf</a> = &#39;hurricane_&#39;nf&#39;.txt&#39;<br>
&gt;<br>
&gt;    ts=2<br>
&gt;    tt=ts<br>
&gt;    while ( tt &lt;= nt )<br>
&gt;      &#39;set t &#39;tt<br>
&gt;      say<br>
&gt;      say &#39;---------&#39;<br>
&gt;      say &#39; tt = &#39;tt<br>
&gt;      say &#39;---------&#39;<br>
&gt;      say<br>
&gt;      &#39;set lon &#39;lon1&#39; &#39;lon2<br>
&gt;      &#39;set lat &#39;lat1&#39; &#39;lat2<br>
&gt;      &#39;set z 1&#39;<br>
&gt; *<br>
&gt; * here is where you put in the name of your sea-level pressure variable:<br>
&gt; *<br>
&gt;      &#39;define psl = name_of_your_sea_level_pressure_variable&#39;<br>
&gt;      &#39;set gxout shaded&#39;<br>
&gt;      &#39;d psl&#39;<br>
&gt; *<br>
&gt; * following interactive &quot;box location&quot; code adapted from &quot;<a href="http://cbar_line_box.gs" target="_blank">cbar_line_box.gs</a><br>
&gt; &quot;...<br>
&gt; *<br>
&gt;      say &#39;Click where you want the left upper corner of the box&#39;<br>
&gt;      &#39;query bpos&#39;<br>
&gt;      xb1o = subwrd(result,3)<br>
&gt;      yb2o = subwrd(result,4)<br>
&gt;      say &#39;Click where you want the right lower corner of the box&#39;<br>
&gt;      &#39;query bpos&#39;<br>
&gt;      xb2o = subwrd(result,3)<br>
&gt;      yb1o = subwrd(result,4)<br>
&gt;<br>
&gt;      say<br>
&gt;      say &#39;...left upper corner of the box at X Y: &#39;xb1o&#39; &#39;yb2o<br>
&gt;      say &#39;...right lower corner of the box at X Y: &#39;xb2o&#39; &#39;yb1o<br>
&gt; *<br>
&gt; * draw box to show specified area for defining hurricane center...<br>
&gt; *<br>
&gt;      &#39;set line 1 1&#39;<br>
&gt;      &#39;draw rec &#39;xb1o&#39; &#39;yb1o&#39; &#39;xb2o&#39; &#39;yb2o<br>
&gt; *<br>
&gt; * convert box &quot;xy&quot; coordinates to &quot;grid&quot; coordinates<br>
&gt; * for finding minimum &quot;pressure&quot; (hurricane center)...<br>
&gt; *<br>
&gt;      &#39;q xy2gr &#39;xb1o&#39; &#39;yb2o<br>
&gt;      x1o = subwrd(result,3) ; x1 = math_int(x1o+0.5)<br>
&gt;      y2o = subwrd(result,6) ; y2 = math_int(y2o+0.5)<br>
&gt;      &#39;q xy2gr &#39;xb2o&#39; &#39;yb1o<br>
&gt;      x2o = subwrd(result,3) ; x2 = math_int(x2o+0.5)<br>
&gt;      y1o = subwrd(result,6) ; y1 = math_int(y1o+0.5)<br>
&gt;      say<br>
&gt;      say &#39;Specified grid coordinates (x1o,x2o), (y1o,y2o) from box,&#39;<br>
&gt;      say &#39;rounded to (x1,x2), (y1,y2) for finding minimum &quot;pressure&quot;...&#39;<br>
&gt;      say<br>
&gt;      say &#39; (x1o,x2o) = (&#39;x1o&#39;,&#39;x2o&#39;)  --&gt;  (x1,x2) = (&#39;x1&#39;,&#39;x2&#39;)&#39;<br>
&gt;      say &#39; (y1o,y2o) = (&#39;y1o&#39;,&#39;y2o&#39;)  --&gt;  (y1,y2) = (&#39;y1&#39;,&#39;y2&#39;)&#39;<br>
&gt; *<br>
&gt; * find minimum &quot;pressure&quot; within box area (x1,x2), (y1,y2)...<br>
&gt; *<br>
&gt;      &#39;set x &#39;x1<br>
&gt;      &#39;set y &#39;y1<br>
&gt;      &#39;set z 1&#39;<br>
&gt;      &#39;set t &#39;tt<br>
&gt;      &#39;set gxout print&#39;<br>
&gt;      nxp = x2-x1+1<br>
&gt; *      say &#39;    nxp = &#39;nxp<br>
&gt;      nyp = y2-y1+1<br>
&gt; *      say &#39;    nyp = &#39;nyp<br>
&gt;      &#39;d minloc(min(psl,y=&#39;y1&#39;,y=&#39;y2&#39;),x=&#39;x1&#39;,x=&#39;x2&#39;)&#39;<br>
&gt;        rec=sublin(result,nxp+3)<br>
&gt;        xc=subwrd(rec,1)<br>
&gt;      &#39;d minloc(min(psl,x=&#39;x1&#39;,x=&#39;x2&#39;),y=&#39;y1&#39;,y=&#39;y2&#39;)&#39;<br>
&gt;        rec=sublin(result,nyp+3)<br>
&gt;        yc=subwrd(rec,1)<br>
&gt;      say<br>
&gt;      say &#39;location of minimum &quot;pressure&quot;...&#39;<br>
&gt;      say<br>
&gt;      say &#39; (xc,yc) = (&#39;xc&#39;,&#39;yc&#39;)&#39;<br>
&gt; *<br>
&gt; * find &quot;world&quot; coordinates of (xc,yc) and convert &quot;world&quot; coordinates<br>
&gt; * to &quot;xy&quot; coordinates for plotting track of hurricane center...<br>
&gt; *<br>
&gt;      &#39;set x &#39;xc<br>
&gt;      lonval = subwrd(result,4)<br>
&gt;      &#39;set y &#39;yc<br>
&gt;      latval = subwrd(result,4)<br>
&gt;      say &#39; (lonval,latval) = (&#39;lonval&#39;,&#39;latval&#39;)&#39;<br>
&gt;      &#39;q w2xy &#39;lonval&#39; &#39;latval<br>
&gt;      xpos = subwrd(result,3)<br>
&gt;      ypos = subwrd(result,6)<br>
&gt;      say &#39; (xpos,ypos) = (&#39;xpos&#39;,&#39;ypos&#39;)&#39;<br>
&gt; *<br>
&gt; * write (xpos,ypos) to output file...<br>
&gt; *<br>
&gt;      res = write (&#39;&#39;<a href="http://hurricane.nf" target="_blank">hurricane.nf</a>&#39;&#39;,&#39;t = &#39;tt&#39; xpos = &#39;xpos&#39; ypos = &#39;ypos&#39;&#39;)<br>
&gt;<br>
&gt;      pause()<br>
&gt; *<br>
&gt; * erase box showing area for defining hurricane center...<br>
&gt; *<br>
&gt; *      &#39;set line 0 1&#39;<br>
&gt; *      &#39;draw rec &#39;xb1o&#39; &#39;yb1o&#39; &#39;xb2o&#39; &#39;yb2o<br>
&gt;      tt=tt+1<br>
&gt;    endwhile<br>
&gt;<br>
&gt;    res = close (&#39;&#39;<a href="http://hurricane.nf" target="_blank">hurricane.nf</a>&#39;&#39;)<br>
&gt;<br>
&gt;    nf=nf+1<br>
&gt;  endwhile<br>
&gt;<br>
&gt; *  Clear frame and re-plot base map...<br>
&gt;<br>
&gt;  plot_base_map(plot_vpage,plot_area,lon1,lon2,lat1,lat2)<br>
&gt;<br>
&gt; *  Plot track(s)...<br>
&gt;<br>
&gt;  nf=1<br>
&gt;  while ( nf &lt;= nmodels )<br>
&gt;<br>
&gt; *  Specify drawing primitives<br>
&gt;<br>
&gt;    marktype = <a href="http://marks.nf" target="_blank">marks.nf</a><br>
&gt;    lcolor = <a href="http://colors.nf" target="_blank">colors.nf</a><br>
&gt;<br>
&gt;    say &#39; marktype, marksize, lcolor, lstyle and lthick:&#39;<br>
&gt;    say &#39; &#39;marktype &#39; &#39; marksize &#39; &#39; lcolor &#39; &#39; lstyle &#39; &#39; lthick<br>
&gt;<br>
&gt;    tt=ts<br>
&gt; *<br>
&gt; *  read (xpos,ypos) from ASCII text input file...<br>
&gt; *<br>
&gt;    res = read(&#39;&#39;<a href="http://hurricane.nf" target="_blank">hurricane.nf</a>&#39;&#39;) ; rc = sublin(res,1)<br>
&gt;    if( rc != 0 )<br>
&gt;      say &#39;&quot;read&quot; status = &#39;rc<br>
&gt;      exit<br>
&gt;    endif<br>
&gt;    rec = sublin(res,2) ; say &#39;&quot;hurricane.&#39;nf&#39;&quot; record: &#39;rec<br>
&gt;    xc = subwrd(rec,6) ; yc = subwrd(rec,9)<br>
&gt;<br>
&gt;    &#39;set line &#39;lcolor&#39; &#39;lstyle&#39; &#39;lthick<br>
&gt;    &#39;draw mark &#39;marktype&#39; &#39;xc&#39; &#39;yc&#39; &#39;marksize<br>
&gt;    tt = tt+1<br>
&gt;    while ( tt &lt;= nt )<br>
&gt;      res = read(&#39;&#39;<a href="http://hurricane.nf" target="_blank">hurricane.nf</a>&#39;&#39;) ; rc = sublin(res,1)<br>
&gt;      if( rc != 0 )<br>
&gt;        say &#39;&quot;read&quot; status = &#39;rc<br>
&gt;        exit<br>
&gt;      endif<br>
&gt;      rec = sublin(res,2) ; say &#39;&quot;hurricane.&#39;nf&#39;&quot; record: &#39;rec<br>
&gt;      xn = subwrd(rec,6) ; yn = subwrd(rec,9)<br>
&gt;      say &#39; &#39;xc&#39; &#39;yc&#39; &#39;xn&#39; &#39;yn<br>
&gt;      &#39;draw line &#39;xc&#39; &#39;yc&#39; &#39;xn&#39; &#39;yn<br>
&gt;      &#39;draw mark &#39;marktype&#39; &#39;xn&#39; &#39;yn&#39; &#39;marksize<br>
&gt;      tt=tt+1<br>
&gt;      xc = xn<br>
&gt;      yc = yn<br>
&gt;    endwhile<br>
&gt;<br>
&gt;    nf=nf+1<br>
&gt;  endwhile<br>
&gt;<br>
&gt;  say<br>
&gt;  say &#39;...plot legend&#39;<br>
&gt;  &#39;&#39;legend_plotter&#39; -x 3 -y 6 &#39;legend_info<br>
&gt;  say<br>
&gt;<br>
&gt;  &#39;print&#39;<br>
&gt;  &#39;disable print&#39;<br>
&gt;<br>
&gt;  say<br>
&gt;  say &#39;************************************************&#39;<br>
&gt;  say &#39; Finished with this script.&#39;<br>
&gt;  say &#39;************************************************&#39;<br>
&gt;  say<br>
&gt;<br>
&gt; return<br>
&gt;<br>
&gt; function plot_base_map(plot_vpage,plot_area,lon1,lon2,lat1,lat2)<br>
&gt;<br>
&gt; *  Plots base map...<br>
&gt;<br>
&gt;  &#39;c&#39;<br>
&gt;  &#39;set vpage &#39;plot_vpage<br>
&gt;  &#39;set parea &#39;plot_area<br>
&gt;  &#39;set lon &#39;lon1&#39; &#39;lon2<br>
&gt;  &#39;set lat &#39;lat1&#39; &#39;lat2<br>
&gt;  &#39;set mproj latlon&#39;<br>
&gt;  &#39;set mpt * 1 1 5&#39;<br>
&gt;  &#39;set mpdset hires&#39;<br>
&gt;  &#39;set poli on&#39;<br>
&gt; *  &#39;set grid on&#39;<br>
&gt;  &#39;set grads off&#39;<br>
&gt;  &#39;set gxout contour&#39;<br>
&gt;  &#39;draw map&#39;<br>
&gt;<br>
&gt; *  add lon/lat lines...<br>
&gt;<br>
&gt;  &#39;set ccolor 1&#39;<br>
&gt;  &#39;set cstyle 5&#39; ; &#39;d lon&#39;<br>
&gt;  &#39;set cstyle 5&#39; ; &#39;d lat&#39;<br>
&gt;<br>
&gt; *  add title, x- and y-labels...<br>
&gt;<br>
&gt;  &#39;draw title \\ Track of Minimum Surface Pressure&#39;<br>
&gt;  &#39;draw xlab Longitude \\&#39;<br>
&gt;  &#39;draw ylab \\ Latitude&#39;<br>
&gt;<br>
&gt; return<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20110429/e2936a21/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20110429/e2936a21/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 14, Issue 62<br>
****************************************<br>
</blockquote></div><br>