Again about plotting rain.ctl example...

[Daniele Gandini]

Marco,

Time 1 must have station count=4 not station
count=1 as you can see when you run stnmap.

Check how you wrote your binary .dat file.

It hope it helps.

Bye

Daniele



  ----- Original Message ----- 
  From: Marco Sartore 
  To: GRADSUSR at LIST.CINECA.IT 
  Sent: Tuesday, April 29, 2008 8:50 AM
  Subject: Again about plotting rain.ctl example...


  Dear Users,
  I have followed the indications received to try plotting the Station data example "rain.dat" created using the C code available in the Grads documentation starting from the following .txt file:

  Year    Month     Stid     Lat     Lon     Rainfall
  1980     1     QQQ     34.3     -85.5     123.3
  1980     1     RRR     44.2     -84.5     87.1
  1980     1     SSS     22.4     -83.5     412.8
  1980     1     TTT     33.4     -82.5     23.3
  1980     2     QQQ     34.3     -85.5     145.1
  1980     2     RRR     44.2     -84.5     871.4
  1980     2     SSS     22.4     -83.5     223.1
  1980     2     TTT     33.4     -82.5     45.5


  I first have created the .map file from inside Grads, getting the following result:

  ga-> !stnmap -i rain.ctl
    Name of binary data set: rain.dat
    Number of times in the data set: 2
    Number of surface variables: 1
    Number of level dependent variables: 0

  Starting scan of station data binary file.
  Binary data file open: rain.dat
  Processing time = 1
       Time = 1 has stn count = 1 
  Processing time = 2
       Time = 2 has stn count = 4 
    Max reports per time:  4 reports at t = 2
    Max data elements in largest report: 1

  Version 2 Station map file created: rain.map

  Then I tried to display either a number or a mark at the station locations, with "set gxout value" and "set gxout stnmark"  without success.
  When I try to print out the station list I get the following:

  ga-> set stat on
  ga-> set stnprint on
  ga-> d p
  Data Type = station
  Dimensions = 0 1
  I Dimension = 0 to 360
  J Dimension = -90 to 90
  Stn count = 1
  Undef value = -999
  Undef count = 0  Valid count = 1 
  Min, Max = 3.89155e-41 3.89155e-41
  Cmin, cmax, cint = -5 5 1
  Stats[sum,sumsqr,root(sumsqr),n]:     3.89155e-41 0 0 0
  Stats[(sum,sumsqr,root(sumsqr))/n)]:     3.89155e-41 0 0
  Stats[(sum,sumsqr,root(sumsqr))/(n-1))]: 3.89155e-41 0 0
  Stats[(sigma,var)(n)]:     0 0
  Stats[(sigma,var)(n-1)]:   0 0
  Printing station values:  #obs = 1
  OB    ID       LON      LAT      LEV      VAL
  1     Stid      La  0        0        -9.99e+33 3.89155e-41


  Can anybody please explain what is going wrong and why I cannot see anything else than the world map ?

  Quite tricky and difficult to me this Grads, I wonder if I will be never able to plot local shaded maps of SPI data as I have to do in my final goal :-(((


  Thanks for your help, regards
  Marco


  _______________________________________________
  Marco Sartore - ElbaTech srl - Via Roma, 10
  57030 Marciana (LI) - Italy - Tel/Fax ++39-0565-90.10.02
  www.elbatech.com - www.spmagic.com 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080429/2afb8f8e/attachment.html