<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'>
Hi All,<br><br>If there are any coders out there that can help, I would greatly appreciate it! I am running the latest grads out of the Scientific Linux 6 (same as CentOS/RHEL) repository and want to set up basic verification for a local Wrf mesoscale model and the 0.5 GFS global models against actuals. <br><br>I create a file called stationtest which has umpteen stations with their latitude and longitude in (OMDB xxxx xxxx) and then want to cron run a batch script which pulls out 48 hours from the GFS at these particular station points. The basic script below runs fine - but I would like to use gr2stn in order to do the bilinear interpolation within the ~50kmX~50km GFS gridspace. The following is my working script and this part runs fine:<br><br>Filename: gfsprinttext<br><br>'open gfsgrads.ctl'<br>returnlinetext = read(stationtest)<br>rtncode = sublin(returnlinetext,1)<br>loclatlon = sublin(returnlinetext,2)<br>while (sublin (rtncode,1) = 0)<br>locationvar = subwrd (loclatlon,1)<br>latvariable = subwrd (loclatlon,2)<br>lonvariable = subwrd (loclatlon,3)<br><br>'set t 1 17'<br>'set lat ' latvariable<br>'set lon ' lonvariable<br>'fprintf (tmp2m) /gfs/grads2/gfsoutput/'locationvar'temp.txt'<br><br>returnlinetext = read (stationtest)<br>rtncode = sublin(returnlinetext,1)<br>
loclatlon = sublin(returnlinetext,2)<br>endwhile<br><br>As I said, this above part works fine - but I have attempted in various ways to incorporate the gr2stn function and am not getting it right, has anybody else done something similar or have suggestions for me?<br>Thanks kindly,<br>Theo<br> </div></body>
</html>