Changeset 8239 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

Timestamp:

2017-06-28T17:55:50+02:00 (7 years ago)

Author:

clem

Message:

merge with v3_6_CMIP6_ice_diagnostics@r8238

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90 (modified) (2 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -734 +734 @@
       END DO  ! End of the do while iterative procedure
 
+      ! MV SIMIP 2016
+      !--- Snow-ice interfacial temperature (diagnostic SIMIP)
+      DO ji = kideb, kiut
+         zfac        = 1. / MAX( epsi10 , rn_cdsn * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji) )
+         t_si_1d(ji) = ( rn_cdsn        * zh_i(ji) * t_s_1d(ji,1) + &
+            &            ztcond_i(ji,1) * zh_s(ji) * t_i_1d(ji,1) ) * zfac
+      END DO
+      WHERE( ( zh_s < 1.0e-3 ) ) ; t_si_1d(:) = t_su_1d(:) ; END WHERE
+      ! END MV SIMIP 2016
+
       IF( ln_limctl .AND. lwp ) THEN
          WRITE(numout,*) ' zerritmax : ', zerritmax
@@ -751 +761 @@
          fc_bo_i(ji)     =  - zkappa_i(ji,nlay_i) * ( zg1*(t_bo_1d(ji) - t_i_1d(ji,nlay_i)) )
       END DO
+
+      ! MV SIMIP 2016
+      !--- Conduction fluxes (positive downwards)
+      diag_fc_bo_1d(:)   = diag_fc_bo_1d(:) + fc_bo_i(:) * a_i_1d(:) / at_i_1d(:)
+      diag_fc_su_1d(:)   = diag_fc_su_1d(:) + fc_su(:)   * a_i_1d(:) / at_i_1d(:)
+      ! END MV SIMIP 2016
 
       ! --- computes sea ice energy of melting compulsory for limthd_dh --- !