ncl文件引用的问题

leos

因为我python的xarray不能直接读取grib，所以我打算通过ncl将grib转换成nc再读取。
但由于我的文件大小太大（7个多g），ncl_covert2nc不能工作，所以就只能手搓一个转换代码。
ncl里面的`=`可以对文件使用，把一个文件赋值给另一个。但我在尝试的过程中发现这样的赋值只是把文件指针给传递了过去。通过copy_VarMeta方法似乎也只能在统一格式下进行转换。报错：fatal:Could not create (/mnt/e/data/200901.nc)

所以有什么办法可以转换大的grib文件吗？

leos

这是用copy_VarMeta方法转换的代码
begin
    setfileoption("nc", "Format", "NetCDF4Classic")
    do ye = 2009, 2023
        y = tostring(ye)
        do mon = 1, 12
            ; 用来生成数据索引
            if(ye .eq. 2023) .and. (mon .eq. 7)
                break
            end if
            if(mon.lt.10)
                m = tostring(0) + tostring(mon)
            else
                m = tostring(mon)
            end if
            path0 = tostring("/mnt/e/data/")+y+m+tostring(".grib")  ; 待转换文件
            path1 = tostring("/mnt/e/data/")+y+m+tostring(".nc")    ; 保存文件
            infile = addfile(path0, "r")
            outfile  = addfile(path1,"c")
            copy_VarMeta(infile, outfile)
            print(outfile)
            exit()
        end do
    end do
end

是冉冉升起的冉

用cdo更快！