Changeset 8486 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90

Timestamp:

2017-09-01T15:49:35+02:00 (7 years ago)

Author:

clem

Message:

changes in style - part1 - (now the code looks better txs to Gurvan's comments)

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90 (modified) (8 diffs)

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

@@ -15 +15 @@
 #if defined key_lim3
    !!----------------------------------------------------------------------
-   !!   'key_lim3'                                      LIM3 sea-ice model
+   !!   'key_lim3'                                       LIM3 sea-ice model
    !!----------------------------------------------------------------------
    !!   ice_thd_ent   : ice redistribution of enthalpy
@@ -36 +36 @@
 
    !!----------------------------------------------------------------------
-   !! NEMO/LIM3 4.0 , UCL - NEMO Consortium (2011)
+   !! NEMO/ICE 4.0 , NEMO Consortium (2017)
    !! $Id: icethd_ent.F90 8420 2017-08-08 12:18:46Z clem $
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -86 +86 @@
             zeh_cum0(ji,jk0) = zeh_cum0(ji,jk0-1) + eh_i_old(ji,jk0-1)
             zh_cum0 (ji,jk0) = zh_cum0 (ji,jk0-1) + h_i_old (ji,jk0-1)
-         ENDDO
-      ENDDO
+         END DO
+      END DO
 
       !------------------------------------
@@ -95 +95 @@
       DO ji = 1, nidx
          zhnew(ji) = SUM( h_i_old(ji,0:nlay_i+1) ) * r1_nlay_i  
-      ENDDO
+      END DO
 
       ! new layers interfaces
@@ -102 +102 @@
          DO ji = 1, nidx
             zh_cum1(ji,jk1) = zh_cum1(ji,jk1-1) + zhnew(ji)
-         ENDDO
-      ENDDO
+         END DO
+      END DO
 
       zeh_cum1(:,0:nlay_i) = 0._wp 
@@ -115 +115 @@
                      &             / ( zh_cum0(ji,jk0) - zh_cum0(ji,jk0-1) )
                ENDIF
-            ENDDO
-         ENDDO
-      ENDDO
+            END DO
+         END DO
+      END DO
       ! to ensure that total heat content is strictly conserved, set:
       zeh_cum1(:,nlay_i) = zeh_cum0(:,nlay_i+2) 
@@ -126 +126 @@
             rswitch      = MAX( 0._wp , SIGN( 1._wp , zhnew(ji) - epsi20 ) ) 
             qnew(ji,jk1) = rswitch * ( zeh_cum1(ji,jk1) - zeh_cum1(ji,jk1-1) ) / MAX( zhnew(ji), epsi20 )
-         ENDDO
-      ENDDO
+         END DO
+      END DO
 
       ! --- diag error on heat remapping --- !
@@ -143 +143 @@
    !!   Default option                               NO  LIM3 sea-ice model
    !!----------------------------------------------------------------------
-CONTAINS
-   SUBROUTINE ice_thd_ent          ! Empty routine
-   END SUBROUTINE ice_thd_ent
 #endif