- 积分
- 126
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-6-4
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 滟宝贝吖! 于 2025-10-1 18:11 编辑
我想将"D:\grads\Contents\Frc_byme_vdt.nc"共v、d、t三个变量(1,20,64,128)和"D:\grads\Contents\Frc_byme_p.nc"共1个变量p(1,1,64,128)写入"D:\grads\Contents\newFrc.grd"。
grads的gs代码如下(参考https://664787022.github.io/2023/12/29/LBM/):
'reinit'
'sdfopen d:\grads\Contents\Frc_byme_vdt.nc'
'sdfopen d:\grads\Contents\Frc_byme_p.nc'
'set gxout fwrite'
* 输出参数sq顺序输出,cl覆写,be大字节序输出,可在谷歌查到
'set fwrite -sq -cl -be d:\grads\Contents\newFrc.grd'
'set lon 0 357.1875'
'set lat -87.864 87.864'
*循环写入顺序: 高度-变量-时间
i=1
while(i=1)
'set t 'i''
'set dfile 1'
j=1
while(j=20)
'set z 'j''
'd v'
j=j+1
endwhile
j=1
while(j=20)
'set z 'j''
'd d'
j=j+1
endwhile
j=1
while(j=20)
'set z 'j''
'd t'
j=j+1
endwhile
'set dfile 2'
j=1
while(j=1)
'set z 'j''
'd p'
j=j+1
endwhile
i=i+1
endwhile
'disable fwrite''
但想使用ctl打开输出的newFrc.grd时,出现错误:
ga- open DgradsContentsnewFrc.ctl
Scanning description file DgradsContentsnewFrc.ctl
Open Error Missing or invalid dimension size.
-- The invalid description file record is
-- ZDEF lev 20 LEVELS 0.99500 0.97999 0.94995 0.89988 0.82977 0.74468
The data file was not opened.
请问是为什么出错呀?又应该如何修改呢
|
|