problem with interpcdo density

for some corrections we had to do to the official omips runs, I am redoing most of the diagnostics for the new formatted omips runs.

I try to interpolate to regular grid the diagnostics density.

The calculation of the density goes smooth, nevertheless the interpolation fails.

I noticed that the generated variable is: float sigma0(time, lev, dim2, dim3) as you can check here:

/esarchive/exp/ecearth/omip1/diags/OMIP/EC-Earth-Consortium/EC-Earth3/omip1/r1i1p1f1/Omon/sigma0/gn/v20200310

while in the past for the old omip, it was: float sigma0(time, lev, j, i) ;, as you can check here:

/esarchive/exp/ecearth/a2l5_sixth/diags/OMIP/EC-Earth-Consortium/EC-Earth3/omip1/r1i1p1f1/Omon/sigma0/gn/v20200316

the error at the interpolation level can be found here:

/esarchive/scratch/Earth/vsicardi/temporal_earth/earthdiags_suite/suite_home/omip1_diagnostics/19480101/Member_0/Chunks_1

thank you