Screamer Project  V3.3.1
Screamer Structure
 All Files Functions Variables
rstpmitl.f
Go to the documentation of this file.
1  subroutine reset_pmitl (index)
2 c
3 c Define passed variables
4 c
5  integer index
6 c
7 c Change log
8 c
9 c 2014-02-06 RBS: Changed real*4 to real
10 c
11 c Resets the conductances in an MITL. Uses the perveance value.
12 c Per LXS and DHM SCEPTRE model.
13 c
14  include 'zdemmax.h'
15  include 'zdemcomm.h'
16  include 'zdemwork.h'
17 c
18  real iave
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)
33  parameter(dep2 = 0.3)
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)
40 c
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)
49 c
50 c Reset the conductances.
51 c
52 c Find the voltage across the capacitor and the average current in the line.
53 c
54  do i = node1, node2
55  vcap = abs(v(i,ibranch)) + 1.0
56  efield = vcap/gap
57  iave = abs(zir(i-1,ibranch))
58  vcapmv = vcap * v2mv
59 c
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
66  else
67  res = alpha4 * vcapmv**beta4
68  end if
69  res2 = res * res
70  cond2 = 1.0 / res2
71 c
72  epi = zline * iave / vcap
73  epi2 = epi * epi
74  if (epi2 .lt. cond2) then
75  u = vcap * remass
76  hrelfu = hrelf(u)
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
84  shuntg = shuntgmin
85  else if (efield .lt. esat) then
86  shuntg = shuntg*(efield-eturnon)/(esat-eturnon)
87  endif
88  else
89  shuntg = shuntgmin
90  end if
91 c
92  g(i,ibranch) = shuntg
93 c
94  end do
95 c
96  return
97  end
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
Definition: zdemfmt.h:84
function hrelf(u)
Definition: relcrrct.f:1
subroutine reset_pmitl(index)
Definition: rstpmitl.f:1