Right now GPEC calculates 5 resonant coupling matrices:
│ ├── C_resonant_flux # [n_rational × numpert_total] coupling matrix
│ ├── C_resonant_current
│ ├── C_island_width_sq
│ ├── C_penetrated_field
│ ├── C_delta_prime
As per #233 all should use coordinate independent conventions. On the RHS of A = C . x it is clear x needs to be the half-area weighted field. All should use the definition of half-area decomposed field divided by sqrt(A_control_surface) so the units are Tesla. On the LHS of the standard coupling matrix, [Park 2008] made it clear the right quantity to use is the flux (not the field), but there is an Area normalization here to make that have units of Tesla as well. It should be clearly documented and annotated which Area is used (local singular surface areas for each mode or a common control surface area for the full vector).
For the others, each should also be explicitly confirmed to be coordinate independent. They should be (as long as psi is always normalized poloidal flux), but I'd like to see plotted proof. In particular, I would like someone to double check if C_penetrated_field has field or flux on the LHS.
Right now GPEC calculates 5 resonant coupling matrices:
As per #233 all should use coordinate independent conventions. On the RHS of
A = C . xit is clear x needs to be the half-area weighted field. All should use the definition of half-area decomposed field divided by sqrt(A_control_surface) so the units are Tesla. On the LHS of the standard coupling matrix, [Park 2008] made it clear the right quantity to use is the flux (not the field), but there is an Area normalization here to make that have units of Tesla as well. It should be clearly documented and annotated which Area is used (local singular surface areas for each mode or a common control surface area for the full vector).For the others, each should also be explicitly confirmed to be coordinate independent. They should be (as long as psi is always normalized poloidal flux), but I'd like to see plotted proof. In particular, I would like someone to double check if
C_penetrated_fieldhas field or flux on the LHS.