WRF | WRF散度计算步骤及Python可视化实现 Image Name 关注我获取更多学习资料,第一时间收到我的Python学习资料,也可获取我的联系方式沟通合作 环境设置 依赖 wrf-python...np import xarray as xr import metpy.calc as mpcalc from metpy.interpolate import cross_section from wrf...matplotlib.pyplot as plt import cartopy.crs as ccrs from netCDF4 import Dataset from metpy.units import units # 读取WRF...coords=u.coords, attrs={'units': 's^-1', 'description': 'divergence'} ) div_da 取剖面 from wrf...Divergence Vertical Cross-Section') plt.show() 细化绘图 小结 本来想基于metpy的插值函数做个两种方法做个对比,但是其函数对于数据的要求比较苛刻,不太适合wrf
这部分包含的时 wrf-python 模块中的API,如果wrf-python提供的函数不能满足你的需求,你也可以根据已有的API重新编写一个处理函数或是其它的诊断函数。...当然,大多数时候我们仅需要wrf-python已经提供的诊断函数和其它可调用函数,可能都不会使用到这些API。...用户API 方法 诊断 从WRF输出中返回基本诊断变量 插值 返回3D场插值到2D平面(指定垂直层)的值 坐标转换 返回经纬度坐标对应的X,Y坐标 网格去栅格 返回去栅格化后的变量,与NCL中的wrf_user_unstagger...Projection 基本类 存储地图投影信息的基类 Projection 基本类方法 返回一个matplotlib.mpl_toolkits.basemap.Basemap 地图投影对象 Projection 子类 wrf.WrfProj
近几年,python在气象领域的发展也越来越快,同时出现了很多用于处理气象数据的python包。比如和NCL中的 WRF_ARWUser库类似的 wrf-python模块。...wrf-python是用于WRF模式后处理的python模块,其中提供了很多有用的函数,下面就来详细说一下其用法: 基本用法 计算诊断变量 wrf.getvar 函数的主要作用是返回需要计算的诊断变量...meta=False) print (type(p_no_meta)) 从DataArray中提取 numpy 数组 如果你需要将 xarray.DataArray 转换为 numpy.ndarray, wrf-python...wrf-python中有算法会对缺省值数组进行检查,但是当你编译模块时,如果模块代码中使用了wrf-python,那么就要小心了,应尽量避免出现上述情况。...lat_lon = xy_to_ll(ncfile, [400,105], [200,205]) x_y = ll_to_xy(ncfile, lat_lon[0,:], lat_lon[1,:]) 添加地图 wrf-python
最近需要针对对WRF应用做一些性能优化;下面简单记录WRF model 的安装过程: 优秀UP主的视频: https://www.bilibili.com/video/av71541427/ UP主提供的资源...; WRF 和 WPS 都要下载; https://www2.mmm.ucar.edu/wrf/users/download/get_source.html 配置编译环境;gfortran,perl,.../users/download/get_sources_wps_geog.html 其他相关优秀教程: https://xg1990.com/blog/archives/tag/wrf 包含wrf安装运行...; 2020/8/29日更新步骤: 工作目录 /home/nfs /home/wrf; 配置共享nfs文件系统; 在jn-01 机器上,使用 /home/nfs/init-WRF-environment.sh...脚本,进行初始化WRF执行环境,并获取用于执行的csh; 运行 run-wps-wrf.sh 脚本,进行wps wrf 样例数据; 根据样例配置,修改运行过程和数据输入输出过程,完成WRF软件的执行;
前言 1.使用os库循环读取文件夹下的wrf数据,并用nc库的dataset读取,可使用wrf_list = [Dataset(f) for f in wrf_files] ,wrf_files是os...读取形成的文件列表 2.使用wrfpython的getvar读取多个wrf文件的RAINC,RAINNC,RAINSH,利用cat将多时次数据合并 例如,RAINC = getvar(wrf_list,...关于三个降水变量的区别可以参考 WRF后处理:降雨量的说明以及降雨的绘制_wrf模拟降水量偏小-CSDN博客 https://blog.csdn.net/islandowner2017/article/...import getvar # WRF数据目录 wrfout_path = '/home/mw/input/typhoon9537/' # 读取 WRF 模拟数据 wrf_file = Dataset...,cartopy_ylim,to_np,latlon_coords # 读取 WRF 模拟数据 wrf_file = Dataset(wrfout_path + 'wrfout_d01_2019-08-
目前,有很多工具可以进行WRF模式后处理,比如NCL,Grads,Python,MATLAB等等,而且每一种语言都有其优势。...NCL中有WRF模式后处理的包,可以非常方便的处理WRF模式结果,而近些年Python在气象上的应用也越来越广泛,各种气象相关包也日渐完善。...MATLAB在WRF模式后处理方面就显得捉襟见肘了,倒不是MATLAB不适合做WRF模式后处理,而是关于这方面的开源包不多。...由于之前主要使用MATLAB处理WRF模式结果(虽然已投入Python的怀抱),而MATLAB中又没有相关函数,便对NCL中的部分函数进行了移植。...转换的源程序包括: 气象参数计算: 1) 露点计算 2) 相对湿度计算 3) 根据 风向风速计算u , v 4) 根据u, v分量 计算 风向风速 根据WRF输出结果计算并绘制: 1) 计算温度
以下文章来源于气海同途 ,作者气海同途 编者按:这是新开的一个系列,有时间会逐步将WRF官方培训的ppt挑选个人认为重要的进行翻译,以及结合个人的使用经验进行一些解释。...本篇内容来源于WRF官网2021年的培训ppt:Application of WRF: How to Get Better Performance。...WRF作为成熟的区域中尺度气象模式,文档齐全且教程详细,对于用户较为友好,但是想要获得一个好的模拟结果,需要注意很多地方, 1 模拟区域domain设置 模拟区域不能太小,否则模拟结果基本为全球模式侧边界的强迫结果...3 侧边界条件 包括WRF在内的所有区域模式,人为引入的侧边界(LBC)是制约模式性能的一个重要因素。...7 物理参数方案的选择 WRF模式中,存在众多次网格过程需要使用物理参数方案进行表达,而每种物理参数化方案都有很多种选项,其不同的排列组合使得选项太多。
WRF模式是数值天气预报和大气模拟系统,其开发目的就是用语研究和实际应用。运行WRF模式时,可以利用多种初始场数据来驱动,然后配置好选项之后便可以模拟天气过程(说的好像很简单的样子==)。...既然涉及气象领域,那么自然就少不了对WRF模式的支持,NCL中有一个脚本库,专门应用于WRF模式的后处理[注1]。...关于python不作过多介绍,这也不是重点。重点是使用python进行WRF后处理。...Python进行WRF模式后处理,主要使用三个库:matplotlib(python中最火的可视化库),netCDF4(处理nc文件),Basemap(处理地图投影)。...中有一个库和NCL 中的 WRF相关库具有几乎相同的功能,就是 wrf-python。
版本:python3.7 参考了气象备忘录的WRF模拟区域绘制和局地放大 改进点: 1. 增加了经纬度刻度 2. 根据现有的资料,针对部分设置进行删改(在绘图方面) 3....对读取文件进行优化 from netCDF4 import Dataset import numpy as np from wrf import getvar, get_cartopy import matplotlib.pyplot
cartopy.crs as crs import cartopy.feature as cfeature from cartopy.feature import NaturalEarthFeature from wrf...smooth2d, CoordPair, GeoBounds,interpline import warnings warnings.filterwarnings('ignore') 读取数据 wrf_file.../input/wrf3880/wrfout_d01_2016-10-07_00_00_00_subset.nc' # Open the NetCDF file ncfile = Dataset(wrf_file...slp) # Create a figure fig = plt.figure(figsize=(12,6)) # Set the GeoAxes to the projection used by WRF...cross_start = CoordPair(lat=43.5, lon=-116.5) cross_end = CoordPair(lat=43.5, lon=-114) # Get the WRF
(): global _WRF_FILES return _WRF_FILES[0] def multiple_wrf_files(): global _WRF_FILES...() wrf_file = Dataset(file_path) slp = getvar(wrf_file, "slp", timeidx=0) td2 = getvar(wrf_file, "td2...() wrf_file = Dataset(file_path) p = getvar(wrf_file, "pressure") z = getvar(wrf_file, "z", units="dm...") ua = getvar(wrf_file, "ua", units="kt") va = getvar(wrf_file, "va", units="kt") wspd = getvar(wrf_file...代码获取在好奇心Log公众号后台回复wrf绘图
“分段错误可能难以追踪。由于通常没有明确的错误消息,因此可能需要反复试验才能找出问题所在。我试了好久(•́へ•́╬)!大致总结了一下,给大家参考,如果还有其他情...
做WRF后处理软件的过程中主要关注了整体流程拉通和产品输出,画的图片比较粗糙,今天静下心来(被迫营业)慢慢修图,经过小组会议的讨论,总结的问题主要可以归纳为以下几点: y轴的纬度标签长短不一 省界过于突出...gsLineThicknessF = 3.0 prov_plot = gsn_add_shapefile_polylines(wks,plot,shp_prov,prov_res) 数据平滑 wrf_smooth...solidstar plat = xxx plon = xxx dum = gsn_add_polymarker(wks,plot,plon,plat,pmres) 裁减图片白边 这里主要是利用python
WPS只是WRF运用实际数据模拟的准备过程,这就好比开车要先上车。讲完前处理过程之后 WRF模式之WPS的配置运行,就要开始模拟了。.../wrf.exe 上述方式为串行运行,如果编译支持并行运算的话,可以执行以下命令: mpirun -np 4 ./wrf.exe 表示使用4个处理器进行运算。...关于模式后处理,之前也提到过一些: 真・WRF模式后处理之Python版 WRF模式后处理 高维数据可视化 后处理过程其实就是处理数据的过程。使用什么工具或是编程语言都可以,没有限制。
WRF散度计算步骤及Python可视化实现2.0 项目概述 在气象家园论坛中,有用户对WRF模式输出的散度计算感到困惑。本项目旨在介绍WRF散度的计算方法,并提供Python可视化实现。...环境设置 依赖库 wrf-python metpy numpy xarray matplotlib cartopy 数据处理 数据读取 import numpy as np import xarray...as xr import metpy.calc as mpcalc from metpy.interpolate import cross_section from wrf import getvar,...import cartopy.crs as ccrs from netCDF4 import Dataset from metpy.units import units import metpy # 读取WRF...特别是处理多层数据时 计算结果可以通过xarray和matplotlib方便地进行可视化 参考资料 气象家园论坛讨论 Stack Overflow上的相关讨论 MetPy官方文档 作者信息 公众号: 气python
那今天就给大家介绍一下当前已经发布的社区版GPU加速的WRF模式-WRFg。...WRFg是利用GPU加速的社区版中小尺度WRF模式,包含了完整的WRF动力核及一些微物理选项,可以利用GPU、OpenACC和CUDA加速WRF模式的运行。...官方给出的性能对比是,GPU加速的WRF模式比CPU版本的速度提高了7倍。 官方网站给出的信息,加速后的WRF模式可以运行1km分辨率,但是没有指明多大范围。...加速后的WRF模式可以同化更多的观测数据,以改善数值模式的初始场。 网站发布的是基于WRF V3.8.1的加速版本。以下是当前已经移植到GPU的一些子模块信息。 ? ?
利用python中的cartopy、wrf-python等库,绘制wrf中的土地利用类型。...主要使用了pcolormesh函数进行绘制,绘制效果如下: type3 原始版本 主要参考了Python气象数据处理与绘图:绘制WRF模式模拟所用的土地利用数据来进行绘制。...具体使用的版本如下: cartopy 0.18.0 matplotlib 3.5.1 wrf-python 1.3.3 其他库如下,一般版本也没啥大的限制,就没有一一列举了。...LATITUDE_FORMATTER, LONGITUDE_FORMATTER import shapely.geometry as sgeom from copy import copy from wrf...如果大家有其他思路和改进方法的,欢迎积极投稿和分享~~ Reference [1] [Python气象数据处理与绘图:绘制WRF模式模拟所用的土地利用数据] (https://mp.weixin.qq.com
Weather Research and Forcasting Model WRF-ARW模式是针对大气研究和业务预报应用设计的中尺度数值天气预报系统,适用于从几十米到几千公里的广泛气象应用。...本文将以WRF-ARW V3.9.1版本为例,介绍其在超算中心上模拟飓风马修登陆。给大家简单运行一下WRF在超算平台的实际应用案例,并建立起数值模式模拟的一个初步印象。...---- ---- 本教程将以官网—2016年10月飓风马修登陆为例,介绍WRF-ARW在北京超算中心上的基础运行操作。...有关namelist.wps的详细参数说明可以参考WRF-ARW官网: https://www2.mmm.ucar.edu/wrf/users/namelist_best_prac_wps.html...和边界场wrfbdy_d01,适配背景超算中心的slurm作业脚本如下: 3.2、wrf.exe:模式运行积分 提交wrf.exe作业,得到模式结果和Restart文件。
其公式为 from wrf import uvmet, to_np, getvar, interplevel, smooth2d, get_cartopy, cartopy_xlim, cartopy_ylim
前言 本项目旨在解决在WRF的ungrib报错GRIB文件类型不一致 错误详情 在链接era5的single数据时后 进行ungrid.exe时出现以下报错 output format is WPS Path...era5的single层数据变量 grib_ls -p editionNumber era5_sf.grib 结果如下 Image Name 原来是一群grib1变量中出了一个grib2变量 对于这种情况,wrf