Changeset 13665


Ignore:
Timestamp:
2020-10-23T15:30:58+02:00 (6 months ago)
Author:
dancopsey
Message:

Make sure ice fraction stays below 1.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_fix_nans/src/OCE/SBC/sbccpl.F90

    r11715 r13665  
    16471647      ! 
    16481648      IF( ln_mixcpl )   zmsk(:,:) = 1. - xcplmask(:,:,0) 
     1649 
     1650      ! Apply a limit to the ice fraction to avoid NANs being generated in ice_thd_snwblow 
     1651      ! (negative number to power of 0.66)  
     1652      WHERE( picefr(:,:) > 1._wp ) picefr(:,:) = 1._wp 
     1653 
    16491654      ziceld(:,:) = 1._wp - picefr(:,:) 
    16501655      zcptn (:,:) = rcp * sst_m(:,:) 
Note: See TracChangeset for help on using the changeset viewer.