!$Author: fredfox $ !$Date: 2004-06-03 06:57:27 $ !$Revision: 1.1 $ !$Source: /weru/cvs/weps/weps.src/soil/depthini.for,v $ subroutine depthini(nlay, bszlyt, bszlyd) integer nlay real bszlyt(*), bszlyd(*) integer idx ! nlay - number of soil layers ! bszlyt - soil layer thickness (mm) ! bszlyd - depth to bottom of soil layer (mm) bszlyd(1) = bszlyt(1) do idx = 2, nlay bszlyd(idx) = bszlyt(idx) + bszlyd(idx-1) end do return end