[gradsusr] stnmap error

[Jennifer Adams]

To plot a time series at an individual station,
'set x 1'
'set y 1'
'set t 1 12'
'd var(stid=XXXXXX)'

Where XXXXXX is the station id, not more than 8 characters. You can  
find the station id by doing a spatial plot and 'set stat on' and 'set  
stnprint on' and then GrADS will list all the stations plotted in the  
command window.

--Jennifer


On Jun 20, 2011, at 2:36 PM, Andre Pattantyus wrote:

>
> 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
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110621/97a05337/attachment-0003.html