首页 > 其他分享 >【如何用Cost733进行天气环流形势分型】

【如何用Cost733进行天气环流形势分型】

时间:2024-09-20 13:20:37浏览次数:8  
标签:cost733 cost733class 分类 PCT ECV Cost733 环流 分型 2018

提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档

文章目录


一、cost733是什么?

简要介绍:cost733是一款由欧洲机构COST开发的天气分型软件,可用于气象工作者进行天气类型划分。

详细介绍
1.谁发明的:cost733是由COST(European Cooperation in Science and Technology)机构开发的一个项目,项目代号为733,它专注于天气分型的研究。COST是一个欧洲的国际研发合作网络,它通过科学研究项目促进欧洲科技的发展。
2.能做什么:cost733是使用Fortran语言编写的软件,它内置了一系列算法,这些算法能够根据历史和实时气象数据,将天气划分为不同的类型,目前cost733在linux系统上使用。
3.处理对象:不同天气类型对降水量、气压分布及空气质量的影响。
4.数学基础:通过cost733进行天气分型,可以采用正交旋转主成分分析法和斜交旋转主成分分析法,能够更准确地反映不同天气状况下的气压和风场特征。

二、cost733完整安装步骤

2.1 安装

下载cost733class−1.4安装包:https://projects.met.no/
-tar xvfj解压安装包 cd进入文件夹

2.2 配置

sh compile_gnu_debug_opengl.sh
./configureFC=gfortranCC=gcc
./configureFC=ifortCC=iccFCFLAGS="−parallel−openmp"
./configure−−disable−grib
./configure−−disable−opengl

2.3 编译

make
src/cost733class
sudo su −
make install
src/cost733class
sudo su−
make install
sudo copy src/cost733class /usr/local/bin

三、cost733使用(举例)

比如我要根据2018年夏季的小时700hpa高度场进行天气环流形势分型(如采用PCT方法;即一种时间模式下的主成分分析法; Hersbach et al, 2020)
3.1 预处理
首先进行数据预处理:
a)小时统一成每日数据,并作出每天的500/700hpa天气环流图,根据环流图剔除掉台风天气(台风会影响分型)
b)时间为列,站点为行,700hpa数据输入到test1.dat文件(同一日的二维网格点也输入到文件的同一维行中)

3.2 簇类计算
需要注意的是,PCT分类法默认将原始数据按照时间先分子集,cost733默认是10个时间单位一组(Hersbach et al, 2020)。也就是说,2018年夏季的90天数据最多分九组,下面对数据进行簇类计算,即选择分别分成2,3,…,9组,计算每组的ECV(解释簇类方差)
计算分2类:

cost733class -dat pth:test1.dat lon:15:90:0.25 lat:95:120:0.25 fdt:2018:6:1:0 ldt:2018:8:31:0 ddt:1d -met PCT -ncl 2 -cla 02.PCT -dcol 3  

#-pth:输入路径,lon,lat:经纬度,fdt,ldt,ddt:起止时间和间隔, -met:分类方法,-ncl:分类数, -cla:输出路径
cost733也可以输出ECV的值,这个值很重要,关系到最终的分类数量选择。
在这里插入图片描述

同理计算分3类:

cost733class -dat pth:test1.dat lon:15:90:0.25 lat:95:120:0.25 fdt:2018:6:1:0 ldt:2018:8:31:0 ddt:1d -met PCT -ncl 3 -cla 03.PCT -dcol 3  

与分2类比起来,分3类改变的只有-ncl 3 -cla 03.PCT这两个参数,此时也应输出分3类时的ECV值

同理计算分4类,5类,6类,7类,8类,9类:得到其他的这几个输出文件和ECV值

3.3 最终选择
得到了分2~9类时的ECV值,接下来我们就要计算ΔECV,即该ECV减去前一个ECV的值。比如分3类的ΔECV为:分3类的ECV减去分2类的ECV;注意分1类的ECV恒为0.

for example: 下图是某篇文章分5~15类时的ΔECV图(蓝色折线),当分类数为7时ΔECV最大,所以最终分类选分类数为7时。
在这里插入图片描述

回到原例子:计算出分2~9类的ΔECV后,我们选择最大的ΔECV对应的分类数为最终的分类数和分类结果,因为此时主成分(分类数)较少但是解释方差较大。比如分5类的ΔECV最大,那么原始数据分5类最好,05.PCT文件就是最终分类的输出结果!!

参考文献
Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., et al. (2020) The ERA5 global reanalysis. Quarterly Journal Of The Royal Meteorological Society, 146, 1999-2049. https://doi.org/10.1002/qj.3803.

标签:cost733,cost733class,分类,PCT,ECV,Cost733,环流,分型,2018
From: https://blog.csdn.net/weixin_53167226/article/details/142178034

相关文章

  • Nature Plants | 基因组所张兴坦团队开发无需参考基因组的单倍体分型挂载工具HapHiC
    “近日,《自然·植物(NaturePlants)》在线发表了中国农业科学院深圳农业基因组研究所(岭南现代农业科学与技术广东省实验室深圳分中心,以下简称“基因组所”)张兴坦团队联合南方科技大学陈国安副教授课题组、湖南农业大学易自力教授团队的研究论文,题为“Chromosome-levelscaffolding......
  • C#使用自己写的海龟绘图类绘制递归分型树
     paython里有个库,叫turtle,俗称海龟作图。自己看《程序员的数学》时,在递归章节看了用海龟绘图画递归树,我想用C#在winform上用gdi+去实现,我却没有在网上找到有适合C#语言的类似的库,我就自己简单写了一个海龟作图帮助类,实现了上图的图案。代码如下:publicclassTurtleHelper......
  • GB | 华中农大焦文标团队开发适用于植物基因组的基于图形的集成式分型算法
    今年4月,华中农业大学焦文标团队在GenomeBiology上发表论文:Acomprehensivebenchmarkofgraph‑basedgeneticvariantgenotypingalgorithmsonplantgenomesforcreatinganaccurateensemblepipeline,主要研究了基于图谱的植物基因组变异基因分型算法,并创建了一个准确的......
  • GB | 华中农大焦文标团队开发适用于植物基因组的基于图形的集成式分型算法
    今年4月,华中农业大学焦文标团队在GenomeBiology上发表论文:Acomprehensivebenchmarkofgraph‑basedgeneticvariantgenotypingalgorithmsonplantgenomesforcreatinganaccurateensemblepipeline,主要研究了基于图谱的植物基因组变异基因分型算法,并创建了一个准确的......
  • 【量化交易】顶底分型策略
    在众多的量化策略中,顶底分型策略因其独特的市场趋势捕捉能力和简洁的实现方式而受到许多投资者的青睐。本文将详细介绍顶底分型策略的原理,并展示如何使用Python在聚宽平台上实现这一策略。感兴趣的朋友,可以在下方公号内回复:001,即可获取源码,共同交流!顶底分型策略概述顶底分型是......
  • 【二分+容斥】【ST表/单调栈】【划分型DP】
    二分+容斥题目链接https://leetcode.cn/problems/kth-smallest-amount-with-single-denomination-combination/description/题目大意题目思路假设有一个x元硬币思考只有一种面额为3的硬币时,3可以组成不超过x的面额的数量有x/3种!思考有两种面额【3,6】,可以组成不......
  • 第二章投资技术《第四节 让走势画地为牢(分型,笔)》
    1.零部件大大小小的零部件组装一起变成机器股市中,对应到其中就是“关节”-》转折,连接的地方,关节两头的表现形式肯定不一样,甚至是相反的找到关节,分析规律与方式,才能游刃有余2.关节故事里的关节,就是分型指数不再创新高或新低,那么最高或者最低这根k与前后两根k组成了顶分型和底......
  • 龙哥量化:缠中说禅(缠论)分型
    声明:看到研究非常细致深入的文章,转载到我的博客园,以便学习和研究。(转载自知乎的DennisWin)分型是缠论K线系统中一个极为关键的概念,源自K线组合的一个完全分类,是一个纯理论的推导。最早提岀分型概念的是诺贝尔奖获得者一一美国的比尔•威廉姆斯博士,他在《混沌操作法》一书中,将分型......
  • 基因分型数据与碱基序列的输入
    基因分型数据和碱基序列的输入都是对DNA信息的编码,但它们的表达方式和所提供的信息不同。为了理解它们之间的联系,让我们首先明确这两者的定义:基因分型数据:基因分型数据通常是在特定的单核苷酸位置上(即SNP位置)对个体的DNA的描述。每个SNP位置可以有三种情况:两种纯合子和一种杂合......
  • 基因分型数据
    基因分型数据是对一个个体在特定的DNA位点或基因座的等位基因组成的记录。换句话说,基因分型是描述特定位置上DNA变化的方法。DNA和变异:DNA由四种碱基:腺嘌呤(A)、胸腺嘧啶(T)、鸟嘌呤(C)和鸟苷酸(G)组成。大部分人类的DNA序列是相同的,但某些位置上存在变异。这些变异点上的不同......