53 USE yomhook
,ONLY : lhook, dr_hook
54 USE parkind1
,ONLY : jprb
62 REAL,
DIMENSION(:,:),
INTENT(IN) :: pz
63 REAL,
DIMENSION(:,:),
INTENT(IN) :: pin
64 REAL,
INTENT(IN) :: ph
66 REAL,
DIMENSION(:) ,
INTENT(OUT) :: pout
71 REAL(KIND=JPRB) :: zhook_handle
78 IF (lhook) CALL dr_hook(
'INTERPOL_SBL',0,zhook_handle)
83 DO WHILE(any(pout(:)==xundef) .AND. ilevel/=
SIZE(pz,2))
87 WHERE(pout(:)==xundef .AND. pz(:,ilevel+1)>=10.)
90 pout(:)=pin(:,ilevel) + &
91 (pin(:,ilevel+1)-pin(:,ilevel)) * &
92 (ph-pz(:,ilevel)) / (pz(:,ilevel+1)-pz(:,ilevel))
97 IF (lhook) CALL dr_hook(
'INTERPOL_SBL',1,zhook_handle)
subroutine interpol_sbl(PZ, PIN, PH, POUT)