首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

数据处理 | python计算经典相当位温2.0

转载本文时,请务必文字注明“来自:和鲸社区:酷炫用户名”,并附带本项目超链接。 前言 本项目旨在通过Python实现经典相当位温的计算方法,帮助大家更好理解位温概念。...内容包括:相当位温的基本概念,计算方法,Python代码示例与简单可视化 相当位温(Equivalent Potential Temperature)是描述大气状态的一个重要指标。...包含冻结过程的广义位温及位涡特征分析 [J]. 大气科学, 44(4): 816−834. ZHOU Kuo, RAN Lingkun, QI Yanbin, et al. 2020....小知识 html中上标的输入格式是数字 数字 latex中分号是frac{{a}}{{b}} 下标格式为 a_b ,有复数字母时需要用{}括起来 1....实际应用 : WRF后处理提取相关变量计算相当位温 还是从老伙计wrfout中提取需要的变量:位温 温度 气压 from wrf import uvmet, to_np, getvar, interplevel

13510
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    如何计算WRF台风模拟的假相当位温

    转载本文时,请务必文字注明“来自:和鲸社区:酷炫用户名”,并附带本项目超链接。...温馨提示 由于可视化代码过长隐藏,可点击如何计算WRF台风模拟的假相当位温运行Fork查看 若没有成功加载可视化图,点击运行可以查看 ps:隐藏代码在【代码已被隐藏】所在行,点击所在行,可以看到该行的最右角...实际应用 : WRF后处理提取相关变量计算假相当位温 还是从老伙计wrfout中提取需要的变量:温度 气压 混合比 等等 设置函数 In [4]: def calculate_theta_se_wrf(...(wrf_list, 'pressure', timeidx=-1).data t = getvar(wrf_list, 'tk', timeidx=-1).data r = (getvar(wrf_list..., 'QVAPOR', timeidx=-1)/1000).data td = getvar(wrf_list, 'td', timeidx=-1).data lons = getvar(wrf_list

    21110

    C++的netCDF读取NC数据:踩坑经历

    针对不同语言,netCDF库也有着对应的不同版本,我们之前的文章也介绍过在Python、C++等代码的开发环境中,配置netCDF库的具体方法,大家直接参考文章Anaconda配置h5py与netCDF4...NcFile file(path, NcFile::read); NcVar var = file.getVar("ssrd");   那么此时,假设我上述的2句代码是在一个函数中运行的,然后我将得到的...NcVar格式的变量var作为这个函数的返回值,返回给调用者;随后,又将这个被返回的NcVar格式变量var作为参数,输入到另一个新的函数中——那么,在这个新的函数中,我们如果还想读取var这个变量(例如用如下的代码来读取...4 时间维度需要放在第一个位置   有时,我们需要按照不同维度,对变量数据加以读取。例如,我这里的.nc格式文件中,变量是1种气象数据,其具有3种维度,包括经度、纬度与时间等。   ...如果是如此,就需要注意将时间维度放在start的第一个元素的位置上;具体代码如下所示。

    36010

    基于WRFOUT计算相对涡度,绝对涡度,位涡并可视化

    版本:python3.7 数据:wrfout模拟数据 核心代码:metpy.calc.vorticity 前言 涡度是流体力学中的一个重要概念,用于描述流体运动中的旋转性质。...=None, latitude=None, longitude=None, crs=None) 因此我们需要从wrfout文件中提取 u,v ,而p是插值需要的参数 # 用 netCDF4 包读取 WRF...=0) v = getvar(wrf_file, 'va', timeidx=0) p = getvar(wrf_file, 'pressure', timeidx=0) # 提取WRF模拟的经纬度数组...,其单位是10-5 s-1,所以在绘图时不需要乘1e5 绘制850hPa绝对涡度分布图 位涡 罗斯贝提出的一个类似位温的用于垂直涡度的概念 ,其公式为 位涡数据获取 pvo = getvar(wrf_file...完整代码与文件在这里,文件在注册社区账号点击左侧文件标识可下载,代码需要右上角在线运行

    52510

    【专业技术】从4行代码看右值引用

    右值引用实际上并没有那么复杂,其实是关于4行代码的故事,通过简单的4行代码我们就能清晰的理解右值引用相关的概念了。...本文希望带领读者通过4行代码来理解右值引用相关的概念,理清他们之间的关系,并最终能透彻地掌握C++11的新特性--右值引用。...四行代码的故事 第1行代码的故事 int i = getVar();   上面的这行代码很简单,从getVar()函数获取一个整形值,然而,这行代码会产生几种类型的值呢?...第2行代码的故事 T&& k = getVar();   第二行代码和第一行代码很像,只是相比第一行代码多了“&&”,他就是右值引用,我们知道左值引用是对左值的引用,那么,对应的,对右值的引用就是右值引用...第3行代码的故事 T(T&& a) : m_val(val){ a.m_val=nullptr; }   这行代码实际上来自于一个类的构造函数,构造函数的一个参数是一个右值引用,为什么将右值引用作为构造函数的参数呢

    1.6K71

    python计算与绘制WRF降水量

    读取形成的文件列表 2.使用wrfpython的getvar读取多个wrf文件的RAINC,RAINNC,RAINSH,利用cat将多时次数据合并 例如,RAINC = getvar(wrf_list,...关于三个降水变量的区别可以参考 WRF后处理:降雨量的说明以及降雨的绘制_wrf模拟降水量偏小-CSDN博客 https://blog.csdn.net/islandowner2017/article/...# 导入数据读取模块 import numpy as np import pandas as pd from netCDF4 import Dataset import xarray as xr #...) # 添加标题 plt.title('08-08 1800 - 08-09 0600(UTC)', loc='left', fontsize=20) plt.show() 每小时降水量组图绘制 为了代码不繁琐直接利用...xarray的plot作图,更多细致的作图敬请自己实现,以下示例小时降水量的组图绘制 此处使用了xarray的data.diff计算每小时的降水量 wrfout中的降水变量都是累计降水量,因此需要根据用后一时次减去前一时次才能得出这小时下了多少

    57811

    经典位温的Python 实现

    内容包括:相当位温的基本概念,计算方法,Python代码示例与简单可视化 相当位温(Equivalent Potential Temperature)是描述大气状态的一个重要指标。...包含冻结过程的广义位温及位涡特征分析 [J]. 大气科学, 44(4): 816−834. ZHOU Kuo, RAN Lingkun, QI Yanbin, et al. 2020....小知识 html中上标的输入格式是数字 数字 latex中分号是frac{{a}}{{b}} 下标格式为 a_b ,有复数字母时需要用{}括起来 1....实际应用 : WRF后处理提取相关变量计算相当位温 还是从老伙计wrfout中提取需要的变量:位温 温度 气压 In [22]: from wrf import uvmet, to_np, getvar...> 4. metpy的相当位温 公式为 from here metpy In [26]: from metpy.calc import equivalent_potential_temperature

    20211

    经典相当位温的Python 实现

    内容包括:相当位温的基本概念,计算方法,Python代码示例与简单可视化 相当位温(Equivalent Potential Temperature)是描述大气状态的一个重要指标。...包含冻结过程的广义位温及位涡特征分析 [J]. 大气科学, 44(4): 816−834. ZHOU Kuo, RAN Lingkun, QI Yanbin, et al. 2020....小知识 html中上标的输入格式是数字 数字 latex中分号是frac{{a}}{{b}} 下标格式为 a_b ,有复数字母时需要用{}括起来 1....实际应用 : WRF后处理提取相关变量计算相当位温 还是从老伙计wrfout中提取需要的变量:位温 温度 气压 In [22]: from wrf import uvmet, to_np, getvar...> 4. metpy的相当位温 公式为 from here metpy In [26]: from metpy.calc import equivalent_potential_temperature

    31510

    WRFOUT风向变量逐时次作差绘图

    同时从“wrfout中提取变量,然后用08:10的风向wdir【ncl函数wind_direction(u,v,0)】减去08:00时刻的风向, 做上循环语句do,就会出现差一个数值对不上的情况。...因为wrfout的变量是xarray格式,想必大家知道要用哪个函数了。 没错,就是xarray.diff() 废话半天了,开始写代码吧。...(wrf_list, 'uvmet10_wspd_wdir', timeidx=ALL_TIMES, method='cat')[1] uvmet10 Out[3]: 以上代码提取了wrfout数小时的...关于风变量的讲解可以看看https://www-k12.atmos.washington.edu/~ovens/wrfwinds.html 简单绘图 In [4]: # #计算逐小时风向变化 hourly_dirchange...matplotlib.cm.register_cmap(name=cname, cmap=cmap) 当然,风电场一般要看离地一百米左右的风,提取方法在风能密度讲过 In [5]: wdir = getvar

    16710

    开工!wrfout 计算台风准地转omega方程右侧项

    在编写代码时,我们需要设置静态稳定度、f0和Rd的常数,并结合MetPy库提供的功能来进行数值计算。...由于可视化代码过长隐藏,可点击运行查看 若没有成功加载可视化图,点击运行可以查看 ps:隐藏代码在【代码已被隐藏】所在行,点击所在行,可以看到该行的最右角,会出现个三角形,点击查看即可 03、导入库与读取文件...units import xarray as xr import os import numpy as np import matplotlib.pyplot as plt from netCDF4...=0) v = getvar(wrf_list[5], 'va', timeidx=0) t = getvar(wrf_list[5], 'tc', timeidx=0) p = getvar(wrf_list...interplevel(u, p, 500) v500 = interplevel(v, p, 500) sigma = 2.0e-6 * units('m^2 Pa^-2 s^-2') f0 = 1e-4

    19910

    更快更强!四种Python并行库批量处理nc数据

    前言 当前镜像:气象分析3.9 资源:4核16g 注意分开运行,不然会爆内存 阅读本文你将学到: 远超循环批量处理nc文件效率的技巧 四种并行库的基本使用与区别 wrf变量极值经纬度索引 Dask...温馨提示 由于可视化代码过长隐藏,可点击 更快更强!...四种Python并行库批量处理nc数据 运行Fork查看 若没有成功加载可视化图,点击运行可以查看 ps:隐藏代码在【代码已被隐藏】所在行,点击所在行,可以看到该行的最右角,会出现个三角形,点击查看即可...for循环 In [10]: %%timeit import os import glob from netCDF4 import Dataset from wrf import getvar, latlon_coords...资源改为4核16g时,并行超越了单循环 当你核数和内存都没困扰时当然是上并行快 ,但是环境不一定能适应多线程 资源匮乏或者无法解决环境问题时还是老实循环或者在列表推导式上做点文章

    66610
    领券