Changeset 13666


Ignore:
Timestamp:
2020-10-23T16:33:46+02:00 (6 months ago)
Author:
dancopsey
Message:

Move the point where I correct the total ice fraction

Location:
NEMO/branches/UKMO/NEMO_4.0.1_fix_nans/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_fix_nans/src/ICE/icesbc.F90

    r11715 r13666  
    138138      alb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:) 
    139139      ! 
     140 
     141      ! Apply a limit to the ice fraction to avoid NANs being generated in ice_thd_snwblow 
     142      ! (negative number to power of 0.66)  
     143      WHERE( at_i_b(:,:) > 1._wp ) at_i_b(:,:) = 1._wp 
     144 
    140145      ! 
    141146      SELECT CASE( ksbc )   !== fluxes over sea ice ==! 
  • NEMO/branches/UKMO/NEMO_4.0.1_fix_nans/src/OCE/SBC/sbccpl.F90

    r13665 r13666  
    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  
    16541649      ziceld(:,:) = 1._wp - picefr(:,:) 
    16551650      zcptn (:,:) = rcp * sst_m(:,:) 
Note: See TracChangeset for help on using the changeset viewer.