循环写入变量

shzljg

用的grads2ascii程序，我下载的是2019-2023小时数据，编写了相应的gs程序。2019年写入正常，到2020年就出错。哪位高手帮帮我看看怎么回事。
'reinit'

j=2019
while(j<=2020)
  'sdfopen e:\byblk\'j'.nc'
  'set lat 43.03'      
  'set lon 84.15'     
  'set lev 750 150'


  if(j=2020)
    i=1
    while(i<=8684)
          'set t 'i''  
      'grads2ascii t-273.15 e:\byblk\'j't.txt a %5.2f 15 2'
      i=i+1
    endwhile
  else
    i=1
    while(i<=8760)
      'set t 'i''   
      'grads2ascii t-273.15 e:\byblk\'j't.txt a %5.2f 15 2'
      i=i+1
    endwhile
  endif
  j=j+1
endwhile

;

lqouc

j=j+1前面加一行close 1，或者单独只运行2020试试看

river

你这个nc文件里面的温度变量名就叫 t ？

shzljg

river 发表于 2024-8-19 14:17
你这个nc文件里面的温度变量名就叫 t ？

是的，温度（t）.不知道能不能把变量名改了。

shzljg

2019,2021,2022,2023都没问题，就是2020（闰年）不行。

shzljg

lqouc 发表于 2024-8-19 09:58
j=j+1前面加一行close 1，或者单独只运行2020试试看

加了close 1，还是运行出错，提示14行语法错误。只算2020年，也不行。

river

你把判断闰年的方式修改一下试试，写成类似下面这样的
if ((math_mod(yy,4)=0&math_mod(yy,100)!=0)|math_mod(yy,400)=0)