<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Hi John, <br>Thanks for the mail. <br>****************************<br>I have tried Tdef 1 as you have told. The data file that I am using is for One time Only (00UTC)<br>Once again the ascii data file looks like ths.<br>YYYY MM DD TT LAT LON RF<br>2010 06 qqq 34.3 85.5 123.3<br>2010 06 qrq 34.6 82.5 023.3<br>2010 06 qql 34.9 83.5 003.3<br>2010 06 qqm 35.3 84.5 000.3<br>2010 06 qqn 36.3 87.5 155.3<br>**********************************************************<br>The program for converting it into the Binary for Grads is taken from the Grads mannual. <br>Character*8 stid<br> open(8,file='rain.ch',status='old')<br> open(10,file='rain.bin',status='unknown',form='unformatted'<br> + ,access='stream')<br><br> iflag=0<br>C read and write<br> 10 read(8,9000,END=90)Iyear,imonth,stid,rlat,rlon,rval<br> 9000 format(i4,3x,i2,2x,a8,3f8.1)<br> print 9000, Iyear,imonth,stid,rlat,rlon,rval<br> if (iflag.eq.0) then<br> iflag=1<br> iyrold=iyear<br> imnold=imonth<br> endif<br>C write this report<br> tim=0.0<br> nlev=1<br> nflag=1<br> write(10)stid,rlat,rlon,tim,nlev,nflag<br> write(10)rval<br> go to 10<br> 90 continue<br>C nlev=0<br>C write(10)stid,rlat,rlon,tim,nlev,nflag<br> stop<br> end<br>********************************************************<br><br>The rain.bin is created <br>-rw-r--r-- 1 root root 160 2010-10-04 20:43 rain.bin<br>***********************************8<br><br>Following is the ctl file structure. <br>ga-> !vi stn.ctl<br>dset ^rain.bin<br>dtype station<br>STNMAP rain.map<br>undef -999.0<br>title rainfall<br>tdef 1 linear 12z18jan1980 12hr<br>vars 1<br>rf 0 99 rainfall<br>endvars<br>**************************************************<br><br>I got the following message from Grads <br><br>Config: v1.8SL11 32-bit little-endian readline printim<br><br>Issue 'q config' command for more information.<br><br>Landscape mode? (no for portrait): <br>GX Package Initialization: Size = 11 8.5 <br>ga-> open stn.ctl<br>Scanning description file: stn.ctl<br>Open Error: Can't open Station/Index map file rain.map<br> --> The invalid description file record is: <br> --> stnmap rain.map <br><br> The data file was not opened. <br>ga-> <br>****************************************************<br>I Further need to convert this Rainfall into the Gridded one using the OACRES Function. <br><br>thanks again<br><br>geeta<br><br> </body>
</html>