cost733 nc运行问题

ゼヅ欢チニ

各位大佬，文字描述较多，请担待
一、core dumped错误
半个月前下载的ERA5数据，nc格式。采用ncatted -a units,time,o,c,"hours since 1900-01-01 00:00:00" infile.nc outfile.nc去掉时间中的小数点，然后使用简单的./cost733class -dat pth:/data/data/ERA/year22.nc slo:60:130:0.25 sla:20:80:0.25 -met PCT -ncl 9 -cla ../MYPCT22.txt可以成功。但前天，重新下载的22年的数据（时间、经纬度范围、时空分辨率均完全相同），同样去掉小数点后，同样用上述cost733运行，结果报错core dumped，但直接用半月前的数据，仍能正确运行。请问如何解决问题
二、PCT和SOM并行运行问题。
昨日尝试了下，在a.sh中写入
../cost733class -dat pth:/data/cost733/data/year22.nc var:u10 wgt:1.D0 -dat pth:/data/cost733/data/year22.nc var:v10 wgt:1.D0 -dat pth:/data/cost733/data/year22.nc var:sp wgt:1.D0 slo:60:130:0.25 sla:20:80:0.25 -met PCT -ncl 9 -cla ./pct2022.txt &../cost733class -dat pth:/data/cost733/data/year23.nc var:u10 wgt:1.D0 -dat pth:/data/cost733/data/year23.nc var:v10 wgt:1.D0 -dat pth:/data/cost733/data/year23.nc var:sp wgt:1.D0 slo:60:130:0.25 sla:20:80:0.25 -met PCT -ncl 9 -cla ./pct2023.txt &然后  csh a.sh发现PCT无法并行运行（top发现初始有2个cost733，几秒钟就只有1个cost733程序了），但SOM可以并行运行（top一直都有2个cost733）。但今日重新运行其余年份，发现PCT和SOM均无法并行运行？请问为啥

SunJiaming

ncrcat -v time_bnds-x infile.nc outfile.nc
把time_bnd删了试试，还有检查一下文件的经纬度范围是不是超出范围，时间步长最好是12/6/3/2/1小时这样，如果不行，那就用cdo求个daymean再试试

ゼヅ欢チニ

并行哪里sh里面是两行代码，给出的两行都是PCT分型，SOM分型就只是-met后改为了SOM，其余未变化
../cost733class -dat pth:/data/cost733/data/year22.nc var:u10 wgt:1.D0 -dat pth:/data/cost733/data/year22.nc var:v10 wgt:1.D0 -dat pth:/data/cost733/data/year22.nc var:sp wgt:1.D0 slo:60:130:0.25 sla:20:80:0.25 -met PCT -ncl 9 -cla ./pct2022.txt &
../cost733class -dat pth:/data/cost733/data/year23.nc var:u10 wgt:1.D0 -dat pth:/data/cost733/data/year23.nc var:v10 wgt:1.D0 -dat pth:/data/cost733/data/year23.nc var:sp wgt:1.D0 slo:60:130:0.25 sla:20:80:0.25 -met PCT -ncl 9 -cla ./pct2023.txt &

SunJiaming

兄弟，第一个问题我怀疑是没有指定变量或者高度层，出现报错。还有可能是pth不是完整的文件路径，最好把根路径也加上试试，第二个问题可能同理。要不咱把文件放上来看看？
还有，SOM现在应该有比较成熟的python程序了，要不是因为不懂PCT的代码，才不和这个软件较劲呢

ゼヅ欢チニ

SunJiaming 发表于 2024-1-16 17:54
兄弟，第一个问题我怀疑是没有指定变量或者高度层，出现报错。还有可能是pth不是完整的文件路径，最好把根 ...

我都试过了     我下载只有地面气压、10米uv风速，都试过了    都不行   以前下载的都行 现在下载的都不行  不清楚为啥

SunJiaming

ゼヅ欢チニ 发表于 2024-1-16 20:11
我都试过了     我下载只有地面气压、10米uv风速，都试过了    都不行   以前下载的都行 现在下 ...

欸，我也不明白了，要不要重新下一遍数据看看，或者用命令行运行试试

ゼヅ欢チニ

SunJiaming 发表于 2024-1-17 11:41
欸，我也不明白了，要不要重新下一遍数据看看，或者用命令行运行试试

都试过了  命令行   服务器直接下数据    或者说电脑下数据上传到服务器  等等    都不行

ゼヅ欢チニ

SunJiaming 发表于 2024-1-18 11:24
ncrcat -v time_bnds-x infile.nc outfile.nc
把time_bnd删了试试，还有检查一下文件的经纬度范围是不是超 ...

哎，还是不行   大佬还有其他可能性吗？  再不行，我就重置服务器试试，但也很奇怪啊，早先下载的数据现在仍然可以的啊，  只是现在下载不行，  老火

SunJiaming

ゼヅ欢チニ 发表于 2024-1-23 15:41
哎，还是不行   大佬还有其他可能性吗？  再不行，我就重置服务器试试，但也很奇怪啊，早先下载的 ...

不至于重置服务器吧。。。要不咱也用anaconda装下cost733试试，或者换其他家的再分析数据做做也成啊

ゼヅ欢チニ

SunJiaming 发表于 2024-1-23 21:23
不至于重置服务器吧。。。要不咱也用anaconda装下cost733试试，或者换其他家的再分析数据做做也成啊

其实最开始我就试过用fnl（grib）来做cost，发现怎么都不对，倒是nc一下就可以了。而且服务器才开的，就一个cost，无所谓了{:lol:}