ch_sedbud_output.f90 Source File

      subroutine ch_sedbud_output

      use time_module
      use hydrograph_module
      use sd_channel_module
      
      implicit none
      
      integer :: ichan
      integer :: iob

      !! loop through and print each use object
      do ichan = 1, sp_ob%chandeg
        iob = sp_ob1%chandeg + ichan - 1
        
        !! sum monthly variables
        ch_morphm(ichan) = ch_morphm(ichan) + ch_morph(ichan)
        
        !! daily print
        if (pco%sd_chan%d == "y") then
          write (3171,*) time%day, time%mo, time%day_mo, time%yrc, ichan, &
            ob(iob)%gis_id, ob(iob)%name, ch_morph(ichan)

          if (pco%csvout == "y") then
          write (3175,'(*(G0.6,:","))') time%day, time%mo, time%day_mo, time%yrc, ichan, &
            ob(iob)%gis_id, ob(iob)%name, ch_morph(ichan)
          end if
        end if
       
        !! zero daily
        ch_morph(ichan) = ch_morphz

        !! monthly print
        if (time%end_mo == 1) then
          !! add into the yearly total
          ch_morphy(ichan) = ch_morphy(ichan) + ch_morphm(ichan)

          if (pco%sd_chan%m == "y") then
          write (3172,*) time%day, time%mo, time%day_mo, time%yrc, ichan, &
            ob(iob)%gis_id, ob(iob)%name, ch_morphm(ichan)
 
          if (pco%csvout == "y") then
          write (3176,'(*(G0.6,:","))') time%day, time%mo, time%day_mo, time%yrc, ichan, &
            ob(iob)%gis_id, ob(iob)%name, ch_morphm(ichan)
          end if
          end if

          !! zero monthly
          ch_morphm(ichan) = ch_morphz

        end if

      !! yearly print
      if (time%end_yr == 1) then
        !! add into the yearly total
        ch_morphy(ichan) =  ch_morphy(ichan) + ch_morphm(ichan)
          
        if (pco%sd_chan%y == "y") then
          write (3173,*) time%day, time%mo, time%day_mo, time%yrc, ichan, &
            ob(iob)%gis_id, ob(iob)%name, ch_morphy(ichan)
  
              if (pco%csvout == "y") then
          write (3177,'(*(G0.6,:","))') time%day, time%mo, time%day_mo, time%yrc, ichan, &
            ob(iob)%gis_id, ob(iob)%name, ch_morphy(ichan)
          end if
        end if

        !! accumulate the year's total into the average-annual accumulator
        !! (fix: yearly write used ch_morpha and zeroed it each year, so AA divided from 0)
        ch_morpha(ichan) = ch_morpha(ichan) + ch_morphy(ichan)

        !! zero yearly
        ch_morphy(ichan) = ch_morphz

      end if

      !! average annual print
      if (time%end_sim == 1) then
        !! convert the accumulated total to an average annual value
        ch_morpha(ichan) = ch_morpha(ichan) / time%yrs_prt

        if (pco%sd_chan%a == "y") then
        write (3174,*) time%day, time%mo, time%day_mo, time%yrc, ichan, &
          ob(iob)%gis_id, ob(iob)%name, ch_morpha(ichan)

        if (pco%csvout == "y") then
        write (3178,'(*(G0.6,:","))') time%day, time%mo, time%day_mo, time%yrc, ichan, &
          ob(iob)%gis_id, ob(iob)%name, ch_morpha(ichan)
        end if
       end if
      end if

      end do    ! do ichan = 1, sp_ob%chandeg
      
      return
      
      end subroutine ch_sedbud_output