Skip to content

Fix EPSOFF logic for post-198183 shots#547

Open
yumouwei wants to merge 1 commit intodevfrom
wei/d3d/fix-epsoff-logic
Open

Fix EPSOFF logic for post-198183 shots#547
yumouwei wants to merge 1 commit intodevfrom
wei/d3d/fix-epsoff-logic

Conversation

@yumouwei
Copy link
Copy Markdown
Contributor

@yumouwei yumouwei commented Apr 9, 2026

Implemented changes in get_ip_parameters and get_rt_ip_parameters as described in Issue 506

@yumouwei yumouwei added the machine: DIII-D Related to the DIII-D tokamak label Apr 9, 2026
@yumouwei yumouwei changed the title Fix EPSOFF logic for shot_id >= 198183 Fix EPSOFF logic for post-198183 shots Apr 9, 2026
@gtrevisan gtrevisan linked an issue Apr 9, 2026 that may be closed by this pull request
Copy link
Copy Markdown
Member

@gtrevisan gtrevisan left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

what's the final effect of this change?
are you able to share a couple of plots for some shot?

also, have you tried executing full tests to see whether we need to mark anything as an expected failure?

-------
- original source: [get_Ip_parameters.m](https://github.com/MIT-PSFC/disruption-py/blob
/matlab/DIII-D/get_Ip_parameters.m)
- issues: #[506](https://github.com/MIT-PSFC/disruption-py/issues/506)
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

you might want to add this PR, as well.

t_epsoff += 0.001
epsoff = interp1(t_epsoff, epsoff, params.times, "linear")
railed_indices = np.where(np.abs(epsoff) > 0.5)
if params.shot_id >= 198183:
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

since we are using this in two places, let's choose a good name and make it a configuration parameter.

- original source: [get_Ip_parameters_RT.m](https://github.com/MIT-PSFC/disruption-py
/blob/matlab/DIII-D/get_Ip_parameters_RT.m)
- pull requests: #[254](https://github.com/MIT-PSFC/disruption-py/pull/254)
- issues: #[506](https://github.com/MIT-PSFC/disruption-py/issues/506)
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

you might want to add this PR, as well.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

machine: DIII-D Related to the DIII-D tokamak

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Fix EPSOFF check logic in D3D physics method

2 participants