Examples - Regenerate the DIII-D-like g-file without boundary truncation

未关闭
#528 0 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看

维护者通常 3 天内回复

还没有人认领这个 Issue。

评估

难度
4/5
预计耗时
3-5 天
新手友好度
25/100
Issue 类型
缺陷
描述清晰度
基本清楚
活跃度
活跃
技术栈
julia
领域
data

调研方向

Start with the header comments of examples/DIIID-like_ideal_example/TkMkr_D3Dlike_Hmode.geqdsk to find the TokaMaker export settings, then check the q column and the psi_N values near the boundary. The kinetic files TkMkr_D3Dlike_Hmode_kinetic.h5 and .gpeckf must be checked for matching psi_N. Done when the g-file is re-exported with truncate_eq=False, the kinetic files agree, and the regression run in test/runtests_coils.jl and test/runtests_rerun_from_h5.jl is explained in the PR.

由索引模型根据 Issue 内容生成。

描述

examples/DIIID-like_ideal_example/TkMkr_D3Dlike_Hmode.geqdsk was written by TokaMaker's save_eqdsk with lcfs_pad = 1e-4 and the default truncate_eq=True. The equilibrium is cut at 1 − 1e-4 and renormalised, so the file's X-point sits at ψ_N = 1.0001 and every rational above the boundary q (8.79) lies outside the file: q = 9, 10, 11 and 12 at ψ_N = 1.0000274, 1.0000848, 1.0000969 and 1.0000994. With truncate_eq=False the flux map, F and p reach the separatrix, and only the last entry of the file's q column is a stand-in, which GPEC does not read.

Proposal: re-export from the same TokaMaker equilibrium with truncate_eq=False, and check that TkMkr_D3Dlike_Hmode_kinetic.h5 and .gpeckf are on the matching ψ_N.

This needs its own PR. ψ_N is renormalised by 1e-4 everywhere, so every deck that reads this file moves: the ideal, multi-n, Riccati, SLAYER, error-field, galerkin-resistive and truncation examples, test/runtests_coils.jl and test/runtests_rerun_from_h5.jl, and benchmark_inputs.jld2. Each deck's psihigh will then refer to a slightly different surface, and the truncation deck's cut should be re-checked.

Done when:

  • The g-file is re-exported (generating script or settings recorded in the example's header comments).
  • The kinetic files are verified or regenerated on the same ψ_N.
  • A full regression-harness run, with the moved quantities explained in the PR.

🤖 Generated with Claude Code

主要语言
Julia
星标
11
派生
7
平均合并
4 天 3 小时
30 天内合并 PR
70

环境准备

  • 没有 Dockerfile 或 Docker Compose 文件
  • 有 Pull Request 模板
  • 没有贡献指南

从这里开始

  1. 先读完整个 Issue,再读项目的贡献指南。
  2. 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
  3. Fork 仓库,在一个分支上完成修改。
  4. 提交 Pull Request,并在描述里引用这个 Issue 编号。

OpenFUSIONToolkit/GPEC 的其他 Issue

查看 OpenFUSIONToolkit/GPEC 的全部 Issue

相似的 Issue

更多 Julia Issue

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。