Changeset 11832

Timestamp:

2019-10-29T18:21:12+01:00 (4 years ago)

Author:

dancopsey

Message:

Convert evap_ice from grid box means to sea ice means

File:

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

NEMO/branches/UKMO/NEMO_4.0_fix_cpl_oce_only/src/OCE/SBC/sbccpl.F90

@@ -1764 +1764 @@
          sprecip (:,:)   = zsprecip (:,:)
          tprecip (:,:)   = ztprecip (:,:)
-         evap_ice(:,:,:) = zevap_ice(:,:,:)
+         IF (ln_scale_ice_fluxes) THEN
+            ! Convert from grid box means to sea ice means
+            WHERE( a_i(:,:,:) > 0.0_wp ) evap_ice(:,:,:) = zevap_ice(:,:,:) / a_i(:,:,:)
+            WHERE( a_i(:,:,:) <= 0.0_wp ) evap_ice(:,:,:) = 0.0
+         ELSE
+            evap_ice(:,:,:) = zevap_ice(:,:,:)
+         ENDIF
          DO jl = 1, jpl
             devap_ice(:,:,jl) = zdevap_ice(:,:)