Changeset 15478 for NEMO/branches/UKMO/NEMO_4.0.4_penetrating_solar_no_neg

Timestamp:

2021-11-08T12:00:43+01:00 (3 years ago)

Author:

timgraham

Message:

Functional change to prevent negative values

File:

• NEMO/branches/UKMO/NEMO_4.0.4_penetrating_solar_no_neg/src/OCE/SBC/sbccpl.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.4_penetrating_solar_no_neg/src/OCE/SBC/sbccpl.F90

@@ -2131 +2131 @@
          CASE ('coupled')
             IF (ln_scale_ice_flux) THEN
-               WHERE( a_i(:,:,:) > 0.0_wp ) zqtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:)
+               WHERE( a_i(:,:,:) > 0.0_wp ) zqtr_ice_top(:,:,:) = MAX(0, frcv(jpr_qtr)%z3(:,:,:)) * a_i_last_couple(:,:,:) / a_i(:,:,:)
                WHERE( a_i(:,:,:) <= 0.0_wp ) zqtr_ice_top(:,:,:) = 0.0_wp
             ELSE
-               zqtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:)
+               zqtr_ice_top(:,:,:) = MAX(0, frcv(jpr_qtr)%z3(:,:,:))
             ENDIF