cost733超简单安装运行

SunJiaming

最近无意间发现有人将cost733class 1.4上传到了anaconda，可以仿照conda安装NCL的方式，安装并运行cost733class原网址：Cost733Class :: Anaconda.org
我的系统是Ubuntu 22.04 LTS

(1) 创建环境，命名为cost733，这步不用管python版本

1. conda create -n cost733

复制代码

(2) 创建好环境后，启动环境cost733并运行如下代码安装cost733

1. conda activate cost733
   
2. conda install -c conda-forge cost733class

复制代码

如果这步出现问题，最好换一下conda的源，不要用镜像源。
(3) 等待安装完成后，我们拿hgt.1999.nc做个测试，空间分辨率2.5度，时间分辨率6小时，算法使用PCT（T-mode PCA）

1. cost733class -dat pth:/home/sunjm/data/ncep_r2/daily/hgt.1999.nc var:hgt sle:500 slo:80:130:2.5 sla:10:60:2.5 -met PCT -ncl 4 -cla  output-PCT.txt  

复制代码

pth换成自己的文件地址，sle代表高度层，slo代表经度范围，sla代表纬度范围，2.5表示文件分辨率，-met表示算法名称（PCT），-ncl表示分型数量，-cla表示输出文件名称，自己起个名字就行，如果不指定输出位置，就会在根目录输出（我的是/home/sunjm/)
(4) output-PCT.txt中长这样，前4列为年月日小时，最后一列代表每个时次对应的天气类型

(5) 再计算一下ECV （explain cluster variance）

1. cost733class -dat pth:/home/sunjm/data/ncep_r2/daily/hgt.1999.nc  -clain pth:/home/sunjm/output-PCT.txt dtc:4 ddt:6h -met ECV -idx clain-ECV
   

复制代码

或者：

1. cost733class -dat pth:/home/sunjm/data/ncep_r2/daily/hgt.1999.nc var:hgt sle:500 slo:80:130:2.5 sla:10:60:2.5 -clain pth:/home/sunjm/output-PCT.txt dtc:4 ddt:6h -met ECV -idx clain-ECV

复制代码

-dat是原始数据的位置，-clain是输出结果的位置，dtc代表文件的日期数据有4列，这个必须要加；ddt代表时间分辨率6小时，-met代表计算ECV
(6) ECV最后会在屏幕输出，我的结果是0.0556

                               
登录/注册后可看大图

2024-03-31 更新
在计算PCT时，最好用var声明下变量名，不然cost733有可能识别不出来参与运算的变量而报错，并且对netcdf文件还有个要求，那就是nc文件的纬度（lat）只能从高纬度排序至低纬度（即从90°~-90°），但是在计算时sla的范围要求从小到大填写（如例子中10:60，而不能是60:10，否则也会报错）

AnneChen

点赞

464731123

太棒了！感谢楼主分享！

南极的七色鹿

你好，博主大大，这个分型数量是根据什么弄出来的哇

SunJiaming

南极的七色鹿 发表于 2024-3-7 12:31
你好，博主大大，这个分型数量是根据什么弄出来的哇

Ning, G., Yim, S.H.L., Wang, S. et al. Synergistic effects of synoptic weather patterns and topography on air quality: a case of the Sichuan Basin of China. Clim Dyn 53, 6729–6744 (2019). https://doi.org/10.1007/s00382-019-04954-3
可以参考下这篇论文，简单来说，每一种分型数都对应一个ECV，取前后两个ECV的差值最大时，对应的分形数就是最佳的分型

南极的七色鹿

SunJiaming 发表于 2024-3-7 14:53
Ning, G., Yim, S.H.L., Wang, S. et al. Synergistic effects of synoptic weather patterns and topogr ...

好嘞博主大大！超级感谢诶

SunJiaming

南极的七色鹿 发表于 2024-3-7 16:35
好嘞博主大大！超级感谢诶

客气客气~一起学习！

南极的七色鹿

SunJiaming 发表于 2024-3-7 22:01
客气客气~一起学习！

博主大大你好~可以麻烦再问您一个问题吗？就是如果我对某个省份进行天气分型，那么用PCT方法进行命令操作时，经纬度的范围是不是要比该省的实际经纬度更大点哇？

SunJiaming

南极的七色鹿 发表于 2024-3-7 23:43
博主大大你好~可以麻烦再问您一个问题吗？就是如果我对某个省份进行天气分型，那么用PCT方法进行命令操作 ...

最好能大一些，因为范围太小的话数据点太少，可能会报错，并且分型出来的结果可能太相似，没法分析

南极的七色鹿

SunJiaming 发表于 2024-3-8 10:16
最好能大一些，因为范围太小的话数据点太少，可能会报错，并且分型出来的结果可能太相似，没法分析

好嘞~太感谢博主大大啦