!
!$Author$
!$Date$
!$Revision$
!$HeadURL$
!
      subroutine openfils(residue)
! ***************************************************************** wjr
! Contains init code from main
!
!       Edit History
!       10-Mar-99       wjr     created

      use weps_interface_defs, ignore_me=>openfils
      use file_io_mod
      use biomaterial, only: biomatter
      use erosion_data_struct_defs, only: am0efl
      use barriers_mod, only: barseas, output_done
      use hydro_data_struct_defs, only: am0hfl, am0hdb
      use soil_data_struct_defs, only: am0sfl, am0sdb
      use manage_data_struct_defs, only: am0tfl, am0tdb
      use crop_data_struct_defs, only: am0cfl, am0cdb
      use decomp_data_struct_defs, only: am0dfl, am0ddb
      use input_run_mod, only: old_run_file

      include 'p1werm.inc'
      include 'wpath.inc'

      include 'm1flag.inc'
      include 'command.inc'
      include 'm1subr.inc'

!     + + + ARGUMENT DECLARATIONS + + +
      type(biomatter), dimension(:,:), intent(out) :: residue

!     + + +   LOCAL VARIABLES + + +
      integer :: nsubr
      integer idx, jdx, alloc_stat, sum_stat
      character*30, dimension(:), allocatable :: subr_text ! subregion subdirectory text string
      character*30 :: dec_text ! decomposition detail age pool output file name text string
      logical :: flag_set

      ! use allocation of residue array for number of subregions
      nsubr = size(residue,2)

      ! allocate the subregion name, number combination text for subregions
      allocate( subr_text(nsubr), stat=alloc_stat)
      if( alloc_stat .gt. 0 ) then
         Write(*,*) 'ERROR: unable to allocate subr_text array'
      end if

      ! create subregion directory names
      do idx = 1, nsubr
         if( old_run_file .and. (nsubr .eq. 1) ) then
            ! create the name
            subr_text(idx) = ''
         else
            ! create the name
            subr_text(idx) = makenamnum( 'subregion', idx, nsubr, '/' )
            ! create the subdirectory
            call makedir(trim(rootp)//trim(subr_text(idx)) )
         end if
      end do

!     these files are opened at all times

      sum_stat = 0
      allocate( luogui1(0:nsubr), stat=alloc_stat )
      sum_stat = sum_stat + alloc_stat
      allocate( luomandate(0:nsubr), stat=alloc_stat )
      sum_stat = sum_stat + alloc_stat
      allocate( luoharvest_si(nsubr), stat=alloc_stat )
      sum_stat = sum_stat + alloc_stat
      allocate( luoharvest_en(nsubr), stat=alloc_stat )
      sum_stat = sum_stat + alloc_stat
      allocate( luohydrobal(nsubr), stat=alloc_stat )
      sum_stat = sum_stat + alloc_stat
      allocate( luoseason(nsubr), stat=alloc_stat )
      sum_stat = sum_stat + alloc_stat
      if( sum_stat .gt. 0 ) then
         Write(*,*) 'ERROR: unable to allocate luomandate, luoharvest_, luohydrobal, luoseason arrays'
      end if
      if( .not. old_run_file .or. (nsubr .gt. 1) ) then
         call fopenk (luogui1(0), trim(rootp) // 'gui1_data.out', 'unknown')
         call fopenk (luomandate(0), trim(rootp) // 'mandate.out', 'unknown')
      end if
      do idx = 1, nsubr
         call fopenk (luogui1(idx), trim(rootp) // trim(subr_text(idx)) // 'gui1_data.out', 'unknown')
         call fopenk (luomandate(idx), trim(rootp) // trim(subr_text(idx)) // 'mandate.out', 'unknown')
         call fopenk (luoharvest_si(idx), trim(rootp) // trim(subr_text(idx)) // 'harvest_si.out', 'unknown')
         call fopenk (luoharvest_en(idx), trim(rootp) // trim(subr_text(idx)) // 'harvest_en.out', 'unknown')
         call fopenk (luohydrobal(idx), trim(rootp) // trim(subr_text(idx)) // 'hydrobal.out', 'unknown')
         call fopenk (luoseason(idx), trim(rootp) // trim(subr_text(idx)) // 'season.out', 'unknown')
      end do

      if (calibrate_crops .gt. 0) then
         sum_stat = 0
         allocate( luoharvest_calib(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luoharvest_calib_parm(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luoharvest_calib, luoharvest_calib_parm arrays'
         end if
         do idx = 1, nsubr
            ! calibration harvest output file
            call fopenk (luoharvest_calib(idx), trim(rootp) // trim(subr_text(idx)) // 'luoharvest_calib.out', 'unknown')
            ! calibration harvest output file for GUI
            call fopenk (luoharvest_calib_parm(idx), trim(rootp) // trim(subr_text(idx)) // 'luoharvest_calib_parm.out', 'unknown')
         end do
      endif

      ! open barrier output file
      if( size(barseas) .gt. 0 ) then
          call fopenk (luo_barr, rootp(1:len_trim(rootp)) // 'barriers.out', 'unknown')
          output_done = .false.
      else
          output_done = .true.
      end if

!     open erosion output files
      if (am0efl.gt.0) then
          call fopenk (luo_subday, rootp(1:len_trim(rootp)) // 'subday.out', 'unknown')
      endif

      if (btest(am0efl,0)) then
       call fopenk (luo_erod, rootp(1:len_trim(rootp)) // 'daily_out.erod', 'unknown')
      endif

!     open plot data file
      if(     (maxval(am0hfl).gt.0) .or. (maxval(am0sfl).gt.0) .or. (maxval(am0tfl).gt.0) &
         .or. (maxval(am0cfl).gt.0) .or. (maxval(am0dfl).gt.0) .or. (am0efl.gt.0)) then
         allocate( luoplt(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luoplt array'
         end if
         do idx = 1, nsubr
            if(     (am0hfl(idx).gt.0) .or. (am0sfl(idx).gt.0) .or. (am0tfl(idx).gt.0) &
               .or. (am0cfl(idx).gt.0) .or. (am0dfl(idx).gt.0) .or. (am0efl.gt.0)) then
               call fopenk (luoplt(idx), trim(rootp) // trim(subr_text(idx)) // 'plot.out', 'unknown')
            end if
         end do
      endif

!     open output file for soil conditioning index
      if( soil_cond .gt. 0 ) then
         sum_stat = 0
         allocate( luosci(0:nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luostir(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luosci, luostir arrays'
         end if
         if( .not. old_run_file .or. (nsubr .gt. 1) ) then
            call fopenk (luosci(0), trim(rootp) // 'sci_energy.out', 'unknown')
         end if
         do idx = 1, nsubr
            call fopenk (luosci(idx), trim(rootp) // trim(subr_text(idx)) // 'sci_energy.out', 'unknown')
            call fopenk (luostir(idx), trim(rootp) // trim(subr_text(idx)) // 'stir_energy.out', 'unknown')
         end do
      end if

!     open detailed output files for hydro
      flag_set = .false.
      do idx = 1, nsubr
         if ((am0hfl(idx) .eq. 1) .or. (am0hfl(idx) .eq. 3) .or. (am0hfl(idx) .eq. 5) .or. (am0hfl(idx) .eq. 7)) then
            flag_set = .true.
         end if
      end do
      if( flag_set ) then
         sum_stat = 0
         allocate( luohydro(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luohlayers(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luohydro, luohlayers arrays'
         end if
         do idx = 1, nsubr
            if ((am0hfl(idx) .eq. 1) .or. (am0hfl(idx) .eq. 3) .or. (am0hfl(idx) .eq. 5) .or. (am0hfl(idx) .eq. 7)) then
               call fopenk (luohydro(idx), trim(rootp) // trim(subr_text(idx)) // 'hydro.out', 'unknown')
               call fopenk (luohlayers(idx), trim(rootp) // trim(subr_text(idx)) // 'hlayers.out', 'unknown')
            end if
         end do
      endif

      flag_set = .false.
      do idx = 1, nsubr
         if ((am0hfl(idx) .eq. 2) .or. (am0hfl(idx) .eq. 6) .or. (am0hfl(idx) .eq. 3) .or. (am0hfl(idx) .eq. 7)) then
            flag_set = .true.
         end if
      end do
      if( flag_set ) then
         allocate( luowater(nsubr), stat=alloc_stat )
         allocate( luosurfwat(nsubr), stat=alloc_stat )
         allocate( luoweather(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luowater array'
         end if
         do idx = 1, nsubr
            if ((am0hfl(idx) .eq. 2) .or. (am0hfl(idx) .eq. 6) .or. (am0hfl(idx) .eq. 3) .or. (am0hfl(idx) .eq. 7)) then
               call fopenk (luowater(idx), trim(rootp) // trim(subr_text(idx)) // 'water.out', 'unknown')
               call fopenk (luosurfwat(idx), trim(rootp) // trim(subr_text(idx)) // 'surfwat.out', 'unknown')
               call fopenk (luoweather(idx), trim(rootp) // trim(subr_text(idx)) // 'weather.out', 'unknown')
            end if
         end do
      end if

      flag_set = .false.
      do idx = 1, nsubr
         if ((am0hfl(idx) .eq. 4) .or. (am0hfl(idx) .eq. 5) .or. (am0hfl(idx) .eq. 6) .or. (am0hfl(idx) .eq. 7)) then
            flag_set = .true.
         end if
      end do
      if( flag_set ) then
         allocate( luotempsoil(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luotempsoil array'
         end if
         do idx = 1, nsubr
            if ((am0hfl(idx) .eq. 4) .or. (am0hfl(idx) .eq. 5) .or. (am0hfl(idx) .eq. 6) .or. (am0hfl(idx) .eq. 7)) then
               call fopenk (luotempsoil(idx), trim(rootp) // trim(subr_text(idx)) // 'temp.out', 'unknown')
            end if
         end do
      end if

! open files for outputing the crop and decomp biomass variables - LEW
      flag_set = .false.
      do idx = 1, nsubr
         if ((am0dfl(idx) .eq. 1).or.(am0dfl(idx).eq.3)) then
            flag_set = .true.
         end if
      end do
      if( flag_set ) then
         sum_stat = 0
         allocate( luocrp1(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luobio1(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luod_above(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luod_above, luocrp1 or luobio1 array'
         end if

         do idx = 1, nsubr
            if ((am0dfl(idx) .eq. 1).or.(am0dfl(idx).eq.3)) then
               call fopenk (luocrp1(idx), trim(rootp) // trim(subr_text(idx)) // 'decomp.out', 'unknown')
               call fopenk (luobio1(idx), trim(rootp) // trim(subr_text(idx)) // 'bio1.btmp', 'unknown')
               call fopenk (luod_above(idx), trim(rootp) // trim(subr_text(idx)) // 'dabove.out', 'unknown')
            end if
         end do
      endif

      flag_set = .false.
      do idx = 1, nsubr
         if ((am0dfl(idx) .eq. 2).or.(am0dfl(idx).eq.3)) then
            flag_set = .true.
         end if
      end do
      if( flag_set ) then
         ! create dbelow.out unit number array for subregions
         allocate( luod_below(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            write(*,*) 'ERROR: unable to allocate luod_below array'
         end if
         do idx = 1, nsubr
            if ((am0dfl(idx) .eq. 2).or.(am0dfl(idx).eq.3)) then
               ! open dbelow.out in each subregion
               call fopenk (luod_below(idx), trim(rootp) // trim(subr_text(idx)) // 'dbelow.out', 'unknown')
               ! open files to match number of biomass pools
               ! create age pool output file names, set unit numbers and open files
               do jdx = 1,mnbpls
                  ! create the name
                  dec_text = makenamnum( 'dec', jdx, mnbpls, '.btmp' )
                  ! display created file
                  !write(*,*) 'File name created: ', dec_text
                  ! assign logical unit number of opening file to array
                  call fopenk (residue(jdx,idx)%luo%dec, trim(rootp) // trim(subr_text(idx)) // trim(dec_text), 'unknown')
               end do
            end if
         end do
      endif

      if( maxval(am0cfl) .gt. 0) then
         sum_stat = 0
         allocate( luocrop(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luoshoot(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luoinpt(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luocrop, luoshoot, luoinpt arrays'
         end if

         do idx = 1, nsubr
            if (am0cfl(idx) .gt. 0) then
               ! daily crop output of most state variables 
               call fopenk (luocrop(idx), trim(rootp) // trim(subr_text(idx)) // 'crop.out', 'unknown')
               call fopenk (luoshoot(idx), trim(rootp) // trim(subr_text(idx)) // 'shoot.out', 'unknown')
               ! echo crop input data - AR
               call fopenk (luoinpt(idx), trim(rootp) // trim(subr_text(idx)) // 'inpt.out', 'unknown')
            end if
         end do
      endif

        ! print headings for crop output files
        ! season.out, crop.out, shoot.out, inpt.out
      do idx = 1, nsubr
        call cpout(idx)
      end do

      if( maxval(am0sfl) .eq. 1 ) then
         ! soil detail output files
         sum_stat = 0
         allocate( luosoilsurf(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luosoillay(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luosoilsurf, luosoillay arrays'
         end if
         do idx = 1, nsubr
            if( am0sfl(idx) .eq. 1 ) then
               ! soil surface
               call fopenk(luosoilsurf(idx), trim(rootp) // trim(subr_text(idx)) // 'soilsurf.out', 'unknown')
               ! soil layers
               call fopenk(luosoillay(idx), trim(rootp) // trim(subr_text(idx)) // 'soillay.out', 'unknown')
            end if
         end do
      endif

      if (maxval(am0tfl) .eq. 1) then
         allocate( luomanage(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luomanage array'
         end if
         do idx = 1, nsubr
            if (am0tfl(idx) .eq. 1) then
               call fopenk (luomanage(idx), trim(rootp) // trim(subr_text(idx)) // 'manage.out', 'unknown')
            end if
         end do
      end if

      if ((calc_confidence .gt. 0)) then
         ! Confidence Interval output file
         call fopenk(luoci, rootp(1:len_trim(rootp)) // 'ci.out', 'unknown')
      endif

         ! create arrays for subregion debug output files
      if (maxval(am0hdb) .eq. 1) then
         allocate( luohdb(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luohdb array'
         end if
         do idx = 1, nsubr
            if (am0hdb(idx) .eq. 1) then
               call fopenk (luohdb(idx), trim(rootp) // trim(subr_text(idx)) // 'hdbug.out', 'unknown')
            end if
         end do
      end if
      if (maxval(am0sdb) .eq. 1) then
         allocate( luosdb(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luosdb array'
         end if
         do idx = 1, nsubr
            if (am0sdb(idx) .eq. 1) then
               call fopenk (luosdb(idx), trim(rootp) // trim(subr_text(idx)) // 'sdbug.out', 'unknown')
            end if
         end do
      end if
      if (maxval(am0tdb) .eq. 1) then
         allocate( luotdb(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luotdb array'
         end if
         do idx = 1, nsubr
            if (am0tdb(idx) .eq. 1) then
               call fopenk (luotdb(idx), trim(rootp) // trim(subr_text(idx)) // 'tdbug.out', 'unknown')
            end if
         end do
      end if
      if (maxval(am0cdb) .eq. 1) then
         allocate( luocdb(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luocdb array'
         end if
         do idx = 1, nsubr
            if (am0cdb(idx) .eq. 1) then
               call fopenk (luocdb(idx), trim(rootp) // trim(subr_text(idx)) // 'cdbug.out', 'unknown')
            end if
         end do
      end if
      if (maxval(am0ddb) .eq. 1) then
         allocate( luoddb(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luoddb array'
         end if
         do idx = 1, nsubr
            if (am0ddb(idx) .eq. 1) then
               call fopenk (luoddb(idx), trim(rootp) // trim(subr_text(idx)) // 'ddbug.out', 'unknown')
            end if
         end do
      end if

!   WEPP Related files
!
       if (wepp_hydro .gt. 1) then
         allocate( luowepphdrive(nsubr), stat=alloc_stat )
         if( alloc_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luowepphdrive array'
         end if
         do idx = 1, nsubr
            call fopenk (luowepphdrive(idx), trim(rootp) // trim(subr_text(idx)) // 'wepp_runoff.out', 'unknown')
            write(luowepphdrive(idx),*) ' WEPP Flow Routing Output'
            write(luowepphdrive(idx),*) ' # day   mon  yr     precip  runoff    peakro  effdrn    effint   effdrr/rainfall excess'
            write(luowepphdrive(idx),*) '                      (mm)   (mm)     (mm/hr)  (min)    (mm/hr)     (min)'
         end do
       endif

       if ((run_erosion.eq.2).or.(run_erosion.eq.3)) then
         sum_stat = 0
         allocate( luowepperod(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luoweppplot(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         allocate( luoweppsum(nsubr), stat=alloc_stat )
         sum_stat = sum_stat + alloc_stat
         if( sum_stat .gt. 0 ) then
            Write(*,*) 'ERROR: unable to allocate luowepperod, luoweppplot, luoweppsum arrays'
         end if
         do idx = 1, nsubr
            call fopenk(luowepperod(idx), trim(rootp) // trim(subr_text(idx)) // 'wepp_eroevents.out','unknown')
            write(luowepperod(idx),*) 'WEPP Erosion Events Output'
            write(luowepperod(idx),*) &
            'day mo  year    Precp  Runoff  IR-det Av-det Mx-det  Point  Av-dep Max-dep  Point Sed.Del    ER'
            write(luowepperod(idx),*) &
            '--- --  ----     (mm)    (mm)  kg/m^2 kg/m^2 kg/m^2    (m)  kg/m^2  kg/m^2    (m)  (kg/m)  ----'

            call fopenk(luoweppplot(idx), trim(rootp) // trim(subr_text(idx)) // 'wepp_eroplot.out','unknown')
     
            call fopenk(luoweppsum(idx), trim(rootp) // trim(subr_text(idx)) // 'wepp_summary.out','unknown')
            write(luoweppsum(idx),*) 'WEPS/WEPP Common Model'
            write(luoweppsum(idx),*) 'March 3, 2009  (2009.3)'
            write(luoweppsum(idx),*) '---------------------------------------'
         end do
       endif

      ! free memory from local subregion text strings
      deallocate( subr_text, stat=alloc_stat)
      if( alloc_stat .gt. 0 ) then
         Write(*,*) 'ERROR: unable to deallocate subr_text array'
      end if


      end