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

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)

File:

• NEMO/branches/2019/dev_r11842_SI3-10_EAP/src/ICE/icedyn_rhg_eap.F90 (modified) (2 diffs)

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