isttok_mag['RPfcPrim'], isttok_mag['ZPfcPrim'], isttok_mag['TurnsPfcPrim']
]).T
radiusHCoil0 = 0.13893
angleHCoil0 = 0.5280 # np.arctan(7./12)
radiusHCoil = 0.136
angleHCoil = -0.8
RHCoil2 = isttok_mag['RM'] + radiusHCoil * np.cos(angleHCoil)
ZHCoil2 = radiusHCoil * np.sin(angleHCoil)
print("Rcoi %g, Zcoil %g" % (RHCoil2, ZHCoil2))
Hcl2 = np.array([[0.58, 0.07, 0.0], [RHCoil2, ZHCoil2,
-4.0]]) # Ideal Coil Position
BpolEst2 = buildIs2Bpol(Vcl1, Hcl2, Pcl)
Hcl11 = np.array([[0.58, 0.07, 4.0], [0.58, -0.07,
0.0]]) # Ideal Coil Position
Hcl11 = np.array([[0.58, 0.07, 4.0], [0.58, -0.07,
0.0]]) # Ideal Coil Position
radiusHCoil = 0.13
angleHCoil = 1.035
RHCoil11 = isttok_mag['RM'] + radiusHCoil * np.cos(angleHCoil)
ZHCoil11 = radiusHCoil * np.sin(angleHCoil)
Hcl11 = np.array([[RHCoil11, ZHCoil11, 4.0], [0.58, -0.07,
0.0]]) # Ideal Coil Position
#Estimation of probe signals
# Alternative Coil Position
Vcl = np.array(
[isttok_mag['RPfcVer'], isttok_mag['ZPfcVer'],
isttok_mag['TurnsPfcVer']]).T
Hcl = np.array(
[isttok_mag['RPfcHor'], isttok_mag['ZPfcHor'],
isttok_mag['TurnsPfcHor']]).T
Pcl = np.array([
isttok_mag['RPfcPrim'], isttok_mag['ZPfcPrim'], isttok_mag['TurnsPfcPrim']
]).T
#print("Rcoi-11 %g, Zcoil %g" %(RHCoil11, ZHCoil11))
BpolEst = buildIs2Bpol(Vcl, Hcl, Pcl)
Pcl1 = np.array([[0.615, 0.14, 14.0], [0.62, -0.13,
0.0]]) # Ideal Coil Position
Pcl2 = np.array([
[0.62, 0.13, 0.0], # Ideal Coil Position
[0.615, -0.14, 14.0]
])
BpolEst1 = buildIs2Bpol(Vcl, Hcl, Pcl1)
BpolEst2 = buildIs2Bpol(Vcl, Hcl, Pcl2)
prbNums = np.arange(1, 13)
Vcl1 = np.array([[0.60, 0.08, -1.], [0.59, 0.07, -1.], [0.58, 0.07, -1.],
[0.59, 0.07, -1.], [0.60, 0.06, -1.], [0.58, -0.07, -5.],
[0.35, 0.07, 5.], [0.35, -0.07, 5.]])
VclOptim = np.array([[0.526368, 0.111782, -3.], [0.554752, -0.0975604, -3.],
[0.380273, 0.113862, 3.], [0.3905, -0.120378, 3.]])
Hcl = np.array(
[isttok_mag['RPfcHor'], isttok_mag['ZPfcHor'],
isttok_mag['TurnsPfcHor']]).T
Pcl = np.array([
isttok_mag['RPfcPrim'], isttok_mag['ZPfcPrim'], isttok_mag['TurnsPfcPrim']
]).T
BpolEstOptim = buildIs2Bpol(VclOptim, Hcl, Pcl)
# Nominal coil Position
radiusHCoil0 = 0.13893
angleHCoil0 = 0.5280 # np.arctan(7./12)
radiusHCoil = 0.136
angleHCoil = -0.8
RHCoil2 = isttok_mag['RM'] + radiusHCoil * np.cos(angleHCoil)
ZHCoil2 = radiusHCoil * np.sin(angleHCoil)
print("Rcoi %g, Zcoil %g" % (RHCoil2, ZHCoil2))
Vcl2 = np.array([
[0.58, 0.07, 0.0], # Nominal Coil Position
])
Hcl3 = np.array([[0.58, 0.07, 0.0], [0.58, -0.07,
-1.0]]) # Ideal Coil Position
Vcl = np.array(
[isttok_mag['RPfcVer'], isttok_mag['ZPfcVer'],
isttok_mag['TurnsPfcVer']]).T
Hcl = np.array(
[isttok_mag['RPfcHor'], isttok_mag['ZPfcHor'],
isttok_mag['TurnsPfcHor']]).T
Pcl = np.array([
isttok_mag['RPfcPrim'], isttok_mag['ZPfcPrim'], isttok_mag['TurnsPfcPrim']
]).T
BpolEst = buildIs2Bpol(Vcl, Hcl3, Pcl)
prbNums = np.arange(1, 13)
client = StartSdas()
#Nshot=44330 # Horizontal Field
#PfcHCurrent = -175.0 # in shot=44330
Nshot = 43066 # Primary Field Current 157
PfcHCurrent = -157.0 # in shot=43066
times, dataH = getMirnovInt(client, Nshot, correctWO='Post')
dataHarr = np.array(dataH)
scale = 1e6 # values in uV.s