c$Author: fredfox $
c$Date: 2001-09-27 22:02:44 $
c$Revision: 1.12 $
c$Source: /weru/cvs/weps/weps.src/soil/callsoil.for,v $
      subroutine callsoil(daysim, isr)
C ***************************************************************** wjr
C Wrapper to call soil
C
C Arguments
      integer daysim
      integer isr                   
      integer ldx
C
C Includes
      include 'p1werm.inc'
      include 'b1glob.inc'
      include 'm1subr.inc'
      include 'm1flag.inc'
      include 'm1dbug.inc'
      include 's1agg.inc'
      include 's1layd.inc'
      include 's1layr.inc'
      include 's1dbc.inc'
      include 's1dbh.inc'
      include 's1phys.inc'
      include 's1sgeo.inc'
      include 's1surf.inc'
      include 'h1hydro.inc'
      include 'h1temp.inc'
      include 'h1db1.inc'
      include 'timer.fi'

      call timer(TIMSOIL,TIMSTART)      
C
            if (am0sdb .eq. 1) call sdbug(isr, nslay(isr))
            call soil (daysim, am0irr(isr), ahzirr(isr), ahzsmt(isr),
     &                 ahtsmx(1,isr), ahtsmn(1,isr),
     &                 ahrwc(1,isr), ahrwca(1,isr), ahrwcw(1,isr),
     &                 asfom(1,isr), aszlyt(1,isr), nslay(isr),
     &                 asfsan(1,isr), asfsil(1,isr), asfcla(1,isr),
     &                 asvroc(1,isr),
     &                 asxrgs(isr), aszrgh(isr), aszrho(isr),
     &                 aslrr(isr), aslrro(isr),
     &                 aszcr(isr), asfcr(isr), asecr(isr),
     &                 asdcr(isr), asmlos(isr), asflos(isr),
     &                 asdsblk(1,isr), asdwblk(1,isr),
     &                 asdblk(0,isr), asdagd(0,isr),
     &                 aslagm(0,isr), aslagn(0,isr),
     &                 as0ags(0,isr), aslagx(0,isr), aseags(0,isr),
     &                 aseagm(1,isr), aseagmn(1,isr), aseagmx(1,isr),
     &                 ask4d(1,isr), aslmin(1,isr), aslmax(1,isr),
     &                 abffcv(isr), abfscv(isr),
     &                 asfcce(1,isr), asfcec(1,isr),
     &                 ahzwid(isr)
     &                )
            if (am0sdb .eq. 1) call sdbug(isr, nslay(isr))

c Depth to layer bottoms need to be updated after SOIL and MANAGE
            aszlyd(1, isr) = aszlyt(1, isr)
            do  ldx = 2, nslay(isr)
               aszlyd(ldx,isr) = aszlyt(ldx,isr) + aszlyd(ldx-1, isr)
            end do

      call timer(TIMSOIL,TIMSTOP)      

C
      end