Screamer Project  V3.3.1
Screamer Structure
 All Files Functions Variables
sphflprm.f
Go to the documentation of this file.
1  subroutine sphfoilparm (parms, nparms)
2 c
3 c Define passed variables
4 c
5  real parms(*)
6  integer nparms
7 c ------------------------------------------------------------------------------
8 c
9 c December 10, 1992; hnw
10 c 2014-02-06 RBS: Changed real*4 to real
11 c 2014-05-04 RBS: Changed integer*4 to integer
12 c
13 c Sets up the parameters needed for the SPHERICAL foil implosion model and
14 c returns them in parms(i).
15 c nparms is the number of parameters in parms.
16 c
17 c parms is sent with the basic parameters needed to rearrange and fill
18 c the actual parms array.
19 c ------------------------------------------------------------------------------
20 c
21  real temp(4)
22 c
23 c Include the common block for plotting so that we can set some
24 c initial plotting values.
25 c
26  include 'zdemmax.h'
27  include 'zdemout.h'
28 c
29 c ------------------------------------------------------------------------------
30 c Fill temp with parms.
31 c
32  do i = 1, nparms
33  temp(i) = parms(i)
34  end do
35 c
36 c ------------------------------------------------------------------------------
37 c In parms we want:
38 c 1: initrad (does not change)
39 c = initial foil radius
40 c 2: angl (does not change)
41 c = foil included angle
42 c 3: mass (does not change)
43 c = foil mass
44 c 4: minrad (does not change)
45 c = minimum foil radius
46 c 5: lt3 (initially = 0)
47 c = inductance from last half time step
48 c 6: radt3 (initially = initrad)
49 c = radius from last half time step
50 c 7: velt3 (initially = 0)
51 c = velocity from last half time step
52 c 8: testimpl (initially = 0, which is imploding, set to 1 when radius
53 c shrinks to minrad)
54 c = value to test to see if foil is still imploding or has stagnated at
55 c the minimum radius
56 c
57 c In temp:
58 c 1: initrad
59 c 2: angl
60 c 3: mass
61 c 4: minrad
62 c
63 c ------------------------------------------------------------------------------
64 c
65  nparms = 8
66  parms(1) = temp(1) ! initrad
67  parms(2) = temp(2) ! angl
68  parms(3) = temp(3) ! mass
69  parms(4) = temp(4) ! minrad
70 c
71  parms(5) = 0.0 ! lt3
72  parms(6) = temp(1) ! radt3
73  parms(7) = 0.0 ! velt3
74  parms(8) = 0.0 ! testimpl
75 c
76 c ------------------------------------------------------------------------------
77 c Set initial plotting values for foil radius, velocity,
78 c acceleration, kinetic energy.
79 c
80  foilrad = parms(1)
81  foilvel = 0.0
82  foilacc = 0.0
83  foilke = 0.0
84 c
85 c ------------------------------------------------------------------------------
86 c
87  return
88  end
subroutine sphfoilparm(parms, nparms)
Definition: sphflprm.f:1