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





########### section 1 choice of differencing scheme ###########
scheme:                         h
scheme_turb:                    h


########### section 2 turbulence models ###########
cmu:                          0.09
kom:                          True
keps:                         False
EARSM:                        False
PINN:                         False
c_omega_1:                    0.556
c_omega_2:                    0.075
prand_k:                      2.0
prand_omega:                  2.0


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


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


########### section 5 relaxation factors ###########
urfvis:                       0.5
urf_vel:                      0.5
urf_p:                        1.0
urf_k:                        0.5


########### section 6 number of iteration and convergence criterira ###########
sormax:                       1e-07
maxit:                        3000
solver_pp:                    lgmres
solver_vel:                   lgmres
solver_turb:                  lgmres
nsweep_vel:                   50
nsweep_pp:                    50
nsweep_turb:                  50
convergence_limit_u:          -1e-06
convergence_limit_v:          -1e-06
convergence_limit_pp:         0.0005
convergence_limit_k:          -1e-06
convergence_limit_om:         -1e-06


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


########### section 8 time-averaging ###########


########### section 9 residual scaling parameters ###########
resnorm_p:                    40.0
resnorm_vel:                  800.0


########### Section 10 grid and boundary conditions ###########
ni:                           3
nj:                           34
cyclic_x:                     True




------boundary conditions for u
     u_bc_west_type:               n
     u_bc_east_type:               n
     u_bc_south_type:              n
     u_bc_north_type:              n
------boundary conditions for v
     v_bc_west_type:               n
     v_bc_east_type:               n
     v_bc_south_type:              d
     v_bc_north_type:              d
     v_bc_south[0]:                0.0
     v_bc_north[0]:                0.0
------boundary conditions for k
     k_bc_west_type:               n
     k_bc_east_type:               n
     k_bc_south_type:              d
     k_bc_north_type:              d
     k_bc_south[0]:                0.0
     k_bc_north[0]:                0.0
------boundary conditions for omega
     om_bc_west_type:              n
     om_bc_east_type:              n
     om_bc_south_type:             d
     om_bc_north_type:             d
     om_bc_south[0]:               711.4
     om_bc_north[0]:               711.4



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



hostname:  tfd152.am.chalmers.se
read_restart_data called
vist_kom called
hybrid scheme, prand= 1
calcu called
south: tauw [0.99996997 0.99996997 0.99996997]
north: tauw [0.99996997 0.99996997 0.99996997]
south: np.sum(su2d[:,0]).np.sum(su2d) -2.823439330550881 0.00018016242765028778
north: np.sum(su2d[:,-1],np.sum(vol) -2.8234393305508805 5.999999999999998
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  1.89e-05              final residual: 1.06e-11                    delta_phi: 4.51e-07
time u: 1.09e-02
calcv called
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  2.38e-14              final residual: 6.96e-21                    delta_phi: 5.57e-16
time v: 2.37e-03
calcp called
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  9.17e-11              final residual: 9.17e-11                    delta_phi: 0.00e+00
correct_u_v_p called
time pp: 1.35e-03
vist_kom called
hybrid scheme, prand= [[2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.
  2. 2. 2. 2. 2. 2. 2. 2. 2. 2.]
 [2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.
  2. 2. 2. 2. 2. 2. 2. 2. 2. 2.]
 [2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.
  2. 2. 2. 2. 2. 2. 2. 2. 2. 2.]]
calck_kom called
south: dy 0.029411764705882353
south wall: dy,ustar_min,max,mean 0.029411764705882353 0.9999850054210294 0.9999850054210335 0.9999849863516577
south wall: kwall_min,max 3.3332333702229873 3.333233370223015
north: dy 0.02941176470588225
north wall: dy,ustar_min,max,mean 0.02941176470588225 1.0 1.0 0.9999849863516577
north wall: kwall_min,max 3.3332333702229873 3.333233370223015
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  3.10e-06              final residual: 2.04e-12                    delta_phi: 1.29e-07
time k: 5.65e-03
hybrid scheme, prand= 2.0
calcom called
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  8.60e-05              final residual: 2.95e-11                    delta_phi: 5.27e-06
time omega: 1.35e-03

--iter:0, max residual:1.08e-07, u:2.36e-08, v:2.97e-17, cont:2.29e-12, k:4.85e-12, om:1.08e-07, eps:0.00e+00


monitor iteration:   0, u: 2.19e+01, v:-2.37e-12, p: 2.73e-11, k: 1.35e+00, om: 1.22e+01, eps: 1.00e-10

kmax 3.3332333702230126

---iter:  0, umax: 2.26e+01, vismax: 1.46e+02, kmin: 7.83e-01, ommin: 7.25e+00, epsmmin: 1.00e-10

time one iteration: 2.18e-02
south: tauw [0.99997001 0.99997001 0.99997001]
north: tauw [0.99997001 0.99997001 0.99997001]
south: np.sum(su2d[:,0]).np.sum(su2d) -2.8234394449654037 0.00017993359860468683
north: np.sum(su2d[:,-1],np.sum(vol) -2.8234394449654037 5.999999999999998
residual history in solve_2d: initial residual:  1.88e-05              final residual: 1.05e-11                    delta_phi: 4.51e-07
time u: 3.38e-03
residual history in solve_2d: initial residual:  2.37e-14              final residual: 2.19e-20                    delta_phi: 5.57e-16
time v: 1.00e-03
residual history in solve_2d: initial residual:  9.17e-11              final residual: 9.17e-11                    delta_phi: 0.00e+00
time pp: 8.58e-04
south: dy 0.029411764705882353
south wall: dy,ustar_min,max,mean 0.029411764705882353 0.999985023382161 0.9999850233821651 0.999985005421031
south wall: kwall_min,max 3.3332334899620704 3.333233489962098
north: dy 0.02941176470588225
north wall: dy,ustar_min,max,mean 0.02941176470588225 1.0 1.0 0.999985005421031
north wall: kwall_min,max 3.3332334899620704 3.333233489962098
fix_k saved
ustar_mean_s 0.999985005421031
residual history in solve_2d: initial residual:  2.98e-06              final residual: 2.33e-12                    delta_phi: 1.20e-07
time k: 2.72e-03
residual history in solve_2d: initial residual:  8.12e-05              final residual: 4.86e-11                    delta_phi: 4.96e-06
time omega: 8.61e-04

--iter:1, max residual:1.01e-07, u:2.36e-08, v:2.97e-17, cont:2.29e-12, k:4.66e-12, om:1.01e-07, eps:0.00e+00


monitor iteration:   1, u: 2.19e+01, v:-2.38e-12, p: 2.73e-11, k: 1.35e+00, om: 1.22e+01, eps: 1.00e-10

kmax 3.3332334899620952

---iter:  1, umax: 2.26e+01, vismax: 1.46e+02, kmin: 7.83e-01, ommin: 7.25e+00, epsmmin: 1.00e-10

time one iteration: 8.98e-03
south: tauw [0.99997005 0.99997005 0.99997005]
north: tauw [0.99997005 0.99997005 0.99997005]
south: np.sum(su2d[:,0]).np.sum(su2d) -2.8234395527305782 0.00017971806825511916
north: np.sum(su2d[:,-1],np.sum(vol) -2.8234395527305782 5.999999999999998
residual history in solve_2d: initial residual:  1.88e-05              final residual: 1.05e-11                    delta_phi: 4.50e-07
time u: 2.82e-03
residual history in solve_2d: initial residual:  2.37e-14              final residual: 1.97e-20                    delta_phi: 5.57e-16
time v: 7.72e-04
residual history in solve_2d: initial residual:  9.17e-11              final residual: 9.17e-11                    delta_phi: 0.00e+00
time pp: 6.41e-04
south: dy 0.029411764705882353
south wall: dy,ustar_min,max,mean 0.029411764705882353 0.9999850409340119 0.999985040934016 0.9999850233821626
south wall: kwall_min,max 3.3332336069726582 3.3332336069726853
north: dy 0.02941176470588225
north wall: dy,ustar_min,max,mean 0.02941176470588225 1.0 1.0 0.9999850233821626
north wall: kwall_min,max 3.3332336069726582 3.3332336069726853
fix_k saved
ustar_mean_s 0.9999850233821626
residual history in solve_2d: initial residual:  2.94e-06              final residual: 9.27e-13                    delta_phi: 1.19e-07
time k: 2.91e-03
residual history in solve_2d: initial residual:  7.94e-05              final residual: 5.83e-11                    delta_phi: 4.85e-06
time omega: 7.24e-04

--iter:2, max residual:9.92e-08, u:2.35e-08, v:2.97e-17, cont:2.29e-12, k:4.59e-12, om:9.92e-08, eps:0.00e+00


monitor iteration:   2, u: 2.19e+01, v:-2.38e-12, p: 2.73e-11, k: 1.35e+00, om: 1.22e+01, eps: 1.00e-10

kmax 3.333233606972684

---iter:  2, umax: 2.26e+01, vismax: 1.46e+02, kmin: 7.83e-01, ommin: 7.25e+00, epsmmin: 1.00e-10

time one iteration: 8.00e-03
save_data called
program reached normal stop
////////////////// Start of input data ////////////////// 





########### section 1 choice of differencing scheme ###########
scheme:                         h
scheme_turb:                    h


########### section 2 turbulence models ###########
cmu:                          0.09
kom:                          True
keps:                         False
EARSM:                        False
PINN:                         False
c_omega_1:                    0.556
c_omega_2:                    0.075
prand_k:                      2.0
prand_omega:                  2.0


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


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


########### section 5 relaxation factors ###########
urfvis:                       0.5
urf_vel:                      0.5
urf_p:                        1.0
urf_k:                        0.5


########### section 6 number of iteration and convergence criterira ###########
sormax:                       1e-07
maxit:                        3000
solver_pp:                    lgmres
solver_vel:                   lgmres
solver_turb:                  lgmres
nsweep_vel:                   50
nsweep_pp:                    50
nsweep_turb:                  50
convergence_limit_u:          -1e-06
convergence_limit_v:          -1e-06
convergence_limit_pp:         0.0005
convergence_limit_k:          -1e-06
convergence_limit_om:         -1e-06


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


########### section 8 time-averaging ###########


########### section 9 residual scaling parameters ###########
resnorm_p:                    40.0
resnorm_vel:                  800.0


########### Section 10 grid and boundary conditions ###########
ni:                           3
nj:                           34
cyclic_x:                     True




------boundary conditions for u
     u_bc_west_type:               n
     u_bc_east_type:               n
     u_bc_south_type:              n
     u_bc_north_type:              n
------boundary conditions for v
     v_bc_west_type:               n
     v_bc_east_type:               n
     v_bc_south_type:              d
     v_bc_north_type:              d
     v_bc_south[0]:                0.0
     v_bc_north[0]:                0.0
------boundary conditions for k
     k_bc_west_type:               n
     k_bc_east_type:               n
     k_bc_south_type:              d
     k_bc_north_type:              d
     k_bc_south[0]:                0.0
     k_bc_north[0]:                0.0
------boundary conditions for omega
     om_bc_west_type:              n
     om_bc_east_type:              n
     om_bc_south_type:             d
     om_bc_north_type:             d
     om_bc_south[0]:               711.4
     om_bc_north[0]:               711.4



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



hostname:  tfd152.am.chalmers.se
read_restart_data called
vist_kom called
hybrid scheme, prand= 1
calcu called
south: tauw [0.99996997 0.99996997 0.99996997]
north: tauw [0.99996997 0.99996997 0.99996997]
south: np.sum(su2d[:,0]).np.sum(su2d) -2.823439330550881 0.00018016242765028778
north: np.sum(su2d[:,-1],np.sum(vol) -2.8234393305508805 5.999999999999998
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  1.89e-05              final residual: 1.06e-11                    delta_phi: 4.51e-07
time u: 1.09e-02
calcv called
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  2.38e-14              final residual: 6.96e-21                    delta_phi: 5.57e-16
time v: 2.37e-03
calcp called
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  9.17e-11              final residual: 9.17e-11                    delta_phi: 0.00e+00
correct_u_v_p called
time pp: 1.35e-03
vist_kom called
hybrid scheme, prand= [[2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.
  2. 2. 2. 2. 2. 2. 2. 2. 2. 2.]
 [2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.
  2. 2. 2. 2. 2. 2. 2. 2. 2. 2.]
 [2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2. 2.
  2. 2. 2. 2. 2. 2. 2. 2. 2. 2.]]
calck_kom called
south: dy 0.029411764705882353
south wall: dy,ustar_min,max,mean 0.029411764705882353 0.9999850054210294 0.9999850054210335 0.9999849863516577
south wall: kwall_min,max 3.3332333702229873 3.333233370223015
north: dy 0.02941176470588225
north wall: dy,ustar_min,max,mean 0.02941176470588225 1.0 1.0 0.9999849863516577
north wall: kwall_min,max 3.3332333702229873 3.333233370223015
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  3.10e-06              final residual: 2.04e-12                    delta_phi: 1.29e-07
time k: 5.65e-03
hybrid scheme, prand= 2.0
calcom called
solve_2d called
solver_local lgmres
nmax 50
solver in solve_2d: lgmres
residual history in solve_2d: initial residual:  8.60e-05              final residual: 2.95e-11                    delta_phi: 5.27e-06
time omega: 1.35e-03

--iter:0, max residual:1.08e-07, u:2.36e-08, v:2.97e-17, cont:2.29e-12, k:4.85e-12, om:1.08e-07, eps:0.00e+00


monitor iteration:   0, u: 2.19e+01, v:-2.37e-12, p: 2.73e-11, k: 1.35e+00, om: 1.22e+01, eps: 1.00e-10

kmax 3.3332333702230126

---iter:  0, umax: 2.26e+01, vismax: 1.46e+02, kmin: 7.83e-01, ommin: 7.25e+00, epsmmin: 1.00e-10

time one iteration: 2.18e-02
south: tauw [0.99997001 0.99997001 0.99997001]
north: tauw [0.99997001 0.99997001 0.99997001]
south: np.sum(su2d[:,0]).np.sum(su2d) -2.8234394449654037 0.00017993359860468683
north: np.sum(su2d[:,-1],np.sum(vol) -2.8234394449654037 5.999999999999998
residual history in solve_2d: initial residual:  1.88e-05              final residual: 1.05e-11                    delta_phi: 4.51e-07
time u: 3.38e-03
residual history in solve_2d: initial residual:  2.37e-14              final residual: 2.19e-20                    delta_phi: 5.57e-16
time v: 1.00e-03
residual history in solve_2d: initial residual:  9.17e-11              final residual: 9.17e-11                    delta_phi: 0.00e+00
time pp: 8.58e-04
south: dy 0.029411764705882353
south wall: dy,ustar_min,max,mean 0.029411764705882353 0.999985023382161 0.9999850233821651 0.999985005421031
south wall: kwall_min,max 3.3332334899620704 3.333233489962098
north: dy 0.02941176470588225
north wall: dy,ustar_min,max,mean 0.02941176470588225 1.0 1.0 0.999985005421031
north wall: kwall_min,max 3.3332334899620704 3.333233489962098
fix_k saved
ustar_mean_s 0.999985005421031
residual history in solve_2d: initial residual:  2.98e-06              final residual: 2.33e-12                    delta_phi: 1.20e-07
time k: 2.72e-03
residual history in solve_2d: initial residual:  8.12e-05              final residual: 4.86e-11                    delta_phi: 4.96e-06
time omega: 8.61e-04

--iter:1, max residual:1.01e-07, u:2.36e-08, v:2.97e-17, cont:2.29e-12, k:4.66e-12, om:1.01e-07, eps:0.00e+00


monitor iteration:   1, u: 2.19e+01, v:-2.38e-12, p: 2.73e-11, k: 1.35e+00, om: 1.22e+01, eps: 1.00e-10

kmax 3.3332334899620952

---iter:  1, umax: 2.26e+01, vismax: 1.46e+02, kmin: 7.83e-01, ommin: 7.25e+00, epsmmin: 1.00e-10

time one iteration: 8.98e-03
south: tauw [0.99997005 0.99997005 0.99997005]
north: tauw [0.99997005 0.99997005 0.99997005]
south: np.sum(su2d[:,0]).np.sum(su2d) -2.8234395527305782 0.00017971806825511916
north: np.sum(su2d[:,-1],np.sum(vol) -2.8234395527305782 5.999999999999998
residual history in solve_2d: initial residual:  1.88e-05              final residual: 1.05e-11                    delta_phi: 4.50e-07
time u: 2.82e-03
residual history in solve_2d: initial residual:  2.37e-14              final residual: 1.97e-20                    delta_phi: 5.57e-16
time v: 7.72e-04
residual history in solve_2d: initial residual:  9.17e-11              final residual: 9.17e-11                    delta_phi: 0.00e+00
time pp: 6.41e-04
south: dy 0.029411764705882353
south wall: dy,ustar_min,max,mean 0.029411764705882353 0.9999850409340119 0.999985040934016 0.9999850233821626
south wall: kwall_min,max 3.3332336069726582 3.3332336069726853
north: dy 0.02941176470588225
north wall: dy,ustar_min,max,mean 0.02941176470588225 1.0 1.0 0.9999850233821626
north wall: kwall_min,max 3.3332336069726582 3.3332336069726853
fix_k saved
ustar_mean_s 0.9999850233821626
residual history in solve_2d: initial residual:  2.94e-06              final residual: 9.27e-13                    delta_phi: 1.19e-07
time k: 2.91e-03
residual history in solve_2d: initial residual:  7.94e-05              final residual: 5.83e-11                    delta_phi: 4.85e-06
time omega: 7.24e-04

--iter:2, max residual:9.92e-08, u:2.35e-08, v:2.97e-17, cont:2.29e-12, k:4.59e-12, om:9.92e-08, eps:0.00e+00


monitor iteration:   2, u: 2.19e+01, v:-2.38e-12, p: 2.73e-11, k: 1.35e+00, om: 1.22e+01, eps: 1.00e-10

kmax 3.333233606972684

---iter:  2, umax: 2.26e+01, vismax: 1.46e+02, kmin: 7.83e-01, ommin: 7.25e+00, epsmmin: 1.00e-10

time one iteration: 8.00e-03
save_data called
program reached normal stop
