6 SUBROUTINE layer_e_budget( PT, PTSTEP, PIMPL, PHC, PTC, PD, PA, PB, PC, PY, PDQS )
73 USE modi_tridiag_ground
75 USE yomhook
,ONLY : lhook, dr_hook
76 USE parkind1
,ONLY : jprb
82 REAL,
DIMENSION(:,:),
INTENT(INOUT) :: pt
83 REAL,
INTENT(IN) :: ptstep
84 REAL,
DIMENSION(:,:),
INTENT(IN) :: phc
85 REAL,
DIMENSION(:,:),
INTENT(IN) :: ptc
87 REAL,
DIMENSION(:,:),
INTENT(IN) :: pd
88 REAL,
DIMENSION(:,:),
INTENT(IN) :: pa
89 REAL,
DIMENSION(:,:),
INTENT(IN) :: pb
90 REAL,
DIMENSION(:,:),
INTENT(IN) :: pc
91 REAL,
DIMENSION(:,:),
INTENT(IN) :: py
92 REAL,
DIMENSION(:),
INTENT(OUT) :: pdqs
93 REAL,
INTENT(IN) :: pimpl
101 REAL,
DIMENSION(SIZE(PT,1),SIZE(PT,2)) :: zx
103 REAL,
DIMENSION(SIZE(PT,1)) :: zei
104 REAL,
DIMENSION(SIZE(PT,1)) :: zpei
109 REAL(KIND=JPRB) :: zhook_handle
111 IF (lhook) CALL dr_hook(
'LAYER_E_BUDGET',0,zhook_handle)
126 zei(jj)=zei(jj) + ( phc(jj,jlayer)*pd(jj,jlayer)*pt(jj,jlayer) )
138 pt(:,jlayer) = zx(:,jlayer)
150 zpei(jj) = zpei(jj)+ ( phc(jj,jlayer)*pd(jj,jlayer)*pt(jj,jlayer) )
157 pdqs(:)=(zpei(:)-zei(:))/ptstep
159 IF (lhook) CALL dr_hook(
'LAYER_E_BUDGET',1,zhook_handle)
subroutine tridiag_ground(PA, PB, PC, PY, PX)
subroutine layer_e_budget(PT, PTSTEP, PIMPL, PHC, PTC, PD, PA, PB, PC, PY, PDQS)