do_write_crossphase Subroutine

private subroutine do_write_crossphase(time)

Uses

Calculate the cross-phase (see get_cross_phase) and write to the phase_unit file

Arguments

Type IntentOptional Attributes Name
real, intent(in) :: time

Current simulation time

Contents

Source Code

do_write_crossphase

Source Code

  subroutine do_write_crossphase(time)
    use mp, only: proc0
    implicit none
    !> Current simulation time
    real, intent(in) :: time
    real :: phase_tot, phase_theta
    if (.not. write_ascii) return
    call get_cross_phase (phase_tot, phase_theta)
    if (proc0) write (unit=phase_unit, fmt="('t= ',e17.10,' phase_tot= ',e11.4,' phase_theta= ',e11.4)") &
         & time, phase_tot, phase_theta
  end subroutine do_write_crossphase