Pencil1: iapp, nprocs, ncpus= 0 4 4 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 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 ------------------------------------------------------------- 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 ----------------------------------- 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.0000000000000000E+000 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.0000000000000000E+000 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.0000000000000000E+000 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.0000000000000000E+000 daa_dt: max(diffus_eta3) = 0.0000000000000000E+000 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] = 9.164E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 2.517081 Maximum used memory per cpu [MBytes] = 54.152 Maximum used memory [MBytes] = 215.234 real 3.78 user 13.63 sys 0.26