- Timestamp:
- 2018-12-18T16:57:34+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icedyn.F90
r10399 r10411 218 218 ! ! -- check h_i -- ! 219 219 ! if h_i is larger than the surrounding 9 pts => reduce h_i and increase a_i 220 zhi = v_i 220 zhi = v_i(ji,jj,jl) / a_i(ji,jj,jl) 221 221 !!clem zdv = v_i(ji,jj,jl) - v_i_b(ji,jj,jl) 222 222 !!clem IF ( ( zdv > 0.0 .AND. zh > phmax(ji,jj,jl) .AND. at_i_b(ji,jj) < 0.80 ) .OR. & … … 228 228 ! ! -- check h_s -- ! 229 229 ! if h_s is larger than the surrounding 9 pts => put the snow excess in the ocean 230 zhs = v_s 230 zhs = v_s(ji,jj,jl) / a_i(ji,jj,jl) 231 231 IF( v_s(ji,jj,jl) > 0._wp .AND. zhs > phs_max(ji,jj,jl) .AND. a_i(ji,jj,jl) < 0.15 ) THEN 232 232 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.