- Timestamp:
- 2020-05-28T16:21:45+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/iceistate.F90
r12890 r12988 85 85 !! ** Steps : 1) Set initial surface and basal temperatures 86 86 !! 2) Recompute or read sea ice state variables 87 !! 3) Fill in the ice thickness distribution using gaussian 88 !! 4) Fill in space-dependent arrays for state variables 89 !! 5) snow-ice mass computation 90 !! 6) store before fields 87 !! 3) Fill in space-dependent arrays for state variables 88 !! 4) snow-ice mass computation 91 89 !! 92 90 !! ** Notes : o_i, t_su, t_s, t_i, sz_i must be filled everywhere, even … … 284 282 ENDIF 285 283 286 !------------- !287 ! fill fields !288 !------------- !284 !----------------! 285 ! 3) fill fields ! 286 !----------------! 289 287 ! select ice covered grid points 290 288 npti = 0 ; nptidx(:) = 0 … … 393 391 t1_ice(:,:,:) = t_i (:,:,1,:) 394 392 ! 393 ! ice concentration should not exceed amax 394 at_i(:,:) = SUM( a_i, dim=3 ) 395 DO jl = 1, jpl 396 WHERE( at_i(:,:) > rn_amax_2d(:,:) ) a_i(:,:,jl) = a_i(:,:,jl) * rn_amax_2d(:,:) / at_i(:,:) 397 END DO 398 at_i(:,:) = SUM( a_i, dim=3 ) 399 ! 395 400 ENDIF ! ln_iceini 396 401 ! 397 at_i(:,:) = SUM( a_i, dim=3 )398 !399 402 !---------------------------------------------- 400 ! 3) Snow-ice mass (case ice is fully embedded)403 ! 4) Snow-ice mass (case ice is fully embedded) 401 404 !---------------------------------------------- 402 405 snwice_mass (:,:) = tmask(:,:,1) * SUM( rhos * v_s(:,:,:) + rhoi * v_i(:,:,:), dim=3 ) ! snow+ice mass … … 448 451 ENDIF 449 452 ENDIF 450 451 !------------------------------------452 ! 4) store fields at before time-step453 !------------------------------------454 ! it is only necessary for the 1st interpolation by Agrif455 a_i_b (:,:,:) = a_i (:,:,:)456 e_i_b (:,:,:,:) = e_i (:,:,:,:)457 v_i_b (:,:,:) = v_i (:,:,:)458 v_s_b (:,:,:) = v_s (:,:,:)459 e_s_b (:,:,:,:) = e_s (:,:,:,:)460 sv_i_b (:,:,:) = sv_i (:,:,:)461 u_ice_b(:,:) = u_ice(:,:)462 v_ice_b(:,:) = v_ice(:,:)463 ! total concentration is needed for Lupkes parameterizations464 at_i_b (:,:) = at_i (:,:)465 453 466 454 !!clem: output of initial state should be written here but it is impossible because
Note: See TracChangeset
for help on using the changeset viewer.