[gradsusr] stnmap error

[Huddleston, John]

Andre

There is another option which many scientists use; that is, to take your area of interest and apply an interpolation technique to grid the station data.

I just finished a two year study of Montana (USA) 30 year precipitation using: ordinary Kriging; inverse distance squared; regression on value and distance; multivariate inverse distance squared; regression on distance, value, and elevation; and universal Kriging on elevation, distance, and value.

You may want to think deeper about your studies on the best approach.


John Huddleston

Cooperative Institute for Research in the Atmosphere


________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On Behalf Of Andre Pattantyus [apattantyus2008 at my.fit.edu]
Sent: Monday, June 20, 2011 12:36 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] stnmap error


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<mailto: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> [mailto: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<mailto: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<tel:%28321%29%20674-8330> | Fax: (321) 674-7212<tel:%28321%29%20674-7212> | Email: apattantyus2008 at fit.edu<mailto:apattantyus2008 at fit.edu>

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto: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<mailto:apattantyus2008 at fit.edu>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110620/fa95a1e1/attachment-0003.html