96 LOGICAL,
DIMENSION(nt,np),
INTENT(in) :: &
98 TYPE(t_mxl),
DIMENSION(np),
INTENT(in) :: &
100 TYPE(t_sit),
DIMENSION(nt,np),
INTENT(inout) :: &
102 TYPE(t_tfl),
DIMENSION(np),
INTENT(inout) :: &
107 REAL,
DIMENSION(np) :: &
109 REAL,
DIMENSION(nt,np) :: &
110 zdssi,zssieq,zqsalt2,zssi
118 zssi(:,:) = tpsit(:,:)%ssi
131 zssieq(:,:) = ssisummer0/ssw0 * spread( tpmxl(:)%sml,1,nt )
132 WHERE( tpsit(:,:)%ssi > zssieq(:,:) .AND. gsmelt(:,:) )
134 ( zssieq(:,:) - tpsit(:,:)%ssi )*dtt / ( ssisummer_ts * xday2sec )
135 tpsit(:,:)%ssi = max( tpsit(:,:)%ssi + zdssi(:,:), zssieq(:,:) )
141 zssieq(:,:) = ssiwinter0/ssw0 * spread( tpmxl(:)%sml,1,nt )
142 WHERE( tpsit(:,:)%ssi > zssieq(:,:) .AND. &
143 tpsit(:,:)%tsf < spread( tpmxl(:)%mlf,1,nt ) )
145 ( zssieq(:,:) - tpsit(:,:)%ssi )*dtt / ( ssiwinter_ts * xday2sec )
146 tpsit(:,:)%ssi = max( tpsit(:,:)%ssi + zdssi(:,:), zssieq(:,:) )
148 zdssi(:,:) = tpsit(:,:)%ssi - zssi(:,:)
157 zqsalt2(:,:) = -rhoice*tpsit(:,:)%fsi*tpsit(:,:)%hsi*zdssi(:,:)
158 zqsalt(:) = sum( zqsalt2(:,:), dim=1 ) / dtt
162 IF ( nleviti==0 )
THEN
163 tptfl%wio = tptfl%wio - 1.e-3*zqsalt(:)
172 tptfl(:)%sio = tptfl(:)%sio + 1.e-3 * zqsalt(:)
subroutine glt_updsal_r(gsmelt, tpmxl, tpsit, tptfl)
subroutine glt_salflx_r(pqsalt, tpmxl, tptfl, pdmass, psalt)