init_vegetation_legacy Subroutine

public subroutine init_vegetation_legacy()

Uses

  • proc~~init_vegetation_legacy~~UsesGraph proc~init_vegetation_legacy init_vegetation_legacy module~modglobal modglobal proc~init_vegetation_legacy->module~modglobal

Arguments

None

Called by

proc~~init_vegetation_legacy~~CalledByGraph proc~init_vegetation_legacy init_vegetation_legacy proc~init_vegetation init_vegetation proc~init_vegetation->proc~init_vegetation_legacy program~udales uDALES program~udales->proc~init_vegetation

Source Code

  subroutine init_vegetation_legacy
    use modglobal, only : ib, ie, jb, je, kb, ke, dzf
    implicit none
    integer :: i, j, k, m
    real, allocatable :: lai_3d(:,:,:)

    allocate(lai_3d(ib:ie, jb:je, kb:ke+1))
    lai_3d = 0.

    do k = ke, kb, -1
      lai_3d(:, :, k) = lai_3d(:, :, k+1) + lad_3d(ib:ie, jb:je, k) * dzf(k)
    end do

    do m = 1, veg%npts
      i = veg%ijk(m,1)
      j = veg%ijk(m,2)
      k = veg%ijk(m,3)
      veg%laiv(m) = lai_3d(i, j, k)
    end do

    deallocate(lai_3d)
  end subroutine init_vegetation_legacy