deseasonalization by substraction of ydaymean -> "Internal error, varIDs of input streams differ!"
Added by Simon Chabrillat over 1 year ago
Hello,
I need to de-seasonalize a daily timeseries by removal of the climatological annual cycle.
Following Adrian Tompkins' advice (found on StackOverflow, "How to remove seasonality from time series data?"),
I try the following on the attached dataset:
cdo -ydaymean -seldate,20050101,20191231 mp04e5_19960104-20200629.01dzm.lev0.nc ydaymean_2005-2019.nc
do -sub mp04e5_19960104-20200629.01dzm.lev0.nc ydaymean_2005-2019.nc deseas_2005-2019.nc
... but I get an internal error message:
cdo sub: Filling up stream2 >mp04e5.01dzm.lev0.ydaymean_2005-2019.nc< by copying all timesteps.
cdo sub (Abort): Internal error, varIDs of input streams differ!
I believe that my input dataset is correctly encoded as it results from many prior (and successful) processing steps by cdo.
I see that this is flagged as an "internal" error... Could I have encountered a bug inside cdo?
Thanks for your help,
Simon.
Replies (3)
RE: deseasonalization by substraction of ydaymean -> "Internal error, varIDs of input streams differ!" - Added by Simon Chabrillat over 1 year ago
Sorry you can discard this question, I think that I have found the issue.
RE: deseasonalization by substraction of ydaymean -> "Internal error, varIDs of input streams differ!" - Added by Karin Meier-Fleischer over 1 year ago
Hi Simon,
the second line of the error message points to the varIDs. I selected the variable ps explicitly before doing the calculations and it works.
cdo -ydaymean -seldate,20050101,20191231 -selname,ps mp04e5_19960104-20200629.01dzm.lev0.nc ps_ydaymean_2005-2019.nc cdo -sub -selname,ps mp04e5_19960104-20200629.01dzm.lev0.nc ps_ydaymean_2005-2019.nc ps_deseas_2005-2019.nc
RE: deseasonalization by substraction of ydaymean -> "Internal error, varIDs of input streams differ!" - Added by Simon Chabrillat over 1 year ago
Thanks a lot Karin, I had found already that I had a problem with variable selection but I had not pin-pointed it yet.