[Fwd: another grads question]

Austin Conaty Austin.L.Conaty at NASA.GOV
Fri Oct 30 16:51:09 EDT 2009


Sim,
Here's the comments from our local station map "expert"...
Caveat: he's still using grads1

1 check your options statement
2 check "endian ness"
3 see if you can stream the data (station map is written in C
and would be happier reading streaming data)

Good luck,
Austin Conaty


   XXXXXX1   18.91675      -51.42435       0.000000               1           0
     500.0000       23.24000       91.25000      -25.13000       3.820000
   XXXXXX2   18.91410      -51.41340       0.000000               1           0
     1000.000       20.18000       91.60000      -25.86500       8.020000
   XXXXXX3   18.91055      -51.40285       0.000000               1           0
     1500.000       18.57000       79.40000      -25.35000       8.525000
   XXXXXX4   18.90680      -51.39260       0.000000               1           0
     2000.000       15.55500       82.45000      -24.31500       9.820000
   XXXXXX5   18.90230      -51.38280       0.000000               1           0
     2500.000       13.33000       82.55000      -25.01000       14.16000
   XXXXXX6   18.89670      -51.37290       0.000000               1           0
     3000.000       10.93500       80.70000      -25.08000       15.18000
   XXXXXX7   18.89175      -51.36465       0.000000               1           0
     3500.000       8.165000       68.85000      -23.95500       16.06000
   XXXXXX7   18.89175      -51.36465       0.000000               0           0

Sim Aberson wrote:
> I am having another frustrating GrADS issue with GrADS2 on my mac.  I am trying
> to make more station plots.  I have code to write out the station plot data.  I
> have done this before with no problem, but get an error now.  I am attaching a
> binary version of that file (fort.10) and a control file (fort.13).  When I run
> that through stnmap, I get an error:
>
> stnmap -v -i fort.13
>     Name of binary data set: fort.10
>     Number of times in the data set: 1
>     Number of surface variables: 0
>     Number of level dependent variables: 5
>
> Starting scan of station data binary file.
> Binary data file open: fort.10
>
> Processing time step 1
>     Invalid station hdr found in station binary file
>     Possible causes:  Invalid level count in hdr
>                       Descriptor file mismatch
>                       File not station data
>                       Invalid relative time
>                       Invalid sequential format
>       levs = 0  flag = 27  time = 1.4013e-45
>
> The contents of the file fort.13, in ascii, are:
>
>    XXXXXX1   18.91675      -51.42435       0.000000               1           0
>      500.0000       23.24000       91.25000      -25.13000       3.820000
>    XXXXXX2   18.91410      -51.41340       0.000000               1           0
>      1000.000       20.18000       91.60000      -25.86500       8.020000
>    XXXXXX3   18.91055      -51.40285       0.000000               1           0
>      1500.000       18.57000       79.40000      -25.35000       8.525000
>    XXXXXX4   18.90680      -51.39260       0.000000               1           0
>      2000.000       15.55500       82.45000      -24.31500       9.820000
>    XXXXXX5   18.90230      -51.38280       0.000000               1           0
>      2500.000       13.33000       82.55000      -25.01000       14.16000
>    XXXXXX6   18.89670      -51.37290       0.000000               1           0
>      3000.000       10.93500       80.70000      -25.08000       15.18000
>    XXXXXX7   18.89175      -51.36465       0.000000               1           0
>      3500.000       8.165000       68.85000      -23.95500       16.06000
>    XXXXXX7   18.89175      -51.36465       0.000000               0           0
>
>
> So, according to the first line, levs should be 1, flag should be zero, and time
> should be zero.  But the above from stnmap has very different values.  Am I
> missing something very fundamental from the documentation?  I'm stumped.
>
> Thanks,
> Sim
>
>
>
> --
> Sim Aberson         AOML/Hurricane Research Division          Miami, FL
>      "Hurricane states erect statues in my honor." -- David Sedaris
>
> ------------------------------------------------------------------------
>
> DSET ^fort.10
> DTYPE station
> STNMAP station_data.map
> OPTIONS sequential
> UNDEF -999.
> TITLE dropwindsondes
> TDEF 1 LINEAR 11:11Z01JAN2000 1mo
> VARS 5
> height    1  99 height
> tmphgt    1  99 temperature
> rhhgt     1  99 relative humidity
> ugrdhgt   1  99 u-component
> vgrdhgt   1  99 v-component
> ENDVARS
>
>


--
Austin Conaty, SAIC        Global Modeling & Assimilation Office
301-614-6149 (ph)          NASA Goddard Space Flight Center
301-614-6297 (fax)         Code 610.1 Greenbelt, MD 20771

Austin.L.Conaty at .nasa.gov

http://gmao.gsfc.nasa.gov



More information about the gradsusr mailing list