Screamer Project  V3.3.1
Screamer Structure
 All Files Functions Variables
relcrrct.f
Go to the documentation of this file.
1  function hrelf (u)
2 c
3 c Define passed variables
4 c
5  real u
6 c
7 c Change log
8 c 2014-02-06 RBS: Changed real*4 to real
9 c
10 c Function from LXS and DHM SCEPTRE model of perveance specified MITLs.
11 c Relativistic correction to electron motion.
12 c
13  parameter(a0 = 9.99992673e-1)
14  parameter(a1 = -1.07000372e-1)
15  parameter(a2 = 2.35920406e-2)
16  parameter(a3 = -5.98692605e-3)
17  parameter(a4 = 1.33650594e-3)
18  parameter(a5 = -2.1510132e-4)
19  parameter(a6 = 2.08370041e-5)
20  parameter(a7 = -8.91792265e-7)
21  parameter(b0 = 1.0)
22  parameter(b1 = -1.0 / 3.0)
23  parameter(b2 = 1.0 / 3.0)
24  parameter(b3 = -4.28571429e-1)
25  parameter(b4 = 6.19047619e-1)
26  parameter(b5 = -9.56709957e-1)
27  parameter(b6 = 1.54545455)
28  parameter(b7 = -2.57575758)
29  parameter(b8 = 4.39393939)
30  parameter(c1 = -0.847213)
31  parameter(c2 = 3.181980515)
32  parameter(uswitch = 5.0)
33 c
34  if (u .le. uswitch) then
35  hrelf = a0 + u *
36  & (a1 + u *
37  & (a2 + u *
38  & (a3 + u *
39  & (a4 + u *
40  & (a5 + u *
41  & (a6 + u * a7))))))
42 c
43  else
44  a = 1.0 / u
45  s = u + b0 + a *
46  & (b1 + a *
47  & (b2 + a *
48  & (b3 + a *
49  & (b4 + a *
50  & (b5 + a *
51  & (b6 + a *
52  & (b7 + a * b8)))))))
53  scon = sqrt(u*u + u + u)
54  s = s / sqrt(scon) + c1
55  hrelf = c2 * s * s / (u * sqrt(u))
56 c
57  end if
58 c
59  return
60  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 c &trapped field c153 c c c Sin c c c remaining calls in rdscrelem c Tabular a10 Cond c c c &used with this model c Exponential model of a resistive c c c Rise model of a resistive c c Z FLOW POS c c c c SW1 a3
Definition: zdemfmt.h:140
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 a6
Definition: zdemfmt.h:36
function hrelf(u)
Definition: relcrrct.f:1