Changeset 8556

Timestamp:

2017-09-22T11:35:47+02:00 (7 years ago)

Author:

clem

Message:

bug fix

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rdgrft.F90 (modified) (6 diffs)

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

@@ -501 +501 @@
       !-----------------------------------------------------------------
 
-      zaksum(:) = apartf(:,0)
       zfac = 1._wp / hi_hrft
+      zaksum(1:nidx) = apartf(1:nidx,0)
       ! Transfer function
       DO jl = 1, jpl !all categories have a specific transfer function
@@ -534 +534 @@
       END WHERE
       
-      DO ji = 1, nidx
-         
+      DO ji = 1, nidx  
          ! correction to closing rate and opening if closing rate is excessive
          !---------------------------------------------------------------------
@@ -586 +585 @@
       REAL(wp), DIMENSION(jpij) ::   ersw             ! enth of water trapped into ridges
       REAL(wp), DIMENSION(jpij) ::   zswitch, fvol    ! new ridge volume going to jl2
-      !
-      REAL(wp), DIMENSION(jpij,jpl) ::   z1_ai        ! 1 / a
+      REAL(wp), DIMENSION(jpij) ::   z1_ai            ! 1 / a
       !
       REAL(wp), DIMENSION(jpij,nlay_i) ::   eirft     ! ice energy of rafting ice
@@ -614 +612 @@
             IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN
 
-               z1_ai(ji,jl1) = 1._wp / a_i_2d(ji,jl1)
+               z1_ai(ji) = 1._wp / a_i_2d(ji,jl1)
 
                !--------------------------------------------------------------------
@@ -628 +626 @@
                ! 2.3) Compute ridging /rafting fractions, make sure afrdg <=1 
                !---------------------------------------------------------------
-               afrdg = airdg1 * z1_ai(ji,jl1)
-               afrft = airft1 * z1_ai(ji,jl1)
+               afrdg = airdg1 * z1_ai(ji)
+               afrft = airft1 * z1_ai(ji)
 
                ! volume and enthalpy (J/m2, >0) of seawater trapped into ridges
@@ -721 +719 @@
                IF( apartf(ji,jl1) > 0._wp .AND. closing_gross(ji) > 0._wp ) THEN
                   ! Compute ridging /rafting fractions
-                  afrdg = aridge(ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji,jl1)
-                  afrft = araft (ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji,jl1)
+                  afrdg = aridge(ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji)
+                  afrft = araft (ji,jl1) * closing_gross(ji) * rdt_ice * z1_ai(ji)
                   ! Compute ridging ice and new ridges for ei
                   eirdg(ji,jk) = ze_i_2d (ji,jk,jl1) * afrdg + ersw(ji) * r1_nlay_i