diff --git a/src/SURFEX/get_vegn.F90 b/src/SURFEX/get_vegn.F90 index da43404278db6d1905c415c67c8f2979221d3785..28d2044a9316bab6a4c9cb4c94af3971610a9af7 100644 --- a/src/SURFEX/get_vegn.F90 +++ b/src/SURFEX/get_vegn.F90 @@ -69,7 +69,6 @@ REAL, DIMENSION(KI), INTENT(OUT) :: PLAI ! INTEGER :: JI,JJ ! loop index over tiles INTEGER :: ILUOUT ! unit numberi -REAL, DIMENSION(U%NSIZE_FULL) :: ZH_TREE_FULL, ZLAI_FULL REAL, DIMENSION(U%NSIZE_NATURE) :: ZH_TREE, ZLAI,ZWORK INTEGER:: IPATCH_TRBE, IPATCH_TRBD, IPATCH_TEBE, IPATCH_TEBD, IPATCH_TENE, & IPATCH_BOBD, IPATCH_BONE, IPATCH_BOND @@ -85,10 +84,6 @@ INTEGER:: IPATCH_TRBE, IPATCH_TRBD, IPATCH_TEBE, IPATCH_TEBD, IPATCH_TENE, & !* 1. Passage dur le masque global ! ------------------------------- - -ZH_TREE_FULL(:) = 0. -ZLAI_FULL (:) = XUNDEF - IPATCH_TRBE = VEGTYPE_TO_PATCH(NVT_TRBE, I%NPATCH) IPATCH_TRBD = VEGTYPE_TO_PATCH(NVT_TRBD, I%NPATCH) IPATCH_TEBE = VEGTYPE_TO_PATCH(NVT_TEBE, I%NPATCH) @@ -131,38 +126,36 @@ DO JJ=1,U%NSIZE_NATURE ( I%XLAI(JJ,IPATCH_BONE) * I%XVEGTYPE(JJ,NVT_BONE) )+ & ( I%XLAI(JJ,IPATCH_BOND) * I%XVEGTYPE(JJ,NVT_BOND) ) - ZH_TREE_FULL(U%NR_NATURE(JJ)) = ZH_TREE(JJ) - ZLAI_FULL (U%NR_NATURE(JJ)) = ZLAI(JJ) ! END IF ! END DO ! -ZLAI_FULL(:) = U%XNATURE(:) * ZLAI_FULL(:) +ZLAI(:) = U%XNATURE(:) * ZLAI(:) ! !* 2. Envoi les variables vers mesonH ! ------------------------------ -IF ( SIZE(PVH) /= SIZE(ZH_TREE_FULL) ) THEN +IF ( SIZE(PVH) /= SIZE(ZH_TREE) ) THEN WRITE(ILUOUT,*) 'try to get VH field from atmospheric model, but size is not correct' WRITE(ILUOUT,*) 'size of field expected by the atmospheric model (PVH) :', SIZE(PVH) - WRITE(ILUOUT,*) 'size of field inthe surface (XVH) :', SIZE(ZH_TREE_FULL) + WRITE(ILUOUT,*) 'size of field inthe surface (XVH) :', SIZE(ZH_TREE) CALL ABOR1_SFX('GET_VHN: VH SIZE NOT CORRECT') ELSE - PVH = ZH_TREE_FULL + PVH = ZH_TREE END IF ! !============================================================================== ! !------------------------------------------------------------------------------- ! -IF ( SIZE(PLAI) /= SIZE(ZLAI_FULL) ) THEN +IF ( SIZE(PLAI) /= SIZE(ZLAI) ) THEN WRITE(ILUOUT,*) 'try to get LAI field from atmospheric model, but size is not correct' WRITE(ILUOUT,*) 'size of field expected by the atmospheric model (PLAI) :', SIZE(PLAI) - WRITE(ILUOUT,*) 'size of field inthe surface (XLAI) :', SIZE(ZLAI_FULL) + WRITE(ILUOUT,*) 'size of field inthe surface (XLAI) :', SIZE(ZLAI) CALL ABOR1_SFX('GET_LAIN: LAI SIZE NOT CORRECT') ELSE - PLAI = ZLAI_FULL + PLAI = ZLAI END IF ! !==============================================================================