text:
!##############################################################################
!# ****************************************************************************
!# <name> genoutput </name>
!# ****************************************************************************
!#
!# <purpose>
!# This module contains several routines for input/output of messages to the
!# terminal and/or to log files.
!#
!# The following routines can be found here:
!#
!# 1.) output_init
!#     -> Initialises the output. Opens a log file and assigns it an
!#        output channel OU_LOG.
!#
!# 2.) output_line
!#     -> Writes a message to the terminal and/or the log file.
!#
!# 3.) output_separator
!#     -> Writes a separator line to the terminal and/or the log file.
!#
!# 4.) output_lbrk
!#     -> Writes an empty line / line break to the terminal and/or the log
!#        file.
!#
!# 5.) output_simple
!#     -> Writes a message to the terminal and/or the log file.
!#        For old FEAT compatibility. Uses a priority identifier to decide on
!#        whether to write only to the log file or to both, log file
!#        and terminal.
!#
!# 6.) output_simple_sep
!#     -> Writes a separator line to the terminal and/or the log file.
!#        For old FEAT compatibility. Uses a priority identifier to decide on
!#        whether to write only to the log file or to both, log file
!#        and terminal.
!#
!# 7.) output_multiline
!#     -> Writes a formatted message to the terminal and/or the log file.
!#
!# 8.) output_done
!#     -> Closes the log file, releases all ressources in use.
!#
!# HOW TO USE:
!# -----------
!# As long as output_init is not called, all output is redirected only to the
!# standard output channel. Directly after program start, output_init should
!# be called so that the log file is written.
!#
!# Proceed as follows:
!#
!# <verb>
!# 1.) call output_init ()
!#     -> Initialise the output system for output on the terminal
!#
!#     Alternatively one can use:
!#
!#     call output_init ('mylogfile.txt')
!#     -> Opens a log file 'mylogfile.txt' for the output additional to
!#        terminal output
!#
!#     Alternatively one can use
!#
!#     call output_init ('mylogfile.txt','myerrorlogfile.txt')
!#     -> Opens a log file 'mylogfile.txt' for the output and
!#        'myerrorlogfile.txt' for error output (both additionally to
!#        terminal output).
!#
!# 2.) call output_line ('This is a message')
!#
!#     -> Writes a message to the terminal. If output_init was used
!#        before, output_line_std will also write the message to the log file.
!#
!#     Alternatively, one can specify where to write output to by using
!#     different variants of output_line:
!#
!#     a) call output_lbrk ()
!#
!#     -> Writes an empty line to the terminal and log file. This is the same
!#        as call output_line ('').
!#
!#     b) call output_line('A message only to the terminal.', &
!#                         OU_CLASS_MSG,OU_MODE_TERM)
!#
!#     -> Writes a message only to the terminal
!#
!#     c) call output_line ('A message only to the log file.', &
!#                          OU_CLASS_MSG,OU_MODE_LOG)
!#
!#     -> Writes a message only to the log file
!#
!#     d) call output_line ('A special debug message.', &
!#                          OU_CLASS_TRACE1,OU_MODE_STD,'mysubroutine')
!#
!#     -> Writes a debug message with '*** (mysubroutine):' in front to the
!#        terminal and the log file. This is usually used for debug purposes.
!#
!#     e) call output_line ('This is an error message.', &
!#                          OU_CLASS_ERROR,OU_MODE_STD,'mysubroutine')
!#
!#     -> Writes a debug message with 'Error (mysubroutine):' in front to the
!#        error log file and error output channel.
!#
!#       or even simpler:
!#
!#        call output_line ('This is an error message.',OU_CLASS_ERROR)
!#        call output_line ('This is an warning message.',OU_CLASS_WARNING)
!#
!#     f) call output_separator (OU_SEP_MINUS)
!#
!#     -> Writes a separation line with '-' signs to the terminal / log file
!#
!#     g) MT = 1
!#        call output_simple (MT,'A log file message.')
!#
!#     -> Writes a message to the log file, not to the terminal. FEAT1.0
!#        compatibility routine for "MT=1"
!#
!#     h) MT = 2
!#        call output_simple (MT,'A log file message')
!#
!#     -> Writes a message to the terminal and to the log file. FEAT1.0
!#        compatibility routine for "MT=2".
!#
!#     i) MT = 2
!#        call output_simple (MT)
!#
!#     -> Writes an empty line to the terminal and to the log file. FEAT1.0
!#        compatibility routine for "MT=2".
!#
!#     i) MT = 2
!#        call output_simple_sep (MT,OU_SEP_MINUS)
!#
!#     -> Writes a separation line with '-' signs to the terminal and
!#        to the log file. FEAT1.0 compatibility routine for "MT=2".
!#
!#     For futher possibilities, consider the documentation out output_line.
!#
!# 3.) call output_done()
!#
!#     -> Closes the output channel(s).
!#
!#
!# Logging benchmark data
!# ----------------------
!# The output library furthermore supports the (semi-)automatic output of
!# deterministic data to a so callen 'benchmark log file'. This is typically
!# used to write out data to a specific file which is compared to reference
!# data in regression tests. Such data is meant to be deterministic,
!# reproducable in every run and usually formatted in a 'nice' style that
!# differences to reference results can easily be checked.
!#
!# To activate the benchmark log file, one has to specify an additional
!# parameter in the call to output_init:
!#
!#
!#    call output_init ('mylogfile.txt','myerrorlogfile.txt','benchmarkresultfile')
!#
!#
!# This opens a file 'benchmarkresultfile' where benchmark data is written to.
!# To write a string to this file, the application has to call output_line
!# with an extended output-mode:
!#
!#
!#   call output_line('A message only to the benchmark log.', &
!#                    OU_CLASS_MSG,OU_MODE_BENCHLOG)
!#
!#
!# This writes a message directly to the benchmark log file. It is also possible
!# to write out data to the standard terminal plus the benchmark log file
!# by combining the output mode constants:
!#
!#
!#   call output_line('A message only to the benchmark log.', &
!#                    OU_CLASS_MSG,OU_MODE_STD+OU_MODE_BENCHLOG)
!#
!#
!# In all cases, the additional constant OU_MODE_BENCHLOG must be manually
!# specified as logging to the benchmark log file is not done automatically.
!#
!# Adding date/time to the output
!# ------------------------------
!# The output module supports automatic adding of the current date and/or
!# time to the output. This helps to keep track of the execution time
!# and allows basic profiling of the application. To enable the output
!# of date/time data, it can be activated using the global
!# cdefaultDateTimeLogPolicy variable. Using cdatetimeLogFormat allows to
!# change the format of the output.
!#
!# Example: The following code enables printing of date and time:
!#
!#
!#    call output_init()
!#    cdefaultDateTimeLogPolicy = OU_DTP_ADDDATETIME
!#
!#
!# Using the optional variable cdateTimeLogPolicy in the output subroutines
!# allows to manually define whether to print date/time or not. This can be
!# used to switch off the output if some lines are printed without line break.
!# Example:
!#
!#
!#    call output_init()
!#    cdefaultDateTimeLogPolicy = OU_DTP_ADDDATETIME
!#
!#    call output_line ("Hello")
!#
!#    call output_line ("This is some output.")
!#
!#    call output_line ("This sentence is incorrectly printed ", bnolinebreak=.true.)
!#    call output_line ("in two statements.")
!#
!#    call output_line ("This sentence is correctly printed ", bnolinebreak=.true.)
!#    call output_line ("in two statements.", cdateTimeLogPolicy = OU_DTP_NONE)
!#
!#
!# This gives the output:
!#
!# <verb>
!#    02-06-2010 16:21:00: Hello.
!#    02-06-2010 16:21:00: This is some output.
!#    02-06-2010 16:21:00: This sentence is incorrectly printed 02-06-2010 16:21:00: in two statements.
!#    02-06-2010 16:21:00: This sentence is correctly printed in two statements.
!# </verb>
!#
!# If cdefaultDateTimeLogPolicy is not set, specifying cdateTimeLogPolicy=OU_DTP_NONE
!# as parameter in the output does not harm since it has no effect.
!#
!# </purpose>
!##############################################################################

module genoutput

!$ use omp_lib
  use fsystem

  implicit none

  private

!<constants>

!<constantblock description="Output mode. Decides on where the output is written to.">

  ! Output mode: Write to main log file / error log file (depending on whether
  ! a message is considered as an error or not by OU_CLASS_ERROR/OU_CLASS_WARNING)
  integer(I32), parameter, public :: OU_MODE_LOG      = 2_I32**0

  ! Output mode: Write to terminal
  integer(I32), parameter, public :: OU_MODE_TERM     = 2_I32**1

  ! Output mode: Write to benchmark log file
  integer(I32), parameter, public :: OU_MODE_BENCHLOG = 2_I32**2

  ! Output mode: Write to both, log file and terminal
  integer(I32), parameter, public :: OU_MODE_STD      = OU_MODE_LOG+OU_MODE_TERM

!</constantblock>

!<constantblock description="Output classification. Prints an additional classification string\\
!                            when writing a string to a log file / terminal">

  ! Output classification: Standard message
  integer, parameter, public :: OU_CLASS_MSG      = 0

  ! Output classification: Trace information, level 1
  integer, parameter, public :: OU_CLASS_TRACE1   = 1

  ! Output classification: Trace information, level 2
  integer, parameter, public :: OU_CLASS_TRACE2   = 2

  ! Output classification: Trace information, level 3
  integer, parameter, public :: OU_CLASS_TRACE3   = 3

  ! Output classification: System/Timer message
  integer, parameter, public :: OU_CLASS_SYSTEM   = 4

  ! Output classification: Error message.
  integer, parameter, public :: OU_CLASS_ERROR    = 5

  ! Output classification: Warning message.
  integer, parameter, public :: OU_CLASS_WARNING  = 6

!</constantblock>

!<constantblock description="Type of separator line">

  ! Separator line: MINUS character
  integer, parameter, public :: OU_SEP_MINUS  = 0

  ! Separator line: STAR character
  integer, parameter, public :: OU_SEP_STAR   = 1

  ! Separator line: EQUAL character
  integer, parameter, public :: OU_SEP_EQUAL  = 2

  ! Separator line: DOLLAR character
  integer, parameter, public :: OU_SEP_DOLLAR = 3

  ! Separator line: @ character
  integer, parameter, public :: OU_SEP_AT     = 4

  ! Separator line: + character
  integer, parameter, public :: OU_SEP_PLUS   = 5

  ! Separator line: ~ character
  integer, parameter, public :: OU_SEP_TILDE  = 6

  ! Separator line: & character
  integer, parameter, public :: OU_SEP_AMPAND = 7

  ! Separator line: % character
  integer, parameter, public :: OU_SEP_PERC   = 8

  ! Separator line: # character
  integer, parameter, public :: OU_SEP_HASH   = 9

!</constantblock>

!<constantblock>

  ! Length of a line on the terminal / in the log file.
  ! Standard value = 80 characters
  integer, public :: OU_LINE_LENGTH         = 80

  ! Global device number for terminal output
  integer, public :: OU_TERMINAL            = 6

  ! Global device number for terminal input
  integer, public :: IN_TERMINAL            = 5

  ! Global device number for errors on terminal
  integer, public :: OU_ERROR               = 6

  ! Global device number for log file. The log file is opened in output_init.
  ! <=0: no log file output.
  integer, public :: OU_LOG                 = 0

  ! Global device number for error log file. The error log file is opened
  ! in output_init and usually coincides with OU_LOG.
  ! <=0: write error messages to standard log file.
  integer, public :: OU_ERRORLOG            = 0

  ! Global device number for benchmark log file. The error log file is opened
  ! in output_init and usually coincides with OU_LOG.
  ! <=0: write error messages to standard log file.
  integer, public :: OU_BENCHLOG            = 0

!</constantblock>

!<constantblock description = "Constants for setting up adding date/time to each message.">

  ! Do not add date/time.
  integer, parameter, public :: OU_DTP_NONE            = 0

  ! Add date.
  integer, parameter, public :: OU_DTP_ADDDATE         = 1

  ! Add time.
  integer, parameter, public :: OU_DTP_ADDTIME         = 2

  ! Add date + time.
  integer, parameter, public :: OU_DTP_ADDDATETIME     = 3

!</constantblock>

!</constants>

!<publicvars>

  ! Current benchmark file logging policy.
  ! =0: None. Do not log anything to the benchmark log file.
  ! =1: Standard. Write data to benchmark log file if OU_BENCHLOG is specified
  !     in coutputMode.
  ! =2: Full. Write everything into the benchmark log file independent on whether
  !     OU_BENCHLOG is specified in coutputMode or not.
  integer, public, save :: cbenchLogPolicy = 1

  ! Default date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  integer, public, save :: cdefaultDateTimeLogPolicy = OU_DTP_NONE

  ! Date/Time format.
  ! =0: english: MM-DD-YYYY HH:MM:SS
  ! =1: german: DD.MM.YYYY HH:MM:SS
  integer, public, save :: cdatetimeLogFormat = 0

  ! Defines an automatic indention by spaces which is put in front of
  ! all messages (After a possible date).
  ! Must be in the range 0..255.
  integer, public, save :: output_iautoOutputIndent = 0

!</publicvars>

  ! Empty string used for indented output.
  ! By Fortran convention, this automatically expands to spaces.
  character(len=256), save :: semptystring = ""

  interface output_line
    module procedure output_line_std
    module procedure output_line_feast
  end interface

  interface output_init
    module procedure output_init_simple
    module procedure output_init_logfile
    module procedure output_init_standard
  end interface

  public :: output_init
  public :: output_line
  public :: output_separator
  public :: output_lbrk
  public :: output_simple
  public :: output_simple_sep
  public :: output_multiline
  public :: output_done

contains

!************************************************************************

!<subroutine>
  subroutine output_openLogfile(sfilename, iunit)

  !<description>
    ! This routine opens a file sfilename for writing. This is
    ! usually used to create a log file. The routine will return a
    ! handle iunit with an output channel where text can be written to.
    ! If the file exists, it is overwritten.
    !
    ! As long as this routine is not called, all output is redirected
    ! only to the standard output channel.
  !</description>

  !<input>

    ! Name of the file to overwrite.
    character(*), intent(in) :: sfilename

  !</input>

  !<output>
    !unit of the opened file
    integer, intent(out) :: iunit
  !</output>

!</subroutine>

    integer :: i
    integer :: istatus ! status variable for opening procedure
    character(len=len(sfilename)) :: sfilepath

    ! use wrapper for C system call stat() in kernel/System/isdirectory.c
    external isdirectory
    ! use wrapper for C system call mkdir() in kernel/System/mkdir_recursive.c
    external mkdir_recursive

    if (len_trim(sfilename) .eq. 0) then
      write (*,'(A)') 'Error: output_openLogfile. sfilename undefined!'
      return
    endif

    iunit = sys_getFreeUnit()
    if (iunit .eq. -1) then
      write (*,'(A)') 'Error: output_openLogfile. No output channel available!'
      return
    endif

    ! Ensure that the path up to the given file name does exist
    ! In short, this could be realised by
    !    call io_pathExtract(sfilename, sfilepath)
    !    if (.not. io_isDirectory(sfilepath)) then
    !      call io_makeDirectory(sfilepath)
    !    end if
    ! but as this would introduce circular module dependencies inline ethe relevant code
    ! here manually.
    i = scan(sfilename, "/\\", .true.)
    if (i .ne. 0) then
      ! Directory ends at position i.
      sfilepath = sfilename(1:i-1)
      call isdirectory(trim(sfilepath) // achar(0), i, istatus)
      if (istatus .ge. 0) then
        if (i .le. 0) then
          ! Path does not exist, create it
          call mkdir_recursive(trim(sfilepath) // achar(0), istatus)
          if (istatus .lt. 0) then
            write (*,'(A)') 'Error: output_openLogfile. Could not (auto-)create path to',&
                            ' output file "', trim(sfilename), '". ***'
          endif
        endif
      end if

    else
      ! No subdirectory contained in log file, no path needs creation
    end if

    open(unit=iunit, file=trim(sfilename), iostat=istatus, status="replace", &
          action="write")

    if (istatus .ne. 0) then
      write (unit=*, fmt=*) &
          'Error: output_openLogfile. Error while opening file "',&
          trim(sfilename), '". ***'
      iunit = -1
    endif

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_init_simple ()

!<description>
  ! Initialises the output system. The output system is configured to show all
  ! messages on the terminal, no log file is opened.
!</description>

!</subroutine>

    call output_init_standard ("","")

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_init_logfile (slogFilename)

!<description>
  ! Initialises the output system.
  ! If sfilename is given, a log file with that name is opened and all log and error
  ! messages are written to it. If not given, all output is directed to the
  ! terminal only.
!</description>

!<input>

  ! Name of a log file for standard messages. If "" is specified,
  ! output messages are written to the standard output.
  character(LEN=*), intent(in) :: slogFilename

!</input>

!</subroutine>

    call output_init_standard (slogFilename,"")

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_init_standard (slogFilename,serrorFilename,sbenchLogfile)

!<description>
  ! Initialises the output system.
  ! If sfilename is given, a log file with that name is opened and all log
  ! messages are written to it. If not given, all output is directed to the
  ! terminal only.
  ! If serrorFilename is given, an error log file with filename serrorFilename
  ! is opened and all error messages are redirected to it. If not given,
  ! the error messages are directed to the log file -- or to the terminal, if
  ! slogFilename does not exist.
!</description>

!<input>

  ! Name of a log file for standard messages. If ""
  ! is specified, output messages are written to the standard output.
  character(LEN=*), intent(in) :: slogFilename

  ! Name of an log file for error messages. If ""
  ! is specified, errors are written to the standard log file. The name of
  ! the file may also coincide with slogFilename.
  character(LEN=*), intent(in) :: serrorFilename

  ! OPTIONAL: Name of an log file for deterministic benchmark messages. If
  ! not present or set to "", benchmark messages are not written out. The
  ! name of the file may also coincide with slogFilename.
  character(LEN=*), intent(in), optional :: sbenchLogfile

!</input>

!</subroutine>

    ! Close previously opened log files.
    call output_done ()

    ! Name of the standard logfile given?
    if (slogFilename .ne. "") then

      ! Open a log file
      call output_openLogfile(slogFilename, OU_LOG)

    end if

    ! Name of the error logfile given?
    if (serrorFilename .ne. "") then

      ! Both filenames the same?
      if (slogFilename .eq. serrorFilename) then
        OU_ERRORLOG = OU_LOG
      else
        ! Open an error log file
        call output_openLogfile(serrorFilename, OU_ERRORLOG)
      end if

    else

      ! Write error messages to standard log file.
      OU_ERRORLOG = OU_LOG

    end if

    ! Name of the benchmark logfile given?
    if (present(sbenchLogfile)) then
      if (sbenchLogfile .ne. "") then

        ! Both filenames the same?
        if (slogFilename .eq. serrorFilename) then
          OU_BENCHLOG = OU_LOG
        else
          ! Open a benchmark log file
          call output_openLogfile(sbenchLogfile, OU_BENCHLOG)
        end if

      else

        ! No benchmark output
        OU_BENCHLOG = 0

      end if

    else

      ! No benchmark output
      OU_BENCHLOG = 0

    end if

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_multiline (smessage, &
                               coutputClass, coutputMode, ssubroutine, &
                               bnolinebreak, bnotrim, cdateTimeLogPolicy)

!<description>
  ! Writes a formatted output message to the terminal, log file or error log
  ! file, depending on the input parameters.
  ! smessage is the message to be written out.
  ! coutputMode decides (if given) about whether the output if written to file,
  ! terminal or both.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
!</description>

!<input>
  ! The message to be written out.
  character(LEN=*), intent(in) :: smessage

  ! OPTIONAL: Output mode. One of the OU_MODE_xxxx constants. If not specified,
  ! OU_MODE_STD is assumed.
  integer(I32), intent(in), optional :: coutputMode

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: When specifying bnolinebreak=TRUE, the output routine will
  ! not perform a line break after printing.
  logical, intent(in), optional :: bnolinebreak

  ! OPTIONAL: When specifying bnotrim=TRUE, the output routine will
  ! not trim the string when printing. This does only work for
  ! coutputClass=OU_CLASS_MSG and coutputClass=OU_CLASS_ERROR
  ! (or if coutputClass is not specified)!
  logical, intent(in), optional :: bnotrim

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy

!</input>

!</subroutine>

  ! local variables
  integer(I32) :: coMode
  integer :: coClass, iofChannel, iotChannel
  integer :: iindent,istart,ilinelen
  logical :: bntrim, bnnewline
  character(LEN=len(smessage)+20+SYS_NAMELEN+10+8+3) :: smsg

    ! Get the actual parameters.

    coMode = OU_MODE_STD
    coClass = OU_CLASS_MSG
    bntrim = .false.
    bnnewline = .false.
    iindent = min(255,output_iautoOutputIndent)

    if (present(coutputMode))  coMode = coutputMode
    if (present(coutputClass)) coClass = coutputClass
    if (present(bnotrim))      bntrim = bnotrim
    if (present(bnolinebreak)) bnnewline = bnolinebreak

    ! Get the file and terminal output channel
    iotChannel = OU_TERMINAL
    if ((coClass .eq. OU_CLASS_ERROR) .or. &
        (coClass .eq. OU_CLASS_WARNING)) iotChannel = OU_ERROR

    iofChannel = OU_LOG
    if ((coClass .eq. OU_CLASS_ERROR) .or. &
        (coClass .eq. OU_CLASS_WARNING)) iofChannel = OU_ERRORLOG

    if (bntrim .and. &
        ((coClass .eq. OU_CLASS_MSG) .or. &
         (coClass .eq. OU_CLASS_WARNING) .or. &
         (coClass .eq. OU_CLASS_ERROR)) ) then

      ! Where to write the message to?
      if ((iand(coMode,OU_MODE_TERM) .ne. 0) .and. (iotChannel .gt. 0)) then
        istart=1; ilinelen=index(smessage,"\n")
        do while(ilinelen.ne.0)
          if (ilinelen.eq.1) then
            write (iotChannel,'(A)') ''
            istart=istart+ilinelen+1
          else
            write (iotChannel,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage(istart:istart+ilinelen-2)
            istart=istart+ilinelen-1
          end if
          ilinelen=index(smessage(istart:),"\n")
        end do
        if (bnnewline) then
          write (iotChannel,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage(istart:)
        else
          write (iotChannel,'(A)') semptystring(1:iindent)//smessage(istart:)
        end if
      end if

      ! Output to log file?
      if ((iand(coMode,OU_MODE_LOG) .ne. 0) .and. (iofChannel .gt. 0)) then
        istart=1; ilinelen=index(smessage,"\n")
        do while(ilinelen.ne.0)
          if (ilinelen.eq.1) then
            write (iofChannel,'(A)') ''
            istart=istart+ilinelen+1
          else
            write (iofChannel,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage(istart:istart+ilinelen-2)
            istart=istart+ilinelen-1
          end if
          ilinelen=index(smessage(istart:),"\n")
        end do
        if (bnnewline) then
          write (iofChannel,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage(istart:)
        else
          write (iofChannel,'(A)') semptystring(1:iindent)//smessage(istart:)
        end if
      end if

      ! Output to the benchmark log file?
      if (OU_BENCHLOG .gt. 0) then
        if ((cbenchLogPolicy .eq. 2) .or. &
            (cbenchLogPolicy .eq. 1) .and. (iand(coMode,OU_MODE_BENCHLOG) .ne. 0)) then
          istart=1; ilinelen=index(smessage,"\n")
          do while(ilinelen.ne.0)
            if (ilinelen.eq.1) then
              write (OU_BENCHLOG,'(A)') ''
              istart=istart+ilinelen+1
            else
              write (OU_BENCHLOG,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage(istart:istart+ilinelen-2)
              istart=istart+ilinelen-1
            end if
            ilinelen=index(smessage(istart:),"\n")
          end do
          if (bnnewline) then
            write (OU_BENCHLOG,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage(istart:)
          else
            write (OU_BENCHLOG,'(A)') semptystring(1:iindent)//smessage(istart:)
          end if
        end if
      end if

    else

      istart=1; ilinelen=index(smessage,"\n")
      do while(ilinelen.ne.0)
        if (ilinelen.eq.1) then
          ! Build the actual error message
          smsg = output_reformatMsg (semptystring(1:iindent), &
              coutputClass, ssubroutine, cdateTimeLogPolicy)

          ! Where to write the new line to?
          if ((iand(coMode,OU_MODE_TERM) .ne. 0) .and. (iotChannel .gt. 0)) then
            write (iotChannel,'(A)') trim(smsg)
          end if

          ! Output to log file?
          if ((iand(coMode,OU_MODE_LOG) .ne. 0) .and. (iofChannel .gt. 0)) then
            write (iofChannel,'(A)') trim(smsg)
          end if

          ! Output to benchmark log file?
          if (OU_BENCHLOG .gt. 0) then
            if ((cbenchLogPolicy .eq. 2) .or. &
                (cbenchLogPolicy .eq. 1) .and. (iand(coMode,OU_MODE_BENCHLOG) .ne. 0)) then
              write (OU_BENCHLOG,'(A)') trim(smsg)
            end if
          end if
          istart=istart+ilinelen+1
        else
          ! Build the actual error message
          smsg = output_reformatMsg (semptystring(1:iindent)//smessage(istart:istart+ilinelen-2), &
              coutputClass, ssubroutine, cdateTimeLogPolicy)

          ! Where to write the message to?
          if ((iand(coMode,OU_MODE_TERM) .ne. 0) .and. (iotChannel .gt. 0)) then
            write (iotChannel,'(A)',ADVANCE='NO') trim(smsg)
          end if

          ! Output to log file?
          if ((iand(coMode,OU_MODE_LOG) .ne. 0) .and. (iofChannel .gt. 0)) then
            write (iofChannel,'(A)',ADVANCE='NO') trim(smsg)
          end if

          ! Output to benchmark log file?
          if (OU_BENCHLOG .gt. 0) then
            if ((cbenchLogPolicy .eq. 2) .or. &
                (cbenchLogPolicy .eq. 1) .and. (iand(coMode,OU_MODE_BENCHLOG) .ne. 0)) then
              write (OU_BENCHLOG,'(A)',ADVANCE='NO') trim(smsg)
            end if
          end if
          istart=istart+ilinelen-1
        end if
        ilinelen=index(smessage(istart:),"\n")
      end do

      ! Build the actual error message
      smsg = output_reformatMsg (semptystring(1:iindent)//smessage(istart:), &
          coutputClass, ssubroutine, cdateTimeLogPolicy)

      ! Where to write the message to?
      if ((iand(coMode,OU_MODE_TERM) .ne. 0) .and. (iotChannel .gt. 0)) then
        if (bnnewline) then
          write (iotChannel,'(A)',ADVANCE='NO') trim(smsg)
        else
          write (iotChannel,'(A)') trim(smsg)
        end if
      end if

      ! Output to log file?
      if ((iand(coMode,OU_MODE_LOG) .ne. 0) .and. (iofChannel .gt. 0)) then
        if (bnnewline) then
          write (iofChannel,'(A)',ADVANCE='NO') trim(smsg)
        else
          write (iofChannel,'(A)') trim(smsg)
        end if
      end if

      ! Output to benchmark log file?
      if (OU_BENCHLOG .gt. 0) then
        if ((cbenchLogPolicy .eq. 2) .or. &
            (cbenchLogPolicy .eq. 1) .and. (iand(coMode,OU_MODE_BENCHLOG) .ne. 0)) then
          if (bnnewline) then
            write (OU_BENCHLOG,'(A)',ADVANCE='NO') trim(smsg)
          else
            write (OU_BENCHLOG,'(A)') trim(smsg)
          end if
        end if
      end if

    end if

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_done ()

!<description>
  ! Closes all log files, cleans up output system.
!</description>

!</subroutine>

    ! Close all open channels.
    if (OU_LOG .gt. 0) close(OU_LOG)
    if (OU_ERRORLOG .gt. 0) close(OU_ERRORLOG)
    if (OU_BENCHLOG .gt. 0) close(OU_BENCHLOG)

  end subroutine

!************************************************************************************

!<function>

  function output_reformatMsg (smessage, coutputClass, ssubroutine, cdateTimeLogPolicy) result (s)

!<description>
  ! Reformats an output message according to an output class.
!</description>

!<input>
  ! The message to reformat.
  character(LEN=*), intent(in) :: smessage

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of a subroutine to include in the message; may be "".
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy
!</input>

!<result>
  ! Reformatted string
!</result>

  character(LEN=len(smessage)+20+SYS_NAMELEN+10+8+3) :: s

!</function>

    logical :: bsub
    integer :: ioc,cdateTime,iindent
    character(LEN=8) :: sdate
    character(LEN=10) :: stime

    iindent = min(255,output_iautoOutputIndent)

    bsub = .false.
    if (present (ssubroutine)) then
      bsub = (ssubroutine .ne. '')
    end if

    cdateTime = cdefaultDateTimeLogPolicy
    if (present(cdateTimeLogPolicy)) cdateTime = cdateTimeLogPolicy

    ioc = OU_CLASS_MSG
    if (present(coutputClass)) ioc = coutputClass

    if (.not. bsub) then

      select case (ioc)
      case (OU_CLASS_TRACE1)
        s = semptystring(1:iindent)//'*** '//adjustl(smessage)
      case (OU_CLASS_TRACE2)
        s = semptystring(1:iindent)//'***** '//adjustl(smessage)
      case (OU_CLASS_TRACE3)
        s = semptystring(1:iindent)//'******* '//adjustl(smessage)
      case (OU_CLASS_SYSTEM)
        s = semptystring(1:iindent)//'* System: '//adjustl(smessage)
      case (OU_CLASS_ERROR)
        s = semptystring(1:iindent)//'* Error: '//adjustl(smessage)
      case (OU_CLASS_WARNING)
        s = semptystring(1:iindent)//'* Warning: '//adjustl(smessage)
      case default
        s = smessage
      end select

    else

      select case (ioc)
      case (OU_CLASS_TRACE1)
        s = semptystring(1:iindent)//'*** '//trim(ssubroutine)//': '//adjustl(smessage)
      case (OU_CLASS_TRACE2)
        s = semptystring(1:iindent)//'***** '//trim(ssubroutine)//': '//adjustl(smessage)
      case (OU_CLASS_TRACE3)
        s = semptystring(1:iindent)//'******* '//trim(ssubroutine)//': '//adjustl(smessage)
      case (OU_CLASS_SYSTEM)
        s = semptystring(1:iindent)//'* System ('//trim(ssubroutine)//'): '//adjustl(smessage)
      case (OU_CLASS_ERROR)
        s = semptystring(1:iindent)//'* Error ('//trim(ssubroutine)//'): '//adjustl(smessage)
      case (OU_CLASS_WARNING)
        s = semptystring(1:iindent)//'* Warning ('//trim(ssubroutine)//'): '//adjustl(smessage)
      case default
        s = smessage
      end select

    end if

    if (cdateTime .ne. OU_DTP_NONE) then

      ! Get date and time.
      call date_and_time(sdate,stime)

      ! Reformat the message.
      select case (cdatetimeLogFormat)
      case (1)
        select case (cdateTime)
        case (OU_DTP_NONE)
        case (OU_DTP_ADDDATE)
          s = sdate(7:8)//"."//sdate(5:6)//"."//sdate(1:4)//": "//s
        case (OU_DTP_ADDTIME)
          s = stime(1:2)//":"//stime(3:4)//":"//stime(5:6)//": "//s
        case (OU_DTP_ADDDATETIME)
          s = sdate(7:8)//"."//sdate(5:6)//"."//sdate(1:4)//" "// &
              stime(1:2)//":"//stime(3:4)//":"//stime(5:6)//": "//s
        end select
      case default
        select case (cdateTime)
        case (OU_DTP_NONE)
        case (OU_DTP_ADDDATE)
          s = sdate(5:6)//"-"//sdate(7:8)//"-"//sdate(1:4)//": "//s
        case (OU_DTP_ADDTIME)
          s = stime(1:2)//":"//stime(3:4)//":"//stime(5:6)//": "//s
        case (OU_DTP_ADDDATETIME)
          s = sdate(5:6)//"-"//sdate(7:8)//"-"//sdate(1:4)//" "// &
              stime(1:2)//":"//stime(3:4)//":"//stime(5:6)//": "//s
        end select
      end select

    end if

  end function

!************************************************************************************

!<subroutine>

  subroutine output_line_std (smessage, &
                              coutputClass, coutputMode, ssubroutine, &
                              bnolinebreak, bnotrim, cdateTimeLogPolicy)

!<description>
  ! Writes an output message to the terminal, log file or error log file,
  ! depending on the input parameters.
  ! smessage is the message to be written out.
  ! coutputMode decides (if given) about whether the output if written to file,
  ! terminal or both.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
!</description>

!<input>
  ! The message to be written out.
  character(LEN=*), intent(in) :: smessage

  ! OPTIONAL: Output mode. One of the OU_MODE_xxxx constants. If not specified,
  ! OU_MODE_STD is assumed.
  integer(I32), intent(in), optional :: coutputMode

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: When specifying bnolinebreak=TRUE, the output routine will
  ! not perform a line break after printing.
  logical, intent(in), optional :: bnolinebreak

  ! OPTIONAL: When specifying bnotrim=TRUE, the output routine will
  ! not trim the string when printing. This does only work for
  ! coutputClass=OU_CLASS_MSG and coutputClass=OU_CLASS_ERROR
  ! (or if coutputClass is not specified)!
  logical, intent(in), optional :: bnotrim

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy

!</input>

!</subroutine>

  ! local variables
  integer(I32) :: coMode
  integer :: coClass, iofChannel, iotChannel
  integer :: iindent
  logical :: bntrim, bnnewline
  character(LEN=len(smessage)+20+SYS_NAMELEN+10+8+3) :: smsg

    ! Get the actual parameters.

    coMode = OU_MODE_STD
    coClass = OU_CLASS_MSG
    bntrim = .false.
    bnnewline = .false.
    iindent = min(255,output_iautoOutputIndent)

    if (present(coutputMode))  coMode = coutputMode
    if (present(coutputClass)) coClass = coutputClass
    if (present(bnotrim))      bntrim = bnotrim
    if (present(bnolinebreak)) bnnewline = bnolinebreak

    ! Get the file and terminal output channel
    iotChannel = OU_TERMINAL
    if ((coClass .eq. OU_CLASS_ERROR) .or. &
        (coClass .eq. OU_CLASS_WARNING)) iotChannel = OU_ERROR

    iofChannel = OU_LOG
    if ((coClass .eq. OU_CLASS_ERROR) .or. &
        (coClass .eq. OU_CLASS_WARNING)) iofChannel = OU_ERRORLOG

    if (bntrim .and. &
        ((coClass .eq. OU_CLASS_MSG) .or. &
         (coClass .eq. OU_CLASS_WARNING) .or. &
         (coClass .eq. OU_CLASS_ERROR)) ) then

      ! Where to write the message to?
      if ((iand(coMode,OU_MODE_TERM) .ne. 0) .and. (iotChannel .gt. 0)) then
        if (bnnewline) then
          write (iotChannel,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage
        else
          write (iotChannel,'(A)') semptystring(1:iindent)//smessage
        end if
      end if

      ! Output to log file?
      if ((iand(coMode,OU_MODE_LOG) .ne. 0) .and. (iofChannel .gt. 0)) then
        if (bnnewline) then
          write (iofChannel,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage
        else
          write (iofChannel,'(A)') semptystring(1:iindent)//smessage
        end if
      end if

      ! Output to the benchmark log file?
      if (OU_BENCHLOG .gt. 0) then
        if ((cbenchLogPolicy .eq. 2) .or. &
            (cbenchLogPolicy .eq. 1) .and. (iand(coMode,OU_MODE_BENCHLOG) .ne. 0)) then
          if (bnnewline) then
            write (OU_BENCHLOG,'(A)',ADVANCE='NO') semptystring(1:iindent)//smessage
          else
            write (OU_BENCHLOG,'(A)') semptystring(1:iindent)//smessage
          end if
        end if
      end if

    else

      ! Build the actual error message
      smsg = output_reformatMsg (semptystring(1:iindent)//smessage, &
          coutputClass, ssubroutine, cdateTimeLogPolicy)

      ! Where to write the message to?
      if ((iand(coMode,OU_MODE_TERM) .ne. 0) .and. (iotChannel .gt. 0)) then
        if (bnnewline) then
          write (iotChannel,'(A)',ADVANCE='NO') trim(smsg)
        else
          write (iotChannel,'(A)') trim(smsg)
        end if
      end if

      ! Output to log file?
      if ((iand(coMode,OU_MODE_LOG) .ne. 0) .and. (iofChannel .gt. 0)) then
        if (bnnewline) then
          write (iofChannel,'(A)',ADVANCE='NO') trim(smsg)
        else
          write (iofChannel,'(A)') trim(smsg)
        end if
      end if

      ! Output to benchmark log file?
      if (OU_BENCHLOG .gt. 0) then
        if ((cbenchLogPolicy .eq. 2) .or. &
            (cbenchLogPolicy .eq. 1) .and. (iand(coMode,OU_MODE_BENCHLOG) .ne. 0)) then
          if (bnnewline) then
            write (OU_BENCHLOG,'(A)',ADVANCE='NO') trim(smsg)
          else
            write (OU_BENCHLOG,'(A)') trim(smsg)
          end if
        end if
      end if
    end if

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_line_feast (coutputClass, ssubroutine, smsg)

!<description>
  ! Writes an output message to the terminal, log file or error log file,
  ! depending on the input parameters.
  ! smsg is the message to be written out.
  ! coutputMode decides (if given) about whether the output if written to file,
  ! terminal or both.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
  !
  ! Compatibility function for FEAST routines.
!</description>

!<input>
  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in) :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in) :: ssubroutine

  ! The message to be written out.
  character(LEN=*), intent(in) :: smsg

!</input>

!</subroutine>

    ! REMARK: Do not rename 'smsg' to 'smessage' -- this does not comply
    ! with the Fortran standard and will give you an ambigious interface
    ! in output_line as Fortran cannot distinguish between output_line_std
    ! and output_line_feast anymore in that case! (as all parameter names
    ! are the same).

    call output_line_std (smsg, coutputClass, OU_MODE_STD, ssubroutine)

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_lbrk (coutputClass, coutputMode, ssubroutine, nlbrk, cdateTimeLogPolicy)

!<description>
  ! Writes a line break to the terminal, log file or error log file,
  ! depending on the input parameters.
  ! coutputMode decides (if given) about whether the output if written to file,
  ! terminal or both.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
!</description>

!<input>
  ! OPTIONAL: Output mode. One of the OU_MODE_xxxx constants. If not specified,
  ! OU_MODE_STD is assumed.
  integer(I32), intent(in), optional :: coutputMode

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: Number of linebreaks
  integer, intent(in), optional :: nlbrk

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy
!</input>

!</subroutine>

    ! local variables
    integer :: ilbrk

    if (present(nlbrk)) then

      do ilbrk = 1, nlbrk
        call output_line_std ('', coutputClass, coutputMode, ssubroutine,&
            cdateTimeLogPolicy=cdateTimeLogPolicy)
      end do

    else

      call output_line_std ('', coutputClass, coutputMode, ssubroutine,&
            cdateTimeLogPolicy=cdateTimeLogPolicy)

    end if

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_separator (csepType, coutputClass, coutputMode, ssubroutine, cdateTimeLogPolicy)

!<description>
  ! Writes a separator line to the terminal, log file or error log file,
  ! depending on the input parameters. The parameter csepType decides
  ! on the type of separator line.
  ! coutputMode decides (if given) about whether the output if written to file,
  ! terminal or both.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
!</description>

!<input>
  ! Type of the separator line. One of the OU_SEP_xxxx constants
  ! (OU_SEP_MINUS, OU_SEP_STAR, OU_SEP_EQUAL,...).
  integer, intent(in) :: csepType

  ! OPTIONAL: Output mode. One of the OU_MODE_xxxx constants. If not specified,
  ! OU_MODE_STD is assumed.
  integer(I32), intent(in), optional :: coutputMode

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy
!</input>

!</subroutine>

    ! local variables
    character(LEN=max(1,OU_LINE_LENGTH-1)) :: cstr
    character(LEN=20) :: saux
    integer :: isub,cdateTime

    cdateTime = cdefaultDateTimeLogPolicy
    if (present(cdateTimeLogPolicy)) cdateTime = cdateTimeLogPolicy

    ! Determine how much to reduce the line length.
    isub = 0
    select case (cdateTime)
    case (OU_DTP_NONE)
    case (OU_DTP_ADDDATE)
      isub = 12
    case (OU_DTP_ADDTIME)
      isub = 10
    case (OU_DTP_ADDDATETIME)
      isub = 21
    end select

    ! Create the separator line
    select case (csepType)
    case (OU_SEP_MINUS)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''-''))'
    case (OU_SEP_PLUS)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''+''))'
    case (OU_SEP_STAR)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''*''))'
    case (OU_SEP_EQUAL)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''=''))'
    case (OU_SEP_DOLLAR)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''$''))'
    case (OU_SEP_AT)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''@''))'
    case (OU_SEP_TILDE)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''~''))'
    case (OU_SEP_AMPAND)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''&''))'
    case (OU_SEP_PERC)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''%''))'
    case (OU_SEP_HASH)
      write (saux,'(A,I3,A)') '(',len(cstr)-isub,'(''#''))'
    case default
      write (unit=*, fmt=*) 'output_separator: Unknown separator type: ',&
                            csepType
      call sys_halt()
    end select

    write (cstr,saux)

    call output_line_std (trim(cstr), coutputClass, coutputMode, ssubroutine)

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_simple (ioutputLevel, smsg, coutputClass, ssubroutine, cdateTimeLogPolicy)

!<description>
  ! Writes an output message to the terminal, log file or error log file,
  ! depending on the input parameters.
  ! smsg is the message to be written out.
  ! ioutputLevel decides on where the message is written to, to the terminal
  ! and/or the log file.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
!</description>

!<input>
  ! Output level. DEcides on where to write the message to.
  ! <=0: No output,
  !  =1: write to log file.
  ! >=2: write to log file and terminal.
  integer, intent(in) :: ioutputLevel

  ! OPTIONAL: The message to be written out.
  character(LEN=*), intent(in), optional :: smsg

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy
!</input>

!</subroutine>

  integer(I32) :: coMode

    select case (ioutputLevel)
    case (:0)
      coMode = 0
    case (1)
      coMode = OU_MODE_LOG
    case (2:)
      coMode = OU_MODE_LOG+OU_MODE_TERM
    end select

    if (present(smsg)) then
      call output_line_std (smsg, coutputClass, coMode, ssubroutine, &
          cdateTimeLogPolicy=cdateTimeLogPolicy)
    else
      call output_line_std ('', coutputClass, coMode, ssubroutine, &
          cdateTimeLogPolicy=cdateTimeLogPolicy)
    end if

  end subroutine

!************************************************************************************

!<subroutine>

  subroutine output_simple_sep (ioutputLevel, csepType, coutputClass, ssubroutine, cdateTimeLogPolicy)

!<description>
  ! Writes a separator line to the terminal, log file or error log file,
  ! depending on the input parameters. The parameter csepType decides
  ! on the type of separator line.
  ! ioutputLevel decides on where the message is written to, to the terminal
  ! and/or the log file.
  ! coutputClass classifies (if given) the message as standard message, trace
  ! or error message.
!</description>

!<input>
  ! Output level. DEcides on where to write the message to.
  ! <=0: No output,
  !  =1: write to log file.
  ! >=2: write to log file and terminal.
  integer, intent(in) :: ioutputLevel

  ! Type of the separator line. One of the OU_SEP_xxxx constants
  ! (OU_SEP_MINUS, OU_SEP_STAR or OU_SEP_EQUAL).
  integer, intent(in) :: csepType

  ! OPTIONAL: Output classification. One of the OU_CLASS_xxxx constants.
  ! If not specified, OU_CLASS_MSG is assumed.
  integer, intent(in), optional :: coutputClass

  ! OPTIONAL: Name of the subroutine that calls this function
  character(LEN=*), intent(in), optional :: ssubroutine

  ! OPTIONAL: Date/time appending flag. Allows to configure the
  ! output submodule to automatically add the current date/time to the output.
  ! =OU_DTP_NONE:    do not add date/time (standard).
  ! =OU_DTP_ADDDATE: add time to the output.
  ! =OU_DTP_ADDTIME:    add date to the output.
  ! =OU_DTP_ADDDATETIME:    add both, date and time to the output.
  ! If the parameter is not present, cdefaultDateTimeLogPolicy will be used
  ! as default parameter.
  integer, intent(in), optional :: cdateTimeLogPolicy
!</input>

!</subroutine>

  integer(I32) :: coMode

    select case (ioutputLevel)
    case (:0)
      coMode = 0
    case (1)
      coMode = OU_MODE_LOG
    case (2:)
      coMode = OU_MODE_LOG+OU_MODE_TERM
    end select

    call output_separator (csepType, coutputClass, coMode, ssubroutine, cdateTimeLogPolicy)

  end subroutine

end module

max_stars_repo_licenses: Intel Unlicense

max_issues_repo_licenses: Intel Unlicense

max_forks_repo_licenses: Intel Unlicense


text:
      SUBROUTINE CMP3OUT(JCMN,IM,UTDELTA,VTDELTA,YM,NDT,FRY)

!     Legacy Mars GCM v24
!     Mars Climate Modeling Center
!     NASA Ames Research Center

C
C  PURPOSE:
C      CMP3OUT ACCUMULATES GLOBAL AND ZONAL AVERAGES, CALCULATES A
C      VARIETY OF OUTPUT QUANTITIES, AND PRODUCES SELECTED PRINTER
C      OUTPUT.
C  AUTHOR
C      STEVE POHORSKY    INFORMATICS     TASK 605    OCT 81
C  FOR
C      JIM POLLACK
C  ENVIRONMENT
C      Cray-2            UNICOS 3.0      FORTRAN
C  REVISION HISTORY
C     4-82  SP   MODIFICATIONS MADE AS IN NOTE OF 3-8-82.
C      Re-written September 1993.
C      Removed references to all commons and created an argument list
C      to pass all variables.
C  OUTPUT PARAMETERS
C      FRY(J,L) - The rate at which heat is being gained by the surface
C                 in exchange with the atmosphere through thermal
C                 radiation and convection.
C  CALLED BY
C      COMP3
C
      use grid_h
      use defines_h
      implicit none
C
C######################################################################

      REAL*8 NDT

      REAL*8 FKSX(L_LEVELS), FKSY(L_LEVELS)
      REAL*8 UTDELTA(L_LAYERS), VTDELTA(L_LAYERS), YM(L_LEVELS)
      REAL*8 FRY(L_JSIZE,L_LAYERS)

      integer :: k, l, jcmn, im

C#=====================================================================

      FKSX(3) = 0.0
      FKSY(3) = 0.0
      FKSX(5) = UTDELTA(1)*YM(4)/NDT
      FKSY(5) = VTDELTA(1)*YM(4)/NDT

C     Velocity change and mass

      DO L = 2, L_LAYERS
        K       = 2*L+3
        FKSX(K) = UTDELTA(L)*YM(K-1)/NDT+FKSX(K-2)
        FKSY(K) = VTDELTA(L)*YM(K-1)/NDT+FKSY(K-2)
      END DO

      DO L=1,L_LAYERS
        FRY(JCMN,L) = FRY(JCMN,L)+(FKSY(2*L+3)-FKSY(2*L+1))/FLOAT(IM)
      END DO

      RETURN
      END

ext: f

max_stars_repo_name: nasa/legacy-mars-global-climate-model

max_issues_repo_name: nasa/legacy-mars-global-climate-model

max_forks_repo_name: nasa/legacy-mars-global-climate-model