Changeset 12915 for NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90

Timestamp:

2020-05-12T19:25:04+02:00 (4 years ago)

Author:

clem

Message:

implement a convergence check on heat diffusion scheme

File:

• NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90 (modified) (1 diff)

NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_zdf_bl99.F90

@@ -771 +771 @@
          
       END DO  ! End of the do while iterative procedure
-      
-      IF( ln_icectl .AND. lwp ) THEN
-         WRITE(numout,*) ' zdti_max : ', zdti_max
-         WRITE(numout,*) ' iconv    : ', iconv
+      !
+      ! convergence tests (only for output)
+      IF( ln_zdf_chkcvg ) THEN
+         tice_cvg_1d(1:npti) = 0._wp
+         DO jk = 1, nlay_s          ! snow temperature
+            DO ji = 1, npti
+               tice_cvg_1d(ji) = MAX( tice_cvg_1d(ji), ABS( t_s_1d(ji,jk) - ztsb(ji,jk) ) )
+            ENDDO
+         ENDDO
+         DO jk = 1, nlay_i          ! ice temperature
+            DO ji = 1, npti
+               tice_cvg_1d(ji) = MAX( tice_cvg_1d(ji), ABS( t_i_1d(ji,jk) - ztib(ji,jk) ) )
+            ENDDO
+         ENDDO
+         IF( k_cnd == np_cnd_OFF .OR. k_cnd == np_cnd_EMU ) THEN
+            DO ji = 1, npti         ! surface temperature
+               tice_cvg_1d(ji) = MAX( tice_cvg_1d(ji), ABS( t_su_1d(ji) - ztsub(ji) ) )
+            ENDDO
+         ENDIF
       ENDIF
-      
       !
       !-----------------------------