<div dir="ltr">I share this example, I hope it serves you, I use it to extract the forecast data of rain and temperatures.<div><br></div><div><br><div><div>'open wrf_00z.ctl'</div><div><br></div><div>* imports a file with the location data of the stations</div><div>input='/home/meteorology/Documents/scripts/stations_pptmaxmin.txt'</div><div><br></div><div>'set t 2'</div><div>'query time'</div><div>say result</div><div>cdate=subwrd(result,3)</div><div>date=substr(cdate,4,12)</div><div><br></div><div>cont=1</div><div>while(cont<=3)</div><div>rec=read(input)</div><div>plat=subwrd(rec,3)</div><div>plon=subwrd(rec,2)</div><div>name_est=subwrd(rec,4)</div><div><br></div><div>'set lat 'plat</div><div>'set lon 'plon</div><div><br></div><div>'d apcpsfc(t=37)-apcpsfc(t=13)'</div><div>respuesta=result</div><div>lluvias=sublin(respuesta,1)</div><div>dpp1=subwrd(lluvias,4)</div><div><br></div><div>'d max(tmax2m-273.15,t=14,t=37)'<br></div><div>respuesta=result</div><div>tmaxs=sublin(respuesta,2)</div><div>tmax1=subwrd(tmaxs,4)</div><div><br></div><div>'d min(tmin2m-273.15,t=14,t=37)'</div><div>respuesta=result</div><div>tmins=sublin(respuesta,2)</div><div>tmin1=subwrd(tmins,4)</div><div><br></div><div>* Write a file</div><div>res=write('/home/meteorology/Documents/data/data_'date'.txt', plon-360'<span style="white-space:pre"> </span>'plat'<span style="white-space:pre"> </span>'name_est'<span style="white-space:pre"> </span>'dpp1'<span style="white-space:pre"> </span>'tmax1'<span style="white-space:pre"> </span>'tmin1)</div><div><br></div><div>cont=cont+1</div><div>endwhile</div></div></div><div><br></div><div><br></div><div><br></div><div>format of file stations</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-08-13 12:22 GMT-06:00 Charles Langton Vanya <span dir="ltr"><<a href="mailto:dolesibeni@gmail.com" target="_blank">dolesibeni@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear,<br>I believe you are meaning extracting data from netcdf file. <br>The best way I do is to average data from a grid box of the area.<br> The reason for doing this is to make sure that your data is a <br>representation of the area you are interested inother than just<br> a single a single point. <br>This can be in so many ways. But the simolest I can remember is like this:<br><br>*****THIS SCRIPT EXTRACT SEA SURFACE TEMPERATURES FOR GRID BOX<br>**LONG 129E to 131E and LAT 20.0N to 24.0N<br><br>'reinit'<br>'sdfopen C:\fortran\sst\<a href="http://2018sst.nc" target="_blank">2018sst.nc</a>'<br>'set x 1'<br>'set y 1'<br>** set your area1 here..If your data has high resoultion you can set your area too closer<br>lon1 =129.0<br>lon2 =131.0<br>lat1 =20.0<br>lat2 =24.0<br><br>** Keep time fixed as below.. Here time is taken as starting from value 1 to 703<br>** you can also set your time by exactly setting 'set t hhddmmyy'<br><br>'set t 1 703'<br>** Here we are finding areal average of the grid.<br><br>'define ssst = aave(sst,lon = 'lon1',lon = 'lon2', lat = 'lat1',lat='lat2')'<br><br>** set file name to pring oout... Here we are saving data as 12 month data table<br>'fprintf ssst C:\fortran\SST_EXTRACT\DJF\<wbr>REG2\REG2\R2DJF_SJPAN.txt %g 12 1'<br><br>The data is save as a text file and in a single column. If you you want <br>the data to be extracted as monthly then you use 1 12 instead of 1 1 meaning <br>that you will have 12 columns representing months of the year.<br>Hope this helps.<br>Thank<br>Regards,<br>Charles Vanya<br></div><br><div class="gmail_quote"><div dir="ltr">On Mon, 13 Aug 2018 at 15:12, Davide Sacchetti <<a href="mailto:davide.sacchetti@arpal.gov.it" target="_blank">davide.sacchetti@arpal.gov.it</a><wbr>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>is file in station data format? if yes use gxout print<br><br><br><div><strong>
From:
</strong>
E Mobarak <<a href="mailto:mobarak_e@yahoo.com" target="_blank">mobarak_e@yahoo.com</a>>
<br>
<strong>
To:
</strong>
<<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>>
<br>
<strong>
Sent:
</strong>
8/11/2018 6:09 PM
<br>
<strong>
Subject:
</strong>
[gradsusr] extract station data
<br><br><blockquote class="m_-8373820826899355839m_-4609455438561977606mori" style="margin:0 0 0 .8ex;border-left:1px solid #ccc;padding-left:1ex"><div style="font-family:times new roman,new york,times,serif;font-size:16px"><div>Hi </div><div>how i can extract data for some station with different latitude and longitude, by loop or function in scripts</div><div>thanks</div><div style="font-family:times new roman,new york,times,serif;font-size:16px"></div></div><br><br>______________________________<wbr>_________________
<br>gradsusr mailing list
<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>
<br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a>
<br></blockquote></div></div>______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><span class="HOEnZb"><font color="#888888"><br>
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><br clear="all"><br>-- <br><div dir="ltr" class="m_-8373820826899355839gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div>Meteorological Services<br>P.O.Box 1808,<br>Blantyre,<br>Malawi<br>Cell:+265-888-980 545</div>
<div> 00265111939610<br> <br></div></div></div></div>
</font></span><br>______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Atte. <div>Walter Arnoldo Bardales Espinoza</div><div><br></div><div> <img src="https://lh5.googleusercontent.com/-54z6MDt7TZs/URJGTsrgxII/AAAAAAAAAb0/TQVdixKgDtA/s512/Logo%2520de%2520la%2520empresa4.png" width="181" height="200"></div></div></div></div>
</div>