Treatment of B0 in mapping of equilibrium to IDS is inconsistent.
In my opinion, since no one expressed any objection to my last remark in meq#41, I think the 0.996
correction should be removed completely.
Should you choose to disagree, then the following errors/approximations in matlab/TCV_IMAS/tcv_get_ids_equilibrium.m
need to be addressed:
-
f_df_dpsi
comes only from the plasma so the correction does not apply there. -
f
is not directly proportional tof0
, here the correction would bef_new = sign(f0)*sqrt((0.996*f0)^2 + (f_old)^2 - (f0)^2)
-
Phi
is also not directly proportional tof0
, you would have to recompute it from the newf
. - [ For
f
andPhi
in the limit of small diamagnetism they are indeed proportional tof0
] -
rho_tor
andq
should also have corrections but I do not see them (maybe they are elsewhere).