////////////////// Start of input data ////////////////// 





########### section 0 section 0 choice of CPU or GPU ##########
GPU:                          False


########### section 1 choice of differencing scheme ###########
scheme:                         c
scheme_turb:                    u
acrank:                         1.0
acrank_conv:                    0.5
acrank_conv_keps:               1
acrank_conv_kom:                1


########### section 2 turbulence models ###########
cmu:                          0.09
pans:                         False
keps:                         False
kom_des:                      True
keps_des:                     False
k_eq_les:                     False
kom:                          False
sst:                          False
smag:                         False
wale:                         False
c_omega_1:                    0.556
c_omega_2:                    0.075
prand_k:                      2.0
prand_omega:                  2.0
cdes:                         0.67


########### section 3 restart/save ###########
restart:                      False
save:                         True


########### section 4 fluid properties ###########
viscos:                       1.92e-04


########### section 5 relaxation factors ###########
urfvis:                       0.5


########### section 6 number of iteration and convergence criterira ###########
sormax:                       0.001
norm_order:                   2
maxit:                        5
min_iter:                     1
solver_vel:                   lgmres
solver_p:                     pyamg
coeff_v:                      False
coeff_w:                      False
embedded:                     False
solver_turb:                  gmres
amg_relax:                    cg
amg_cycle:                    V
nsweep_vel:                   50
nsweep_keps:                  50
nsweep_kom:                   50
convergence_limit_u:          1e-05
convergence_limit_v:          1e-05
convergence_limit_w:          1e-05
convergence_limit_p:          0.0005
convergence_limit_k:          0.0001
convergence_limit_om:         1e-06


########### section 7 all variables are printed during the iteration at node ###########
imon:                         0
jmon:                         0
kmon:                         0


########### section 8 time-averaging ###########
ntstep:                       15000
dt[0]:                        1.17e-03
itstep_start:                 7500
itstep_save:                  2000
itstep_stats:                 1


########### section 9 residual scaling parameters ###########
resnorm_p:                    64.0
resnorm_vel:                  1280.0


########### Section 10 grid and boundary conditions ###########
ni:                           96
nj:                           96
nk:                           32
dz:                           5.00e-02


cyclic_x:                     False
cyclic_z:                     True


L_t_synt:                     0.2
nmodes_synt:                  150
dmin_synt:                    0.025
jmirror_synt:                 48


------boundary conditions for u
     u_bc_west_type:               d
     u_bc_east_type:               n
     u_bc_west[0,0]:               0.0
     u_bc_south_type:              d
     u_bc_north_type:              d
     u_bc_south[0,0]:              0.0
     u_bc_north[0,0]:              0.0
------boundary conditions for v
     v_bc_west_type:               d
     v_bc_east_type:               n
     v_bc_west[0,0]:               0.0
     v_bc_south_type:              d
     v_bc_north_type:              d
     v_bc_south[0,0]:              0.0
     v_bc_north[0,0]:              0.0
------boundary conditions for w
     w_bc_west_type:               d
     w_bc_east_type:               n
     w_bc_west[0,0]:               0.0
     w_bc_south_type:              d
     w_bc_north_type:              d
     w_bc_south[0,0]:              0.0
     w_bc_north[0,0]:              0.0
------boundary conditions for p
     p_bc_west_type:               n
     p_bc_east_type:               n
     p_bc_south_type:              n
     p_bc_north_type:              n
------boundary conditions for k
     k_bc_west_type:               d
     k_bc_east_type:               n
     k_bc_west[0,0]:               0.0
     k_bc_south_type:              d
     k_bc_north_type:              d
     k_bc_south[0,0]:              0.0
     k_bc_north[0,0]:              0.0
------boundary conditions for omega
     om_bc_west_type:              d
     om_bc_east_type:              n
     om_bc_west[0,0]:              0.0
     om_bc_south_type:             d
     om_bc_north_type:             d
     om_bc_south[0,0]:             18318761.7
     om_bc_north[0,0]:             18318761.7



 ////////////////// End of input data //////////////////



hostname:  tfd152.am.chalmers.se

sli:  0.2, visc: 1.92e-04, nmodes: 150, dxmin: 2.500e-02, dkn: 1.663e+00, dmin_synt: 2.500e-02

wnre:  3.73e+00, wnr1: 1.87e+00, epsm: 2.818e+01, wnrn: 2.513e+02

eigenvalue 1, 2 and 3: 1.848e+00, 3.368e-01, 8.151e-01

eigenvector R11, R12 and R13: -9.792e-01, -2.031e-01, 0.000e+00

eigenvector R21, R22 and R23: 2.031e-01, -9.792e-01, 0.000e+00

eigenvector R31, R32 and R33: 0.000e+00, 0.000e+00, 1.000e+00

time step no, 0
uvmean_time 1.6772932689883184
synt: uvmean [1.62980671] uv_rans_max= 0.959787802649839 at j= (array([21]),)
uin,a_synt,b_synt,tturb 23.576600425090536 0.9941577576591378 0.10793680042578163 0.2

sli:  0.2, visc: 1.92e-04, nmodes: 150, dxmin: 2.500e-02, dkn: 1.663e+00, dmin_synt: 2.500e-02

wnre:  3.73e+00, wnr1: 1.87e+00, epsm: 2.818e+01, wnrn: 2.513e+02

eigenvalue 1, 2 and 3: 1.848e+00, 3.368e-01, 8.151e-01

eigenvector R11, R12 and R13: -9.792e-01, -2.031e-01, 0.000e+00

eigenvector R21, R22 and R23: 2.031e-01, -9.792e-01, 0.000e+00

eigenvector R31, R32 and R33: 0.000e+00, 0.000e+00, 1.000e+00

time step no, 0
uvmean_time 1.6772932689883184
synt: uvmean [1.62980671] uv_rans_max= 0.959787802649839 at j= (array([21]),)
uvmean_aver at max 10.157114489238658

sli:  0.2, visc: 1.92e-04, nmodes: 150, dxmin: 2.500e-02, dkn: 1.663e+00, dmin_synt: 2.500e-02

wnre:  3.73e+00, wnr1: 1.87e+00, epsm: 2.818e+01, wnrn: 2.513e+02

eigenvalue 1, 2 and 3: 1.848e+00, 3.368e-01, 8.151e-01

eigenvector R11, R12 and R13: -9.792e-01, -2.031e-01, 0.000e+00

eigenvector R21, R22 and R23: 2.031e-01, -9.792e-01, 0.000e+00

eigenvector R31, R32 and R33: 0.000e+00, 0.000e+00, 1.000e+00

time step no, 0
uvmean_time 1.6772932689883184
synt: uvmean [1.62980671] uv_rans_max= 0.959787802649839 at j= (array([21]),)
uin,a_synt,b_synt,tturb 23.576600425090536 0.9941577576591378 0.10793680042578163 0.2

sli:  0.2, visc: 1.92e-04, nmodes: 150, dxmin: 2.500e-02, dkn: 1.663e+00, dmin_synt: 2.500e-02

wnre:  3.73e+00, wnr1: 1.87e+00, epsm: 2.818e+01, wnrn: 2.513e+02

eigenvalue 1, 2 and 3: 1.848e+00, 3.368e-01, 8.151e-01

eigenvector R11, R12 and R13: -9.792e-01, -2.031e-01, 0.000e+00

eigenvector R21, R22 and R23: 2.031e-01, -9.792e-01, 0.000e+00

eigenvector R31, R32 and R33: 0.000e+00, 0.000e+00, 1.000e+00

time step no, 0
uvmean_time 1.6772932689883184
synt: uvmean [1.62980671] uv_rans_max= 0.959787802649839 at j= (array([21]),)
uvmean_aver at max 10.157114489238658
conv called
vist_kom called

sli:  0.2, visc: 1.92e-04, nmodes: 150, dxmin: 2.500e-02, dkn: 1.663e+00, dmin_synt: 2.500e-02

wnre:  3.73e+00, wnr1: 1.87e+00, epsm: 2.818e+01, wnrn: 2.513e+02

eigenvalue 1, 2 and 3: 1.848e+00, 3.368e-01, 8.151e-01

eigenvector R11, R12 and R13: -9.792e-01, -2.031e-01, 0.000e+00

eigenvector R21, R22 and R23: 2.031e-01, -9.792e-01, 0.000e+00

eigenvector R31, R32 and R33: 0.000e+00, 0.000e+00, 1.000e+00

time step no, 0
uvmean_time 1.6772932689883184
synt: uvmean [1.62980671] uv_rans_max= 0.959787802649839 at j= (array([21]),)
uin,a_synt,b_synt,tturb 23.576600425090536 0.9941577576591378 0.10793680042578163 0.2

sli:  0.2, visc: 1.92e-04, nmodes: 150, dxmin: 2.500e-02, dkn: 1.663e+00, dmin_synt: 2.500e-02

wnre:  3.73e+00, wnr1: 1.87e+00, epsm: 2.818e+01, wnrn: 2.513e+02

eigenvalue 1, 2 and 3: 1.848e+00, 3.368e-01, 8.151e-01

eigenvector R11, R12 and R13: -9.792e-01, -2.031e-01, 0.000e+00

eigenvector R21, R22 and R23: 2.031e-01, -9.792e-01, 0.000e+00

eigenvector R31, R32 and R33: 0.000e+00, 0.000e+00, 1.000e+00

time step no, 0
uvmean_time 1.6772932689883184
synt: uvmean [1.62980671] uv_rans_max= 0.959787802649839 at j= (array([21]),)
uvmean_aver at max 10.157114489238658
CDS scheme, prand= 1
calcu called
solve_3d called
nmax,acrank_conv_local 50 0.5
solver in solve_3d: lgmres,tol,atol 1e-05 1e-05
variable= u
residual history in solve_3d. Variable u: initial residual:  2.75e+00    final residual: 2.00e-02          delta_phi: 1.66e+00
time solve_3d: 1.68e-01
time u: 2.39e-01
calcv called
solve_3d called
nmax,acrank_conv_local 50 0.5
solver in solve_3d: lgmres,tol,atol 1e-05 1e-05
variable= v
residual history in solve_3d. Variable v: initial residual:  1.42e+00    final residual: 1.03e-05          delta_phi: 7.48e-01
time solve_3d: 3.55e-01
time v: 3.62e-01
calcw called
solve_3d called
nmax,acrank_conv_local 50 0.5
solver in solve_3d: lgmres,tol,atol 1e-05 1e-05
variable= w
residual history in solve_3d. Variable w: initial residual:  1.90e+00    final residual: 1.48e-05          delta_phi: 9.34e-01
time solve_3d: 3.56e-01
time w: 3.63e-01
conv called
area_out 3.1999999999999993
flow_in 75.4451213602897 flow_out 75.43948412024857 area_out 3.1999999999999993 flow_out_new 75.4451213602897 uinc: 0.001761637512855963
calcp called
solve_p called
relation method= cg
A and M computed,tol_conv= 0.0005
in solve_p
Residual history in solve_p ['3.4091e+02', '5.5841e+01', '1.7096e+01', '8.5016e+00', '3.3331e+00', '2.9617e+00', '2.7296e+00', '2.3942e+00', '1.7574e+00', '2.7977e+00', '5.8571e+00', '7.2721e+00', '2.9530e+00', '1.0898e+00', '3.8000e-01', '1.7987e-01', '8.2950e-02']
time solve_p: 9.15e-01
time p: 9.51e-01
vist_kom called
upwind scheme, prand= 2.0
calck_kom called
solve_3d called
nmax,acrank_conv_local 50 1
solver in solve_3d: gmres
variable= k
residual history in solve_3d. Variable k: initial residual:  2.81e+00    final residual: 2.41e-03          delta_phi: 3.42e+00
time solve_3d: 4.13e-01
time k: 4.48e-01
upwind scheme, prand= 2.0
calcom called
solve_3d called
nmax,acrank_conv_local 50 1
solver in solve_3d: gmres
variable= om
residual history in solve_3d. Variable om: initial residual:  1.41e+07    final residual: 1.36e+01          delta_phi: 1.46e+07
time solve_3d: 7.33e-01
time omega: 7.52e-01
time Smag or Wale: 7.15e-07
compute_fk called

--time step:   0, iter:0, max residual:5.46e-05, u:1.56e-05, v:8.02e-09, w:1.15e-08, p:5.46e-05, k:1.47e-09, eps:0.00e+00, om:1.06e-02


monitor time step:   0, iter:0, u: 6.39e-01, v: 1.48e-03, w:-4.40e-04, p:-2.24e-13, k: 4.68e-01, eps: 1.00e+00, om: 1.46e+07


-- cfl_max_x: 3.24e-01, cfl_max_y: 8.34e-01, at i=  0, j=  8


No of points cfl_x > 1:  0, No of points cfl_y > 1:  0


--time step:  1.17e-03, iter:  0, umax: 2.59e+01, vismax: 4.87e+00, kmin: 1.56e-01, epsmin: 1.00e+00, ommin: 2.73e+02, fkmax: 1.00e+00

time one iteration: 3.12e+00
vismean_mean 5.1999999999999994e-17
time step no, 1
uvmean_time 1.1085700593403693
synt: uvmean [0.96777094] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  3.83e+00    final residual: 1.97e-02          delta_phi: 1.98e+00
time solve_3d: 1.49e-01
time u: 1.82e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  3.53e+00    final residual: 3.79e-05          delta_phi: 1.41e+00
time solve_3d: 2.85e-01
time v: 2.93e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  2.78e+00    final residual: 3.86e-05          delta_phi: 1.93e+00
time solve_3d: 3.73e-01
time w: 3.80e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.43423227880156 area_out 3.1999999999999993 flow_out_new 75.4451213602897 uinc: 0.0034028379650496996
in solve_p
Residual history in solve_p ['7.3263e+01', '1.1265e+01', '3.7216e+00', '2.9876e+00', '2.1888e+00', '1.4602e+00', '9.1993e-01', '4.6134e-01', '2.9798e-01', '3.6898e-01', '6.1083e-01', '9.0294e-01', '1.2738e+00', '9.5602e-01', '3.7963e-01', '1.3385e-01', '4.7618e-02', '1.4672e-02']
time solve_p: 6.43e-01
time p: 6.76e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  1.89e+02    final residual: 2.33e-03          delta_phi: 3.58e+00
time solve_3d: 2.08e-01
time k: 2.43e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  1.23e+09    final residual: 1.35e+01          delta_phi: 1.44e+07
time solve_3d: 2.38e-01
time omega: 2.57e-01
time Smag or Wale: 7.15e-07

--time step:   1, iter:0, max residual:1.54e-05, u:1.54e-05, v:2.96e-08, w:3.02e-08, p:1.05e-05, k:1.42e-09, eps:0.00e+00, om:1.06e-02


monitor time step:   1, iter:0, u: 6.57e-01, v: 5.27e-03, w: 2.91e-01, p: 6.06e-14, k: 6.74e-04, eps: 1.00e+00, om: 2.00e+05


--time step:  1.17e-03, iter:  0, umax: 2.66e+01, vismax: 6.22e+00, kmin: 5.69e-04, epsmin: 1.00e+00, ommin: 2.57e+02, fkmax: 1.00e+00

time one iteration: 2.03e+00
vismean_mean 2.5999999999999997e-17
time step no, 2
uvmean_time 0.9588594643903332
synt: uvmean [0.61778771] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  3.55e+00    final residual: 1.97e-02          delta_phi: 1.68e+00
time solve_3d: 1.72e-01
time u: 2.04e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  3.38e+00    final residual: 6.64e-05          delta_phi: 1.19e+00
time solve_3d: 2.96e-01
time v: 3.03e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  2.52e+00    final residual: 6.11e-05          delta_phi: 1.82e+00
time solve_3d: 3.65e-01
time w: 3.72e-01
area_out 3.1999999999999993
flow_in 75.44512136028973 flow_out 75.43328592978774 area_out 3.1999999999999993 flow_out_new 75.44512136028973 uinc: 0.0036985720318716284
in solve_p
Residual history in solve_p ['2.7981e+02', '5.0557e+01', '1.5602e+01', '7.4129e+00', '2.0556e+00', '9.3650e-01', '4.7509e-01', '5.0057e-01', '6.2844e-01', '5.4928e-01', '5.9713e-01', '1.1679e+00', '1.4147e+00', '1.1909e+00', '7.3584e-01', '3.7795e-01', '1.5313e-01', '6.5135e-02']
time solve_p: 6.44e-01
time p: 6.74e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  1.51e+00    final residual: 2.25e-03          delta_phi: 5.52e-01
time solve_3d: 2.49e-01
time k: 2.84e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  1.37e+07    final residual: 1.37e+01          delta_phi: 6.87e+06
time solve_3d: 4.88e-01
time omega: 5.08e-01
time Smag or Wale: 7.15e-07

--time step:   2, iter:0, max residual:5.17e-05, u:1.54e-05, v:5.19e-08, w:4.77e-08, p:5.17e-05, k:1.38e-09, eps:0.00e+00, om:1.07e-02


monitor time step:   2, iter:0, u: 4.86e-01, v: 5.47e-03, w: 5.29e-01, p:-6.66e-14, k: 1.17e-02, eps: 1.00e+00, om: 6.72e+06


--time step:  1.17e-03, iter:  0, umax: 2.73e+01, vismax: 7.85e+00, kmin: 5.80e-03, epsmin: 1.00e+00, ommin: 2.48e+02, fkmax: 1.00e+00

time one iteration: 2.35e+00
vismean_mean 1.733333333333333e-17
time step no, 3
uvmean_time 1.0835414414252784
synt: uvmean [0.77269712] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  2.12e+00    final residual: 2.05e-02          delta_phi: 1.27e+00
time solve_3d: 1.49e-01
time u: 1.82e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  1.52e+00    final residual: 7.17e-05          delta_phi: 4.87e-01
time solve_3d: 2.92e-01
time v: 3.00e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  1.38e+00    final residual: 6.63e-05          delta_phi: 7.16e-01
time solve_3d: 3.42e-01
time w: 3.49e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.43606377756578 area_out 3.1999999999999993 flow_out_new 75.44512136028972 uinc: 0.002830494601231415
in solve_p
Residual history in solve_p ['2.8469e+02', '5.3388e+01', '1.6138e+01', '7.9662e+00', '2.8904e+00', '2.3163e+00', '1.8983e+00', '1.4552e+00', '7.3125e-01', '7.7135e-01', '1.5253e+00', '3.4954e+00', '2.4347e+00', '9.9827e-01', '3.4510e-01', '1.6156e-01', '7.4107e-02']
time solve_p: 6.07e-01
time p: 6.39e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  7.54e+00    final residual: 2.06e-03          delta_phi: 5.47e-01
time solve_3d: 1.89e-01
time k: 2.24e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  2.61e+08    final residual: 1.37e+01          delta_phi: 6.68e+06
time solve_3d: 2.52e-01
time omega: 2.70e-01
time Smag or Wale: 9.54e-07

--time step:   3, iter:0, max residual:5.03e-05, u:1.60e-05, v:5.60e-08, w:5.18e-08, p:5.03e-05, k:1.26e-09, eps:0.00e+00, om:1.07e-02


monitor time step:   3, iter:0, u: 3.71e-01, v: 1.49e-03, w: 4.48e-01, p:-4.41e-13, k: 1.38e-04, eps: 1.00e+00, om: 4.13e+05


--time step:  1.17e-03, iter:  0, umax: 2.76e+01, vismax: 9.35e+00, kmin: 6.81e-05, epsmin: 1.00e+00, ommin: 2.42e+02, fkmax: 1.00e+00

time one iteration: 1.97e+00
vismean_mean 1.2999999999999999e-17
time step no, 4
uvmean_time 1.0120990287784957
synt: uvmean [0.59409282] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  2.09e+00    final residual: 2.00e-02          delta_phi: 8.75e-01
time solve_3d: 1.45e-01
time u: 1.77e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  3.64e+00    final residual: 6.18e-05          delta_phi: 1.12e+00
time solve_3d: 2.84e-01
time v: 2.91e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  2.66e+00    final residual: 6.67e-05          delta_phi: 1.59e+00
time solve_3d: 3.71e-01
time w: 3.78e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.43955410222148 area_out 3.1999999999999993 flow_out_new 75.4451213602897 uinc: 0.0017397681463249275
in solve_p
Residual history in solve_p ['8.6385e+01', '1.1668e+01', '3.9883e+00', '3.6834e+00', '2.4854e+00', '1.8770e+00', '1.0296e+00', '7.1987e-01', '8.0841e-01', '1.5840e+00', '2.3979e+00', '2.7232e+00', '1.9450e+00', '1.0401e+00', '3.9611e-01', '1.6419e-01', '5.8944e-02', '2.0271e-02']
time solve_p: 6.47e-01
time p: 6.78e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  1.10e+00    final residual: 2.13e-03          delta_phi: 3.10e-01
time solve_3d: 1.91e-01
time k: 2.26e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  1.29e+07    final residual: 1.35e+01          delta_phi: 4.34e+06
time solve_3d: 4.34e-01
time omega: 4.52e-01
time Smag or Wale: 9.54e-07

--time step:   4, iter:0, max residual:1.56e-05, u:1.56e-05, v:4.83e-08, w:5.21e-08, p:1.40e-05, k:1.30e-09, eps:0.00e+00, om:1.05e-02


monitor time step:   4, iter:0, u: 3.15e-01, v:-1.75e-03, w: 1.81e-01, p:-6.08e-14, k: 2.86e-03, eps: 1.00e+00, om: 4.57e+06


--time step:  1.17e-03, iter:  0, umax: 2.78e+01, vismax: 1.05e+01, kmin: 1.95e-03, epsmin: 1.00e+00, ommin: 2.38e+02, fkmax: 1.00e+00

time one iteration: 2.20e+00
vismean_mean 1.0399999999999999e-17
time step no, 5
uvmean_time 1.0293939534340686
synt: uvmean [0.7497456] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  2.02e+00    final residual: 1.95e-02          delta_phi: 9.06e-01
time solve_3d: 1.49e-01
time u: 1.81e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  3.61e+00    final residual: 6.24e-05          delta_phi: 1.27e+00
time solve_3d: 2.96e-01
time v: 3.03e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  2.69e+00    final residual: 7.01e-05          delta_phi: 1.57e+00
time solve_3d: 3.67e-01
time w: 3.74e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.44064056594779 area_out 3.1999999999999993 flow_out_new 75.44512136028972 uinc: 0.0014002482318531408
in solve_p
Residual history in solve_p ['2.1917e+02', '4.8237e+01', '1.8291e+01', '7.5322e+00', '2.3645e+00', '1.2476e+00', '6.8514e-01', '7.2465e-01', '5.9316e-01', '5.4417e-01', '6.2899e-01', '8.7355e-01', '1.1875e+00', '1.6920e+00', '1.3003e+00', '5.6445e-01', '2.4982e-01', '9.3188e-02']
time solve_p: 6.41e-01
time p: 6.72e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  1.49e+00    final residual: 2.03e-03          delta_phi: 6.92e-01
time solve_3d: 1.77e-01
time k: 2.11e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  1.10e+08    final residual: 1.41e+01          delta_phi: 4.16e+06
time solve_3d: 2.70e-01
time omega: 2.88e-01
time Smag or Wale: 7.15e-07

--time step:   5, iter:0, max residual:7.70e-05, u:1.53e-05, v:4.88e-08, w:5.48e-08, p:7.70e-05, k:1.24e-09, eps:0.00e+00, om:1.10e-02


monitor time step:   5, iter:0, u: 2.88e-01, v:-1.78e-03, w:-4.85e-02, p:-4.42e-14, k: 1.99e-05, eps: 1.00e+00, om: 5.93e+05


--time step:  1.17e-03, iter:  0, umax: 2.82e+01, vismax: 1.14e+01, kmin: 1.99e-05, epsmin: 1.00e+00, ommin: 2.35e+02, fkmax: 1.00e+00

time one iteration: 2.03e+00
vismean_mean 8.666666666666665e-18
time step no, 6
uvmean_time 1.0823992887306013
synt: uvmean [0.9482409] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  1.52e+00    final residual: 2.04e-02          delta_phi: 5.93e-01
time solve_3d: 1.27e-01
time u: 1.59e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  1.50e+00    final residual: 6.26e-05          delta_phi: 5.50e-01
time solve_3d: 3.09e-01
time v: 3.17e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  1.35e+00    final residual: 7.40e-05          delta_phi: 6.36e-01
time solve_3d: 3.37e-01
time w: 3.43e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.43849998164862 area_out 3.1999999999999993 flow_out_new 75.4451213602897 uinc: 0.0020691808253436865
in solve_p
Residual history in solve_p ['2.5077e+02', '5.1509e+01', '1.7379e+01', '7.8033e+00', '2.4187e+00', '1.4406e+00', '9.9337e-01', '9.5570e-01', '4.8031e-01', '2.8088e-01', '2.8637e-01', '5.2010e-01', '7.4029e-01', '1.0359e+00', '9.7110e-01', '4.7115e-01', '1.8371e-01', '6.6900e-02']
time solve_p: 6.45e-01
time p: 6.75e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  8.83e-01    final residual: 1.93e-03          delta_phi: 2.67e-01
time solve_3d: 1.78e-01
time k: 2.12e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  1.18e+07    final residual: 1.40e+01          delta_phi: 3.05e+06
time solve_3d: 4.14e-01
time omega: 4.33e-01
time Smag or Wale: 7.15e-07

--time step:   6, iter:0, max residual:5.59e-05, u:1.60e-05, v:4.89e-08, w:5.78e-08, p:5.59e-05, k:1.18e-09, eps:0.00e+00, om:1.09e-02


monitor time step:   6, iter:0, u: 2.62e-01, v: 1.23e-03, w:-8.23e-02, p: 1.05e-13, k: 3.18e-04, eps: 1.00e+00, om: 3.53e+06


--time step:  1.17e-03, iter:  0, umax: 2.85e+01, vismax: 1.20e+01, kmin: 3.18e-04, epsmin: 1.00e+00, ommin: 2.33e+02, fkmax: 1.00e+00

time one iteration: 2.14e+00
vismean_mean 7.428571428571428e-18
time step no, 7
uvmean_time 1.1260120221771313
synt: uvmean [0.94914884] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  2.36e+00    final residual: 1.96e-02          delta_phi: 1.07e+00
time solve_3d: 1.43e-01
time u: 1.75e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  2.95e+00    final residual: 6.53e-05          delta_phi: 5.31e-01
time solve_3d: 2.75e-01
time v: 2.83e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  1.86e+00    final residual: 7.59e-05          delta_phi: 9.49e-01
time solve_3d: 3.56e-01
time w: 3.63e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.43593491583994 area_out 3.1999999999999993 flow_out_new 75.44512136028972 uinc: 0.0028707638905567383
in solve_p
Residual history in solve_p ['1.1020e+02', '2.1971e+01', '9.1902e+00', '5.0324e+00', '3.5762e+00', '3.1627e+00', '2.2764e+00', '1.6235e+00', '1.6556e+00', '2.6293e+00', '4.1269e+00', '5.3936e+00', '4.8074e+00', '2.4835e+00', '1.0598e+00', '4.1288e-01', '1.7151e-01', '5.9945e-02', '2.2513e-02']
time solve_p: 6.79e-01
time p: 7.10e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  8.29e-01    final residual: 1.78e-03          delta_phi: 3.01e-01
time solve_3d: 1.64e-01
time k: 1.98e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  5.90e+07    final residual: 1.34e+01          delta_phi: 2.89e+06
time solve_3d: 2.63e-01
time omega: 2.81e-01
time Smag or Wale: 4.77e-07

--time step:   7, iter:0, max residual:1.77e-05, u:1.53e-05, v:5.10e-08, w:5.93e-08, p:1.77e-05, k:1.08e-09, eps:0.00e+00, om:1.05e-02


monitor time step:   7, iter:0, u: 2.32e-01, v: 4.26e-03, w: 5.14e-02, p: 1.31e-14, k: 4.19e-06, eps: 1.00e+00, om: 7.54e+05


--time step:  1.17e-03, iter:  0, umax: 2.86e+01, vismax: 1.24e+01, kmin: 4.19e-06, epsmin: 1.00e+00, ommin: 2.31e+02, fkmax: 1.00e+00

time one iteration: 2.01e+00
vismean_mean 6.499999999999999e-18
time step no, 8
uvmean_time 1.150172054502432
synt: uvmean [0.9207401] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  2.43e+00    final residual: 2.00e-02          delta_phi: 1.36e+00
time solve_3d: 1.44e-01
time u: 1.76e-01
variable= v
residual history in solve_3d. Variable v: initial residual:  2.73e+00    final residual: 7.87e-05          delta_phi: 7.51e-01
time solve_3d: 2.91e-01
time v: 2.99e-01
variable= w
residual history in solve_3d. Variable w: initial residual:  2.01e+00    final residual: 8.59e-05          delta_phi: 1.28e+00
time solve_3d: 3.67e-01
time w: 3.74e-01
area_out 3.1999999999999993
flow_in 75.44512136028972 flow_out 75.43571719883337 area_out 3.1999999999999993 flow_out_new 75.4451213602897 uinc: 0.0029388004551078603
in solve_p
Residual history in solve_p ['1.8529e+02', '4.5405e+01', '1.8543e+01', '7.0704e+00', '2.5135e+00', '1.4407e+00', '1.0670e+00', '1.3178e+00', '1.2328e+00', '1.3489e+00', '1.8785e+00', '2.3713e+00', '3.0161e+00', '2.8515e+00', '1.5033e+00', '5.9184e-01', '2.5970e-01', '9.4313e-02', '3.9748e-02']
time solve_p: 6.79e-01
time p: 7.09e-01
variable= k
residual history in solve_3d. Variable k: initial residual:  7.37e-01    final residual: 1.81e-03          delta_phi: 2.37e-01
time solve_3d: 1.64e-01
time k: 1.98e-01
variable= om
residual history in solve_3d. Variable om: initial residual:  1.06e+07    final residual: 1.31e+01          delta_phi: 2.25e+06
time solve_3d: 4.02e-01
time omega: 4.21e-01
time Smag or Wale: 4.77e-07

--time step:   8, iter:0, max residual:3.58e-05, u:1.56e-05, v:6.15e-08, w:6.71e-08, p:3.58e-05, k:1.11e-09, eps:0.00e+00, om:1.02e-02


monitor time step:   8, iter:0, u: 1.98e-01, v: 4.80e-03, w: 2.20e-01, p: 1.67e-13, k: 9.00e-05, eps: 1.00e+00, om: 2.93e+06


--time step:  1.17e-03, iter:  0, umax: 2.86e+01, vismax: 1.27e+01, kmin: 5.00e-05, epsmin: 1.00e+00, ommin: 2.27e+02, fkmax: 1.00e+00

time one iteration: 2.18e+00
vismean_mean 5.777777777777777e-18
time step no, 9
uvmean_time 1.0630265995817036
synt: uvmean [0.75420779] uv_rans_max= 0.959787802649839 at j= (array([21]),)
variable= u
residual history in solve_3d. Variable u: initial residual:  1.84e+00    final residual: 1.92e-02          delta_phi: 1.03e+00
time solve_3d: 1.31e-01
time u: 1.63e-01
