关于fwrite命令

hanyunxiang

为了从某日某时我国的卫星观测数据（.AWX文件，文件全名为“FY2D_AMV_IR1_OTG_20120317_0300.AWX”）中提取风的三维矢量信息和温度信息：
（1）用Grads对AWX配置好相关的ctl文件，文件命名为“sat1.ctl”，内容如下
dset d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/FY2D_AMV_IR1_OTG_20120317_0300.AWX
title 20120317 fy2d_wind
undef 32766
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 7 levels 1000 600 500 400 300 200 100
tdef 1 linear 03z17MAR2012 6hr
vars  4
u         7   99 Eastward wind [m/s]
v         7   99 Northward wind [m/s]
w         7   99 Vertical wind [m/s]
T         7   99 Air Temperature [K]
endvars
该ctl文档能正常运行！
（2）在脚本gs中用fwrite命令转化为grd格式，gs文件命名为“wf1.gs”
'reinit'
****************** u分量 ***********************
'open d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/sat1.ctl'
'set fwrite d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/wind_u.grd'
'set gxout fwrite'
tt=1
while(tt=1)
'set t 'tt''
zz=7
while(zz>=1)
'set z 'zz''
'set lon 45 165'
'set lat -60 60'
'd u'
zz=zz-1
endwhile
tt=tt-1
endwhile
'disable fwrite'
'reinit'
*******************************************
***************** v分量 ************************
'open d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/sat1.ctl'
'set fwrite d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/wind_V.grd'
'set gxout fwrite'
tt=1
while(tt=1)
'set t 'tt''
zz=7
while(zz>=1)
'set z 'zz''
'set lon 45 165'
'set lat -60 60'
'd v'
zz=zz-1
endwhile
tt=tt-1
endwhile
'disable fwrite'
'reinit'
*******************************************
***************** w分量 ************************
'open d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/sat1.ctl'
'set fwrite d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/wind_W.grd'
'set gxout fwrite'
tt=1
while(tt=1)
'set t 'tt''
zz=7
while(zz>=1)
'set z 'zz''
'set lon 45 165'
'set lat -60 60'
'd w'
zz=zz-1
endwhile
tt=tt-1
endwhile
'disable fwrite'
'reinit'
*******************************************
***************** T分量 ************************
'open d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/sat1.ctl'
'set fwrite d:/Program Files/grads-2.0/OpenGrADS/Contents/Resources/SampleDatasets/wind_T.grd'
'set gxout fwrite'
tt=1
while(tt=1)
'set t 'tt''
zz=7
while(zz>=1)
'set z 'zz''
'set lon 45 165'
'set lat -60 60'
'd T'
zz=zz-1
endwhile
tt=tt-1
endwhile
'disable fwrite'
*******************************************
启动 grads软件，输入
run wf1.gs
出现以下结果：






请问是哪里出问题了？？谢谢大家

兰溪之水

路径不要有空格

兰溪之水

路径不要有空格

godenflame135

你确定你的ctl能够被正确打开？

hanyunxiang

godenflame135 发表于 2012-4-7 21:13

                               
登录/注册后可看大图

你确定你的ctl能够被正确打开？

ctl能够被正确打开   

纳川

gs文件中open语句的目录符应该用“\”，而不是“/”。

godenflame135

纳川 发表于 2012-4-7 22:35

                               
登录/注册后可看大图

gs文件中open语句的目录符应该用“\”，而不是“/”。

对头，考虑win系统下斜杠与反斜杠的区别。

hanyunxiang

按照楼上各位朋友的指导 进行了以下两方面的更改，
(1)在d盘新建了文件夹-名字为“qx”，来存储相关文档以消除“program files”的空格影响；
（2）更改“/”为“\”

还是不行啊  


附ctl文档：
dset d:\qx\FY2D_AMV_IR1_OTG_20120317_0300.AWX
title 20120317 fy2d_wind
undef 32766
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 7 levels 700 600 500 400 300 200 100
tdef 1 linear 03z17MAR2012 6hr
vars  4
u         7   99 Eastward wind [m\s]
v         7   99 Northward wind [m\s]
w         7   99 Vertical wind [m\s]
T         7   99 Air Temperature [K]
endvars
附gs文档--以提取u分量为例：
'reinit'
****************** u ***********************
'open d:\qx\sat1.ctl'
'set fwrite d:\qx\wind_u.grd'
'set gxout fwrite'
tt=1
while(tt=1)
'set t 'tt''
zz=7
while(zz>=1)
'set z 'zz''
'set lon 45 165'
'set lat -60 60'
'd u'
zz=zz-1
endwhile
tt=tt-1
endwhile
'disable fwrite'
'reinit'
*******************************************

圣诞2011

我有个问题：你的AWX文件的描述文件时怎么来的？是grads自己生成的，还是你自己写的呢？

冬日情愫+

我自己试了一下，你这个ctl 不能正常运行呀？