NCL-有沒有辦法從索引信息中減去相同的變量、不同的值（相同的文件）？

Question

我有一個 3D nc 變量 PP = (time, lat, lon) 一整年的累積降水量，我想計算每日值。

我想在示例中做一些鏈接，但找不到任何示例。

我試圖做一個循環，其中減法應該在 [index+1] 減去值 [index] 的值之間（如下所示）

t = f->time(:)   ;size = 365

pp = f->TOTALRAIN(:,:,:) ; time, lat, lon dimensions


do i = 0,dimsizes(t)-1

DailyPp = pp(i,:,:) - pp(i-1,:,:)

end do

但這種方式會產生錯誤。

我需要每天的值並且只將累積的降水量作為輸入信息，並且認為這樣的循環應該是最好的方法。

我怎樣才能從變量中進行計算？？？

我正在與 ncl 合作，但如果您有 cdo 或 nco 提議，歡迎他們。

Answer 1

您可以使用兩個 NCO shell 命令執行此操作

下面的命令將記錄上移一位

ncks -v TOTALRAIN --msa_usr_rdr -d time,1,364 -d time,364 in.nc shift.nc

現在使用 ncbo 找出差異

ncbo -v TOTALRAIN shift.nc in.nc diff.nc

現在 netcdf 文件 diff 應該包含您的每日 pp 此文件中 TOTALRAIN 的最終記錄將為零

...亨利

Answer 2

在 cdo 中，您可以使用分解累積變量

cdo deltat in.nc diff.nc 

這與此問題有關，請參閱那里的答案以獲取更多詳細信息，包括插入初始時間步長。