爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 102|回复: 0

[脚本编辑] grads小白求助,nc转grd层数错误

[复制链接]
发表于 前天 18:09 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

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.
请问是为什么出错呀?又应该如何修改呢

Frc_byme_vdt的ctl.png
Frc_byme_p的ctl.png
newFrc.ctl.png
ctl错误.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表