Changeset 13792 for NEMO

Timestamp:

2020-11-13T13:22:04+01:00 (3 years ago)

Author:

dancopsey

Message:

Merge in changes from NEMO4.0.1_fix_nans branch

File:

• NEMO/branches/UKMO/NEMO_4.0.3_fix_nans/src/ICE/icesbc.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.3_fix_nans/src/ICE/icesbc.F90

@@ -133 +133 @@
       CALL ice_alb( t_su, h_i, h_s, ln_pnd_alb, a_ip_eff, h_ip, cloud_fra, alb_ice )
 
+      ! Apply a limit to the ice fraction to avoid NANs being generated in ice_thd_snwblow
+      ! (negative number to power of 0.66) 
+      WHERE( at_i_b(:,:) > 1._wp ) at_i_b(:,:) = 1._wp
+
       !
       SELECT CASE( ksbc )   !== fluxes over sea ice ==!