cost733class

sorrymaker · 发表于 2022-9-30 23:14:50

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 sorrymaker 于 2022-10-10 16:32 编辑

cost733class天气分型软件安装及使用记录
本科是环境工程，去年上岸了，进了大气环境的组，但是学校没有大气的专业，研究很痛苦，由于做了天气分型这个方向，但是cost的攻略实在太少了，断断续续弄了一年，才勉强了解，所以放了一点自己的研究过程记录一下。
一.安装这个软件是在linux环境下运行，我是用的centos7，其他的不了解，安装包应该还是比较好找，Http://cost733class.geo.uni-augsburg.de/cost733class-1.2，官网是很早就停止运行和维护了，所以不好进。
需要C和Fortran

yum install gcc-c++

yum install gcc-gfortran

tar -zxvf cost733clas-1.4_rev23.tg正常的解压

./configure FC=gfortran CC=gcc--disable-jpeg --disable-grib（因为我报错了，所以禁用了grib）

make

make install

我还报了个错，需要automake1.15，我的版本过低，而且不好更新，因此我先给原先的删了，重新装的新版本automake（rpm -e automake如果提示有依赖，或者冲突，后面加上--nodeps --force即可强制卸载rpm -e --nodeps automake）

后面就没啥问题了。

sorrymaker · 发表于 2022-9-30 23:19:04

本帖最后由 sorrymaker 于 2022-9-30 23:27 编辑

二.运行
cost文件夹里面有一个slp.dat的文件可以测试装成功没有
2.1 数据
我先用的ncep的再分析数据（https://www.psl.noaa.gov/data/gridded/data.ncep.reanalysis2.html），每天四个时刻的数据，然后分型会报错，错误是找不到文件格式，我把他用ncl_convert2nc把ncep再转换一次，输出的文件就是????.nc.nc这个时候就可以读了。
命令
cost733class -v 4 -dat pth:mslp.????.nc.nc var:mslp fmt:netcdfslo:75:150:2.5 sla:0:60:2.5 fdt:2015 ldt:2020 ddt:1d -met PCT -cla ncep.cla -dcol3
2.2 数据处理
由于ncep的网格分辨率太大，我改用了ERA5的再分析数据，然后报错是时间格式出错。
5(BMX`NJ7L)6C6~{FE(R[YF.png

解决方案
多次尝试各种时间格式，直接上命令！利用nco进行ERA5文件的时间格式修改了，我尝试后是可以读的
ncatted -a units,time,o,c,"hours since 1900-01-01 00:00:00" infile.nc outfile.nc
我此时使用的是每四个时刻的文件，也就是分出来的型每天有四种
LLK$A9OQPQ7TAWKPOLO1V}N.png

前三列是时间，最后一列是类别
我发现单时刻的文件可读，由于分型的对象是海平面气压，我就下载了逐时的数据也就是24小时的，然后算出日均值，对日均值进行分型，ERA5是世界时需要和我的北京时对应，因此我下载了需要时间段的前一天的最后8小时以及最后一天的前16个小时数据。具体处理方法是用的cdo进行处理。
cdo timselmean,24 mslp_all.nc mslp_all_daily.nc
这样就对每24小时数据进行了均值计算，但是它同样对时间变量也进行了均值计算，相当于把24个时刻的数据对应到了一个时刻。这个时候报了错。
)W165B}]MB3FMGQ%T56PKJ3.png