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

李文帅

practice.gs (2.36 KB, 下载次数: 1)

2020-7-3 12:49 上传

点击文件名下载附件

李文帅

相关脚本在附件中，为方便查看，这里再贴上：
'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
;

苏拉苏拉

是你gs里有'set dfile 3'和 'set dfile 4'这两句，但实际上你只打开了两个文件，所以那两个语句报错了。

李文帅

苏拉苏拉 发表于 2020-7-3 13:32
是你gs里有'set dfile 3'和 'set dfile 4'这两句，但实际上你只打开了两个文件，所以那两个语句报错了。

那这两句应该都写成'set dfile 2'是吗？

苏拉苏拉

李文帅 发表于 2020-7-4 08:47
那这两句应该都写成'set dfile 2'是吗？

可以的，或者直接去掉，也不影响结果。你的变量使用了.2下标，所以set dfile可用可不用

李文帅

苏拉苏拉 发表于 2020-7-5 15:35
可以的，或者直接去掉，也不影响结果。你的变量使用了.2下标，所以set dfile可用可不用

好的，十分感谢！！！！！