<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Dear Dai,<BR>
<BR>
If you want to plot contour using station data, you need to open a empty grid file and interpolate to the grid first. It seemed that you missed something.<BR>
<BR>
Here are a sample<BR>
<BR>
'open grid.ctl'<BR>'open stn.ctl'<BR>
........<BR>
<BR>
'define a=oacres(g,prep.2,4)'<BR>
'd a'<BR>
<BR>
<BR>
Hope help.<BR>
<BR>
<BR>
<HR>
Date: Tue, 22 Jan 2008 18:00:00 -0700<BR>From: qdai@MATH.UALBERTA.CA<BR>Subject: Re: [GrADS] How to plot stations climatology contour in a map?<BR>To: GRADSUSR@LIST.CINECA.IT<BR><BR>
<META content="Microsoft SafeHTML" name=Generator><BR>
<BLOCKQUOTE>
<FONT size=2>Dear Jayakrishnan,</FONT><BR>
<FONT size=2>In order to plot the station climatology <STRONG>contour</STRONG>, I used the C++ program, ctl file and gs file similar to yours. But it is still not working. The attachment is my TXT file (station ID, latitude, longitude, climatology), format is "A7, F6.2 F8.2, F8.1". So could you please use my TXT data file and your C++ program, ctl file and gs file to test again?? Thanks a lot.</FONT><BR>
<FONT size=2>Best,</FONT><BR>
<FONT size=2>Qingfang</FONT><BR>
<FONT size=2></FONT><BR>
<FONT size=2><B>On Tue, 22 Jan 2008 12:21:51 +0530, Jayakrishnan PR wrote</B> <BR>> Dear Qingfang Dai, <BR>> I have followed the same method as given by you in plotting station data and was successful in plotting the shaded plots of the required value. I herewith attaches the ctl file which I used and the c program. The program looks almost same with some modifications. May be you can compare it with your purposes because It is working well for me. <BR>> <BR>> I have two suggestions....One is that it is better to create station map from out side the gs file and incorperate it in the ctl file. Next suggestion is that before running in the gs file kindly try the procedures in the command terminal. That will be better to know where the error has occured. <BR>> Hope this helps......best wishes.... <BR>> <BR>> <BR>> <BR>> <BR>> <SPAN class=EC_gmail_quote>On 1/22/08, <B class=EC_gmail_sendername>Qingfang D
ai</B> <<A href="mailto:qdai@math.ualberta.ca">qdai@math.ualberta.ca</A>> wrote:</SPAN> </FONT><BR>
<BLOCKQUOTE class=EC_gmail_quote style="PADDING-LEFT: 1ex; BORDER-LEFT: #ccc 1px solid"><BR>> <FONT size=2><BR>> <FONT face=verdana,arial,helvetica,sans-serif size=2>Dear GrADS Users,</FONT> <BR>> <FONT face=verdana,arial,helvetica,sans-serif size=2>I want to plot the stations' climatology contour. My TXT data file is the format of (stnid, lat, lon, climatology value). So I did as followings but it doesn't work. I don't how I shuld do. </FONT><BR>> <FONT face=verdana,arial,helvetica,sans-serif size=2><STRONG>First, I used the following C++ program to write the TXT file into DAT file:</STRONG></FONT> <BR>> #include "stdafx.h" <BR>> #include <stdio.h> <BR>> int main(int argc, char* argv[]) <BR>> { <BR>> /* Structure that describes a report header in a stn file */ <BR>> struct rpthdr <BR>> { <BR>> char id[8]; /* Station ID */ <BR>> float lat; 
; /* Latitude of Station */ <BR>> float lon; /* Longitude of Station */ <BR>> float t; /* Time in grid-relative units */ <BR>> int nlev; /* Number of levels following */ <BR>> int flag; /* Level independent var set flag */ <BR>> } hdr; <BR>> FILE *ifile, *ofile; <BR>> char rec[80]; <BR>> //int flag,year,month,yrsav,mnsav; <BR>> int i; <BR>> float val; <BR>> /* Open files */ <BR>> ifile = fopen ("clim.txt","r"); <BR>> ofile = fopen ("clim.dat","wb"); <BR>> if (ifile==NULL || ofile==NULL) { <BR>> printf("Error opening files\n"); <BR>> r
eturn 1; <BR>> } <BR>> /* Read, write loop */ <BR>> val = 0.0; <BR>> while (fgets(rec,79,ifile)!=NULL) { <BR>> /* Format conversion */ <BR>> sscanf (rec+8," %g %g %g",&hdr.lat,&hdr.lon,&val); <BR>> for (i=0; i<8; i++) hdr.id[i] = rec[i]; <BR>> /* Write this report */ <BR>> hdr.nlev = 1; <BR>> hdr.flag = 1; <BR>> hdr.t = 0.0; <BR>> fwrite (&hdr,sizeof(struct rpthdr), 1, ofile); <BR>> fwrite (&val,sizeof(float), 1, ofile); <BR>> } <BR>> /* Time group terminator */ <BR>> hdr.nlev = 0; <BR>> fwrite (&hdr,sizeof(struct rpthdr), 1, ofile); <BR>> fclose(ifile); <BR>> fclose(ofile); <BR>> return 0; <BR>> } <BR>> <BR>> <STRO
NG>Second, I write the ctl file 'clim.ctl' as</STRONG> <BR>> DSET ^clim.dat <BR>> DTYPE station <BR>> STNMAP ^clim.map <BR>> UNDEF -999.9 <BR>> TITLE 61-90 climatology <BR>> TDEF 1 LINEAR Jan1960 1mo <BR>> VARS 1 <BR>> prep 1 99 precipitation climatology of station <BR>> ENDVARS <BR>> <STRONG>Finally, I write the gs file '<A href="http://clim.gs/" target=_blank>clim.gs</A>' as</STRONG> <BR>> '!stnmap -i ../stnmap/AB_elev.ctl' <BR>> 'reinit' <BR>> 'open ../stnmap/AB_elev.ctl' <BR>> 'set gxout contour' <BR>> 'set grads off' <BR>> 'set mproj nps' <BR>> 'set grid off' <BR>> 'set frame off' <BR>> 'set cstyle 3' <BR>> 'set csmooth off' <BR>> 'set t 1 1' <BR>> 'd prep' <BR>> 'set font 1' <BR>> 'set string 1 tc 3' <BR>> 'set strsiz 0.15 0.15' <BR>> 'draw title 1961-1990 Climatology in January' <BR>> 'printim ../Clim.gif gif w
hite' <BR>> Thanks. <BR>> Qingfang <BR>> </FONT></BLOCKQUOTE><BR>> <BR clear=all>Sincerely <BR>> -- <BR>> *********************************************** <BR>> Jayakrishnan.P.R <BR>> Junior Research Fellow <BR>> Department of Atmospheric Sciences <BR>> Cochin University of Science And Technology (CUSAT), Cochin-16 <BR>> Kerala, India. <BR>> Mob: 09895417565 <BR>> <BR>> *********************************************** <BR><BR><BR>-- <BR>Department of Math and Stat Sciences (<A href="http://www.math.ualberta.ca/" target=_blank>http://www.math.ualberta.ca</A>) <BR><BR></BLOCKQUOTE><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>