Changeset 13666 for NEMO/branches/UKMO/NEMO_4.0.1_fix_nans/src/ICE/icesbc.F90

Timestamp:

2020-10-23T16:33:46+02:00 (4 years ago)

Author:

dancopsey

Message:

Move the point where I correct the total ice fraction

File:

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

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

@@ -138 +138 @@
       alb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:)
       !
+
+      ! 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 ==!