3 & (
time, dt, rdt, ilt2,
9 character(len=80)::buffer
83 real implode, end_implode
84 parameter(implode = 0.0, end_implode = 1.0)
145 if (abs(ilt2) .gt. peakcur)
then
151 if (testimpl .eq. implode)
then
155 acct2 = (ilt2 * ilt2 * aconst) / (radt3)
159 velt1 = velt3 + acct2*dt
166 radt1 = radt3 + ((velt1+velt3) * dt * 0.5)
167 if (radt1 .lt. minrad)
then
169 testimpl = end_implode
172 foilke = 0.5 * mass * velt1 * velt1
177 &
'(/A/,A/,A,1pe10.3/,A,1pe10.3/,A,1pe10.3/,A,1pe10.3/,A,
179 &
' ------------------------------------------------------------',
180 &
' Foil radius has reached the minimum value.',
182 &
' minimum radius: ',minrad,
183 &
' maximum KE: ',foilke,
184 &
' peak velocity: ',velt1,
185 &
' peak current: ',peakcur,
186 &
' ------------------------------------------------------------'
189 &
' ------------------------------------------------'
190 call writebuffer(buffer)
192 &
' Foil radius has reached the minimum value.'
193 call writebuffer(buffer)
194 write(buffer,
'(A,1pe10.3)')
' time: ',
time
195 call writebuffer(buffer)
196 write(buffer,
'(A,1pe10.3)')
' minimum radius: ',minrad
197 call writebuffer(buffer)
198 write(buffer,
'(A,1pe10.3)')
' maximum KE: ',foilke
199 call writebuffer(buffer)
200 write(buffer,
'(A,1pe10.3)')
' peak velocity: ',velt1
201 call writebuffer(buffer)
202 write(buffer,
'(A,1pe10.3)')
' peak current: ',peakcur
203 call writebuffer(buffer)
205 &
' ------------------------------------------------'
206 call writebuffer(buffer)
211 lt1 = lconst * log(initrad/radt1)
215 ldott2 = (lt1 - lt3) * rdt
220 foilke = 0.5 * mass * velt1 * velt1
236 call
yield(13.,mass,radt1,
length,foilke,eta,t,yw_al,ym_al)
237 call
yield(18.,mass,radt1,
length,foilke,eta,t,yw_ar,ym_ar)
238 call
yield(29.,mass,radt1,
length,foilke,eta,t,yw_cu,ym_cu)
239 call
yield(36.,mass,radt1,
length,foilke,eta,t,yw_kr,ym_kr)
240 call
yield(54.,mass,radt1,
length,foilke,eta,t,yw_xe,ym_xe)
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 time
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 length
subroutine yield(Z, massm, rfm, Lzm, Kin, eta, T, Yw, Ym)