z=getvar(ncfile, "height_agl") 获得位势高度后,进行插值,即可获得目标高度的要素值,完整参考代码如下: from netCDF4 import Dataset import numpy as np from wrf import getvar,interpz3d,destagger # 打开NetCDF文件 file_path = r'H:\wrfout_d01_2023-04-08_00_00_00' ncfile = Dataset(file...
ds=xr.open_dataset(flnm_wrfout)data_nc=nc.Dataset(flnm_wrfout)wrfin=data_ncifinterp_type=='ght_agl':units='Km'name_dim0='height'else:units='hPa'name_dim0='pressure'### va=getvar(wrfin,"va")ua=getvar(wrfin,"ua")### rh=getvar(wrfin,"rh")# RH_interp=wrf.vinterp(data_n...
最近,又开始了那个无休止的WRF实验循环。每次跑完都得费尽心思地检查结果,一遍又一遍地盯着那几个变量。偏偏,老板有时候突发奇想想看看其他要素的变化,然后我就得重新翻阅代码,捋一遍。 为了避免这种反复的无聊事情,我决定写个代码来提取任意数量的变量。我这里只测试了我常用的一些二维和三维的变量,并将它们插值到...
另外,WRF 模式中不同边界层参数化方案对边界层 高度(Planetary Boundary Layer Height,PBLH)的 计算方法也不同,例如 ACM2 与 YSU 使用临界理 查德森数计算边界层高度 , 而 MYJ 则由 TKE (Turbulence kinetic energy)廓线诊断得出.这里,实 况边界层高度采用 Nielsen-Gammon et al.(2010) ...
height_agl: 'PH', 'PHB', 'HGT' air_temperature: 'T', 'P', 'PB', 'T2' wind: 'U', 'V', 'U10', 'V10' relative_humidity: 'T', 'P', 'PB', 'QVAPOR', 'T2', 'PSFC', 'Q2' dew_temperature: 'P', 'PB', 'QVAPOR', 'PSFC', 'Q2' air_pressure: 'P', 'PRE...
气候模 known_lat = -89.99583 known_lon = -179.995832 S; j+ P5 e y6 U wordsize = 2 tile_x = 1200 tile_y = 1200 tile_z = 1 tile_bdr=3气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模 units=meters MSL description=Topography height气象,数值,模式,气象数值模式,数值模式,...
虽然提供了众多数据集的Vtables,但是ungrib不可能以GRIB格式预测每种可能的气象数据来源。当ungrib.exe处理新的数据源时,用户可以从头开始创建新的Vtable,也可以使用现有的Vtable作为示例。 每个Vtable包含7个或11个字段,具体取决于Vtable是用于GRIB1的数据还是GRIB2的数据。该Vtable的字段分为三类:描述如何在GRIB文件中识...
wrf.g_geoht.get_height wrf.g_geoht.get_height_agl wrf.g_geoht.get_stag_geopt wrf.g_geoht.get_stag_height wrf.g_helicity.get_srh wrf.g_helicity.get_uh wrf.g_omega.get_omega wrf.g_pressure.get_pressure wrf.g_pressure.get_pressure_hpa wrf.g_pw.get_pw wrf.g_rh.get_rh wrf.g...
z=getvar(ncfile,"height_agl") 获得位势高度后,进行插值,即可获得目标高度的要素值,完整参考代码如下: fromnetCDF4importDataset importnumpyasnp fromwrfimportgetvar,interpz3d,destagger # 打开NetCDF文件 file_path =r'H:\wrfout_d01_2023-04-08_00_00_00' ...
The profiles variations often appear between 400 m and 1200 m AGL. Within this range of height, the nonlocal PBL schemes always predict higher Qvs than the local PBL schemes. Moreover, the difference in Qv profiles over ocean is often greater than that over land. [56] To investigate the ...