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(nx,ny),
INTENT(in) :: &
98 TYPE(t_tfl),
DIMENSION(nx,ny),
INTENT(in) :: &
100 TYPE(t_sit),
DIMENSION(nt,nx,ny),
INTENT(in) :: &
102 REAL,
INTENT(inout) :: &
104 REAL,
DIMENSION(nx,ny),
INTENT(in),
OPTIONAL :: &
107 REAL,
DIMENSION(nx,ny) :: &
110 zsnow_a, zemp_a, zdemp, zdmsnow, zaddterm_a
114 zsnow(:,:) = sum( tpsit(:,:,:)%fsi*tpsit(:,:,:)%rsn*tpsit(:,:,:)%hsn, dim=1 )
115 zsnow_a =
glt_avg(tpdom, zsnow(:,:), 1)
116 zemp(:,:) = tptfl(:,:)%wlo
117 zemp_a =
glt_avg(tpdom, zemp(:,:), 1)
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(tpdom, paddterm(:,:), 1)
131 IF(lwg)
WRITE(noutlu,*)
' Bilan snow content :', zdmsnow + zaddterm_a
132 zaddterm_a =
glt_avg(tpdom, paddterm2(:,:), 1)
133 IF(lwg)
WRITE(noutlu,*)
' Bilan emp :', zdemp + zaddterm_a
134 ELSE IF ( present(paddterm) )
THEN
135 zaddterm_a =
glt_avg(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(kinit, omsg, tpdom, tptfl, tpsit, psnow_a, pemp_a, paddterm, paddterm2)