37 USE modi_vegtype_to_patch
40 USE yomhook
,ONLY : lhook, dr_hook
41 USE parkind1
,ONLY : jprb
48 INTEGER,
INTENT(IN) :: kpatch
49 REAL,
DIMENSION(:,:,:),
INTENT(IN) :: pvegtype_patch
50 REAL,
DIMENSION(:,:),
INTENT(IN) :: ppatch
51 REAL,
DIMENSION(:,:,:),
INTENT(IN) :: pfieldout
52 REAL,
DIMENSION(:,:,:),
INTENT(OUT) :: pw
60 REAL(KIND=JPRB) :: zhook_handle
65 IF (lhook) CALL dr_hook(
'VEGTYPE_GRID_TO_PATCH_GRID',0,zhook_handle)
67 DO jvegtype=1,nvegtype
69 DO jlayer=1,
SIZE(pw,2)
70 pw(:,jlayer,jpatch) = pw(:,jlayer,jpatch) &
71 + pvegtype_patch(:,jvegtype,jpatch) * pfieldout(:,jlayer,jvegtype)
77 DO jlayer=1,
SIZE(pw,2)
78 WHERE(ppatch(:,jpatch)==0.) pw(:,jlayer,jpatch) = xundef
81 WHERE( abs(pw-xundef)/xundef < 1.e-6 ) pw = xundef
82 IF (lhook) CALL dr_hook(
'VEGTYPE_GRID_TO_PATCH_GRID',1,zhook_handle)
integer function vegtype_to_patch(IVEGTYPE, INPATCH)
subroutine vegtype_grid_to_patch_grid(KPATCH, PVEGTYPE_PATCH, PPATCH, PFIELDOUT, PW)