output Subroutine

  1. module_self_force.f90
  2. self_force_base
  3. output

public subroutine output(this)

Uses

  • proc~~output~~UsesGraph proc~output output module~time_info time_info proc~output->module~time_info module~output_base output_base proc~output->module~output_base module~kinds kinds module~time_info->module~kinds module~output_base->module~kinds

Routine that performs output of the self-force.

Arguments

Type IntentOptional AttributesName
class(self_force), intent(inout) :: this

The routine is called on this self_force object.

Calls

proc~~output~~CallsGraph proc~output output proc~next_available_io_id next_available_io_id proc~output->proc~next_available_io_id proc~get_current_time get_current_time proc~output->proc~get_current_time

Contents

Source Code

output

Source Code

  subroutine output ( this )
  !! Routine that performs output of the self-force.

    use output_base
    use time_info, only : get_current_time

    implicit none

    class(self_force), intent(inout) :: this
    !! The routine is called on this [[self_force]] object.
    integer(ip) :: ioo_id, tmp_id
    character(len=14) :: filename = 'self_force.asc'

    ioo_id = this%ioo_id
    if (ioo_id<0) then
      ioo_id = next_available_io_id ()
      this%ioo_id = ioo_id
      print*,'Opening ', filename, ' with id ', ioo_id
      open(ioo_id, file=filename, status='replace', action='write')
    end if

    write(ioo_id, '(*(es23.15e3,1x))') get_current_time ( ), this%f
  end subroutine output