!$Author$
!$Date$
!$Revision$
!$HeadURL$

      subroutine submodels (isr, cd, cm, cy)

      include 'p1werm.inc'
      include 'm1flag.inc'      !am0cgf
      include 'main/main.inc'   !daysim, lopday, lopmon, lopyr, iy

      include 'c1glob.inc'
      include 'c1gen.inc'
      include 'c1db1.inc'
      include 'c1db2.inc'
      include 'h1hydro.inc'
      include 's1sgeo.inc'
      include 's1layr.inc'

! Arguments
      integer isr, cd, cm, cy

!        write(*,*) "Start manage"      !MANAGEment (tillage) submodel
        call manage (isr, cd, cm, cy,iy,lopday,lopmon,lopyr)

        if( am0cropupfl.gt.0 ) then
        ! update all derived globals for crop global variables
          call cropupdate(                                              &
     &      acmstandstem(isr), acmstandleaf(isr), acmstandstore(isr),   &
     &      acmflatstem(isr), acmflatleaf(isr), acmflatstore(isr),      &
     &      acmshoot(isr), acmbgstemz(1,isr),                           &
     &      acmrootstorez(1,isr),acmrootfiberz(1,isr),                  &
     &      aczht(isr), acdstm(isr), aczrtd(isr),                       &
     &      acthucum(isr), aczgrowpt(isr), acmbgstem(isr),              &
     &      acmrootstore(isr), acmrootfiber(isr), acxstmrep(isr),       &
     &      acm(isr), acmst(isr), acmf(isr), acmrt(isr), acmrtz(1,isr), &
     &      acrcd(isr), aszrgh(isr), asxrgs(isr), asargo(isr),          &
     &      acrsai(isr), acrlai(isr), acrsaz(1,isr), acrlaz(1,isr),     &
     &      acffcv(isr), acfscv(isr), acftcv(isr), acfcancov(isr),      &
     &      ac0rg(isr), acxrow(isr),                                    &
     &      nslay(isr), ac0ssa(isr), ac0ssb(isr), ac0sla(isr),          &
     &      accovfact(isr), ac0ck(isr), acxstm(isr), acdpop(isr),       &
     &      ahztranspdepth(isr), ahzfurcut(isr),                        &
     &      ahztransprtmin(isr), ahztransprtmax(isr)  )

          ! dependent variables have been updated
          am0cropupfl = 0
        end if

!        write(*,*) "Start updres"
        call updres(isr)                 !update decomp residue pools

!        write(*,*) "Start callhydr"
        call callhydr(daysim, isr)      !call HYDROLOGY submodel
        ! do not change order. Hydro may set irrigation amounts that
        ! will affect soil.

!        write(*,*) "Start callsoil"
        call callsoil(daysim, isr)       !SOIL submodel

!        write(*,*) "Start callcrop"     !CROP submodel
        ! Crop growth flag indicates growing crop
        ! Harvest flag indicates that harvest occured today. Crop is called
        ! to generate end of growth period report from values retained in
        ! the previous day crop data registers even though growth flag is
        ! turned off.
        if( am0cgf ) then
            call callcrop(daysim, isr)
        end if

!        write(*,*) "Start decomp"
        call decomp(isr)                 !DECOMPosition submodel

!        write(*,*) "Start updres"
        call updres(isr)

!        write(*,*) "Start sumbio"
        call sumbio(isr) ! sum live and dead biomass

      return
      end