<!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.2800.1609" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2>Marco,</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Time 1 must have station count=4 not
station</FONT></DIV>
<DIV><FONT face=Arial size=2>count=1 as you can see when you run
stnmap.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Check how you wrote your binary .dat
file.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>It hope it helps.</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Bye</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Daniele</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<BLOCKQUOTE dir=ltr
style="PADDING-RIGHT: 0px; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #000000 2px solid; MARGIN-RIGHT: 0px">
<DIV style="FONT: 10pt arial">----- Original Message ----- </DIV>
<DIV
style="BACKGROUND: #e4e4e4; FONT: 10pt arial; font-color: black"><B>From:</B>
<A title=sartore@ELBATECH.COM href="mailto:sartore@ELBATECH.COM">Marco
Sartore</A> </DIV>
<DIV style="FONT: 10pt arial"><B>To:</B> <A title=GRADSUSR@LIST.CINECA.IT
href="mailto:GRADSUSR@LIST.CINECA.IT">GRADSUSR@LIST.CINECA.IT</A> </DIV>
<DIV style="FONT: 10pt arial"><B>Sent:</B> Tuesday, April 29, 2008 8:50
AM</DIV>
<DIV style="FONT: 10pt arial"><B>Subject:</B> Again about plotting rain.ctl
example...</DIV>
<DIV><FONT face=Arial size=2></FONT><FONT face=Arial size=2></FONT><FONT
face=Arial size=2></FONT><FONT face=Arial size=2></FONT><FONT face=Arial
size=2></FONT><FONT face=Arial size=2></FONT><FONT face=Arial
size=2></FONT><BR></DIV>
<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><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Starting scan of station data binary
file.<BR>Binary data file open: rain.dat<BR><STRONG>Processing time =
1<BR> Time = 1 has stn count = 1
<BR></STRONG>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><FONT face=Arial size=2></FONT> </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></BLOCKQUOTE></FONT></BODY></HTML>