Hi Martin, I don't see clearly how balanceEmPmR & balanceQnet are working: In do_oceanic_phys.F, the calls that remove the mean are after CALL EXTERNAL_FORCING_SURF (which compute surfaceForcingT & surfaceForcingS). Am I missing something ? Jean-Michel