c$Author: fredfox $
c$Date: 2002-05-06 22:53:17 $
c$Revision: 1.5 $
c$Source: /weru/cvs/weps/weps.src/manage/report_harvest.for,v $

      subroutine report_harvest( sr, mass_rem, mass_left )

c     + + + ARGUMENT DECLARATIONS + + +
      integer sr
      real mass_rem, mass_left

c     + + + ARGUMENT DEFINITIONS + + +
c     sr    - subregion number
c     mass_rem - mass removed by the harvest process
c     mass_left - mass left behind by the harvest process

c     + + + PARAMETERS AND COMMON BLOCKS + + +
      include 'p1werm.inc'
      include 'file.fi'
      include 'main/main.inc'
      include 'c1info.inc'
      include 'c1gen.inc'
      include 'c1report.inc'
      include 'manage/oper.inc'

c     + + + FUNCTION DECLARATIONS + + +
      integer difdat

c     + + + LOCAL DECLARATIONS + + +
      real tot_mass, harvest_index

      if( difdat(prevday(sr),prevmon(sr),prevyr(sr),
     &           lopday,lopmon,lopyr).le.0 ) then
c         write newline
          write(unit=luoharvest,fmt=1001) ''
      end if

      tot_mass = mass_rem + mass_left
      if( tot_mass .le. 0.0 ) then
          harvest_index = 0.0
      else
          harvest_index = mass_rem/tot_mass
      end if

      write(unit=luoharvest,fmt=1000,advance='NO')
     &      lopday, lopmon, lopyr,
     &      ac0nam(sr)(1:len_trim(ac0nam(sr))),
     &      mass_rem, 'kg/m^2', harvest_index,
     &      "Harvest Index", 
     &      mass_rem * acycon(sr) / ( 1.0-acywct(sr)/100.0 ),
     &      acynmu(sr)(1:len_trim(acynmu(sr))),
     &      acywct(sr), 'percent water'
c           the conversion is from dry mass to wet weight
c           and from kg/m^2 to acynmu units

 1000 format(1x,i2,'/',i2,'/',i2,'|',a,'|',f7.3,'|',a,'|',f6.3,
     &       '|',a,'|',f9.2,'|',a,'|',f5.1,'|',a,'|')
 1001 format(a)

      prevday(sr) = lopday
      prevmon(sr) = lopmon
      prevyr(sr) = lopyr

      return
      end