NCL读取ERA5数据的时间解析

一大碗年糕

NCL利用NCL读取ERA5数据时，无法对时间使用calendar_decode类的函数，报错显示：

1. fatal:Argument type mismatch on argument (0) of (calendar_decode2) can not coerce
   
2. fatal:["Execute.c":8637]:Execute: Error occurred at or near line 13 in file test.ncl

复制代码

原因在于，ERA5时间数据类型为int64，需要转为double才能被解析
https://mailman.ucar.edu/pipermail/ncl-talk/2024-August/019416.html
time_o = f1->valid_time
printVarSummary(time_o)
time_new = todouble(time_o)
copy_VarMeta(time_o,time_new)
printVarSummary(time_new)

重拾NCL中泪目...

一大碗年糕

fly100516 发表于 2025-1-25 10:30
我也遇到了这个问题，请问博主目前有没有好的解决方法，谢谢！

我目前也就是重新把时间time转化成float或double这种办法

再见惜别

我用的cd_calendar函数做时间解析，在调用函数之前，得转换成double型，然后还得给时间添加单位
time            = todouble(f4->valid_time)
time@units="seconds since 1970-01-01" ; 单位可以通过ncl_filedump命令查询
printVarSummary(time)

fly100516

我也遇到了这个问题，请问博主目前有没有好的解决方法，谢谢！

Riakira

一大碗年糕 发表于 2025-1-25 18:09
我目前也就是重新把时间time转化成float或double这种办法

为什么转化完后还是一样的报错呀

fly100516

一大碗年糕 发表于 2025-1-25 18:09
我目前也就是重新把时间time转化成float或double这种办法

好的，谢谢！我是用cdo转换，感觉era5更新了，更加麻烦了

风叶ele

fly100516 发表于 2025-2-9 00:08
好的，谢谢！我是用cdo转换，感觉era5更新了，更加麻烦了

啊还有这种事情？前两个月用NCL处理ERA5数据还一切正常...

jianglan

fly100516 发表于 2025-2-9 00:08
好的，谢谢！我是用cdo转换，感觉era5更新了，更加麻烦了

你好，想向您请教一下如何用cdo进行转换呀？

jianglan

再见惜别 发表于 2025-1-26 04:26
我用的cd_calendar函数做时间解析，在调用函数之前，得转换成double型，然后还得给时间添加单位
time      ...

那如果后续要处理某一个时间段的平均，是要把time重新写进去吗？

再见惜别

jianglan 发表于 2025-2-12 18:21
那如果后续要处理某一个时间段的平均，是要把time重新写进去吗？

不需要把。后续的程序，你之前怎么处理就怎么处理