wrf-python后处理之HDF error？？？

kikichan

问题：在利用wrf-python中的getvar函数提取slp的时候报出错误，但很疑惑前后时次都正常，在遍历的时候就仅仅03这个时刻出错，而且该时刻提取其他变量如10m风速都正常，就只有slp提取的时候报错，我不禁怀疑是不是数据有问题了。。。脑壳疼希望有大神能帮忙指点一下，感激不尽{:soso__4953025489317195156_4:}

原代码如下

1. from wrf import getvar
   
2. from netCDF4 import Dataset
   
3. fplist='i:/wrfout_d01_2015-09-29_03_00_00'
   
4. data = Dataset(fplist)
   
5. slp= getvar(data,"slp",units="hPa")

复制代码

出错如下

1. File "E:\getv.py", line 11, in <module>
   
2.     slp= getvar(data,"slp",units="hPa")
   
3. 
   
4.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\routines.py", line 355, in getvar
   
5.     meta, _key, **kwargs)
   
6. 
   
7.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\metadecorators.py", line 129, in func_wrapper
   
8.     result = wrapped(*new_args)
   
9. 
   
10.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\decorators.py", line 45, in func_wrapper
    
11.     return do_conversion(wrapped(*args, **kwargs), unit_type,
    
12. 
    
13.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\g_slp.py", line 79, in get_slp
    
14.     meta=False, _key=_key)
    
15. 
    
16.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\util.py", line 2291, in extract_vars
    
17.     for var in varlist}
    
18. 
    
19.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\util.py", line 2291, in <dictcomp>
    
20.     for var in varlist}
    
21. 
    
22.   File "D:\Anaconda3\envs\py37\lib\site-packages\wrf\util.py", line 2220, in _extract_var
    
23.     result = wrfin.variables[varname][timeidx, :]
    
24. 
    
25.   File "src\netCDF4\_netCDF4.pyx", line 4406, in netCDF4._netCDF4.Variable.__getitem__
    
26. 
    
27.   File "src\netCDF4\_netCDF4.pyx", line 5350, in netCDF4._netCDF4.Variable._get
    
28. 
    
29.   File "src\netCDF4\_netCDF4.pyx", line 1927, in netCDF4._netCDF4._ensure_nc_success
    
30. 
    
31. RuntimeError: NetCDF: HDF error

复制代码

1099221723

就是数据文件有问题

kikichan

1099221723 发表于 2022-3-14 10:00
就是数据文件有问题

可是为什么同一个文件提取其他变量（如10m风速）是正常的呢

1099221723

kikichan 发表于 2022-3-14 11:39
可是为什么同一个文件提取其他变量（如10m风速）是正常的呢

可能的原因很多，最大的可能是多线程的问题，而且这种情况是小概率事件，很难排查，看看wrf的日志吧，里面或许有记录