<div dir="ltr"><div>dear all</div>
<div> </div>
<div>I want to extract point data from a netcdf file for 9 stations ie for 9 different latitudes longitudes .</div>
<div>my script is like this:</div>
<div>'reinit'</div>
<div>'sdfopen <a href="http://air.2009.nc">air.2009.nc</a>'<br>n=1<br>t=1<br>ln1=74.9<br>ln2=76.03<br>ln3=68.9<br>ln4=73.25<br>ln5=71.78<br>ln6=71.68<br>ln7=73.35<br>ln8=70.1<br>ln9=69.72<br>lt1=35.34<br>lt2=35.12<br>
lt3=24.63<br>lt4=29.95<br>lt5=29.4<br>lt6=29.38<br>lt7=34.38<br>lt8=33<br>lt9=29.88<br>while (n<=9)<br>'set lat lt'n<br>'set lon ln'n<br>while(t<3)<br>'set t 't<br>'q time'<br>say result<br>
res=subwrd(result,3)<br>dd=substr(res,4,2)<br>mm=substr(res,6,3)<br>yy=substr(res,9,4)<br>'set gxout value'<br>* Temperature </div>
<div>'d air-273.15'</div>
<div>say result<br>tmp=subwrd(result,4)<br>aph=write('output/lahore.csv',' 'yy' 'mm' 'tmp,append)<br>'c'</div>
<div>t=t+1<br>endwhile<br>n=n+1<br>endwhile</div>
<div> </div>
<div> </div>
<div> </div>
<div>but it results in error that </div>
<div>set error lat /lon missing or invalid arguments.</div>
<div> </div>
<div>regards</div>
<div><br><br clear="all"><br>-- <br>Waheed Iqbal<br><br></div></div>