爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 2156|回复: 4

读取wrfout文件中变量存为nc文件报错

[复制链接]
发表于 2023-7-6 11:14:43 | 显示全部楼层 |阅读模式

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

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

x
报错信息:fatal:FileAddVar: an error occurred while adding variable <times> to file <wrfout2021100100>, check to make sure data type is supported by the output format,求助,为什么新的nc文件写不了times变量,是因为她是string型吗?要如何实现把它写入新的nc文件(别的变量都可以正常写入)
这是我写的代码:
    a1 = addfiles(FILES1+".nc","r")
    numFILES1 = dimsizes(FILES1)
;*******read data******************
    it = -1 ;-1 is all times
    TIMES = wrf_user_getvar(a1,"times",it)  ; get all times in the file   
    Lat   = wrf_user_getvar(a1,"lat",it)
    Lon   = wrf_user_getvar(a1,"lon",it)
    Tk    = wrf_user_getvar(a1,"tk",it) ;Temperature [K]
    U     = wrf_user_getvar(a1,"ua",it)
    V     = wrf_user_getvar(a1,"va",it)
    P     = wrf_user_getvar(a1,"pressure",it) ;Full model pressure [hPa]
    RH    = wrf_user_getvar(a1,"rh",it) ;Relative Humidity [%]
    SLP   = wrf_user_getvar(a1,"slp",it) ;Sea level pressure [hPa]
    U10   = wrf_user_getvar(a1,"U10",it)  ; u at 10 m
    V10   = wrf_user_getvar(a1,"V10",it)  ; v at 10 m
    T2    = wrf_user_getvar(a1,"T2",it)
    TD2   = wrf_user_getvar(a1,"td2",it);2m dew point temperature [C]
    rain_exp = wrf_user_getvar(a1,"RAINNC",-1)  ;累计格点降水
    rain_con = wrf_user_getvar(a1,"RAINC",-1)    ;累计积云对流降水

;*******creat new file******************
    outfile = outdir + "wrfout"+DATE + ".nc"
    print(outfile)
    if (isfilepresent(outfile)) then
      system("rm -rf "+outfile)
    end if
    setfileoption(outfile,"Format","LargeFile")
    fout = addfile(outfile, "c")

    fout->times = TIMES
    fout->lat = Lat
    fout->lon = Lon
    fout->tk = Tk
    fout->u = U
    fout->v = V
    fout->p = P
    fout->rh = RH
    fout->slp = SLP
    fout->u10 = U10
    fout->v10 = V10
    fout->t2 = T2
    fout->td2 = TD2
    fout->RAINNC = rain_exp
    fout->RAINC = rain_con

密码修改失败请联系微信:mofangbao
发表于 2023-7-7 00:15:20 | 显示全部楼层
如果印象没错的话,wrfout里面的时间是单个字符(S1),然后有个字符长度的维度,字符串的话估计得先转换下。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-9-20 09:51:44 | 显示全部楼层
请问楼主解决了吗~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-3-20 16:32:19 | 显示全部楼层
我也出现这样的报错,改了这句命令就好了setfileoption("nc","Format","NetCDF4"),主要是最后一个属性
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2025-4-11 20:17:50 | 显示全部楼层
Sherozbq 发表于 2024-3-20 16:32
我也出现这样的报错,改了这句命令就好了setfileoption("nc","Format","NetCDF4"),主要是最后一个属性

请问,您还有转换脚本吗,我最近也在做这个,如果有的话,非常感谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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