Changeset 13574 for NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE

Timestamp:

2020-10-07T18:02:40+02:00 (4 years ago)

Author:

stefryn

Message:

corrected bug in shear calculation, updated testcase (example output figs for EVP and EAP included)

Location:

NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE

Files:

• icedyn_rdgrft.F90 (modified) (1 diff)
• icedyn_rhg_eap.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rdgrft.F90

@@ -760 +760 @@
       !                              !--------------------------------------------------!
          strength(:,:) = rn_pstar * SUM( v_i(:,:,:), dim=3 ) * EXP( -rn_crhg * ( 1._wp - SUM( a_i(:,:,:), dim=3 ) ) )
-         ismooth = 1
+!         ismooth = 1    ! original code
+         ismooth = 0    ! try for EAP stability
          !                           !--------------------------------------------------!
       ELSE                           ! Zero strength                                    !

NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg_eap.F90

@@ -391 +391 @@
                zdsT = ( zds(ji,jj  ) * e1e2f(ji,jj  ) + zds(ji-1,jj  ) * e1e2f(ji-1,jj  )  &
                   &   + zds(ji,jj-1) * e1e2f(ji,jj-1) + zds(ji-1,jj-1) * e1e2f(ji-1,jj-1)  &
-                  &   ) * 0.25_wp * r1_e1e2t(ji,jj)
+!                  &   ) * 0.25_wp * r1_e1e2t(ji,jj)       .25 was an error
+                  &   ) * r1_e1e2t(ji,jj)
                !WRITE(numout,*) 'shear output', ji, jj, zdsT
                
@@ -972 +973 @@
       ! Factor to maintain the same stress as in EVP (see Section 3)
       ! Can be set to 1 otherwise
-      zinvstressconviso = 1._wp/(1._wp+kfriction*kfriction)
+!      zinvstressconviso = 1._wp/(1._wp+kfriction*kfriction)
+      zinvstressconviso = 1._wp
  
       zinvsin = 1._wp/sin(2._wp*pphi) * zinvstressconviso