如何利用cdo求r95p这个指标

drx

请问如何用cdo极端r95p这个极端指标，手册里的解释我实在是没看懂，cdo eca_r95p rrfile rrn95file ofile这个rrfile和rrn95file该怎么得到呢？我现在有的数据是1971-2010年的日尺度降水数据。 求解 提前感谢各位小伙伴[td]

[/td]

AuuuTi

rrfile 是你的降水数据(根据定义是大于等于1mm的降水需要被参与计算，需要预处理下 可以参考 setrtomiss,0,999 infile outfile)
rrn95file 是你的95分位值数据 利用 ydaypctl,95 来计算(手册里这一块也写了)
cdo ydaymin infile minfile
cdo ydaymax infile maxfile
cdo ydaypctl,95 infile minfile maxfile rrn95file

drx

AuuuTi 发表于 2023-1-11 21:06
rrfile 是你的降水数据(根据定义是大于等于1mm的降水需要被参与计算，需要预处理下 可以参考 setrtomiss,0, ...

我今天下午也是这么做的，但是不知道为什么最后算r95p的时候提示有错

AuuuTi

我记得ECA 和 ETCCDI对于“R95p” 这个指数定义有点不同，你可以算一下etccdi_r95p 再试下

drx

AuuuTi 发表于 2023-1-11 21:06
rrfile 是你的降水数据(根据定义是大于等于1mm的降水需要被参与计算，需要预处理下 可以参考 setrtomiss,0, ...

我详细描述下整个过程吧：
数据：01-pr_1971-2010_Jun-Aug.nc，也就是1971到2010年6月至8月的日尺度降水
第一步：cdo setrtomiss,0,1  01-pr_1971-2010_Jun-Aug.nc  02-pr_1971-2010_Jun-Aug_nan.nc
第二步：cdo ydaypctl,95  02-pr_1971-2010_Jun-Aug_nan.nc -ydaymin  02-pr_1971-2010_Jun-Aug_nan.nc  -ydaymax 02-pr_1971-2010_Jun-Aug_nan.nc  03-rrn95file.nc
第三步：cdo eca_r95p 02-pr_1971-2010_Jun-Aug_nan.nc 03-rrn95file.nc r95p.nc
显示有报错，不知道哪一步出了问题

AuuuTi

drx 发表于 2023-1-11 21:37
我详细描述下整个过程吧：
数据：01-pr_1971-2010_Jun-Aug.nc，也就是1971到2010年6月至8月的日尺度降水 ...

第二部分加个lock试试呢   cdo -L -ydaypctl,95  02-pr_1971-2010_Jun-Aug_nan.nc -ydaymin  02-pr_1971-2010_Jun-Aug_nan.nc  -ydaymax 02-pr_1971-2010_Jun-Aug_nan.nc  03-rrn95file.nc

drx

AuuuTi 发表于 2023-1-12 09:36
第二部分加个lock试试呢   cdo -L -ydaypctl,95  02-pr_1971-2010_Jun-Aug_nan.nc -ydaymin  02-pr_1971- ...

还是不行，最后提示我 “Input streams have different number of time steps!”，我把原始降水数据贴上了，您方便的话可否试一下，万分感谢

AuuuTi

drx 发表于 2023-1-13 15:51
还是不行，最后提示我 “Input streams have different number of time steps!”，我把原始降水数据贴上 ...

第二步：cdo ydaypctl,95  02-pr_1971-2010_Jun-Aug_nan.nc -ydaymin  02-pr_1971-2010_Jun-Aug_nan.nc  -ydaymax 02-pr_1971-2010_Jun-Aug_nan.nc  03-rrn95file.nc

把这个拆成三行来写吧，我看报的错是infile 时间步长不一致。你可以cdo showdate 或者 cdo sinfo 看看timestep是不是有偏差了

wwy1448909180

楼主怎么样解决了吗、我也是这样