- Timestamp:
- 2017-09-27T11:11:01+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icedyn_adv_pra.F90
r8534 r8564 41 41 42 42 SUBROUTINE ice_dyn_adv_pra( kt, pu_ice, pv_ice, & 43 & pato_i, pv_i, pv_s, ps mv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i )43 & pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, pe_s, pe_i ) 44 44 !!---------------------------------------------------------------------- 45 45 !! ** routine ice_dyn_adv_pra ** … … 59 59 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_i ! ice volume 60 60 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pv_s ! snw volume 61 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: ps mv_i! salt content61 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: psv_i ! salt content 62 62 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: poa_i ! age content 63 63 REAL(wp), DIMENSION(:,:,:) , INTENT(inout) :: pa_i ! ice concentration … … 99 99 z0opw(:,:,1) = pato_i(:,:) * e1e2t(:,:) ! Open water area 100 100 DO jl = 1, jpl 101 z0snw(:,:,jl) = pv_s 102 z0ice(:,:,jl) = pv_i 103 z0ai (:,:,jl) = pa_i 104 z0smi(:,:,jl) = ps mv_i(:,:, jl) * e1e2t(:,:) ! Salt content105 z0oi (:,:,jl) = poa_i 106 z0es (:,:,jl) = pe_s 101 z0snw(:,:,jl) = pv_s (:,:, jl) * e1e2t(:,:) ! Snow volume 102 z0ice(:,:,jl) = pv_i (:,:, jl) * e1e2t(:,:) ! Ice volume 103 z0ai (:,:,jl) = pa_i (:,:, jl) * e1e2t(:,:) ! Ice area 104 z0smi(:,:,jl) = psv_i(:,:, jl) * e1e2t(:,:) ! Salt content 105 z0oi (:,:,jl) = poa_i(:,:, jl) * e1e2t(:,:) ! Age content 106 z0es (:,:,jl) = pe_s (:,:,1,jl) * e1e2t(:,:) ! Snow heat content 107 107 DO jk = 1, nlay_i 108 108 z0ei(:,:,jk,jl) = pe_i(:,:,jk,jl) * e1e2t(:,:) ! Ice heat content … … 227 227 pato_i(:,:) = z0opw(:,:,1) * r1_e1e2t(:,:) 228 228 DO jl = 1, jpl 229 pv_i 230 pv_s 231 ps mv_i(:,:, jl) = z0smi(:,:,jl) * r1_e1e2t(:,:)232 poa_i 233 pa_i 234 pe_s 229 pv_i (:,:, jl) = z0ice(:,:,jl) * r1_e1e2t(:,:) 230 pv_s (:,:, jl) = z0snw(:,:,jl) * r1_e1e2t(:,:) 231 psv_i(:,:, jl) = z0smi(:,:,jl) * r1_e1e2t(:,:) 232 poa_i(:,:, jl) = z0oi (:,:,jl) * r1_e1e2t(:,:) 233 pa_i (:,:, jl) = z0ai (:,:,jl) * r1_e1e2t(:,:) 234 pe_s (:,:,1,jl) = z0es (:,:,jl) * r1_e1e2t(:,:) 235 235 DO jk = 1, nlay_i 236 236 pe_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e1e2t(:,:)
Note: See TracChangeset
for help on using the changeset viewer.