关于fwrite命令

hanyunxiang

Fwrite命令疑问为了从某日某时我国的卫星观测数据（.AWX文件，文件全名为“FY2D_AMV_IR1_OTG_20120317_0300.AWX”）中提取风的三维矢量信息和温度信息：
（1）用Grads对AWX配置好相关的ctl文件，文件命名为“sat1.ctl”，内容如下
附ctl文档：
dset d:/qx/FY2D_AMV_IR1_OTG_20120317_0300.AWX
title 20120317 fy2d_wind
undef 0.10000E+16
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 1 levels 250
tdef 1 linear 03z17MAR2012 6hr
vars 1
u    1   99 Eastward wind [m\s]
endvars
该ctl文档可以正常运行

（2）在脚本gs中用fwrite命令转化为grd格式，gs文件命名为“wf1.gs”
附gs文档--以提取u分量为例：
'reinit'
****************** u ***********************
'open d:\qx\sat1.ctl'
'set gxout fwrite'
'set fwrite d:\qx\wind_u.grd'
tt=1
while(tt>=1)
'set t 'tt''
zz=1
while(zz>=1)
'set z 'zz''
'set x 1 50'
'set y 1 50'
'd u'
zz=zz-1
endwhile
tt=tt-1
endwhile
'disable fwrite'
'reinit'
*******************************************

此时能正常运行，生成grd格式的文件
但将'set x 1 50'
'set y 1 50'
改为
'set x 1 500'
'set y 1 500'
就又提示错误了，什么原因呢？？谢谢

兰溪之水

刚刚发错了。。。

兰溪之水

具体错误贴出来。。。。。。。。。。。

hanyunxiang

错误的提示截图如下,thank you!!!!!!!!

兰溪之水

'reinit'
****************** u ***********************
'open d:\qx\sat1.ctl'
'set gxout fwrite'
'set fwrite d:\qx\wind_u.grd'
'set t 1'
'set z 1'
'set x 1 50'
'set y 1 50'
'd u'
'disable fwrite'
'reinit'
*******************************************

hanyunxiang

按照你的建议进行了修改 还是有错误（和贴图提示的错误一样）

兰溪之水

hanyunxiang 发表于 2012-4-10 10:32

                               
登录/注册后可看大图

按照你的建议进行了修改 还是有错误（和贴图提示的错误一样）

图显示你读数据进来就错了，检查一下你的ctl是不是能用的，画图试试~

hanyunxiang

嗯 我运行了一下 提示如下图示 不能正常 “d u”      难道我的ctl文件有问题 但是它能正常打开啊

兰溪之水

hanyunxiang 发表于 2012-4-10 10:42

                               
登录/注册后可看大图

嗯 我运行了一下 提示如下图示 不能正常 “d u”      难道我的ctl文件有问题 但是它能正常打开啊

数据或ctl有问题，自个检查~

fish

2.0版本fwrite好像是有问题的，我写wrfout数据都用1.8，2.0写出来的总是好几个字节~也不晓得是怎么回事~