Write estimated frequency and growth rates to out_unit for an individual point
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(in) | :: | time |
Simulation time |
||
integer, | intent(in) | :: | it |
- and -indices |
||
integer, | intent(in) | :: | ik |
- and -indices |
||
real, | intent(in) | :: | phitot |
Simulation time |
subroutine do_write_line(time,it,ik,phitot)
use kt_grids, only: aky, akx, theta0
use gs2_time, only: woutunits
implicit none
!> Simulation time
real, intent(in) :: time, phitot
!> \(k_y\)- and \(\theta\)-indices
integer, intent(in) :: ik, it
if (.not. write_ascii) return
write (out_unit, "('t= ',e17.10,' aky= ',1p,e12.4, ' akx= ',1p,e12.4, &
&' om= ',1p,2e12.4,' omav= ', 1p,2e12.4,' phtot= ',1p,e12.4,' theta0= ',1p,e12.4)") &
time, aky(ik), akx(it), &
real( omega(it,ik)*woutunits(ik)), &
aimag(omega(it,ik)*woutunits(ik)), &
real( omegaavg(it,ik)*woutunits(ik)), &
aimag(omegaavg(it,ik)*woutunits(ik)), &
phitot,theta0(it,ik)
end subroutine do_write_line