SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.1 (2018/06/17 13:33:10) joern
SVN: run.in,v v. 1.1 (2018/06/17 13:33:10) joern
initialize_mpicomm: enabled MPI
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
The verbose level is ip= 14 (ldebug= F )
This is a 3-D run
nxgrid, nygrid, nzgrid= 64 64 64
Lx, Ly, Lz= 50000.000000000000 50000.000000000000 50000.000000000000
Vbox= 125000000000000.00
rsnap: read snapshot var.dat in 3.0353730000000001E-003 seconds
setup_slices: slice_position = m
setup_slices: ix_loc,iz2_loc, (video files) = 35 4
units_general: unit_velocity= 1000.0000000000000
units_general: unit_density= 1.0000000000000000E-004
units_general: unit_length= 1000.0000000000000
units_general: unit_magnetic= 3.5449077018110318
units_eos: unit_temperature= 10000.000000000000
units_eos: cp, lnTT0, cs0, pp0, Rgas= 311.63690737267581 -0.50454608698564718 11.199999999999999 150.60326399999997 83.144726887029918
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
-------------------------------------------------------------
Parameters to be set in run.in:
Kpara= 2.00E-05
Kperp= 2.76E-03
eta <= 1.00E+01
hcond_grad= 5.00E+12
Assume Umax=150 km/s and R=1 => diff= 1.19E+05
nu_spitzer= 2.21E-09
-------------------------------------------------------------
initialize_gravity: no x-gravity
initialize_gravity: no y-gravity
initialize_gravity: constant gravz= -0.27400000000000002
diffusion: div(D*grad(rho))
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
select_eos_variable: Using lnrho and lnTT
initialize_energy: heat conduction: constant chi.
resistivity: z-dependent
initialize_shock: prenormalised shock_factor sum= 1.0000000000000000
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
viscous force: nu_shock*(XXXXXXXXXXX)
Kspitzer_para= 2.00E-05
| solar_corona: settings for granules
-----------------------------------
| radius [Mm]: 2.3437500000000000
| lifetime [min] 5.0000000000000000
-----------------------------------
sh: 1: Syntax error: Bad fd number
pde: ENTER
SVN: ------- v. ( ) $Id$
No lists found, creating points for level: 1
No lists found, creating points for level: 2
No lists found, creating points for level: 3
calc_pencils_energy: max(advec_cs2) = 9.0200546319888766E-004
duu_dt: SOLVE
Bcs for ux, x: < p>, y: < p>, z: < s:a>
Bcs for uy, x: < p>, y: < p>, z: < s:a>
Bcs for uz, x: < p>, y: < p>, z: < div:a>
duu_dt: max(advec_uu) = 0.0000000000000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < p>, z: < fg:s>
dlnrho_dt: diffrho= 250.00000000000000
dlnrho_dt: max(diffus_diffrho ) = 1.2161000000000001E-003
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
SOLVE dlnTT_dt
denergy_dt: lnTT,cs2= -0.11369644490665429 185.42995296416569
Bcs for lnTT, x: < p>, y: < p>, z: < fg:s>
denergy_dt: cs20= 185.42995296416569
WARNING: calc_heat_cool_RTV: cool acts everywhere!
denergy_dt: max(diffus_chi ) = 0.20268333333333338
denergy_dt: max(diffus_chi3) = 0.0000000000000000
daa_dt: SOLVE
Bcs for Ax, x: < p>, y: < p>, z: < pfe>
Bcs for Ay, x: < p>, y: < p>, z: < nil>
Bcs for Az, x: < p>, y: < p>, z: < nil>
daa_dt: iresistivity=zdep
daa_dt: max(diffus_eta) = 9.7288000000000013E-002
daa_dt: max(diffus_eta2) = 0.0000000000000000
daa_dt: max(diffus_eta3) = 0.0000000000000000
dheatflux_dt: SOLVE dheatflux_dt
Bcs for heatfluxx, x: < p>, y: < p>, z: < a>
Bcs for heatfluxy, x: < p>, y: < p>, z: < a>
Bcs for heatfluxz, x: < p>, y: < p>, z: < a>
dspecial_dt: SOLVE dspecial_dt
---it---------t----------dt-------dtb------dteta------dtnu-----dtrad-----dtchi2--dtspitzer----dtu-------dtc-------urms-----urmsz------umax-----TTmax-----gTmax-----rhomin--rhomzmask-meshRemax---Mamax-----vAmax------j2m-------jm2-------brms-----brmsz-----dtchi------qmax------qrms------dtq-------dtq2----qsatmin---qsatrms---tauqmax--
1 4.6684E-02 4.67E-02 6.36E-01 1.82E-02 7.86E-01 8.99E-02 1.86E-02 2.38E+00 1.17E-03 1.65E-02 4.53E-01 4.53E-01 1.58E+01 1.00E+02 2.44E-02 4.08E-09 4.98E-02 2.27E-01 1.40E+00 1.69E+04 9.02E-02 1.29E+02 4.32E-03 4.32E-03 3.78E-02 9.52E+05 1.02E+05 3.10E-01 1.00E+00 9.71E-01 9.99E-01 1.00E+03
2 9.3369E-02 4.67E-02 6.36E-01 1.82E-02 7.80E-01 8.05E-02 2.15E-02 2.38E+00 1.17E-03 1.65E-02 4.56E-01 4.56E-01 1.58E+01 1.00E+02 2.45E-02 4.09E-09 4.98E-02 1.78E-01 1.40E+00 1.69E+04 9.02E-02 1.29E+02 4.32E-03 4.32E-03 3.78E-02 1.27E+06 1.43E+05 3.10E-01 1.00E+00 9.71E-01 9.99E-01 1.00E+03
3 1.4005E-01 4.67E-02 6.36E-01 1.82E-02 7.75E-01 7.20E-02 2.41E-02 2.37E+00 1.17E-03 1.65E-02 4.60E-01 4.60E-01 1.58E+01 1.00E+02 2.46E-02 4.09E-09 4.98E-02 1.48E-01 1.40E+00 1.69E+04 9.02E-02 1.29E+02 4.32E-03 4.32E-03 3.78E-02 1.36E+06 1.60E+05 3.10E-01 1.00E+00 9.71E-01 9.99E-01 1.00E+03
4 1.8674E-01 4.67E-02 6.36E-01 1.82E-02 7.70E-01 5.87E-02 2.62E-02 2.37E+00 1.17E-03 1.65E-02 4.64E-01 4.64E-01 1.58E+01 1.00E+02 2.47E-02 4.09E-09 4.98E-02 1.30E-01 1.40E+00 1.68E+04 9.02E-02 1.29E+02 4.32E-03 4.32E-03 3.78E-02 1.39E+06 1.67E+05 3.10E-01 1.00E+00 9.71E-01 9.99E-01 1.00E+03
Simulation finished after 5 time-steps
Writing final snapshot at time t = 0.23342542270180064
Wall clock time [hours] = 5.903E-04 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 1.621274
Wall clock time/timestep/local meshpoint [microsec] = 6.485095
Rhs wall clock time/timestep/local meshpoint [microsec] = 5.802945
Maximum used memory per cpu [MBytes] = 53.770
Maximum used memory [MBytes] = 213.082
real 2.54
user 8.84
sys 0.20