While trying to understand the part of the code that manages surface variables (in processing_surf_vars.py), I noticed a potentially inoperative if statement in get_old_4d on line 125 :
if sur_var.ttr.values[i_t, 2, 2] != sur_var.ttr.values[i_t, 2, 2]: