!------------------------------------------------------------------------------
!          Harvard University Atmospheric Chemistry Modeling Group            !
!------------------------------------------------------------------------------
!BOP
!
! !IROUTINE: getFuelWood
!
! !DESCRIPTION: This subroutine ...
!\\
!\\
! !INTERFACE:
!
SUBROUTINE getFuelWood
!
! !USES:
!
  USE defineConstants
  USE loadCASAinput
  USE defineArrays
  
  implicit none
!
! !REVISION HISTORY:
!  09 July 2010 - C. Carouge  - Parallelization
!EOP
!------------------------------------------------------------------------------
!BOC
!
! !LOCAL VARIABLES
!
  REAL*8  :: perc_treea(n_veg, 1)

!$OMP PARALLEL WORKSHARE  &
!$OMP DEFAULT(SHARED)
  perc_treea(:,1)=perc_tree1(:,1)
      
  WHERE (perc_tree1(:,1) == 0d0)
     perc_treea(:,1)=1.000d0
  END WHERE
  
  fuelwooddemand(:,1)=(1d0/perc_treea(:,1))
  fuelwooddemand(:,1)=fuelwooddemand(:,1)*popdens1(:,1)*&
                      fuelneed1(:,1)/12.000d0

  WHERE (perc_tree1(:,1) == 0d0)
     fuelwooddemand(:,1)=0d0
  END WHERE
!$OMP END PARALLEL WORKSHARE

END SUBROUTINE getFuelWood
!EOC