<div>Hi, </div>
<div>I have used "oacres", it works very well. One thing you can try is, open the grid file (<a href="http://cal.controla.1dom.nc">cal.controla.1dom.nc</a>) firstly, then open the station file (a.ctl) (don't forget to change the OACRES command!).
</div>
<div>Good Luck</div>
<div>Qiong<br><br> </div>
<div><span class="gmail_quote">On 11/21/05, <b class="gmail_sendername">Hui Du</b> <<a href="mailto:huidu@ucdavis.edu">huidu@ucdavis.edu</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hello, grads users<br> I am trying to convert the station data to grid data using the function 'oacres'.<br>
I created a station file for one variable. <br> I also successfully got the map file using 'stnmap'<br> Then I opened two files to try to get the grid data, the following are my command and error information.<br>*************************************************************************************
<br> ga-> open a.ctl<br>Scanning description file: a.ctl<br>Data file controla.dat is open as file 1<br>LON set to 0 360<br>LAT set to -90 90<br>LEV set to 500 500<br>Time values set: 1995:8:1:0 1995:8:1:0<br>ga-> sdfopen ../cal.controla.1dom.nc
<br>Scanning self-describing file: ../cal.controla.1dom.nc<br>SDF file ../cal.controla.1dom.nc is open as file 2<br>ga-> set lon 230 250<br>LON set to 230 250<br>ga-> set lat 30 45<br>LAT set to 30 45<br>ga-> set t 100
<br>Time values set: 1995:8:25:18 1995:8:25:18<br>ga-> d oacres(t2.2,t2)<br><font color="#ff0000">Segmentation fault (core dumped)</font><br>*******************************************************************<br>then grads broke, I have no idea what made this happen.
<br><br>For the a.clt, its script is the following<br>***********************************************<br>dset controla.dat<br>title t2<br>dtype station<br>stnmap t2.map<br>undef -999<br>tdef 1697 linear 00z1aug1995 6hr<br>
vars 1<br>t2 0 99 t2 data<br>endvars<br>*****************************************************<br>I don't think it is a good idea to convert my station data because with only one surface variable t2, the size of station data is about 580M.
<br><br>Has anyone have any program to do this outside grads? Or can you tell me where the source code of the function 'oacres' is<br><br>Thanks a lot <br><br></blockquote></div><br>