Half to full mesh interpolation and calculation of df/ds
Closed this issue · 1 comments
lazersos commented
When dealing with VMEC quantities we often need values both at arbitrary points in space, radial derivatives and moved from the half to the full grid. Care must be taken as the even m-modes can be linearly interpolated on the flux grid but the odd m-modes need to be linearly interpolated on the rho (sqrt(s)) grid. This issue is to help track the proper implementation of interpolation in various codes:
- BEAMS3D (fixed in Lazerson branch)
- BNORM
- BOOTSJ (BOOZ_XFORM)
- BOOZ_XFORM
- COBRAVMEC
- COILOPT (fixed in Lazerson branch)
- DIAGNO
- DKES (BOOZ_XFORM)
- FIELDLINES (fixed in Lazerson branch)
- J_INVARIANT (BOOZ_XFORM)
- LIBSTELL
- MAKEGRID
- NEO (BOOZ_XFORM)
- NESCOIL
- PENTA (fixed in Lazerson branch)
- STELLOPTV2 (fixed but may need to revisit ns=1 quantities in the future)
- TORLINES
- VMEC2000
- VMEC2PIES
- VMEC2SPEC
- VMEC2STEL
- VMEC2V690
- VMEC2XGC
- WALL_ACCELERATE
lazersos commented
This has now been checked in all the selected codes. Recently the fix for COBRAVMEC was implemented and merged into the Lazerson branch. A pull request of Lazerson into develop will come soon.