Station data-Stnmap problem

[Xiucheng Wang]

Hi,
   why do you use two flag variables nflag and nflag2 in your programe? delete nflag2 and try again. 
If it does not work again, you can send me your fortran program, .ctl file and source data.   

xiucheng wang
  ----- Original Message ----- 
  From: Ioannis Koletsis 
  To: GRADSUSR at LIST.CINECA.IT 
  Sent: Wednesday, April 16, 2008 5:26 PM
  Subject: Station data-Stnmap problem


  Dear grads users,
  following the example for the file of rain (http://www.iges.org/grads/gadoc/gadocindex.html), I created the following program for 2 variables u and v 

  character*8 stid
         open (8,name='clear_data.txt')
         open (10,name='data.dat',form='unformatted',access='direct',
       + recl=4)
   10    read(8,9000,end=90) iyear,imonth,stid,rlat,rlon,rval,rval2
   9000  format(i4,1x,i1,1x,i2,1x,f9.5,1x,f9.5,1x,f7.2,1x,f7.2)
         if (iflag.eq.0) then
            iflag=1
            iyrold=iyear
            imnold=imonth
         endif
         if (iyrold.ne.iyear.or.imnold.ne.imonth) then
            nlev=0
         write(10) stid,rlat,rlon,tim,nlev,nflag,nflag2
         endif
         iyrold=iyear
         imnold=imonth
         tim=0.0
         nlev=1
         nflag=1
         nflag2=1
         write(10) stid,rlat,rlon,tim,nlev,nflag,nflag2
         write(10) rval,rval2
         goto 10
   90    continue
         nlev=0
         write(10) stid,rlat,rlon,tim,nlev,nflag,nflag2
         stop
         end

  The compilation was completed successfully and a fortran_pr.e file was created. Afterwards running the file a data. dat file was created, also.
  Then using the following station_data.ctl file

  DSET   ^data.dat
  DTYPE  station 
  STNMAP station_data.map
  UNDEF   -99999
  TITLE  Station Data Sample
  TDEF   1 linear 06z27jul2003 12hr
  VARS 2
  u     0  99  surface u-wind
  v     0  99  surface v-wind  
  ENDVARS

  I execute the command stnmap -i station_data.ctl, I received the following 

  Name of binary data set: data.dat
  Number of times in the data set: 1
  Number of surface variables: 2
  Number of level dependent variables: 0

  Starting scan of station data binary file.
  Binary data file open: data.dat
  Processing time=1
  Invalid station hdr found in station binary file
  Possible causes: Invalid level count in hdr
  Descriptor file mismatch
  File not station data
  Invalid relative time
  levs=-1073751116 flag=0 time=32.625


  Do you have an idea for what is wrong?
  Thanks in advance......



   

  -----------------------------------------------------------              
  Ioannis Koletsis
  National Observatory of Athens
  Institute for Environmental Research
  Lofos Koufou
  15236-Penteli
  Athens-Greece
  Tel: (+30) 210 8109140
  e-mail: koletsis at meteo.noa.gr
             ikoletsis at gmail.com
  -----------------------------------------------------------       
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080417/1709ef1b/attachment.html