Screamer Project  V3.3.1
Screamer Structure
 All Files Functions Variables
abdioprm.f
Go to the documentation of this file.
1  subroutine abdiodeparm (parms, nparms)
2  real parms(*)
3  integer nparms
4 c December 10, 1992; hnw
5 c 2014-02-06 RBS: removed real*4 replaced with real
6 c
7 c Sets up the parameters needed for the applied-B diode model and
8 c returns them in parms(i).
9 c nparms is the number of parameters in parms.
10 c
11 c parms is sent with the basic parameters needed to rearrange and fill
12 c the actual parms array.
13 c
14  real temp(10)
15 c
16 c Internal parameters
17 c
18  real pi
19  parameter(pi = 3.14159265)
20  real eps0
21  parameter(eps0 = 8.8542e-12)
22  real cspeed
23  parameter(cspeed = 2.9979e+8)
24  real qelectron
25  parameter(qelectron = 1.602e-19)
26  real massproton
27  parameter(massproton = 1.67262e-27)
28  real q2m
29  parameter(q2m = qelectron / massproton)
30  real enhancemin
31  parameter(enhancemin = 9.0 * pi * pi / 16.0)
32  real s0
33  parameter(s0 = 1.6806)
34  real const
35  parameter(const = 4.0 * eps0 / 9.0)
36 c
37 c Fill temp with parms, then rearrange parms and fill it in.
38 c
39  do i = 1, nparms
40  temp(i) = parms(i)
41  end do
42 c
43 c In parms we want:
44 c 1: gap (does not change)
45 c = diode gap
46 c 2: mconst (does not change)
47 c = constant for simplifying calculations
48 c 3: sconst (does not change)
49 c = constant for simplifying calculations
50 c 4: iclconst (does not change)
51 c = constant for simplifying calculations
52 c 5: vt4
53 c = voltage at the time step t=i-2
54 c 6: gt5
55 c = conductance at the time step t=i-3/2
56 c 7: x0 (does not change)
57 c = distance from cathode corner to gas cell foil
58 c 8: gmax (does not change)
59 c = maximum possible conductance of diode
60 c 9: gmin (does not change)
61 c = minimum possible conductance of diode
62 c 10: vstar (does not change)
63 c = voltage at which current diverges theoretically
64 c 11: ionfrac (does not change)
65 c = ion current as a fraction of the total current
66 c 12: derivt6
67 c = the value of gt5-gt7
68 c 13: ratemax (does not change)
69 c = maximum growth rate for the conductance per time step
70 c
71 c In temp:
72 c 1: gap
73 c 2: area
74 c 3: units of electron charge on ion (basically an integer)
75 c 4: amu of ion (basically an integer)
76 c 5. b0
77 c 6. x0
78 c 7. gmax
79 c 8. gmin
80 c 9. ionfrac
81 c 10. ratemax
82 c
83 c total number of parameters for model
84  nparms = 13
85 c
86 c --gap
87  parms(1) = temp(1)
88 c
89 c --mconst
90  parms(2) = enhancemin / ((temp(6)+ temp(1)) **2)
91 c
92 c --sconst
93  parms(3) = cspeed * temp(5) * temp(1)
94 c
95 c --iclconst
96  parms(4) = const * temp(2) * sqrt(2.0*q2m*temp(3)/temp(4))
97  & / (temp(1)**2)
98 c
99 c --vt4
100  parms(5) = 0.0
101 c
102 c --gt5
103  parms(6) = temp(8)
104 c
105 c --x0
106  parms(7) = temp(6)
107 c
108 c --gmax
109  parms(8) = temp(7)
110 c
111 c --gmin
112  parms(9) = temp(8)
113 c
114 c --vstar
115  parms(10) = parms(3) / s0
116 c
117 c --ionfrac
118  parms(11) = temp(9)
119 c
120 c --derivt6
121  parms(12) = 0.0
122 c
123 c --ratemax
124  parms(13) = temp(10)
125 c
126  return
127  end
subroutine abdiodeparm(parms, nparms)
Definition: abdioprm.f:1