! source file: /Users/mtivig/REPOS/uvic_with_news/RUNS/RUN_UVic_2015_NEWS_BD_BP_V2020/updates/MTivig_2020_BP/source/mtlm/afforestation.F
      subroutine AFFORESTATION (LAND_PTS, LAND_INDEX, DZ_SOIL, M
     &,                         FRAC)

! This is a climate engineering subroutine for irrigating N. Africa
! and Australia.

      implicit none

      include "size.h"
      include "csbc.h"

! LAND_PTS   = IN Number of land points.
! LAND_INDEX = IN Indices of landpoints.
! GREEN_THETA_SAT = Saturated volumetric soil moisture
!                   concentration resulting from irrigation
! DZ_SOIL = IN Soil layer thickness (m).
! FRAC_VS = fraction of grid cell containing veg

      integer LAND_PTS, LAND_INDEX(POINTS)
      integer I, J, L, N
      real GREEN_THETA_SAT, M(POINTS) , DZ_SOIL
      real FRAC(POINTS,NTYPE)

      GREEN_THETA_SAT = 0.40

      return
      end