diff --git a/biogeochem/EDPhysiologyMod.F90 b/biogeochem/EDPhysiologyMod.F90
index e315e8cfe2..67c7cb327d 100644
--- a/biogeochem/EDPhysiologyMod.F90
+++ b/biogeochem/EDPhysiologyMod.F90
@@ -244,7 +244,7 @@ subroutine trim_canopy( currentSite )
                    ! Use slatop for sla value at all leaf levels 
                    sla_levleaf = EDPftvarcon_inst%slatop(ipft)
                 
-                elseif(if(use_slaprofile == 1)then
+                elseif(use_slaprofile == 1)then
                    
                    ! Calculate sla_levleaf following the sla profile with overlying leaf area
                    ! Scale for leaf nitrogen profile