[gradsusr] stnmap error
Andre Pattantyus
apattantyus2008 at my.fit.edu
Mon Jun 20 14:36:41 EDT 2011
Thank you John for the files. I don't think it is compiler dependent because
I tried on a different machine yesterday with the same result. The files you
created work however I cannot plot out time series values for the individual
stations. I can only plot station values for the fields in file. Does anyone
know how this can be resolved?
On Mon, Jun 20, 2011 at 9:56 AM, Huddleston, John <
Huddleston at cira.colostate.edu> wrote:
> 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****
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
--
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110620/1496bacb/attachment-0003.html
More information about the gradsusr
mailing list