请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 474|回复: 5

如何解决"file not open"的问题

[复制链接]

新浪微博达人勋

发表于 2020-7-3 12:52:45 | 显示全部楼层 |阅读模式
GrADS
系统平台: Grads
问题截图:
问题概况: 参照家园中的脚本,编写了一个画地形叠加流线和散度的脚本,但是grads面板显示“file 3 not open”, “file 4 not open”。但是图也能画出来,请问该如何解决"file not open"的问题呢?
我看过提问的智慧: 看过
自己思考时长(天): 1

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

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

x
12.8_2_00.png practice.gs (2.36 KB, 下载次数: 0)
Untitled picture.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2020-7-3 12:56:32 | 显示全部楼层
相关脚本在附件中,为方便查看,这里再贴上:
'reinit'
'open D:\Research_Room\program\meoteorology\surface_elevation\dz.ctl'
'sdfopen D:\Research_Room\Data\ThirdGroup\nc\multiple\ERA5_2019.12.nc'
* time resolution is 3-h(248 times)
* spacial resolution is 0.25(60E-150E, 0-70N)
* z => Geopotential, t => Temp
* r => RH, q => SH
* d => Divergence
* pv => Potential Vorticity, vo => Vortivity
* u,v,w => u,v,w component of wind
* 10-levels 1000 925 850 800 750 700 600 500 300 200

'set grid off'
'set grads off'
'set parea 1.5 8.5 1.0 7.5'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'

*=======================================
*              cycle part  
*=======================================

i = 8
* i is date
* j is hour

while(i<=8)
   j = 2
   while(j<=2)
      'c'
      'set grads off'
      'set font 1'
       a = (8*i)-(10-(j+1)/3)

      'set dfile 1'
      'set lat 30'
      'set lon 110 125'
      'set lev 1000 200'
      'set t 1'
      'define zz=dz.1'
      
      'set dfile 2'
      'set t 'a
      'set lat 30'
      'set lon 110 125'
      'set lev 1000 200'
      'define uu=u.2'

      'set dfile 3'
      'set t 'a
      'set lat 30'
      'set lon 110 125'
      'set lev 1000 200'
      'define ww=-5*w.2'

      'set dfile 4'
      'set t 'a
      'set lat 30'
      'set lon 110 125'
      'set lev 1000 200'
      'define dd=d.2'

      'set zlog on'
      'set gxout shaded'
      'set xlabs 110E|113E|116E|119E|122E|125E'
      'D:\Research_Room\Data\ThirdGroup\gs\colorbars\30colors_rainbow.gs'
      'set clevs -8e-5 -6e-5 -4e-5 -2e-5 0 2e-5 4e-5 6e-5 8e-5'
      'd dd'
      'xcbar 9.0 9.2 1.5 7.3 -fw 0.13 -fh 0.15 -ft 5 -fs 1 -fo 0 -line on'
*      xcbar left-edge, right-edge, bottom-edge, top-edge, -fw font-width,
*      -fh font-height, -ft font-thichness, -fs lable-interval,
*      -fo lable-mask-number, -line on/off(colorbar-boarder on/off)

      'set gxout stream'
      'set strmden 5 0.4 0.1 1'
      'set cthick 5'
      'set ccolor 1'
      'set cint 1'
      'd uu;ww'

      'set gxout shaded'
      'set csmooth on'
      'set rbcols 1'
      'set cmax 0'
      'd zz'

      'drawfigstr TL (a) 2.5'
      'draw title Dec.'%i%'th '%j%':00'
      'printim D:\Mine\Work\Program\My_Program\program\script\cross_section\12.'%i%'_'%j%'_00.png white'
      
       j=j+6
   endwhile
i=i+1
endwhile
;









密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-7-3 13:32:12 | 显示全部楼层
是你gs里有'set dfile 3'和 'set dfile 4'这两句,但实际上你只打开了两个文件,所以那两个语句报错了。

评分

参与人数 1金钱 +10 贡献 +1 收起 理由
李文帅 + 10 + 1 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-7-4 08:47:39 | 显示全部楼层
苏拉苏拉 发表于 2020-7-3 13:32
是你gs里有'set dfile 3'和 'set dfile 4'这两句,但实际上你只打开了两个文件,所以那两个语句报错了。

那这两句应该都写成'set dfile 2'是吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-7-5 15:35:03 | 显示全部楼层
李文帅 发表于 2020-7-4 08:47
那这两句应该都写成'set dfile 2'是吗?

可以的,或者直接去掉,也不影响结果。你的变量使用了.2下标,所以set dfile可用可不用
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-7-7 09:17:41 | 显示全部楼层
苏拉苏拉 发表于 2020-7-5 15:35
可以的,或者直接去掉,也不影响结果。你的变量使用了.2下标,所以set dfile可用可不用

好的,十分感谢!!!!!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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