Update stand_bal.f

correction

source_code/version2.2_windows/stand_bal.f

@@ -104,7 +104,12 @@ do
             sumbio  = sumbio  + ntr * zeig%coh%totBio
             sumNPP  = sumNPP  + ntr * zeig%coh%NPP
             Ndem    = Ndem + ntr * zeig%coh%Ndemc_c
-            autresp = autresp + ntr * zeig%coh%maintres
+            select case (flag_dis)
+	         case (0,1)
+             autresp = autresp + ntr * zeig%coh%maintres
+             case (2)
+             autresp = autresp + ntr * (zeig%coh%maintres+zeig%coh%biocost_all)
+            end select
             totfol  = totfol  + ntr * zeig%coh%x_fol
             totsap  = totsap  + ntr * zeig%coh%x_sap
             totfrt  = totfrt  + ntr * zeig%coh%x_frt
@@ -498,10 +503,6 @@ svar%frt        = 0.
                       END SELECT 
                     Enddo 
 
-
-
-
-
                     IF(spar(spec_new)%phmodel==4) THEN
                         svar(spec_new)%daybb = svar(spec_new)%ext_daybb
                     ELSE
@@ -543,7 +544,7 @@ do i=1,nspecies
           svar(i)%mean_diam   = svar(i)%mean_diam / ntr
           svar(i)%mean_jrb    = svar(i)%mean_jrb / ntr
           
-          svar(i)%basal_area = pi*ntr*(svar(i)%med_diam*svar(i)%med_diam/40000)*10000/kpatchsize
+          svar(i)%basal_area = pi*(ntr-helpdiam(i))*(svar(i)%med_diam*svar(i)%med_diam/40000)*10000/kpatchsize
        else
           svar(i)%sum_ntreea = 0.
        endif
@@ -1083,8 +1084,14 @@ IMPLICIT NONE
   DO
     IF(.not.associated(zeig)) exit
      if (zeig%coh%species.ne.nspec_tree+2) then                          ! exclude mistletoe from senescence
-      zeig%coh%sfol = spar(zeig%coh%species)%psf * zeig%coh%x_fol
-      zeig%coh%sfrt = spar(zeig%coh%species)%psr * zeig%coh%x_frt
+         select case (flag_dis)
+	     ! case (1,2)
+	     !  zeig%coh%sfol = spar(zeig%coh%species)%psf * zeig%coh%x_fol + zeig%coh%x_fol_loss
+         !  zeig%coh%sfrt = spar(zeig%coh%species)%psr * zeig%coh%x_frt + zeig%coh%x_frt_loss
+	      case (0,1,2)
+	       zeig%coh%sfol = spar(zeig%coh%species)%psf * zeig%coh%x_fol
+           zeig%coh%sfrt = spar(zeig%coh%species)%psr * zeig%coh%x_frt
+	     end select
       IF (zeig%coh%height.ge.thr_height .and.zeig%coh%species.LE. nspec_tree) THEN
          zeig%coh%ssap = spar(zeig%coh%species)%pss * zeig%coh%x_sap
       ELSE