Skip to content
Snippets Groups Projects
Commit 7732cb09 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Jean W. 18/07/2024: SURFEX cosmetics, if statements expanded

parent 03ca7421
No related branches found
No related tags found
No related merge requests found
......@@ -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
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment