<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.2900.3132" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2>Dear Users,</FONT></DIV>
<DIV><FONT face=Arial size=2>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:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Year
Month Stid
Lat Lon
Rainfall<BR>1980 1
QQQ 34.3
-85.5 123.3<BR>1980
1 RRR
44.2 -84.5
87.1<BR>1980 1
SSS 22.4
-83.5 412.8<BR>1980
1 TTT
33.4 -82.5
23.3<BR>1980 2
QQQ 34.3
-85.5 145.1<BR>1980
2 RRR
44.2 -84.5
871.4<BR>1980 2
SSS 22.4
-83.5 223.1<BR>1980
2 TTT
33.4 -82.5 45.5<BR></FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>I first have created the .map file from inside
Grads, getting the following result:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>ga-> !stnmap -i rain.ctl<BR> Name of
binary data set: rain.dat<BR> Number of times in the data set: 2<BR>
Number of surface variables: 1<BR> Number of level dependent variables:
0</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial size=2>Starting scan of station data binary
file.<BR>Binary data file open: rain.dat<BR>Processing time =
1<BR> Time = 1 has stn count = 1 <BR>Processing time =
2<BR> Time = 2 has stn count = 4 <BR> Max reports
per time: 4 reports at t = 2<BR> Max data elements in largest
report: 1</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial size=2>Version 2 Station map file created:
rain.map<BR></FONT></DIV>
<DIV><FONT face=Arial size=2>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.</FONT></DIV>
<DIV><FONT face=Arial size=2>When I try to print out the station list I get the
following:</FONT></DIV>
<DIV><FONT face=Arial size=2><BR>ga-> set stat on<BR>ga-> set stnprint
on<BR>ga-> d p<BR>Data Type = station<BR>Dimensions = 0 1<BR>I Dimension = 0
to 360<BR>J Dimension = -90 to 90<BR>Stn count = 1<BR>Undef value =
-999<BR>Undef count = 0 Valid count = 1 <BR>Min, Max = 3.89155e-41
3.89155e-41<BR>Cmin, cmax, cint = -5 5
1<BR>Stats[sum,sumsqr,root(sumsqr),n]: 3.89155e-41 0 0
0<BR>Stats[(sum,sumsqr,root(sumsqr))/n)]: 3.89155e-41 0
0<BR>Stats[(sum,sumsqr,root(sumsqr))/(n-1))]: 3.89155e-41 0
0<BR>Stats[(sigma,var)(n)]: 0
0<BR>Stats[(sigma,var)(n-1)]: 0 0<BR>Printing station values:
#obs = 1<BR>OB ID
LON LAT
LEV VAL<BR>1
Stid La
0
0 -9.99e+33
3.89155e-41<BR></FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Can anybody please explain what is going wrong and
why I cannot see anything else than the world map ?</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>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 :-(((</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Thanks for your help, regards</FONT></DIV>
<DIV><FONT face=Arial size=2>Marco</FONT></DIV>
<DIV><FONT face=Arial size=2> </DIV>
<DIV><BR>_______________________________________________<BR>Marco Sartore -
ElbaTech srl - Via Roma, 10<BR>57030 Marciana (LI) – Italy - Tel/Fax
++39-0565-90.10.02<BR><A href="http://www.elbatech.com">www.elbatech.com</A> -
<A href="http://www.spmagic.com">www.spmagic.com</A>
<BR></DIV></FONT></BODY></HTML>