95 REAL,
DIMENSION(nilay,nt,np),
INTENT(in) :: &
97 REAL,
DIMENSION(nt,np),
INTENT(inout) :: &
99 TYPE(t_vtp),
DIMENSION(nl,nt,np),
INTENT(inout) :: &
104 REAL :: zavti,zavtf,zdavt
105 REAL,
DIMENSION(nilay) :: &
107 REAL,
DIMENSION(nilay+1) :: &
109 REAL,
DIMENSION(nt,np) :: &
118 zhsi(:,:) = sum( pdhi(:,:,:),dim=1 )
122 IF ( zhsi(jk,jp)>0. .AND. zhsi(jk,jp)<phsi(jk,jp) )
THEN
124 zsf3tinvo=pdhi(:,jk,jp) / zhsi(jk,jp)
127 zlevo(jl)=zlevo(jl-1)+pdhi(jl-1,jk,jp)
129 zlevo(:)=zlevo/zlevo(nilay+1)
131 zentn=
glt_interpz( height,tpsil(1:nilay,jk,jp)%ent,zlevo )
134 zavti=sum( zsf3tinvo*tpsil(1:nilay,jk,jp)%ent )
135 zavtf=sum( sf3tinv*zentn )
136 zdavt=abs(zavtf-zavti)
137 IF ( zdavt>epsil5 )
THEN
138 print*,
'PB in mltvtp at jp,jk =',jp,jk
139 print*,
'Average temperature is not conserved.'
140 print*,
' Initial =',zavti
141 print*,
' Final =',zavtf
144 tpsil(1:nilay,jk,jp)%ent = zentn
151 phsi(:,:) = zhsi(:,:)
subroutine glt_mltvtp_r(pdhi, phsi, tpsil)