[gradsusr] stnmap error

[Huddleston, John]

Andre

Attached is a ZIP container with the MAP and DAT files computed from your FORTRAN and using the stnmap application. i.e.

gfortran -o makestations makestations.f
stnmap -I grads.ctl

It works for the GrADS Windows application.

Since your DAT file differs from the one I created, I wonder if you are running this on a big endian system?

Other than that, I cannot think of a reason why your FORTRAN compiled version would create a different DAT file.

John Huddleston
Cooperative Institute for Research in the Atmosphere

From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Andre Pattantyus
Sent: Sunday, June 19, 2011 9:39 AM
To: gradsusr at gradsusr.org
Subject: [gradsusr] stnmap error

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?

This is the fortran script to convert data:

   CHARACTER*8 STID
    REAL  TD, TT,RH,WSP,WDIR,PP
       OPEN (8,file='stns_new.txt')
      OPEN (10,file='stn.dat',FORM='UNFORMATTED'
     + ,status='unknown')

       IFLAG = 0

10     READ (8,9000,END=90) ITIME, STID, RLAT, RLON, TT,Td,rh,WSP,
     +   wdir, pp

9000   FORMAT( I2, 1X, A4, 1X, F7.4,1X,F8.3,1X,
     + F4.1,1X,F5.1,1X,F3.0,1x,F4.1,1X,F3.0,1X,F4.1)
!
       IF (IFLAG.EQ.0) THEN
          IFLAG = 1
!         IYROLD = IYEAR
          ITIMEOLD = ITIME
       ENDIF

!    print*,IYEAR,IMONTH,STID,RLAT,RLON,RVAL,TD,UU,VV
!  If new time group, write time group terminator.
!  Assuming no empty time groups.
!
       IF (ITIMEOLD.NE.ITIME) THEN
          NLEV = 0
          WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG
          ENDIF
!         IYROLD = IYEAR
          ITIMEOLD = ITIME
!  Write this report
       TIM = 0.0
       NLEV = 1
       NFLAG = 1
       WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG

!print*,ITIME,STID,RLAT,RLON,TT,TD,UU,VV
!      print*,STID,RLAT,RLON,TIM,NLEV,NFLAG
       WRITE (10) TT,TD,RH,WSP,WDIR,PP
       print*, TT ,TD, RH,WSP, WDIR, PP
       GO TO 10
!  On end of file write last time group terminator.
90     CONTINUE

       NLEV = 0
       WRITE (10) STID,RLAT,RLON,TIM,NLEV,NFLAG
       print*,STID,RLAT,RLON,TIM,NLEV,NFLAG
       STOP
       END

Here is the .ctl file:

DSET   ^stn.dat
DTYPE  station
STNMAP station.map
OPTIONS sequential
UNDEF  -999.00
TITLE  Station Data Sample
TDEF  24 linear 07Z30Nov2007 60mn
VARS 6
ts    0 99  Surface Temperature
td    0 99  Surface Dewpoint Temperature
rh    0 99  Relative Humidity
wsp   0 99  Wind speed
wdir  0 99  wind dir
pp    0 99  accumulated precipitation
ENDVARS


--
Andre Pattantyus, Graduate Student Research Assistant
Marine and Environmental Systems, Florida Institute of Technology
150 W. University Blvd, Melbourne, FL 32901

Phone: (321) 674-8330 | Fax: (321) 674-7212 | Email: apattantyus2008 at fit.edu<mailto:apattantyus2008 at fit.edu>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110620/723ccda9/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: andre.zip
Type: application/x-zip-compressed
Size: 1897 bytes
Desc: andre.zip
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20110620/723ccda9/attachment-0003.bin