Dear Jayakrishnan,
In order to plot the station climatology contour, 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.
Best,
Qingfang
On Tue, 22 Jan 2008 12:21:51 +0530, Jayakrishnan PR wrote
> Dear Qingfang Dai,
> 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.
>
> 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.
> Hope this helps......best wishes....
>
>
>
>
> On 1/22/08, Qingfang Dai <qdai at math.ualberta.ca> wrote:
>
> Dear GrADS Users,
> 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.
> First, I used the following C++ program to write the TXT file into DAT file:
> #include "stdafx.h"
> #include <stdio.h>
> int main(int argc, char* argv[])
> {
> /* Structure that describes a report header in a stn file */
> struct rpthdr
> {
> char id[8]; /* Station ID */
> float lat; /* Latitude of Station */
> float lon; /* Longitude of Station */
> float t; /* Time in grid-relative units */
> int nlev; /* Number of levels following */
> int flag; /* Level independent var set flag */
> } hdr;
> FILE *ifile, *ofile;
> char rec[80];
> //int flag,year,month,yrsav,mnsav;
> int i;
> float val;
> /* Open files */
> ifile = fopen ("clim.txt","r");
> ofile = fopen ("clim.dat","wb");
> if (ifile==NULL || ofile==NULL) {
> printf("Error opening files\n");
> return 1;
> }
> /* Read, write loop */
> val = 0.0;
> while (fgets(rec,79,ifile)!=NULL) {
> /* Format conversion */
> sscanf (rec+8," %g %g %g",&hdr.lat,&hdr.lon,&val);
> for (i=0; i<8; i++) hdr.id[i] = rec[i];
> /* Write this report */
> hdr.nlev = 1;
> hdr.flag = 1;
> hdr.t = 0.0;
> fwrite (&hdr,sizeof(struct rpthdr), 1, ofile);
> fwrite (&val,sizeof(float), 1, ofile);
> }
> /* Time group terminator */
> hdr.nlev = 0;
> fwrite (&hdr,sizeof(struct rpthdr), 1, ofile);
> fclose(ifile);
> fclose(ofile);
> return 0;
> }
>
> Second, I write the ctl file 'clim.ctl' as
> DSET ^clim.dat
> DTYPE station
> STNMAP ^clim.map
> UNDEF -999.9
> TITLE 61-90 climatology
> TDEF 1 LINEAR Jan1960 1mo
> VARS 1
> prep 1 99 precipitation climatology of station
> ENDVARS
> Finally, I write the gs file 'clim.gs' as
> '!stnmap -i ../stnmap/AB_elev.ctl'
> 'reinit'
> 'open ../stnmap/AB_elev.ctl'
> 'set gxout contour'
> 'set grads off'
> 'set mproj nps'
> 'set grid off'
> 'set frame off'
> 'set cstyle 3'
> 'set csmooth off'
> 'set t 1 1'
> 'd prep'
> 'set font 1'
> 'set string 1 tc 3'
> 'set strsiz 0.15 0.15'
> 'draw title 1961-1990 Climatology in January'
> 'printim ../Clim.gif gif white'
> Thanks.
> Qingfang
>
> Sincerely
> --
> ***********************************************
> Jayakrishnan.P.R
> Junior Research Fellow
> Department of Atmospheric Sciences
> Cochin University of Science And Technology (CUSAT), Cochin-16
> Kerala, India.
> Mob: 09895417565
>
> ***********************************************
--
Department of Math and Stat Sciences (http://www.math.ualberta.ca)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080122/2ed6a181/attachment.html
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: clim.txt
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080122/2ed6a181/attachment.txt