8 pwg, pwgi, pwr, psnowswe, pdg, pdzg, &
9 pwg_ini, pwgi_ini, pwr_ini, pswe_ini )
41 USE yomhook
,ONLY : lhook, dr_hook
42 USE parkind1
,ONLY : jprb
51 CHARACTER(LEN=*),
INTENT(IN) :: hisba
55 CHARACTER(LEN=*),
INTENT(IN) :: hsnow_isba
61 REAL,
DIMENSION(:,:),
INTENT(IN) :: pwg
62 REAL,
DIMENSION(:,:),
INTENT(IN) :: pwgi
63 REAL,
DIMENSION(:),
INTENT(IN) :: pwr
64 REAL,
DIMENSION(:,:),
INTENT(IN) :: psnowswe
65 REAL,
DIMENSION(:,:),
INTENT(IN) :: pdg
66 REAL,
DIMENSION(:,:),
INTENT(IN) :: pdzg
68 REAL,
DIMENSION(:),
INTENT(OUT) :: pwg_ini
69 REAL,
DIMENSION(:),
INTENT(OUT) :: pwgi_ini
70 REAL,
DIMENSION(:),
INTENT(OUT) :: pwr_ini
71 REAL,
DIMENSION(:),
INTENT(OUT) :: pswe_ini
75 INTEGER :: ini, inl, inls
77 REAL(KIND=JPRB) :: zhook_handle
81 IF (lhook) CALL dr_hook(
'ISBA_BUDGET_INIT',0,zhook_handle)
96 IF(dgei%LWATER_BUDGET)
THEN
105 pswe_ini(ji)=pswe_ini(ji)+psnowswe(ji,jl)
115 IF(pwg(ji,jl)/=xundef)
THEN
116 pwg_ini(ji)=pwg_ini(ji)+pwg(ji,jl)*pdzg(ji,jl)*xrholw
117 pwgi_ini(ji)=pwgi_ini(ji)+pwgi(ji,jl)*pdzg(ji,jl)*xrholw
122 pwg_ini(:)=pwg(:,2)*pdg(:,2)*xrholw
123 pwgi_ini(:)=pwgi(:,2)*pdg(:,2)*xrholw
125 pwg_ini(:)=pwg_ini(:)+pwg(:,3)*(pdg(:,3)-pdg(:,2))*xrholw
137 IF (lhook) CALL dr_hook(
'ISBA_BUDGET_INIT',1,zhook_handle)
subroutine isba_budget_init(DGEI, HISBA, HSNOW_ISBA, PWG, PWGI, PWR, PSNOWSWE, PDG, PDZG, PWG_INI, PWGI_INI, PWR_INI, PSWE_INI)