Changeset 13601 for NEMO/trunk/src/ICE/icestp.F90

Timestamp:

2020-10-14T17:59:34+02:00 (4 years ago)

Author:

clem

Message:

trunk: rewrite heat budget of sea ice to make it perfectly conservative by construction. Also, activating ln_icediachk now gives an ascii file (icedrift_diagnostics.ascii) containing mass, salt and heat global conservation issues (if any). In addition, 2D drift files can be outputed (icedrift_mass…) and field_def is changed accordingly. Note that advection schemes are not yet commited since there is still a restartability issue that I do not understand.

File:

• NEMO/trunk/src/ICE/icestp.F90 (modified) (4 diffs)

NEMO/trunk/src/ICE/icestp.F90

@@ -197 +197 @@
          IF( ln_icediahsb )             CALL ice_dia( kt )            ! -- Diagnostics outputs 
          !
+         IF( ln_icediachk )             CALL ice_drift_wri( kt )      ! -- Diagnostics outputs for conservation 
+         !
                                         CALL ice_wri( kt )            ! -- Ice outputs 
          !
@@ -281 +283 @@
       !
       CALL ice_dia_init                ! initialization for diags
+      !
+      CALL ice_drift_init              ! initialization for diags of conservation
       !
       fr_i  (:,:)   = at_i(:,:)        ! initialisation of sea-ice fraction
@@ -341 +345 @@
       ENDIF
       !
-      IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('par_init: online conservation check does not work with BDY')
-      !
       rDt_ice   = REAL(nn_fsbc) * rn_Dt          !--- sea-ice timestep and its inverse
       r1_Dt_ice = 1._wp / rDt_ice
@@ -438 +440 @@
       diag_trp_ei(:,:) = 0._wp   ;   diag_trp_es(:,:) = 0._wp
       diag_trp_sv(:,:) = 0._wp
+      diag_adv_mass(:,:) = 0._wp
+      diag_adv_salt(:,:) = 0._wp
+      diag_adv_heat(:,:) = 0._wp
       
    END SUBROUTINE diag_set0