Changeset 15340

Timestamp:

2021-10-07T09:02:23+02:00 (3 years ago)

Author:

clem

Message:

little optimization of ice ridging/rafting

File:

• NEMO/trunk/src/ICE/icedyn_rdgrft.F90 (modified) (8 diffs)

NEMO/trunk/src/ICE/icedyn_rdgrft.F90

@@ -270 +270 @@
 
       ! controls
-      IF( sn_cfctl%l_prtctl )   CALL ice_prt3D   ('icedyn_rdgrft')                                                        ! prints
+      IF( sn_cfctl%l_prtctl )   CALL ice_prt3D('icedyn_rdgrft')                                                           ! prints
       IF( ln_icectl    )   CALL ice_prt     (kt, iiceprt, jiceprt,-1, ' - ice dyn rdgrft - ')                             ! prints
       IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icedyn_rdgrft', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation
@@ -518 +518 @@
       !
       INTEGER , DIMENSION(jpij) ::   itest_rdg, itest_rft   ! test for conservation
+      LOGICAL , DIMENSION(jpij) ::   ll_shift         ! logical for doing calculation or not
       !!-------------------------------------------------------------------
       !
@@ -538 +539 @@
          DO ji = 1, npti
 
-            IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN   ! only if ice is ridging
+            ! set logical to true when ridging
+            IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN   ;   ll_shift(ji) = .TRUE.
+            ELSE                                                                ;   ll_shift(ji) = .FALSE.
+            ENDIF
+            
+            IF( ll_shift(ji) ) THEN   ! only if ice is ridging
 
                IF( a_i_2d(ji,jl1) > epsi10 ) THEN   ;   z1_ai(ji) = 1._wp / a_i_2d(ji,jl1)
@@ -628 +634 @@
          DO jk = 1, nlay_s
             DO ji = 1, npti
-               IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN
+               IF( ll_shift(ji) ) THEN
                   ! Compute ridging /rafting fractions
                   afrdg = aridge(ji,jl1) * closing_gross(ji) * rDt_ice * z1_ai(ji)
@@ -649 +655 @@
          DO jk = 1, nlay_i
             DO ji = 1, npti
-               IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN
+               IF( ll_shift(ji) ) THEN
                   ! Compute ridging /rafting fractions
                   afrdg = aridge(ji,jl1) * closing_gross(ji) * rDt_ice * z1_ai(ji)
@@ -672 +678 @@
             DO ji = 1, npti
 
-               IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN
+               IF( ll_shift(ji) ) THEN
 
                   ! Compute the fraction of ridged ice area and volume going to thickness category jl2
@@ -729 +735 @@
             DO jk = 1, nlay_s
                DO ji = 1, npti
-                  IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp )   &
+                  IF( ll_shift(ji) )   &
                      &   ze_s_2d(ji,jk,jl2) = ze_s_2d(ji,jk,jl2) + ( esrdg(ji,jk) * rn_fsnwrdg * fvol(ji)  +  &
                      &                                               esrft(ji,jk) * rn_fsnwrft * zswitch(ji) )
@@ -738 +744 @@
             DO jk = 1, nlay_i
                DO ji = 1, npti
-                  IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp )   &
+                  IF( ll_shift(ji) )   &
                      &   ze_i_2d(ji,jk,jl2) = ze_i_2d(ji,jk,jl2) + eirdg(ji,jk) * fvol(ji) + eirft(ji,jk) * zswitch(ji)
                END DO