Build EAST physics methods

[yumouwei]

Status

efit.py

Function	Build	Validate
get_efit_parameters	✅
get_pefit_parameters	✅

physics.py

Function	Build	Validate	Note
get_time_until_disrupt	✅
get_ip_parameters	✅
get_v_loop	✅
get_z_error	✅
get_density_parameters	✅
get_power	✅
get_p_ohm	✅
(prad_bulk_xuv2014_2016)	⚠️		Script missing in MATLAB branch
get_n_equal_1_data			Skip, implement at the end
(get_rmp_and_saddle_signals)
(calculate_rmp_saddle_coeff_matrix.m)	⚠️		Both this script and the precomputed mp_saddle_coeff_matrix.mat are missing!
get_efit_gaps	✅		Extremely messy function! Needs to compare line by line with MATLAB!
get_pefit_gaps			Skip, copy of get_efit_gaps
get_kappa_area	✅
get_pkappa_area	✅
get_pcs	➖	➖	Broken into get_pcs_parameters and get_pcs_powers.
get_pcs_parameters	✅
get_pcs_powers	✅		p_ecrh_rt, p_oh_rt, p_icrf_rt, rad_input_frac, & rad_loss_frac weren't implemented in the source MATLAB script.
get_prad_peaking	✅
get_mirnov_std	✅
get_n1rms_n2rms	✅
get_h98	✅

• (...): auxiliary functions

TODO prior to the validation stage

• [ ] get_n_equal_1_data, get_mirnov_std, and get_n1rms_n2rms all use a similar method to compute btor. Consider making it into a standalone physics method.

TODO at the validation stage

• [ ] Verify MDS nodes
• [ ] Verify units
• [ ] Check similar EFIT addresses (e.g. \aout vs \aminor)
• [ ] Check, add, or remove try-except logic.
• [ ] interp1 extrapolation (e.g. in get_power)
• [ ] Restructure repetitive code blocks to make things more efficient (e.g. get icrf power).
• [ ] Divide by 0 error (e.g. greenwald_fraction, rad_input_frac)
• [ ] get_p_ohm: replace 1.85 with r0 signal
• [ ] n1rms, n2rms: database might have Mirnov_std stored in that column.
• [ ] n1rms_normalized, n2rms_normalized: database might have Mirnov_std_normalized stored in that column
• [ ] h98: database might have prad_peaking stored in that column

TODO after the validation stage

• [ ] Modify get_heating_power() to make it compatible with both get_power() and get_pcs_power.