Changeset 12154 for NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/SBC/sbcisf.F90

Timestamp:

2019-12-10T15:44:23+01:00 (4 years ago)

Author:

cetlod

Message:

commit

File:

• NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/SBC/sbcisf.F90 (modified) (11 diffs)

NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/SBC/sbcisf.F90

@@ -303 +303 @@
       !
       ! Allocate public variable
-      IF ( sbc_isf_alloc()  /= 0 )         CALL ctl_stop( 'STOP', 'sbc_isf : unable to allocate arrays' )
+      IF( sbc_isf_alloc()  /= 0 )         CALL ctl_stop( 'STOP', 'sbc_isf : unable to allocate arrays' )
       !
       ! initialisation
@@ -440 +440 @@
             !! Initialize arrays to 0 (each step)
             zt_sum = 0.e0_wp
-            IF ( ik > 1 ) THEN
+            IF( ik > 1 ) THEN
                ! 1. -----------the average temperature between 200m and 600m ---------------------
                DO jk = misfkt(ji,jj),misfkb(ji,jj)
@@ -459 +459 @@
             ELSE
                qisf(ji,jj) = 0._wp   ;   fwfisf(ji,jj) = 0._wp
-            END IF
+            ENDIF
          END DO
       END DO
@@ -496 +496 @@
       ! coeficient for linearisation of potential tfreez
       ! Crude approximation for pressure (but commonly used)
-      IF ( l_useCT ) THEN   ! linearisation from Jourdain et al. (2017)
+      IF( l_useCT ) THEN   ! linearisation from Jourdain et al. (2017)
          zlamb1 =-0.0564_wp
          zlamb2 = 0.0773_wp
@@ -558 +558 @@
                   ! compute s freeze
                   zsfrz=(-zbqe-SQRT(zdis))*zaqer
-                  IF ( zsfrz < 0.0_wp ) zsfrz=(-zbqe+SQRT(zdis))*zaqer
+                  IF( zsfrz < 0.0_wp ) zsfrz=(-zbqe+SQRT(zdis))*zaqer
 
                   ! compute t freeze (eq. 22)
@@ -578 +578 @@
 
          ! define if we need to iterate (nn_gammablk 0/1 do not need iteration)
-         IF ( nn_gammablk <  2 ) THEN ; lit = .FALSE.
+         IF( nn_gammablk <  2 ) THEN ; lit = .FALSE.
          ELSE                           
             ! check total number of iteration
-            IF (nit >= 100) THEN ; CALL ctl_stop( 'STOP', 'sbc_isf_hol99 : too many iteration ...' )
+            IF(nit >= 100) THEN ; CALL ctl_stop( 'STOP', 'sbc_isf_hol99 : too many iteration ...' )
             ELSE                 ; nit = nit + 1
-            END IF
+            ENDIF
 
             ! compute error between 2 iterations
             ! if needed save gammat and compute zhtflx_b for next iteration
             zerr = MAXVAL(ABS(zhtflx-zhtflx_b))
-            IF ( zerr <= 0.01_wp ) THEN ; lit = .FALSE.
+            IF( zerr <= 0.01_wp ) THEN ; lit = .FALSE.
             ELSE                        ; zhtflx_b(:,:) = zhtflx(:,:)
-            END IF
-         END IF
+            ENDIF
+         ENDIF
       END DO
       !
@@ -718 +718 @@
                   pgt(ji,jj) = zustar(ji,jj) / (zgturb + zgmolet)
                   pgs(ji,jj) = zustar(ji,jj) / (zgturb + zgmoles)
-               END IF
+               ENDIF
             END DO
          END DO
@@ -757 +757 @@
                ! determine the deepest level influenced by the boundary layer
                DO jk = ikt+1, mbku(ji,jj)
-                  IF ( (SUM(e3u_n(ji,jj,ikt:jk-1)) < zhisf_tbl(ji,jj)) .AND. (umask(ji,jj,jk) == 1) ) ikb = jk
+                  IF( (SUM(e3u_n(ji,jj,ikt:jk-1)) < zhisf_tbl(ji,jj)) .AND. (umask(ji,jj,jk) == 1) ) ikb = jk
                END DO
                zhisf_tbl(ji,jj) = MIN(zhisf_tbl(ji,jj), SUM(e3u_n(ji,jj,ikt:ikb)))  ! limit the tbl to water thickness.
@@ -789 +789 @@
                ! determine the deepest level influenced by the boundary layer
                DO jk = ikt+1, mbkv(ji,jj)
-                  IF ( (SUM(e3v_n(ji,jj,ikt:jk-1)) < zhisf_tbl(ji,jj)) .AND. (vmask(ji,jj,jk) == 1) ) ikb = jk
+                  IF( (SUM(e3v_n(ji,jj,ikt:jk-1)) < zhisf_tbl(ji,jj)) .AND. (vmask(ji,jj,jk) == 1) ) ikb = jk
                END DO
                zhisf_tbl(ji,jj) = MIN(zhisf_tbl(ji,jj), SUM(e3v_n(ji,jj,ikt:ikb)))  ! limit the tbl to water thickness.
@@ -869 +869 @@
                ! determine the deepest level influenced by the boundary layer
                DO jk = ikt, mbkt(ji,jj)
-                  IF ( (SUM(e3t_n(ji,jj,ikt:jk-1)) .LT. rhisf_tbl(ji,jj)) .AND. (tmask(ji,jj,jk) == 1) ) ikb = jk
+                  IF( (SUM(e3t_n(ji,jj,ikt:jk-1)) .LT. rhisf_tbl(ji,jj)) .AND. (tmask(ji,jj,jk) == 1) ) ikb = jk
                END DO
                rhisf_tbl(ji,jj) = MIN(rhisf_tbl(ji,jj), SUM(e3t_n(ji,jj,ikt:ikb)))  ! limit the tbl to water thickness.
@@ -879 +879 @@
             END DO
          END DO
-      END IF 
+      ENDIF 
       !
       !==   ice shelf melting distributed over several levels   ==!