83 ( kinit,omsg,tpdom,tptfl,tpsit,psnow_a,pemp_a,paddterm,paddterm2)
92 INTEGER,
INTENT(in) :: &
94 CHARACTER(*),
INTENT(in) :: &
96 TYPE(t_dom),
DIMENSION(np),
INTENT(in) :: &
98 TYPE(t_tfl),
DIMENSION(np),
INTENT(in) :: &
100 TYPE(t_sit),
DIMENSION(nt,np),
INTENT(in) :: &
102 REAL,
INTENT(inout) :: &
104 REAL,
DIMENSION(np),
INTENT(in),
OPTIONAL :: &
107 REAL,
DIMENSION(np) :: &
110 zsnow_a, zemp_a, zdemp, zdmsnow, zaddterm_a
114 zsnow(:) = sum( tpsit(:,:)%fsi*tpsit(:,:)%rsn*tpsit(:,:)%hsn, dim=1 )
116 zemp(:) = tptfl(:)%wlo + tptfl(:)%wio
119 zdmsnow = ( zsnow_a - psnow_a) / dtt
120 zdemp = zemp_a - pemp_a
123 '--------------------------------------------------------------------'
124 WRITE(noutlu,*) omsg ,
' Snow Content :', zsnow_a
125 WRITE(noutlu,*)
' Change in snow content :', zdmsnow
126 WRITE(noutlu,*)
' Change in emp :', zdemp
127 WRITE(noutlu,*)
' BILAN DMSNOW-EMP :', zdmsnow + zdemp
129 IF ( present(paddterm) .AND. present(paddterm2) )
THEN
130 zaddterm_a =
glt_avg_r(tpdom, paddterm(:), 1)
131 IF (lwg)
WRITE(noutlu,*)
' Bilan snow content :', zdmsnow + zaddterm_a
132 zaddterm_a =
glt_avg_r(tpdom, paddterm2(:), 1)
133 IF(lwg)
WRITE(noutlu,*)
' Bilan emp :', zdemp + zaddterm_a
134 ELSE IF ( present(paddterm) )
THEN
135 zaddterm_a =
glt_avg_r(tpdom, paddterm(:), 1)
137 WRITE(noutlu,*)
' Terme additif :', zaddterm_a
138 WRITE(noutlu,*)
' BILAN DMSNOW-EMP+terme additif :', zdmsnow + zdemp + zaddterm_a
subroutine glt_updsnow_r(kinit, omsg, tpdom, tptfl, tpsit, psnow_a, pemp_a, paddterm, paddterm2)