FIXME : Add documentation
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | unit |
subroutine wnml_dist_fn(unit)
use species, only: spec, has_electron_species
implicit none
integer, intent(in) :: unit
if (dfexist) then
write (unit, *)
write (unit, fmt="(' &',a)") "dist_fn_knobs"
select case (boundary_option_switch)
case (boundary_option_default)
write (unit, fmt="(' boundary_option = ',a)") '"default"'
case (boundary_option_zero)
write (unit, fmt="(' boundary_option = ',a)") '"zero"'
case (boundary_option_self_periodic)
write (unit, fmt="(' boundary_option = ',a)") '"periodic"'
case (boundary_option_linked)
write (unit, fmt="(' boundary_option = ',a)") '"linked"'
case (boundary_option_alternate_zero)
write (unit, fmt="(' boundary_option = ',a)") '"alternate-zero"'
end select
write (unit, fmt="(' nonad_zero = ',L1)") nonad_zero
write (unit, fmt="(' gridfac = ',e17.10)") gridfac
write (unit, fmt="(' esv = ',L1)") esv
write (unit, fmt="(' opt_source = ',L1)") opt_source
if (.not. has_electron_species(spec)) then
select case (adiabatic_option_switch)
case (adiabatic_option_default)
write (unit, *)
write (unit, fmt="(' adiabatic_option = ',a)") &
& '"no-field-line-average-term"'
case (adiabatic_option_fieldlineavg)
write (unit, fmt="(' adiabatic_option = ',a)") '"field-line-average-term"'
case (adiabatic_option_yavg)
write (unit, fmt="(' adiabatic_option = ',a)") '"iphi00=3"'
end select
end if
if (apfac /= 1.) write (unit, fmt="(' apfac = ',e17.10)") apfac
if (driftknob /= 1.) write (unit, fmt="(' driftknob = ',e17.10)") driftknob
if (tpdriftknob /= 1.) write (unit, fmt="(' tpdriftknob = ',e17.10)") tpdriftknob
if (poisfac /= 0.) write (unit, fmt="(' poisfac = ',e17.10)") poisfac
if (afilter /= 0.) write (unit, fmt="(' afilter = ',e17.10)") afilter
if (mult_imp) write (unit, fmt="(' mult_imp = ',L1)") mult_imp
if (def_parity) then
write (unit, fmt="(' def_parity = ',L1)") def_parity
if (even) write (unit, fmt="(' even = ',L1)") even
end if
write (unit, fmt="(' /')")
endif
if (skexist) then
write (unit, *)
write (unit, fmt="(' &',a)") "source_knobs"
select case (source_option_switch)
case (source_option_full)
write (unit, fmt="(' source_option = ',a)") '"full"'
case(source_option_phiext_full)
write (unit, fmt="(' source_option = ',a)") '"phiext_full"'
write (unit, fmt="(' source0 = ',e17.10)") source0
write (unit, fmt="(' omega0 = ',e17.10)") omega0
write (unit, fmt="(' gamma0 = ',e17.10)") gamma0
write (unit, fmt="(' t0 = ',e17.10)") t0
write (unit, fmt="(' phi_ext = ',e17.10)") phi_ext
case(source_option_homogeneous)
write (unit, fmt="(' source_option = ',a)") '"homogeneous"'
end select
write (unit, fmt="(' /')")
endif
end subroutine wnml_dist_fn