From 7732cb0955664ed795292636eb99d97ec60c961b Mon Sep 17 00:00:00 2001 From: Quentin Rodier <quentin.rodier@meteo.fr> Date: Thu, 18 Jul 2024 14:55:22 +0200 Subject: [PATCH] Jean W. 18/07/2024: SURFEX cosmetics, if statements expanded --- src/SURFEX/veg_from_lai.F90 | 21 ++++++++++++++++----- 1 file changed, 16 insertions(+), 5 deletions(-) diff --git a/src/SURFEX/veg_from_lai.F90 b/src/SURFEX/veg_from_lai.F90 index dd055fe88..59ea4d367 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 -- GitLab