[gradsusr] Query Reg the Station data plotting and gridding

[Geeta Geeta]

Hi John, 
Thanks for the mail. 
****************************
I have tried Tdef 1 as you have told. The data file that I am using is for One time Only (00UTC)
Once again the ascii data file looks like ths.
YYYY MM DD TT LAT      LON        RF
2010   06       qqq    34.3    85.5   123.3
2010   06       qrq    34.6    82.5   023.3
2010   06       qql    34.9    83.5   003.3
2010   06       qqm    35.3    84.5   000.3
2010   06       qqn    36.3    87.5   155.3
**********************************************************
The program for converting it into the Binary for Grads is taken from the Grads mannual. 
Character*8 stid
       open(8,file='rain.ch',status='old')
       open(10,file='rain.bin',status='unknown',form='unformatted'
     +  ,access='stream')

      iflag=0
C read and write
 10    read(8,9000,END=90)Iyear,imonth,stid,rlat,rlon,rval
 9000  format(i4,3x,i2,2x,a8,3f8.1)
       print 9000, Iyear,imonth,stid,rlat,rlon,rval
       if (iflag.eq.0) then
           iflag=1
           iyrold=iyear
           imnold=imonth
       endif
C     write this report
        tim=0.0
        nlev=1
        nflag=1
        write(10)stid,rlat,rlon,tim,nlev,nflag
        write(10)rval
        go to 10
 90     continue
C        nlev=0
C        write(10)stid,rlat,rlon,tim,nlev,nflag
        stop
        end
********************************************************

The rain.bin is created 
-rw-r--r-- 1 root root    160 2010-10-04 20:43 rain.bin
***********************************8

Following is the ctl file structure. 
ga-> !vi stn.ctl
dset ^rain.bin
dtype station
STNMAP  rain.map
undef -999.0
title rainfall
tdef 1 linear 12z18jan1980 12hr
vars 1
rf 0 99 rainfall
endvars
**************************************************

I got the following message from Grads 

Config: v1.8SL11 32-bit little-endian readline printim

Issue 'q config' command for more information.

Landscape mode? (no for portrait):  
GX Package Initialization: Size = 11 8.5 
ga-> open stn.ctl
Scanning description file:  stn.ctl
Open Error:  Can't open Station/Index map file rain.map
  --> The invalid description file record is: 
  --> stnmap  rain.map 

  The data file was not opened. 
ga-> 
****************************************************
I Further need to convert this Rainfall into the Gridded one using the OACRES Function. 

thanks again

geeta

 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101006/06dd5e63/attachment-0003.html