KineticForces.calculate_clar (src/KineticForces/Torque.jl) is currently a @warn+return-zero stub. The CLAR (Cylindrical Large Aspect Ratio) method is needed for low-aspect-ratio tokamak NTV calculations.
Required port
- Fortran reference:
~/Code/gpec/pentrc/torque.F90:423-521 (~100 lines: λ-grid setup, kappadjsum action sum, v_∥(θ) bounce-point root-finding, cubic_interp(fbnce(Λ)), lambdaintgrl_lsode, dT/dψ assembly).
- Paper reference: Logan & Park, PoP 20, 122507 (2013), Sec. IV; Park et al. PRL 102 065002 (2009), Eq. 12.
Acceptance
- Reproduces Fortran CLAR runtime-verified local profile from the
kinetic-clean branch at src/Pentrc/tmp/runtime_clar_case/clar_local_profile_compare.svg (within ε ≤ 5% — note that branch also documents a known cumulative-error issue near the outer edge; resolving that is part of this issue).
- The existing
@warn stub is removed and a regression baseline is added.
KineticForces.calculate_clar(src/KineticForces/Torque.jl) is currently a@warn+return-zero stub. The CLAR (Cylindrical Large Aspect Ratio) method is needed for low-aspect-ratio tokamak NTV calculations.Required port
~/Code/gpec/pentrc/torque.F90:423-521(~100 lines: λ-grid setup,kappadjsumaction sum, v_∥(θ) bounce-point root-finding,cubic_interp(fbnce(Λ)),lambdaintgrl_lsode, dT/dψ assembly).Acceptance
kinetic-cleanbranch atsrc/Pentrc/tmp/runtime_clar_case/clar_local_profile_compare.svg(within ε ≤ 5% — note that branch also documents a known cumulative-error issue near the outer edge; resolving that is part of this issue).@warnstub is removed and a regression baseline is added.