19 parameter(alpha1 = 0.6026)
20 parameter(alpha2 = 0.537)
21 parameter(alpha3 = 0.537)
22 parameter(alpha4 = 0.5755)
23 parameter(beta1 = 0.392)
24 parameter(beta2 = 0.297)
25 parameter(beta3 = 0.203)
26 parameter(beta4 = 0.14)
27 parameter(emass = 5.11e+5)
28 parameter(remass = 1.0 / emass)
29 parameter(v2mv = 1.0e-6)
30 parameter(vcapmv1 = 0.3)
31 parameter(vcapmv2 = 1.0)
32 parameter(vcapmv3 = 3.0)
34 parameter(rdep2 = 1.0 / dep2)
35 parameter(qqq = 0.96432601)
36 parameter(rqqq = 1.0 / qqq)
37 parameter(shuntgmin = 1.0e-6)
38 parameter(eturnon = 2.0e7)
39 parameter(esat = 4.0e7)
41 ibranch = indexmitl(1,index)
42 iblock = indexmitl(2,index)
43 node1 = indexmitl(3,index)
44 node2 = indexmitl(4,index)
45 delperv = pin(1,iblock,ibranch)
46 delt = pin(2,iblock,ibranch)
47 zline = pin(3,iblock,ibranch)
48 gap = pin(5,iblock,ibranch)
55 vcap = abs(v(i,ibranch)) + 1.0
57 iave = abs(zir(i-1,ibranch))
60 if (vcapmv .lt. vcapmv1)
then
61 res = alpha1 * vcapmv**beta1
62 else if (vcapmv .lt. vcapmv2)
then
63 res = alpha2 * vcapmv**beta2
64 else if (vcapmv .lt. vcapmv3)
then
65 res = alpha3 * vcapmv**beta3
67 res = alpha4 * vcapmv**beta4
72 epi = zline * iave / vcap
74 if (epi2 .lt. cond2)
then
77 gc = delperv * sqrt(vcap) * hrelfu
78 epi2res2 = epi2 * res2
79 shuntg = gc * (1.0 - 0.15*epi2res2)
80 shuntg = shuntg * (1.0 - exp((epi2res2 - 1.0) * rdep2))
81 shuntg = shuntg * rqqq
82 shuntg = amax1(shuntg, shuntgmin)
83 if (efield .lt. eturnon)
then
85 else if (efield .lt. esat)
then
86 shuntg = shuntg*(efield-eturnon)/(esat-eturnon)
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_pmitl(index)