28 parameter(pi = 3.1415927)
29 parameter(zmuo = 4.0e-7 * pi)
30 parameter(epso = 8.8541878e-12)
31 parameter(cspeed = 2.99792458e+8)
33 parameter(esat = 6.0e+7)
34 parameter(con1 = 0.0625 * zmuo * zmuo)
36 parameter(xy = 0.9 * bm)
37 dimension gg(max_nodes)
39 ibranch = indexmitl(1,index)
40 iblock = indexmitl(2,index)
41 node1 = indexmitl(3,index)
42 node2 = indexmitl(4,index)
43 circum = pin(1,iblock,ibranch)
44 gap = pin(2,iblock,ibranch)
45 tau = pin(3,iblock,ibranch)
46 rcircum2 = pin(6,iblock,ibranch)
47 rgap = pin(7,iblock,ibranch)
48 rgap2 = pin(8,iblock,ibranch)
49 eturnon = pin(9,iblock,ibranch)
51 elec_len = tau * cspeed
52 elec_area = elec_len * circum
53 if (eturnon .lt. esat)
then
54 edenom = 1.0 / (esat - eturnon)
61 do node = node1, node2
63 v1 = 0.5 * abs(v(node,ibranch)+vold(node,ibranch))
67 if (e1 .ge. eturnon)
then
68 if (e1 .gt. esat)
then
71 emp1 = (e1-eturnon) * edenom
73 v1mv = 0.00001 + 1.0e-6*v1
74 if (v1mv .gt. 0.09855)
then
75 a1 = 1.38 * (v1mv+0.0001)**(-0.036)
79 zir1 = zir(node-1,ibranch) + zir(node,ibranch)
80 zir2 = zirn(node-1,ibranch) + zirn(node,ibranch)
83 con2 = 1.113e-5 * rgap2
84 b1s = con1 * zir1 * zir1 * rcircum2 + 1.0e-6
85 bc1s = con2 * (1.022*v1mv + v1mv*v1mv + 0.0001)
94 arg1 = 1.21 * b1s * rbc1s
99 if (g1x .gt. 40.0)
then
101 else if (g1x .lt. -40.0)
then
106 xx = bm * (b1s * rbc1s)
108 if ((xx-xy).gt. 4.0)
then
111 ftail = exp(-10.0*(xx-xy))
116 fung1 = fung1 * ftail
117 g1y = 1.99e-3 * (v1mv**(a1-1.0)) * argp * fung1
118 g(node,ibranch) = g1y * emp1
120 g(node,ibranch) = 0.0
129 cup = zir(node1-1,ibranch)
130 cdn = zir(node2,ibranch)
132 alossden = anodeloss/(area*float(node2-node1+1))
133 pin(10,iblock,ibranch) = alossden
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line corrected spelling errors in format added statement for Zflow Plasma Loss Model added format for CSV output type fixed more lines longer than characters added format for Measure Zflow Block and forward c reverse current directions in Zflow plasma loss c and Zflow POS models added format for SFC output type c removed from all code calls c c c c c a80 c i10 c No grids on plots c Do not write files containing the plotted points c Execute only one cycle c Do not echo the setup parameters and indicies c c &exitting c a13 c c c c c c c c102 c &described as a function of time c c c &described as a function of time c c c &function of time c c shell c &min A K gap
subroutine reset_mitl(index)