Patrick, why do you have to store HEFFM in seaice_model.F (this is the land mask for the seaice model and should really be replaced by maskC(i,i, 1,bi,bj)). Doesn't this point to something strange? I just hope that is wasn't too difficult the fix the adjoint ... Martin