Changeset 10581

Timestamp:

2019-01-25T13:45:45+01:00 (5 years ago)

Author:

clem

Message:

add a comment

File:

• NEMO/trunk/src/ICE/icectl.F90 (modified) (2 diffs)

NEMO/trunk/src/ICE/icectl.F90

@@ -186 +186 @@
       !!-------------------------------------------------------------------
       CHARACTER(len=*), INTENT(in)    :: cd_routine    ! name of the routine
-      REAL(wp)                        :: zhfx, zsfx, zvfx
+      REAL(wp)                        :: zqmass, zhfx, zsfx, zvfx
       REAL(wp)                        :: zarea, zv_sill, zs_sill, zt_sill
       REAL(wp), PARAMETER             :: zconv = 1.e-9 ! convert W to GW and kg to Mt
@@ -199 +199 @@
       ! heat flux
       ! clem: not the good formulation
-      !!zhfx  = glob_sum( 'icectl', ( qt_atm_oi - qt_oce_ai - diag_heat ) * e1e2t ) * zconv
+      !!zhfx  = glob_sum( 'icectl', ( qt_oce_ai - qt_atm_oi + diag_heat + hfx_thd + hfx_dyn + hfx_res + hfx_sub + hfx_spr  &
+      !!   &                        ) * e1e2t ) * zconv
 
       ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice)