3 & (
time, dt, rdt, ilt2,
9 character(len=80)::buffer
85 real implode, end_implode
86 parameter(implode = 0.0, end_implode = 1.0)
160 if (testimpl .eq. implode)
then
164 if (radt3 .gt. innerrad)
then
165 accilt2 = (ilt2 * ilt2 * aconst) / (radt3 * masst3)
166 dragt2 = (dconst * radt3 * velt3 * velt3 ) / (masst3)
167 acct2 = accilt2 + dragt2
169 accilt2 = (ilt2 * ilt2 * aconst) / (radt3 * masst3)
175 velt1 = velt3 + acct2*dt
182 radt1 = radt3 + ((velt1+velt3) * dt * 0.5)
183 if (radt1 .lt. minrad)
then
185 testimpl = end_implode
190 lt1 = lconst * log(initrad/radt1)
194 ldott2 = (lt1 - lt3) * rdt
200 if (radt1 .gt. innerrad)
then
201 masst1 = mconst * (initrad * initrad - radt1 * radt1)
210 gaske = 0.5 * masst1 * velt1 * velt1
217 if (testimpl .eq. end_implode)
then
219 &
'(/A/,A/,A,1pe10.3/,A,1pe10.3/,A,1pe10.3/,A,1pe10.3/,A/)')
220 &
' ------------------------------------------------------------',
221 &
' Foil radius has reached the minimum value.',
223 &
' minimum radius: ',minrad,
224 &
' KE at min-radius: ',gaske,
225 &
' velocity at min-radius: ',velt1,
226 &
' ------------------------------------------------------------'
229 &
' ------------------------------------------------'
230 call writebuffer(buffer)
232 &
' Foil radius has reached the minimum value.'
233 call writebuffer(buffer)
234 write(buffer,
'(A,1pe10.3)')
' time: ',
time
235 call writebuffer(buffer)
236 write(buffer,
'(A,1pe10.3)')
' minimum radius: ',minrad
237 call writebuffer(buffer)
238 write(buffer,
'(A,1pe10.3)')
' KE at min-radius: ',gaske
239 call writebuffer(buffer)
240 write(buffer,
'(A,1pe10.3)')
' velocity at min-radius: ',velt1
241 call writebuffer(buffer)
243 &
' ------------------------------------------------'
244 call writebuffer(buffer)
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