No subject
Wed Apr 14 09:33:43 EDT 2010
e
not linear scaling. When I used "oacres", I generated a grid data file and
ctl file, like
DSET emefs_fc_grid.dat
UNDEF -999.9
XDEF 30 LINEAR -100 1.0
YDEF 24 LINEAR 30 1.0
ZDEF 1 levels 1000
TDEF 1 LINEAR 12Z01aug1988 1dy
It is not important what kind of data in the data file (can be the
constant), but the dimensions should be consistent with the ctl file.
If you still want to use "cal.controla.1dom.nc" as the grid file, check
(1) Are the grid settings are linear?
(2) Try to plot grid file and station file respectively (such as, "d t2" an=
d
"d t2.2") before using "oacres", to see whether both of files are opened
correctly.
Hope those can be somewhat helpful
Qiong
On 11/22/05, Hui Du <huidu at ucdavis.edu> wrote:
>
> Hi, Qiong
> I have tried your method and failed.
> I got different error message
> I wonder whether I can use 'sdfopen' for function 'oacres'
>
> *****************************************
> ga-> sdfopen ../cal.controla.1dom.nc
> Scanning self-describing file: ../cal.controla.1dom.nc
> SDF file ../cal.controla.1dom.nc is open as file 1
> LON set to 223.674 253.663
> LAT set to 28.4976 46.9737
> LEV set to 0 0
> Time values set: 1995:8:1:0 1995:8:1:0
> ga-> open a.ctl
> Scanning description file: a.ctl
> Data file controla.dat is open as file 2
> ga-> set lon 230 250
> LON set to 230 250
> ga-> set lat 30 45
> LAT set to 30 45
> ga-> d oacres(t2,t2.2)
> Error from OACRES: Grid must have linear scaling
> Operation Error: Error from oacres function
> Error ocurred at column 1
> DISPLAY error: Invalid expression
> Expression =3D oacres(t2,t2.2)
> **********************************************************
>
> I also created a 'ctl' file for 'cal.controla.1dom.nc' and got the
> following error message
>
> ****************************************************
> Warning from OACRES: Less than two stations
> Grid will be all missing values
> Cannot contour grid - all undefined values
>
> ****************************************************
> If you have use 'oacres' a lot, can you give me some hint
> Anyway, I still want to use fortran program to do this outside grads.
> Thanks a lot
>
------=_Part_9615_19198477.1132757866768
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
Content-Disposition: inline
<div>Hi,</div>
<div>From error message, it looks the grid settings in "<a href=3D"htt=
p://cal.controla.1dom.nc">cal.controla.1dom.nc</a>" are not linear sca=
ling. When I used "oacres", I generated a grid data file and=
ctl file, like
</div>
<div> </div>
<div>DSET emefs_fc_grid.dat<br>UNDEF -999.9<br>XDEF =
30 LINEAR -100 1.0<br>YDEF 24 LINEAR 30 1.0<=
br>ZDEF 1 levels 1000<br>TDEF 1 LINEAR 12Z01aug1988=
1dy<br> </div>
<div>It is not important what kind of data in the data file (can =
be the constant), but the dimensions should be consistent with the ctl=
file.</div>
<div> </div>
<div>If you still want to use "<a href=3D"http://cal.controla.1dom.nc"=
>cal.controla.1dom.nc</a>" as the grid file, check</div>
<div>(1) Are the grid settings are linear?</div>
<div>(2) Try to plot grid file and station file respectively (such as, &quo=
t;d t2" and "d t2.2") before using "oacres", to se=
e whether both of files are opened correctly.</div>
<div> </div>
<div>Hope those can be somewhat helpful</div>
<div>Qiong<br> </div>
<div><span class=3D"gmail_quote">On 11/22/05, <b class=3D"gmail_sendername"=
>Hui Du</b> <<a href=3D"mailto:huidu at ucdavis.edu">huidu at ucdavis.edu</a>&=
gt; wrote:</span>
<blockquote class=3D"gmail_quote" style=3D"PADDING-LEFT: 1ex; MARGIN: 0px 0=
px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi, Qiong<br>I have tried your m=
ethod and failed.<br>I got different error message<br>I wonder whether I ca=
n use 'sdfopen' for function 'oacres'
<br><br>*****************************************<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 1<br>LON set to 223.67=
4
253.663<br>LAT set to 28.4976 46.9737<br>LEV set to 0 0<br>Time values set=
: 1995:8:1:0 1995:8:1:0<br>ga-> open a.ctl<br>Scanning description file:=
a.ctl<br>Data file controla.dat is open as file 2<br>ga-> set lon 230 2=
50
<br>LON set to 230 250<br>ga-> set lat 30 45<br>LAT set to 30 45<br>ga-&=
gt; d oacres(t2,t2.2)<br>Error from OACRES: Grid must have linear scaling<b=
r>Operation Error: Error from oacres function<br>Error ocurred at column 1
<br>DISPLAY error: Invalid expression<br>Expression =3D oacres(t2,t2.2)<br>=
**********************************************************<br><br>I also cr=
eated a 'ctl' file for '<a href=3D"http://cal.controla.1dom.nc">cal.control=
a.1dom.nc
</a>' and got the<br>following error message<br><br>***********************=
*****************************<br>Warning from OACRES: Less than two station=
s<br>Grid will be all missing values<br>Cannot contour grid - all undefined=
values
<br><br>****************************************************<br>If you have=
use 'oacres' a lot, can you give me some hint<br>Anyway, I still want to u=
se fortran program to do this outside grads.<br>Thanks a lot<br></blockquot=
e>
</div><br>
------=_Part_9615_19198477.1132757866768--
More information about the gradsusr
mailing list