MPAS(the Model for Prediction Across Scales)最早是美国阿拉莫斯实验室(LANL)和美国大气研究中心(NCAR)的一个合作项目，旨在发展大气、海洋和其他地球系统分量模式用于研究全球或区域的天气和气候。LANL主要负责海洋分量模式MPAS-O的开发，而大气分量模式MPAS-A主要由NCAR负责，MPAS的框架、应用程序和工具则由两者共同 ...

最近气象大模型比较火，分享一下在气象大模型论文中常见的几个评估指标(RMSE,ACC,MAE)的计算方法。 1.计算方法 2计算示例 参考：https://github.com/xiazh18/WeatherBench/blob/master/src/score.py 文章来源于微信公众号：自学气象人 ...

本文作者：王静瑶 联系方式：krisscy@163.com 欢迎大家来到WRF新手村特辑（四） 接上节，我们已经学习了namelist.wps的相关知识，此处将直接依据教学官网（WRF-ARW Online Tutorial）给出的参数进行设置，重合知识点不再赘述。 nano指令进入namelist.wps进行设置: 按照官网步骤，设置完时间就可以运行ungrib程序（也 ...

WRF模式模拟中，前处理WPS有许多有用的小工具，本篇文章列出几个常用的小工具，对一些特殊用途或者诊断错误或有帮助。 1.查看GRIB文件 ungrib.exe解压的是GRIB格式的数据，实际使用中，可能会遇到下载变量不全的情况，如果想要检验，可以使用以下工具打印grib文件的气象信息、变量、层次和日期，进行检验。 位于WPS/uti ...

大气海洋的特点，决定了我们无法做一些真实的实验，因此开展数值模拟，是其重要手段。业务预报中，现在气象预报员基本离不开模式的结果，甚至许多预报员毫不避讳，直言预报结论基本照搬模式结果。科研中，众多领域也是要需要使用数值模式，哪怕不使用数值模式，也需用到模式运行得到的再分析资料。因此对于大气和海洋科学 ...

文章概述 01 文章题目 利用降尺度网格气候数据提升气象指数保险的抗风险保障能力 02 摘要译文 极端气候事件是导致农业生产波动的重要原因，尤其是干旱事件会直接导致农作物产量的降低。因此，有效地应对干旱事件，例如种植耐旱作物、改善灌溉系统等，有助于减轻干旱风险，降低对农业生产的影响，保障农民的收入稳定。但是 ...

集合预报是一种气象预报方法，其中计算机模型不是只运行一次预测，而是多次从略有不同的起始条件开始运行。这种方法的目的是为了应对大气作为一个混沌系统所带来的挑战，即初始状态中的微小错误可能导致预报中的大错误。 图一 红色表示单一预报的不准确性，蓝色代表集合预报的多次采样。集合预报通过采样起始条件的不确定 ...

本文作者：王静瑶 联系方式：krisscy@163.com 欢迎大家来到WRF新手村特辑（七） 欢迎大家来到第七节课，上节已经学习了namelist.input的相关内容，这里将直接根据官网单域个例马修教学的参数进行设置： 设置完成，运行real.exe： 成功后你将得到wrfinput_d01和wrfbdy_d01文件。 之后进行最后一步，运行wrf.exe： 当 ...

WRF模式广泛应用于台风的模拟和预报，但是其并不能输出台风路径信息（除非编译的时候nesting选择vortex following选项，运行时会输出ATCF数据，不过这个不常用）。 因此使用WRF模拟台风，后处理需要解决如何从wrfout数据提取出台风路径的问题。NCL官网提供的脚本中，是以每个时刻海平面气压(SLP)最小处作为台风的中心位置 ...

先说结论 1、目前最推荐的方法是xesmf转换，插值方法多，自定义高，入门快 ，而且可导出权重文件进行复用 2、各种插值方式最常用的是线性或双线性插值，其他插值方法例如最邻近插值对边缘的处理一眼假，cubic慢 3、pyproj加scipy的griddata是第二推荐，进行pyproj投影转换后三种插值方法差别不明显，比之直接插值效果好 ...

使用数值模式往往会涉及到spin-up的问题，在气海同途的数值模式交流群中，有不少模式初学者经常感到疑惑。本篇谈谈自己的一些理解以及整理的总结。 由于难以得到理想的准确的初始场，模式初始场各变量之间，或者与外部条件之间存在着不平衡或者不匹配。因此模式启动后，需要运行一段时间，以达到所需的“平衡状态”。而达 ...

以1713号台风“天鸽”为例，绘制WRF模拟的台风路径和上海台风所的最佳路径做对比。为了更好的对比台风的移动路径和移动快慢，制作动图。 以下为基于Python的实现，重点说明如何绘制台风符号和制作动图。 一. 绘制台风符号 台风符号代码如下： 实际效果： 二. 制作动图 python制作动图使用imageio库，安装可使用如 ...

目前WRF自带的地形数据分辨率最高为30弧秒(约900m)，对于大部分研究和应用，此分辨率的地形数据基本能够满足。但是对于一些特殊应用，如超高分辨率的模拟或者大涡模拟，WRF 静态数据目前提供的地形数据分辨率稍显不足，因此此类应用需要解决如何添加更高分辨率的地形数据进入 WRF 的静态数据中。 SRTM(Shuttle Radar Topo ...

前言 实际应用中探空图可以分析所在区域的动热力特征，是预报员的好朋友 而在WRF应用中可以将其作为模式是否准确的检验工具 下面进行WRFOUT数据的探空图绘制 导入库 # 读数据 简单绘图 Text(0, 0.5, '压强 (hPa)') 缝合绘图 层结曲线（红色），露点曲线（绿色），LCL(黑点) CIN 蓝色阴影 CAPE 红色阴影 干绝热线（ ...

本文作者：王静瑶 联系方式：krisscy@163.com 欢迎大家来到WRF新手村特辑 我将带领大家使用WRF-ARW实现一次单域个例模拟，以WRF官网给出的单域个例飓风马修（指路链接：Single Model Domain case）为实训对象，带领大家转通一次WRF个例的同时把基础要点教给大家。我在刚接触WRF的时候也走了很多弯路（抹泪），希望游览完 ...

NCSS GFS 0.25 snowdepth 上次测试的是0.5分辨率的GFS数据 这次试试更高分辨率0.25的效果 实际上TDSCatalog还有很多产品，可以在这里看看 数据获取 In [8]: ['Best GFS Quarter Degree Forecast Time Series'] In [9]: 看看有啥变量 In [10]: 这得有上百个变量了 看看可获取的时间范围 In [4]: ------------------ ...

本文作者：王静瑶 联系方式：krisscy@163.com 欢迎大家来到WRF新手村特辑（二） 欢迎大家来到第二课，这节开始我们将进入实操部分。 首先大家需要学习一些linux基础命令： 接下来进入到WRF单域个例飓风马修的教学部分： （图为WRF-ARW Online Tutorial单域案例教学） 首先下载好教学官网中给出的飓风马修的 ...

这些NCEP FNL（最终）作战全球分析数据位于每六小时作战准备的1度乘1度网格上。该产品来自全球数据同化系统（GDAS），该系统持续收集来自全球电信系统（GTS）和其他来源的观测数据，用于许多分析。FNL是用NCEP在全球预报系统（GFS）中使用的相同模型制作的，但FNL是在GFS初始化后约一小时左右准备的。FNL被延迟，以便可以 ...

前言 很多时候我们需要拿模拟数据和站点图作对比，那就需要把模拟数据插值到站点 今天来尝试两种WRF数据插值到站点的方法并使用meteva进行简单绘图 方法一：xesmf库重插值后使用meteva进行双线性插值到站点 方法二：proj+scipy重插值后使用meteva进行最临近插值到站点 方法一 xesmf库重插值 In [3]: Out[3]: In [ ...

NCL在WRF模式的后处理中，有很多强大的函数以及现成易用的脚本，因此即使NCL官方不再继续更新，但NCL在WRF模式后处理中仍然不失为最合适的语言之一。有关NCL的安装，可以参考前面的文章，基于conda的安装。本文以NCL为例，介绍如何将WRF模拟结果插值到站点，包括特定的高度层和气压层。 一 .将u10, v10风插值到站点 二 ...