diff --git a/src/SURFEX/veg_from_lai.F90 b/src/SURFEX/veg_from_lai.F90 index dd055fe88de3893ea16ea8a409eae1b527e0f771..59ea4d36737ad6e5353149f0c21abdaa4b568c18 100644 --- a/src/SURFEX/veg_from_lai.F90 +++ b/src/SURFEX/veg_from_lai.F90 @@ -577,15 +577,26 @@ REAL(KIND=JPRB) :: ZHOOK_HANDLE IF (LHOOK) CALL DR_HOOK('MODI_VEG_FROM_LAI:VEG_FROM_LAI_VEGTYPE_1D',0,ZHOOK_HANDLE) PVEG(:) = XUNDEF ! -IF (PLAI(NVT_C4)/=XUNDEF) PVEG(NVT_C4 )= 1. - EXP( -0.6 * PLAI(NVT_C4 ) ) + +IF (PLAI(NVT_C4)/=XUNDEF) THEN + PVEG(NVT_C4 )= 1. - EXP( -0.6 * PLAI(NVT_C4 ) ) +END IF IF (NVT_IRR>0) THEN - IF (PLAI(NVT_IRR)/=XUNDEF) PVEG(NVT_IRR )= 1. - EXP( -0.6 * PLAI(NVT_IRR ) ) + IF (PLAI(NVT_IRR)/=XUNDEF) THEN + PVEG(NVT_IRR )= 1. - EXP( -0.6 * PLAI(NVT_IRR ) ) + END IF ENDIF IF (NVT_C3>0) THEN - IF (PLAI(NVT_C3)/=XUNDEF) PVEG(NVT_C3 )= 1. - EXP( -0.6 * PLAI(NVT_C3 ) ) + IF (PLAI(NVT_C3)/=XUNDEF) THEN + PVEG(NVT_C3 )= 1. - EXP( -0.6 * PLAI(NVT_C3 ) ) + END IF ELSEIF (NVT_C3W>0 .AND. NVT_C3S>0) THEN - IF (PLAI(NVT_C3W)/=XUNDEF) PVEG(NVT_C3W )= 1. - EXP( -0.6 * PLAI(NVT_C3W ) ) - IF (PLAI(NVT_C3S)/=XUNDEF) PVEG(NVT_C3S )= 1. - EXP( -0.6 * PLAI(NVT_C3S ) ) + IF (PLAI(NVT_C3W)/=XUNDEF) THEN + PVEG(NVT_C3W )= 1. - EXP( -0.6 * PLAI(NVT_C3W ) ) + END IF + IF (PLAI(NVT_C3S)/=XUNDEF) THEN + PVEG(NVT_C3S )= 1. - EXP( -0.6 * PLAI(NVT_C3S ) ) + END IF ENDIF ! PVEG(NVT_TEBD)= 0.95