Changeset 15763

Timestamp:

2022-03-18T15:03:54+01:00 (2 years ago)

Author:

mikebell

Message:

hpg_djr bug cffy = MAX (..) inserted

File:

• NEMO/branches/UKMO/NEMO_4.0-TRUNK_r14960_HPG/src/OCE/DYN/dynhpg.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0-TRUNK_r14960_HPG/src/OCE/DYN/dynhpg.F90

@@ -1698 +1698 @@
             zdz_i(ji,jj,jk)   = cffx / SIGN( MAX( ABS(z1_cff), zep ), z1_cff )
 
-            cffy = 2._wp * zdzy  (ji  ,jj-1) * zdzy  (ji,jj  )
+            cffy = MAX( 2._wp * zdzy  (ji  ,jj-1) * zdzy  (ji,jj), 0._wp)
             z1_cff = zdzy(ji,jj-1)   + zdzy(ji,jj)
             zdz_j(ji,jj,jk)   = cffy / SIGN( MAX( ABS(z1_cff), zep ), z1_cff )