sdfwrite函数 怎么在循环语句中使用？

leeol

原始数据：NCEP再分析资料：prate.sfc.day.1981-2010.ltm.nc

使用sdfwrite读取输出NC格式数据（GS如下）：
'reinit'
'sdfopen d:\data\prate.sfc.day.1981-2010.ltm.nc'
'set x 1 192'
'set y 1 94'
'set z 1'
*'set t 1'
'set sdfwrite d:\prate.sfc.pentad.1981-2010.ltm.nc'
day=1
while(day<=365)
'define pp=ave(prate,t='day',t='day'+4)'
'sdfwrite pp'
day=day+5
endwhile
;

grads运行没有报错，但是所输出文件的变量只有1个时次（经检验，为最后一个时次的资料）：file:///C:\Users\Edward Lee\AppData\Roaming\Tencent\Users\330493972\QQ\WinTemp\RichOle\UGQTHNVQ6CXI9RB(I_0ZP3D.jpg


我试着在GS中设置 set t 1 73 （GS中红字部分），结果出来了73个时次的资料，但是每个都相同，为最后一个时次资料。

PS：换用fwrite输出grd文件来经检验，这个用sdfwrite出来的数据是最后一个时次的资料。

求大家指导下，sdfwrite函数 怎么在循环语句中使用？

最新自我回复：

我今天又再次使用到了sdfwrite指令，说明页:http://www.iges.org/grads/gadoc/gradcomdsdfwrite.html 。为什么在while循环中无效？因为sdfwrite不同于fwrite可以多次 d 变量，与netCDF格式有关，在循环体中后一个sdfwrite出来的内容会将前面的内容进行覆盖，想要sdfwrite输出多个间断维度（或者变量）的话，需要先生成多个文件，再用ncks 或者ncrcat指令进行合并，具体可以参见上面贴出来的说明网址

悟空不长角

今天也是写nc时看到了楼主的帖子，也是试了一下sdfwrite,也出现了楼主的问题。所以用了lats4d,觉得这个用起来简单方便。

leeol

在手册中查询sdfwrite:
http://iges.org/grads/gadoc/gradcomdsdfwrite.html
无果

river

是不是你的define那句有问题啊，应该是'define pp=ave(prate,t='day',t='day+4')'

又是那隻貓

看看手册啊

leeol

river 发表于 2013-8-11 20:45

                               
登录/注册后可看大图

是不是你的define那句有问题啊，应该是'define pp=ave(prate,t='day',t='day+4')'

谢谢答复，可惜不是，不影响结果。
我的错误在这个GS中每次循环sdfwrite写出的物理量pp都将前次输出的覆盖，不能接在后面输出。

river

leeol 发表于 2013-8-11 23:59

                               
登录/注册后可看大图

谢谢答复，可惜不是，不影响结果。
我的错误在这个GS中每次循环sdfwrite写出的物理量pp都将前次输出的覆 ...

改一下命令的顺序看看。
'reinit'
'sdfopen d:\data\prate.sfc.day.1981-2010.ltm.nc'
'set x 1 192'
'set y 1 94'
'set z 1'
*'set t 1'
day=1
while(day<=365)
'define pp=ave(prate,t='day',t='day+4')'
'set sdfwrite d:\prate.sfc.pentad.1981-2010.ltm.nc'
'sdfwrite pp'
day=day+5
endwhile
;

mofangbao

如果只是生成一个文件放循环外面去

leeol

river 发表于 2013-8-12 08:22

                               
登录/注册后可看大图

改一下命令的顺序看看。
'reinit'
'sdfopen d:\data\prate.sfc.day.1981-2010.ltm.nc'

已测试，无效

lqouc

leeol 发表于 2013-8-11 23:59

                               
登录/注册后可看大图

谢谢答复，可惜不是，不影响结果。
我的错误在这个GS中每次循环sdfwrite写出的物理量pp都将前次输出的覆 ...

虽然不能得到你想要的结果，不过define那句最好按照river说的写，他说的方法是正确的，你那样写搞不好会出错。

风儿飘飘

你估计没有设置时间维的缘故，你试一下