[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