Merge pull request #360 from ckoven/laicanfix

laican bugfix

Author: rgknox

biogeophys/FatesPlantRespPhotosynthMod.F90

@@ -324,10 +324,12 @@ subroutine FatesPlantRespPhotosynthDrive (nsites, sites,bc_in,bc_out,dtime)
                      ! are there any leaves of this pft in this layer?
                      if(currentPatch%canopy_mask(cl,ft) == 1)then 
 
-                        if(cl==NCL_p)then !are we in the top canopy layer or a shaded layer?
+                        if(cl==1)then !are we in the top canopy layer or a shaded layer?
                            laican = 0._r8
                         else
-                           laican = sum(currentPatch%canopy_layer_tai(cl+1:NCL_p)) 
+
+                           laican = sum(currentPatch%canopy_layer_tai(1:cl-1)) 
+
                         end if
 
                         ! Loop over leaf-layers
@@ -366,7 +368,7 @@ subroutine FatesPlantRespPhotosynthDrive (nsites, sites,bc_in,bc_out,dtime)
                                  laican = laican + 0.5_r8 * vai
                               else
                                  laican = laican + 0.5_r8 * (currentPatch%elai_profile(cl,ft,iv-1)+ &
-                                       currentPatch%esai_profile(cl,ft,iv-1))+vai
+                                       currentPatch%esai_profile(cl,ft,iv-1)+vai)
                               end if
 
                               ! Scale for leaf nitrogen profile