! source file: /gxfs_work/geomar/smomw113/UVic_ESCM/2.9_dk_0814/updates/11_valley/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

! last endif: of _green_africa

! last endif: of 1 (this was missing in Davids code)
      return
      end