bmin and bmax are taken over ntheta2 poloidal points for every run:
! eqfor.f90:401
CALL bextrema (modb, bmin(1,js), bmax(1,js), nzeta, ntheta2)For lasym = T the field is stored on ntheta3 = ntheta1 points and the second poloidal half is not the mirror of the first, so the extrema of |B| along the points in [pi, 2 pi) are never examined, and the threed1 line printed at eqfor.f90:404 reports half-surface values.
Measured on input.cth_like_fixed_bdy_asym (VMEC++ test data) from the boundary-surface bmnc, bmns: max |B| over the stored half is 0.73987, over the full surface 0.78388, and the per-theta maxima on the second half differ from their mirrors by up to 10%. The 2017 VMEC2000 8.52 source and PARVMEC (eqfor.f:375) carry the same call.
bminandbmaxare taken overntheta2poloidal points for every run:For
lasym = Tthe field is stored onntheta3 = ntheta1points and the second poloidal half is not the mirror of the first, so the extrema of |B| along the points in[pi, 2 pi)are never examined, and the threed1 line printed at eqfor.f90:404 reports half-surface values.Measured on
input.cth_like_fixed_bdy_asym(VMEC++ test data) from the boundary-surfacebmnc,bmns: max |B| over the stored half is 0.73987, over the full surface 0.78388, and the per-theta maxima on the second half differ from their mirrors by up to 10%. The 2017 VMEC2000 8.52 source and PARVMEC (eqfor.f:375) carry the same call.