38 USE yomhook
,ONLY : lhook, dr_hook
39 USE parkind1
,ONLY : jprb
49 TYPE(isba_canopy_t
),
INTENT(INOUT) :: icp
50 TYPE(isba_grid_t
),
INTENT(INOUT) :: ig
54 REAL,
DIMENSION(:,:),
ALLOCATABLE :: zzf
55 REAL(KIND=JPRB) :: zhook_handle
63 IF (lhook) CALL dr_hook(
'PREP_ISBA_CANOPY',0,zhook_handle)
70 ALLOCATE(zzf(ig%NDIM,icp%NLVL))
78 ALLOCATE(icp%XZ(ig%NDIM,icp%NLVL))
79 DO jlayer=1,icp%NLVL-1
80 icp%XZ(:,jlayer) = 0.5 * (zzf(:,jlayer)+zzf(:,jlayer+1))
82 icp%XZ(:,icp%NLVL) = 1.5 * zzf(:,icp%NLVL) - 0.5 * zzf(:,icp%NLVL-1)
86 IF (lhook) CALL dr_hook(
'PREP_ISBA_CANOPY',1,zhook_handle)
subroutine prep_isba_canopy(ICP, IG)