<html dir="ltr"><head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="GENERATOR" content="MSHTML 9.00.8112.16430">
<style id="owaTempEditStyle"></style><style title="owaParaStyle"><!--P {
MARGIN-TOP: 0px; MARGIN-BOTTOM: 0px
}
--></style>
</head>
<body ocsi="x">
<div style="FONT-FAMILY: Tahoma; DIRECTION: ltr; COLOR: #000000; FONT-SIZE: 13px">
<div>Andre</div>
<div><font size="2" face="tahoma"></font> </div>
<div><font size="2" face="tahoma">There is another option which many scientists use; that is, to take your area of interest and apply an interpolation technique to grid the station data.
</font></div>
<div><font size="2" face="tahoma"></font> </div>
<div><font size="2" face="tahoma">I just finished a two year study of Montana (USA) 30 year precipitation using: ordinary Kriging<a></a>; inverse distance squared; regression on value and distance; multivariate inverse distance squared; regression on distance,
value, and elevation; and universal Kriging<a></a> on elevation, distance, and value.
</font></div>
<div><font size="2" face="tahoma"></font> </div>
<div><font size="2" face="tahoma">You may want to think deeper about your studies on the best approach.</font></div>
<div><font size="2" face="tahoma"></font> </div>
<font size="2" face="tahoma">
<p>John Huddleston</font></p>
<div><font size="2" face="tahoma">Cooperative Institute for Research in the Atmosphere</font></div>
<div><font size="2" face="tahoma"></font> </div>
<div dir="ltr"><font color="#000000" size="2" face="Tahoma"></font> </div>
<div style="DIRECTION: ltr" id="divRpF385906">
<hr tabindex="-1">
<font color="#000000" size="2" face="Tahoma"><b>From:</b> gradsusr-bounces@gradsusr.org [gradsusr-bounces@gradsusr.org] On Behalf Of Andre Pattantyus [apattantyus2008@my.fit.edu]<br>
<b>Sent:</b> Monday, June 20, 2011 12:36 PM<br>
<b>To:</b> GrADS Users Forum<br>
<b>Subject:</b> Re: [gradsusr] stnmap error<br>
</font><br>
</div>
<div></div>
<div>
<div> </div>
<div>Thank you John for the files. I don't think it is compiler dependent because I tried on a different machine yesterday with the same result. The files you created work however I cannot plot out time series values for the individual stations. I can only
plot station values for the fields in file. Does anyone know how this can be resolved?<br>
<br>
</div>
<div class="gmail_quote">On Mon, Jun 20, 2011 at 9:56 AM, Huddleston, John <span dir="ltr">
<<a href="mailto:Huddleston@cira.colostate.edu">Huddleston@cira.colostate.edu</a>></span> wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<div lang="EN-US">
<div>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">Andre<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">Attached is a ZIP container with the MAP and DAT files computed from your FORTRAN and using the stnmap application. i.e.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">gfortran –o makestations makestations.f<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">stnmap –I grads.ctl<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">It works for the GrADS Windows application.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">Since your DAT file differs from the one I created, I wonder if you are running this on a big endian system?<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">Other than that, I cannot think of a reason why your FORTRAN compiled version would create a different DAT file.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">John Huddleston<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt">Cooperative Institute for Research in the Atmosphere<u></u><u></u></span></p>
<p class="MsoNormal"><span style="COLOR: #1f497d; FONT-SIZE: 11pt"><u></u><u></u></span> </p>
<p class="MsoNormal"><b><span style="FONT-SIZE: 10pt">From:</span></b><span style="FONT-SIZE: 10pt">
<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>]
<b>On Behalf Of </b>Andre Pattantyus<br>
<b>Sent:</b> Sunday, June 19, 2011 9:39 AM<br>
<b>To:</b> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<b>Subject:</b> [gradsusr] stnmap error<u></u><u></u></span></p>
<div>
<div></div>
<div class="h5">
<p class="MsoNormal"><u></u><u></u> </p>
<p class="MsoNormal">Attached is the station data I am trying to use to convert to grads format. It looks like the script to turn into grads binary is a problem or the options sequential needs to be something else in the .ctl file, but what should it be?<br>
<br>
This is the fortran script to convert data:<br>
<br>
CHARACTER*8 STID <br>
REAL TD, TT,RH,WSP,WDIR,PP<br>
OPEN (8,file='stns_new.txt') <br>
OPEN (10,file='stn.dat',FORM='UNFORMATTED'<br>
+ ,status='unknown')<br>
<br>
IFLAG = 0 <br>
<br>
10 READ (8,9000,END=90) ITIME, STID, RLAT, RLON, TT,Td,rh,WSP, <br>
+ wdir, pp<br>
<br>
9000 FORMAT( I2, 1X, A4, 1X, F7.4,1X,F8.3,1X, <br>
+ F4.1,1X,F5.1,1X,F3.0,1x,F4.1,1X,F3.0,1X,F4.1)<br>
!<br>
IF (IFLAG.EQ.0) THEN<br>
IFLAG = 1 <br>
! IYROLD = IYEAR<br>
ITIMEOLD = ITIME<br>
ENDIF <br>
<br>
! print*,IYEAR,IMONTH,STID,RLAT,RLON,RVAL,TD,UU,VV<br>
! If new time group, write time group terminator. <br>
! Assuming no empty time groups. <br>
!<br>
IF (ITIMEOLD.NE.ITIME) THEN <br>
NLEV = 0 <br>
WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG <br>
ENDIF <br>
! IYROLD = IYEAR <br>
ITIMEOLD = ITIME<br>
! Write this report <br>
TIM = 0.0 <br>
NLEV = 1 <br>
NFLAG = 1<br>
WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG <br>
<br>
!print*,ITIME,STID,RLAT,RLON,TT,TD,UU,VV<br>
! print*,STID,RLAT,RLON,TIM,NLEV,NFLAG<br>
WRITE (10) TT,TD,RH,WSP,WDIR,PP<br>
print*, TT ,TD, RH,WSP, WDIR, PP<br>
GO TO 10<br>
! On end of file write last time group terminator. <br>
90 CONTINUE <br>
<br>
NLEV = 0 <br>
WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG <br>
print*,STID,RLAT,RLON,TIM,NLEV,NFLAG<br>
STOP<br>
END<br>
<br>
Here is the .ctl file:<br clear="all">
<br>
DSET ^stn.dat<br>
DTYPE station<br>
STNMAP station.map<br>
OPTIONS sequential<br>
UNDEF -999.00<br>
TITLE Station Data Sample<br>
TDEF 24 linear 07Z30Nov2007 60mn<br>
VARS 6<br>
ts 0 99 Surface Temperature<br>
td 0 99 Surface Dewpoint Temperature<br>
rh 0 99 Relative Humidity<br>
wsp 0 99 Wind speed<br>
wdir 0 99 wind dir<br>
pp 0 99 accumulated precipitation<br>
ENDVARS<br>
<br>
<br>
-- <br>
Andre Pattantyus, Graduate Student Research Assistant<br>
Marine and Environmental Systems, Florida Institute of Technology<br>
150 W. University Blvd, Melbourne, FL 32901<br>
<br>
Phone: <a href="tel:%28321%29%20674-8330" target="_blank" value="+13216748330">(321) 674-8330</a> | Fax:
<a href="tel:%28321%29%20674-7212" target="_blank" value="+13216747212">(321) 674-7212</a> | Email:
<a href="mailto:apattantyus2008@fit.edu">apattantyus2008@fit.edu</a><u></u><u></u></p>
</div>
</div>
</div>
</div>
<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>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
Andre Pattantyus, Graduate Student Research Assistant<br>
Marine and Environmental Systems, Florida Institute of Technology<br>
150 W. University Blvd, Melbourne, FL 32901<br>
<br>
Phone: (321) 674-8330 | Fax: (321) 674-7212 | Email: <a href="mailto:apattantyus2008@fit.edu">
apattantyus2008@fit.edu</a><br>
</div>
</div>
</body>
</html>