1 subroutine zdem(ifile,clen)
2 c -----------------------------------------------------------------------------
4 c ssss cccc rrrrr eeeeee aa m m eeeeee rrrrr
5 c s s
c c r r e a a mm mm e r r
6 c s
c r r e a a m mm m e r r
7 c ssss
c rrrrr eeeee aaaaaa m mm m eeeee rrrrr
8 c s
c r r e a a m m e r r
9 c s s
c c r r e a a m m e r r
10 c ssss cccc r r eeeeee a a m m eeeeee r r
12 c -----------------------------------------------------------------------------
16 c sandia national laboratories
19 c albuquerque nm 87185-1186
21 c mlkiefe@sandia.gov(e-mail)
26 c under the terms of contract de-ac04-94al85000,
27 c there is a non-exclusive license for
use of this work by
28 c or on behalf of the u.s. government
30 c -----------------------------------------------------------------------------
32 c series matrix load and solution
34 c single
line code for pulse power design
35 c array parameters - nr mesh points - 1 at top
36 c v,vn = old, new node voltage
37 c zib = current entering
branch segment
38 c zir,zirn = old, new long. current in r - l segment
39 c zlr = series inductance,
c = shunt capacitance,
40 c g = shunt conductance, rr = series resistance
41 c shunt capacitance at the half
time step.
42 c inductance(vacuum) = impedance/speed of light
43 c capacitance(vacuum) = 1/(impedance*speed of light)
44 c capacitance(vacuum) = epso*width*
length/height
45 c inductance(vacuum) = zmuo*height*
length/width
47 c -----------------------------------------------------------------------------
49 c include the files with the various keywords and
integer flags.
54 c include the files specifying the array dimensions and the
common blocks.
64 c define passed variables
69 c define internal variables
75 c---------------------------------------------------------
84 c get environment variables
86 c Read in the screamer input
file(assigned to
logical unit 4)
87 c this
subroutine also gives a
"nice" print-out of the data.
89 if (status .ne. 0)
then
92 c echo the set-up showing all values(slightly cryptic, but useful
94 if (echoset .eq. yes_echo)
then
98 c initialize the problem
100 if (status .ne. 0)
then
105 c---------------------------------------------------------
107 c first, put out a
banner page for the status cycle
printing for t > 0.
110 c ***** loop over all
time steps. *****
111 c tim is the problem
time and timehalf is the current half-step
time.
112 c icycle is a counter for printout of current circuit status.
115 timehalf = tim - htd2
120 c---------------------------------------------------------
123 c ***** now
do the output requests
if we executed all cycles.
subroutine init_problem(status)
c *****************************************************************************c Parameters for output selections c Added MFI CB parameters c Added callouts for the radiation yield in the c cylindrical foil model c Added line
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 branch
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
subroutine close_outfile(iunit, ierr)
subroutine banner_start_run(bgin)
subroutine get_env_var(inpfile, clen)
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 readscreamerdata(status)
subroutine banner_end_run(ifile, clen, bgin)
c *****************************************************************************c Common blocks for SCREAMER output printing
subroutine zdem(ifile, clen)
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line c
c This is a Fortran header file