This makes the following changes to summary.7.txt: 1. Preserve the sign of J_φ_peak J_φ_max, i.e. if the J_φ is negative, those will be negative. Previously they would always be unsigned. 2. If the sign of J_φ is alternating significantly, specifically if (∫|J_φ|dA - |∫J_φ⋅dA|) / ∫|J_φ|dA > 10% then sign of ρ_avg_J, Δρ_avg_J, ρ_max_J, Δρ_max_J will be negative. Previously they were always positive. 3. If there is no current J_φ=0 (power dPdV=0), then ρ_avg_J=1, rho_max_J=1 (ρ_avg_P=1, rho_max_P=1). Previosuly they would be set to zero. This helps producing continuous results as a beam is gradually focused into the plasma.
24 lines
1.5 KiB
Plaintext
24 lines
1.5 KiB
Plaintext
# Run date/time: 2025/02/21 08:27:23.299
|
||
# GRAY Git revision: e86d293-dirty
|
||
# EQL input: ./F4E_Equil_10002_case001_24VQV5.eqdsk
|
||
# EQL B0 R0 aminor Rax zax: 0 0 0 0 0
|
||
# EQL iequil sgnb sgni factb: 2 -1 -1 0.500
|
||
# EQL icocos ipsinorm idesc ifreefmt: 11 F T F
|
||
# EQL ssplps ssplf ixp: 0.500E-2 0.100E-1 0
|
||
# PRF input: ./profiles.txt
|
||
# PRF iprof iscal factne factte: 0 2 1.000 1.000
|
||
# PRF irho psnbnd sspld: N/A (analytical)
|
||
# PRF Te0 ne0 Zeff0: 0 0 0
|
||
# ANT input: ./beamUSM.data
|
||
# ANT ibeam iox psi chi: 0 2 0.000 0.000
|
||
# ANT alpha beta power fghz: 39.0 20.0 1.00 170.
|
||
# ANT x0 y0 z0: 700. 0.00 441.
|
||
# ANT wx wy Rcix Rciy psiw psir: 5.0470 5.0470 -0.31389E-2 -0.31389E-2 0.0000 0.0000
|
||
# RFL rwall: 0.40000E+01
|
||
# COD igrad idst ipass ipol: 1 1 1 F
|
||
# COD nrayr nrayth nstep rwmax dst: 11 16 8000 1.000 0.1000
|
||
# COD iwarm ilarm imx ieccd: 2 4 -20 3
|
||
# COD ipec nrho istpr istpl: 1 501 5 5
|
||
# I_cd P_abs J_φ_peak dPdV_peak ρ_max_J ρ_avg_J ρ_max_P ρ_avg_P Δρ_avg_J Δρ_avg_P ratio_Ja_max ratio_Jb_max s_max ψ χ index_rt J_φ_max dPdV_max Δρ_J Δρ_P P0 cpl1 cpl2
|
||
-3.90425628E-005 3.84383056E-001 4.66399221E-004 3.00793382E-002 -3.70133716E-001 -3.55972396E-001 3.53038507E-001 3.53578490E-001 -6.45085061E-002 4.10148056E-002 1.02377231E+000 1.00001666E+000 8.00000000E+002 -5.63635279E+001 3.41970443E+001 2 5.70581310E-004 2.98529841E-002 -2.96939387E-002 4.15444773E-002 1.00000000E+000 0.00000000E+000 0.00000000E+000
|