write_boundary_config_summary_to_unit Subroutine

   subroutine write_boundary_config_summary_to_unit(unit_id, params)
      integer, intent(in) :: unit_id
      type(case_params_t), intent(in) :: params

      integer :: p

      write(unit_id,'(a)') ''
      write(unit_id,'(a)') 'Boundary conditions:'
      write(unit_id,'(a,i0)') '  number of patches            : ', params%n_patches
      write(unit_id,'(a)') '  mass-flux velocity convention:'
      write(unit_id,'(a)') '    mass_flux/signed_mass_flux : signed rho*dot(U,n_out) [kg/m^2/s]'
      write(unit_id,'(a)') '    inlet_mass_flux            : positive magnitude directed into the domain'
      write(unit_id,'(a)') '    outlet_mass_flux           : positive magnitude directed out of the domain'

      if (params%n_patches <= 0) return

      do p = 1, params%n_patches
         write(unit_id,'(a,i0,a,a)') '  patch[', p, '] ', trim(params%patch_name(p))
         write(unit_id,'(a,a)') '    patch type                 : ', trim(params%patch_type(p))
         write(unit_id,'(a,a)') '    velocity type              : ', trim(params%patch_velocity_type(p))
         write(unit_id,'(a,3(es12.5,1x))') '    velocity [m/s]             : ', &
            params%patch_u(p), params%patch_v(p), params%patch_w(p)
         write(unit_id,'(a,es12.5)') '    mass flux [kg/m^2/s]       : ', params%patch_mass_flux(p)
         write(unit_id,'(a,a)') '    pressure type              : ', trim(params%patch_pressure_type(p))
         write(unit_id,'(a,es12.5)') '    pressure value             : ', params%patch_p(p)
         write(unit_id,'(a,es12.5)') '    pressure gradient          : ', params%patch_dpdn(p)
         write(unit_id,'(a,a)') '    temperature type           : ', trim(params%patch_temperature_type(p))
         write(unit_id,'(a,es12.5)') '    temperature [K]            : ', params%patch_T(p)
         write(unit_id,'(a,a)') '    species type               : ', trim(params%patch_species_type(p))
         if (len_trim(params%patch_composition(p)) > 0) then
            write(unit_id,'(a,a)') '    composition                : ', trim(params%patch_composition(p))
         end if
      end do
   end subroutine write_boundary_config_summary_to_unit