- 积分
- 2005
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-4-7
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 LemonTree007 于 2021-3-9 11:41 编辑
四维数据:time,level,lat,lon,所以需要循环写入,但nc文件中只输出了一层数据,输出k值的话是能看到1-27的,所以问题应该出在持续添加数据到文件中,请问是哪里需要改一下,才能让他连续写入nc文件而不被覆盖掉?
'open G:/month/q1.ctl'
'set sdfwrite -flt G:/month/q1.nc'
k=1
while(k<=27)
'set gxout fwrite'
'set x 1 240'
'set y 1 119'
'set z 'k''
'set t 1 384'
'define Q1=q1'
k=k+1
'sdfwrite Q1'
endwhile
'disable fwrite'
'reinit'解决办法:不用while循环直接转nc,之前之所以放弃这种办法是因为闪退了,后经人指点发现缩小时间选取范围可以执行,测试显示330-340之间为闪退点,那就分两段转就好了,再拼接起来。gs脚本如下:
'open G:/q1.ctl'
'set gxout fwrite'
'set sdfwrite -flt G:/q1_period1.nc'
'set x 1 240'
'set y 1 119'
'set z 1 27'
'set t 1 324'**另写脚本'set t 325 384'
'define Q1=q1'
'sdfwrite Q1'
'disable fwrite'
'reinit'
|
|