<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=GB2312" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
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>
</body>
</html>