Ignore:
Timestamp:
2018-04-03T17:33:37+02:00 (3 years ago)
Author:
clem
Message:

make sure that agrif is working with sea ice in any conditions

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/NST_SRC/agrif_lim3_update.F90

    r9167 r9454  
    2727   USE ice 
    2828   USE agrif_ice  
     29   USE phycst , ONLY: rt0 
    2930 
    3031   IMPLICIT NONE 
    … …  
    158159         ! 
    159160         ato_i(i1:i2,j1:j2) = 1._wp - at_i(i1:i2,j1:j2)         ! open water fraction   
     161 
     162         DO jl = 1, jpl 
     163            WHERE( tmask(i1:i2,j1:j2,1) == 0._wp )   t_su(i1:i2,j1:j2,jl) = rt0   ! to avoid a division by 0 in sbcblk.F90 
     164         END DO 
    160165          
    161166      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.