2
### BL test COMPRESSIBLE
3
### First part - before first *** exit *** statement
4
### is for scalar variables that will set up the right
7
# Problem type: options -2, -1, 0, and 1
11
# total number of elements
15
# nlev --> controls u_snloc
17
# Number of velocity nodes
19
# Number of spatial nodes associated w velocity
21
# number of nodes associated with CV
23
# Number of spatial nodes
25
# Number of pressure nodes
27
# Number of surface nodes associated with CV
29
# Number of surface nodes associated w velocity
31
# Number of surface pressure nodes
33
# Number of surface spatial nodes
35
# Total number of surface elements
37
# Number of coeffs for polynomial representation for EOS
39
# Number of coeffs for polynomial representation of permeability/absorption
41
# Option for element type: Velocity mesh
43
# Option for element type: Pressure mesh
45
# Option for element type: Material mesh
47
# Option for element type: CV mesh
49
# Option for element type: CV surface mesh
51
# Option for element type: Velocity surface mesh
53
# Total number of time dumps
55
# Total number of non-linear iterations
57
# Option for modified CMC color algorithm(/=0), if <0==>ndpset=cv_nonods
61
# Reference pressure (not in use)
63
# Initial pressure (not in use)
65
# Initial temperature (or any scalar field, not in use)
67
# Conservative(1.)/non-conservative(0.) flag for CV discretisation
69
# Conservative(1.)/non-conservative(0.) flag for Volfra discretisation
71
# Time-stepping discretisation parameter
73
# Time-stepping discretisation parameter
75
# Time-stepping discretisation parameter (not in use)
77
# Disopt: discretisation option in space and time (see cv-adv-dif.f90)
78
# = 8: Finite elements in space (Theta = specified; Limiting: Downwind+)
80
# Disopt: discretisation option in space and time (see cv-adv-dif.f90)
81
# = 0: 1st order in space (Theta=specified; Limiting: universal)
83
# Disopt: discretisation option in space and time (see cv-adv-dif.f90)
84
# = 0: 1st order in space (Theta=specified; Limiting: universal)
86
# t_dg_vel_int_opt: interface scalar field calculation option between elements
88
# u_dg_vel_int_opt: interface velocity calculation option between elements
90
# v_dg_vel_int_opt: interface velocitycalculation option between elements
92
# w_dg_vel_int_opt: interface velocitycalculation option between elements
94
# Lump_eqns: Lump multiphase flow equations
99
### end of scalars variables necessary for the arrays read after
105
# If value < -1000, then a pre-defined function may be used.
107
# wic_vol_bc( stotel * nphase )
108
wic_vol_bc -1001 input_wic_vol_bc_BL1_fcn
110
# wic_d_bc( stotel * nphase )
111
wic_d_bc -1001 input_wic_d_bc_BL1_fcn
113
# wic_u_bc( stotel * nphase )
114
wic_u_bc -1001 input_wic_u_bc_BL1_fcn
116
# wic_p_bc( stotel * nphase )
117
wic_p_bc -1001 input_wic_p_bc_BL1_fcn
119
# wic_t_bc( stotel * nphase )
122
# suf_vol_bc( stotel * cv_snloc * nphase )
123
suf_vol_bc -1001. input_suf_vol_bc_BL1_fcn
125
# suf_d_bc( stotel * cv_snloc * nphase )
126
suf_d_bc -1001. input_suf_d_bc_BL1_fcn
128
# suf_cpd_bc( stotel * cv_snloc * nphase )
131
# suf_t_bc( stotel * cv_snloc * nphase )
132
suf_t_bc -1001 input_suf_t_bc_BL1_fcn
134
# suf_p_bc ( stotel * p_snloc * nphase )
137
# suf_u_bc( stotel * u_snloc * nphase )
138
suf_u_bc -1001. input_suf_u_bc_BL1_fcn
140
# suf_v_bc( stotel * u_snloc * nphase )
143
# suf_w_bc( stotel * u_snloc * nphase )
146
# suf_one_bc( stotel * cv_snloc * nphase )
149
# suf_u_bc_rob1( stotel * u_snloc * nphase )
152
# suf_u_bc_rob2( stotel * u_snloc * nphase )
155
# suf_v_bc_rob1( stotel * u_snloc * nphase )
158
# suf_v_bc_rob2( stotel * u_snloc * nphase )
161
# suf_w_bc_rob1( stotel * u_snloc * nphase )
164
# suf_w_bc_rob2( stotel * u_snloc * nphase )
167
# suf_t_bc_rob1( stotel * cv_snloc * nphase )
170
# suf_t_bc_rob2( stotel * cv_snloc * nphase )
191
# nu( u_nonods * nphase )
194
# nv( u_nonods * nphase )
197
# nw( u_nonods * nphase )
200
# ug( u_nonods * nphase )
203
# vg( u_nonods * nphase )
206
# wg( u_nonods * nphase )
209
# uabs_option( nphase )
212
# uabs_coefs( nphase, nuabs_coefs )
215
# u_abs_stab ( mat_nonods, ndim * nphase, ndim * nphase )
218
# u_absorb ( mat_nonods, ndim * nphase, ndim * nphase )
221
# t_absorb ( cv_pha_nonods, nphase, nphase )
224
# v_absorb ( cv_pha_nonods, nphase, nphase )
227
# u_source( u_pha_nonods )
230
# t_source( cv_pha_nonods )
233
# v_source( cv_pha_nonods )
236
# udiffusion( mat_nonods, ndim, ndim, nphase )
239
# tdiffusion( mat_nonods, ndim, ndim, nphase )
251
# den ( cv_pha_nonods )
254
# satura ( cv_pha_nonods )
255
satura -1001. input_satura_BL1_func
257
# volfra ( cv_pha_nonods )
260
# t ( cv_pha_nonods )
263
# cv_one( cv_pha_nonods )
272
# volfra_pore ( totele )
275
# perm ( totele, ndim, ndim )
276
perm -1001. input_perm_BL1_func
278
# eos_option( nphase )
281
# cp_option( nphase )
284
# eos_coefs( nphase, ncoef )
285
eos_coefs -1001. input_eos_coef_BL1_func
287
# cp_coefs( nphase, ncp_coefs )
290
# end of the input file