Skip to content

GPEC v1.2.0

Compare
Choose a tag to compare
@logan-nc logan-nc released this 17 Apr 16:11
· 403 commits to master since this release

Adds

  • RDCON & RMATCH - Adds resistive DCON packages
    • Calculates inner layer model and performs matching with ideal outer layer
    • Cite [A.H. Glasser, Z.R. Wang, and J.-K. Park, Physics of Plasmas 23, 112506 (2016)]
    • Includes new resistive examples
  • DCON - Adds ability to calculate bounce harmonics in parallel when forming kinetic matrices for Euler-Lagrange Equation
  • DCON - Adds ability to start the Euler-Lagrange ODE integration at or above a minimum safety factor qlow
  • DCON - Adds ability to include electron kinetic terms in Euler-Lagrange equation (controlled by the new electron_flag)
  • VACUUM - Adds ability to read prescribed wall geometry file
  • GPEC - Adds new singfld and signcoup calculations and includes all singcoup and singfld outputs in netcdf
    • Delta: is the published resonant coupling metric [Park, Phys. Plasmas 2007] normalized by B.grad(theta) instead of B.grad(phi).
    • This is similar to what is sometimes called external Delta' in tearing stability theory need to divide by vacuum
    • Penetrated resonant flux: interpolated across singspot, and physically meaningful for kinetic MHD equilibria only
  • DCON, GPEC, & PENTRC - Updates the version based on the compile-time git commit
  • DCON, GPEC & PENTRC - Can use a classical spline coefficient solution for "extrap" boundary condition splines, avoiding a suspected (minor) bug in the original tri-diagonal solution that resulted in large grad-shafranov errors in poorer quality equilibrium (especially inverse or modified equilibrium).
    • Previous tri-diagonal spline solutions can be recovered by setting use_classic_splines to false

Changes

  • DCON - Improves clarity of singular surface search messages
  • GPEC - Improves clarity and consistency of singular coupling outputs
    • Uses iszinv to invert hermitian fldflxmat
    • Uses area normalization of penetrated flux for consistency with effective flux
    • Adds unique names for the singcoup mat and svd ascii outputs (enables python reading)
  • PYPEC - Improves automatic selection of partitions and threads in job submission and adds rdcon to exe options

Fixes

  • DCON - Fixes only the the plasma energy matrix written to dcon.out to include full matrix (previously only 2 columns)
  • GPEC - Improves clarity and consistency of singular coupling outputs
    • Corrects units of Phi_res in netcdf (area normalized, so T not Wb)
    • Corrects units and calculation of island width in netcdf (unitless width in psi_n, required a sqrt)
  • GPEC - Fixes bug in iszinv for m/=mpert matrices (no impact on previous results, which all used m=mpert)
  • GPEC - Fixes bug in the normalization of singular coupling islandwidths (singdfld unchanged)
  • GPEC - Fixes poor formatting in response file header
  • PENTRC - Corrects the sign of the charge when calculating NTV torque and kinetic delta-W for electrons
  • VACUUM - Makes vacuum code robust to namelists without a header line

Removes

  • ALL - Removes official support for all compilers other than intel
    • Parallel openmpi calls unique to intel
    • Move is consistent with RDCON development path