[gradsusr] stnmap error

Tue Jun 21 14:04:26 EDT 2011

thanks. that works great. never would have figured that out. Is it in
documentation somewhere?

On Tue, Jun 21, 2011 at 6:52 AM, Jennifer Adams <jma at cola.iges.org> wrote:

> 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
>
>
>
>
> _______________________________________________
> 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/20110621/d0e69851/attachment-0003.html 