<div dir="ltr">Hi Antonella,<br><div><br></div><div>The script I sent you earlier reads coordinates of the trajectory from text files. An example of such a text file is attached with this mail (longitude in first column and latitude in second column). The script uses functions subline and subword to read data from the text files. Such functions are not required read data when the data is in binary station data format. </div>
<div><br></div>
<div>The script I sent you earlier will not require extensive modification if the trajectory data can be written in text file similar to the attached one (you may have to add a third column for time). I can help you If you find it very hard to tweak the script.</div>
<div><br></div><div style>regards</div><div style>martin</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 14 November 2013 19:45, <span dir="ltr"><<a href="mailto:antonia.dimaio@cnr.it" target="_blank">antonia.dimaio@cnr.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi<br>
I d0n' know whi but the command 'printim' doesn't work. Can be a lib problem?<br>
Anyway, I tried to mofify a script sent to me by Martin to plot track<br>
cyclone and so I could produce a image output.eps.<br>
<br>
The question now is that this layout is just the final image of the<br>
tracking while I want produce an image for eache timestep.<br>
Martin open a netcdf file while I have a ctl file (map and dat,<br>
too). How can I modify the following part of the script to read my ctl<br>
file? I don't know how use subline and subwrd functions to find lat<br>
and lon using a ctl file<br>
<br>
Thank you very much!!<br>
Antonella<br>
<br>
*===========================================================================<br>
* 9) Plotting cyclone tracks<br>
*===========================================================================<br>
dta1=trkNo<br>
dta=0<br>
while(dta1>0)<br>
dta=dta+1<br>
*dtaFile=InFile2Prt%trkFile.dta<br>
dtaFile=datum1.ctl<br>
say 'reading from file 'dtaFile<br>
res=read(dtaFile) ;* skipping first line of cyclone track data file<br>
say 'title line 'res<br>
rows=trk1.dta<br>
res=read (dtaFile);* Reading first lat lon points<br>
anslin2=sublin(res,2)<br>
say 'second row ' res<br>
lonpt1=subwrd (anslin2,1)<br>
latpt1=subwrd (anslin2,2)<br>
'q w2xy 'lonpt1' 'latpt1<br>
<br>
<br>
*'q w2xy 100 28' will give result as 'X = 6.22321 Y = 8.08214'<br>
* 'x1 = subwrd(result,3)', will take '6.22321', which is in 3rd<br>
position in 'result'<br>
* and save it in 'x1'<br>
x1=subwrd(result,3)<br>
y1=subwrd(result,6)<br>
<br>
while(rows>1)<br>
<br>
res=read (dtaFile)<br>
anslin2=sublin(res,2)<br>
lonpt2=subwrd (anslin2,1)<br>
latpt2=subwrd (anslin2,2)<br>
<br>
'q w2xy 'lonpt2' 'latpt2<br>
x2=subwrd(result,3)<br>
y2=subwrd(result,6)<br>
<br>
'set line 'LnClr' 'LnStyl;*set line color# <style> <thickness><br>
'draw line 'x1' 'y1' 'x2' 'y2<br>
<br>
<br>
My ctl file is following (I want plot the variabel 'age'):<br>
<br>
dset datum1.dat<br>
dtype station<br>
stnmap datum1.map<br>
undef -999.0<br>
title PIW-6<br>
tdef 85 linear 23:57z11jul2013 30mn<br>
vars 3<br>
sta 0 99 State<br>
age 0 99 Age<br>
ort 0 99 Orientation<br>
endvars<br>
<br>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<br>
<br>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Martin V. Mathew<br>Research scholar<br>Centre for Oceans, Rivers, Atmosphere and Land Sciences (CORAL)<br>IIT Kharagpur<br>India.<br>
</div>