Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
!**********************************************************************
!
!
SUBROUTINE null
!
!
!***********************************************************************
!
!
! Subroutine to inizialize arrays and variables
!
!
!=======================================================================
USE fly_h
implicit none
INTEGER :: i_c1
!=======================================================================
nbodies=0
one=1
tol=0.
tol2inv=0.
eps=0.
usequad=.TRUE.
rmt_acc=.FALSE.
rsize=0.
rmin=0.
incells=0
root=nb_loc+1
subp=0.
mass=0.
mass_read=0.
pos=0.
pos_rmt=0.0
acc_rmt=0.0
mass_cell=0.
pos_cell=0.
vel=0.
acc=0.
quad=0.
tnow=0.
tpos=0.
dtime=0.
dtime2=0.
tstep=0.
phsmooth=0.
acsmooth=0.
rk=0.
fc=0.
znow=0.
inasc='B'
ouasc='B'
z32_ou=0.
z32_end=0.
z32_last=0
next_out=0
max_time=0
hubble=0.0
wr_out32=.FALSE.
dt_var=.TRUE.
halt_sim=0
PE0=0
!ori lpes=LEADZ(0) -LEADZ(npes) -1
!ori nb_res_loc=0
i_c1=0
DO WHILE(.TRUE.)
IF(2**i_c1.ge.npes) EXIT
i_c1=i_c1+1
END DO
IF(2**i_c1.gt.npes) i_c1=i_c1-1
lpes=i_c1
!-------------------------------------------------------------------------------
! Definitions specific to input/output filenames (fly_names).
!-------------------------------------------------------------------------------
! IF(me.eq.0) THEN
fly_fnames='fly_fnames'
OPEN(UNIT=ufnames,FILE=fly_fnames,STATUS='OLD')
READ(ufnames,30) parsfile
READ(ufnames,30) par1file
READ(ufnames,30) fewgrid
READ(ufnames,30) fewtab
READ(ufnames,30) fpar35
READ(ufnames,30) ql_tab
CLOSE(UNIT=ufnames)
parsfile=TRIM(parsfile)
par1file=TRIM(par1file)
fewgrid=TRIM(fewgrid)
fewtab=TRIM(fewtab)
fpar35=TRIM(fpar35)
ql_tab=TRIM(ql_tab)
IF(me.eq.0) THEN
WRITE(uterm,*) 'INPUT PARAMETER FILENAMES'
WRITE(uterm,*) '-------------------------------------------'
WRITE(uterm,130) 'STAT_PARS =',parsfile(1:LEN_TRIM(parsfile))
WRITE(uterm,130) 'DYN_PARS =',par1file(1:LEN_TRIM(par1file))
WRITE(uterm,130) 'EW_GRID =',fewgrid(1:LEN_TRIM(fewgrid))
WRITE(uterm,130) 'EW_TAB =',fewtab(1:LEN_TRIM(fewtab))
WRITE(uterm,130) 'OUT32_TAB =',fpar35(1:LEN_TRIM(fpar35))
WRITE(uterm,130) 'QL_TAB =',ql_tab(1:LEN_TRIM(ql_tab))
WRITE(uterm,*) '-------------------------------------------'
call flush(uterm)
ENDIF
30 FORMAT(12x,A256)
130 FORMAT(A12,A)
RETURN
END