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

!     file: decout.for

      subroutine  decout(isr, residue)

!       + + +  PURPOSE + + + +
!      This subroutine writes decomposition output

      use datetime_mod, only: get_simdate
      use file_io_mod, only: luod_above, luod_below
      use biomaterial, only: biomatter
      use decomp_data_struct_defs, only: am0dfl

!  +  + +  COMMON  BLOCKS + + +

      include 'p1werm.inc'
      include 's1layr.inc'

!     + + + ARGUMENT DECLARATIONS + + +
      integer    isr
      type(biomatter), dimension(:), intent(in) :: residue

! + + + DECLARE VARIABLES + + +

!      integer outday(13)
      integer cd, cm, cy
      integer :: isz

! + + + FORMATS + + +

 2001 format (4x,i2,'/',i2,'/',i4,'  1',1x,f5.0,1x,f7.2,1x,f7.2,1x,     &
     &        f7.2,1x,f7.2,2x,f6.2,3x,2(f6.2,3x),1x,f6.2,3x,            &
     &        /,15x,' 2',1x,f5.0,1x,f7.2,1x,f7.2,1x,f7.2,1x,f7.2,/,     &
     &        15x,' 3',23x,f7.2,1x,f7.2)
 2005 format (7x,i2,'/',i2,'/',i4,4x,'cumddg',13x,'admbgz',             &
     & 13x,'admrtz',/,6x,'layer',7x,3(' 1',6x,' 2',8x))
 2010 format (6x,i3,4x,3(f7.2,1x,f7.2,3x))

! + + + END SPECIFICATIONS + + +

       call get_simdate(cd, cm, cy)

! write output for subregion 1 only
! standing and surface residues

!  test code  use goto to get output for specific days only

!  and remove 'c' for lines following 25 and 'end if' at 2025

!  buried = goto 10
!  stemcts = goto 15
!  sir residues = goto 20
!  no outdays = goto 25

!      goto 25

!  Buried residue South Section experiment outdays

!10    outday(1)=1
!      outday(2)=34
!      outday(3)=65
!      outday(4)=99
!      outday(5)=126
!      outday(6)=161
!      outday(7)=191
!      outday(8)=216
!
!      goto 25

!    STEMCT experiment  OUTDAYS

!15    outday(1)=18
!      outday(2)=98
!      outday(3)=158
!      outday(4)=223
!      outday(5)=289
!      outday(6)=379
!
!      goto 25
!
!c     surface sir residues output days
!
!20    outday(1)=1
!      outday(2)=34
!      outday(3)=65
!      outday(4)=99
!      outday(5)=126
!      outday(6)=161
!      outday(7)=191
!      outday(8)=216
!      outday(9)=251
!      outday(10)=279
!      outday(11)=307
!      outday(12)=336
!      outday(13)=370

!25    continue

!      if(daysim.eq.outday(1)) i=1
!      if(daysim.eq.outday(i)) then
!     output above ground data on a daily basis
      if ((am0dfl(isr) .eq. 1) .or. (am0dfl(isr) .eq. 3)) then
         write (luod_above(isr),2001) cd, cm, cy,                    &
     &   residue(1)%geometry%dstm, residue(1)%decomp%cumdds, residue(1)%deriv%mst,           &
     &   residue(1)%decomp%cumddf, residue(1)%deriv%mf, residue(1)%deriv%fscv,            &
     &   residue(1)%deriv%ffcv, residue(1)%deriv%rsai,                            &
     &   residue(1)%deriv%rsaz(1),                                            &
     &   residue(2)%geometry%dstm, residue(2)%decomp%cumdds, residue(2)%deriv%mst,           &
     &   residue(2)%decomp%cumddf, residue(2)%deriv%mf, residue(3)%decomp%cumddf,            &
     &   residue(3)%deriv%mf
      end if

!     output below ground residues
      if ((am0dfl(isr) .eq. 2) .or. (am0dfl(isr) .eq. 3)) then
         write(luod_below(isr),2005) cd, cm, cy
         do 100 isz = 1, nslay(isr)
         write (luod_below(isr),2010)                                        &
     &         isz,residue(1)%decomp%cumddg(isz),residue(2)%decomp%cumddg(isz),           &
     &         residue(1)%deriv%mbgz(isz),residue(2)%deriv%mbgz(isz),               &
     &         residue(1)%deriv%mrtz(isz),residue(2)%deriv%mrtz(isz)
  100    continue
      end if

!     end if

      return
      end