求教读取GODAS数据相关的问题

MgLily

报错如图，代码如下
我下载了1982-2020的3月到10月初的GODAS海洋数据(grib格式)，在读取的时候报错，不知道有没有人之前也遇到过类似问题，求解

begin
target_var_num=7
varname=(/"POT","DZDT","U","V","SALI","UF","VF","ML","IL","SSH","SLTFL","THFLX"/)

year=1982

file_list=systemfunc("ls /sky1/home/zhangmy43/data/GODAS/"+year+"/*.grb")
;从第一个文件里提取一些文件共有的信息(文件大小，文件名)
f_in=addfile(file_list(0),"r")
var_names=getfilevarnames(f_in)
print(var_names)
x5_0=f_in->$var_names(target_var_num)$
printVarSummary(x5_0)
dim_x5_0 = dimsizes( x5_0 )

x5=new((/39,44,dim_x5_0(0),dim_x5_0(1)/),typeof(x5_0),x5_0@_FillValue);new((/nyear,nday,nlat,nlon/),float)
do j=0,38;
  year=1982+j
  file_list=systemfunc("ls /sky1/home/zhangmy43/data/GODAS/"+year+"/*.grb")

  do i=0,43;
    f_in=addfile(file_list(i),"r")
    ;printVarSummary(f_in)
    var_names=getfilevarnames(f_in);
    x5(j,i,:,:)=f_in->$var_names(target_var_num)$;(！！！！！报错的位置！！！！！)
    delete(var_names)
    delete(f_in)
  end do

  delete(file_list)
end do

x5!1="ncl_join"
printVarSummary(x5)

exit

MgLily

自顶