39 parameter(pi = 3.14159)
41 parameter(enhancemin = 9.0 * pi * pi / 16.0)
43 parameter(s0 = 1.6806)
45 parameter(nu = 0.6416)
47 parameter(mu = 0.2982)
49 parameter(condvneg = 1.0e-6)
51 parameter(smallv = 0.01)
53 parameter(nearly1 = 0.99)
82 iclconst = var_parms(4)
89 ionfrac = var_parms(11)
90 derivt6 = var_parms(12)
91 ratemax = var_parms(13)
103 if (vt3 .lt. 0.0)
then
104 conductance = condvneg
109 else if (vt3 .lt. vstar*smallv)
then
111 conductance = enhance * iclconst * sqrt(vt3) / ionfrac
112 if (conductance .lt.
gmin)
then
114 else if (conductance .gt.
gmax)
then
125 gs = (gs - 1.0) / (gs - mu)
126 enhance = (x0 + gs*
gap) / gs
127 enhance = mconst * enhance * enhance
128 conductance = enhance * iclconst * sqrt(vt3) / ionfrac
132 rate = (conductance - gt5) / abs(0.5*(conductance + gt5))
133 if (abs(rate) .gt. ratemax)
then
134 conductance = gt5 * (1.0 + ratemax) * sign(1.0,rate)
136 write (88,*)
'growth-rate: ', rate
140 if (conductance .lt.
gmin)
then
142 else if (conductance .gt.
gmax)
then
152 derivt4 = conductance - gt5
153 if (derivt6 .ne. 0.0)
then
154 if (derivt4*derivt6 .lt. 0.0)
then
159 write (88,*)
'gt3,gt5,dgt4: ',
160 & conductance, gt5, derivt4
165 var_parms(6) = conductance
166 var_parms(12) = derivt4
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 c c c158 MFI Insulator CB c c & gmax
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
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 c c c158 MFI Insulator CB c & gmin
subroutine sabdiode_model