subroutine callhydr(daysim, isr)
C ***************************************************************** wjr
C Wrapper to call hydro
C
C Arguments
      integer daysim
      integer isr                   
C
C Includes
      include 'p1werm.inc'
      include 'c1glob.inc'
      include 'd1glob.inc'
      include 'm1subr.inc'
      include 'm1flag.inc'
      include 'm1dbug.inc'
      include 's1layd.inc'
      include 's1layr.inc'
      include 's1dbh.inc'
      include 's1phys.inc'
      include 's1surf.inc'
      include 'h1hydro.inc'
      include 'h1temp.inc'
      include 'h1db1.inc'
      include 'h1scs.inc'
      include 'w1wind.inc'
	  include 'w1clig.inc'
      include 'timer.fi'

      real bdmft
	  integer i

      call timer(TIMHYDR,TIMSTART)      

            if (am0hdb .eq. 1) call hdbug(isr, nslay(isr))

C     Sum up the flat residue from all decomp pools
      bdmft = 0.0
      do i=1,mnbpls
        bdmft = bdmft + admf(i,isr)
      enddo

      call hydro( nslay(isr), amrslp(isr),
     &                  acftcv(isr), acrlai(isr),
     &                  bdmft, aczrtd(isr), ahfwsf(isr),
     &                  aszlyd(1, isr), asdblk(1, isr),
     &                  ahrwc(1, isr), ahrwcs(1, isr),
     &                  ahrwcf(1, isr), ahrwcw(1, isr),
     &                  ah0cb(1,isr), aheaep(1,isr),
     &                  asfsan(1,isr), asfsil(1,isr), asfcla(1,isr),
     &                  ah0cng(isr), ah0cnp(isr),
     &                  ahzper(isr), ahzirr(isr), ahzrun(isr),
     &                  awudav, ahrsk(1, isr),
     &                  ahtsmx(1, isr), ahtsmn(1, isr),
     &                  ahrwc0(1, isr), daysim,
     &                  asfald(isr), asfalw(isr), aszlyt(1,isr),
     *                  awzdpt, awtdmx, awtdmn, ahzwid(isr) )
C

			ahzirr(isr) = 0.0	! zero out irrig amount from previous day
            if (am0hdb .eq. 1) call hdbug(isr, nslay(isr))
      call timer(TIMHYDR,TIMSTOP)      

      end