New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10411 for NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icedyn.F90 – NEMO

Ignore:
Timestamp:
2018-12-18T16:57:34+01:00 (5 years ago)
Author:
clem
Message:

make all the sette tests successful

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icedyn.F90

    r10399 r10411  
    218218                  !                               ! -- check h_i -- ! 
    219219                  ! if h_i is larger than the surrounding 9 pts => reduce h_i and increase a_i 
    220                   zhi = v_i (ji,jj,jl) / a_i(ji,jj,jl) 
     220                  zhi = v_i(ji,jj,jl) / a_i(ji,jj,jl) 
    221221!!clem                  zdv = v_i(ji,jj,jl) - v_i_b(ji,jj,jl)   
    222222!!clem                  IF ( ( zdv >  0.0 .AND. zh > phmax(ji,jj,jl) .AND. at_i_b(ji,jj) < 0.80 ) .OR. & 
     
    228228                  !                               ! -- check h_s -- ! 
    229229                  ! if h_s is larger than the surrounding 9 pts => put the snow excess in the ocean 
    230                   zhs = v_s (ji,jj,jl) / a_i(ji,jj,jl) 
     230                  zhs = v_s(ji,jj,jl) / a_i(ji,jj,jl) 
    231231                  IF( v_s(ji,jj,jl) > 0._wp .AND. zhs > phs_max(ji,jj,jl) .AND. a_i(ji,jj,jl) < 0.15 ) THEN 
    232232                     zfra = a_i(ji,jj,jl) * phs_max(ji,jj,jl) / MAX( v_s(ji,jj,jl), epsi20 ) 
Note: See TracChangeset for help on using the changeset viewer.