Pencil1: iapp, nprocs, ncpus= 0 8 8 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2015/05/11 19:37:40) wlyra SVN: run.in,v v. 1.1 (2015/05/11 19:37:40) wlyra initialize_mpicomm: enabled MPI SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: gravity_r.f90,v v. 1.1 (2018/08/24 15:48:10) wlyra 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: newton_cooling.f90 v. 1.2 (2015/12/11 02:56:37) wlyra SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 16 78 Lx, Ly, Lz= 2.1000000000000001 0.56000000000000005 6.2831853999999998 Vbox= 18.013349034979019 units_general: unit_velocity= 1306000.0000000000 units_general: unit_density= 1.9999999999999999E-011 units_general: unit_length= 77850000000000.000 units_general: unit_magnetic= 3.5449077018110318 units_eos: unit_temperature= 1.0000000000000000 units_eos: cp, lnTT0, cs0, pp0, Rgas= 7.2912257354367427E-005 5.1240244518003530 7.0000000000000007E-002 8.1900000000000011E-003 4.8747052059777068E-005 WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.! initialize_gravity: non-smoothed newtonian gravity diffusion: Dhyper/pi^4 *(Delta(rho))^6/Deltaq^2 initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using lnrho and ss initialize_energy: nheatc_max,iheatcond= 4 hyper3-sph nothing nothing nothing heat conduction: hyperdiffusivity of ss viscous force: nu_hyper3/pi^4 *(Deltav)^6/Deltaq^2 viscous force: nu_shock*(XXXXXXXXXXX) pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 93.204484805507050 WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_polar! duu_dt: SOLVE Bcs for ux, x: < out>, y: < s>, z: < p> Bcs for uy, x: < s>, y: < out>, z: < p> Bcs for uz, x: < a2>, y: < a2>, z: < p> duu_dt: max(advec_uu) = 38.284019983021700 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < str>, z: < p> dlnrho_dt: diffrho_hyper3= 1.0000000000000000E-003 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 4.9518350650971085E-002 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < a2>, y: < a2>, z: < p> denergy_dt: lnTT,cs2,cp1= 6.0800383749066729 1.2746403179318876E-002 0.0000000000000000E+000 calc_heatcond_hyper3_polar: chi_hyper3= 1.0000000000000000E-003 calc_heatcond_hyper3: added thdiff dspecial_dt: SOLVE dspecial_dt ---it-------t---------dt---------dtv------rhomin------rhom------rhomax-------TTm-------TTmin------TTmax-------uz2m-------taucoolm----taucoolmax---taucoolmin-- 0 0.0000 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.370E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 5 0.0050 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 10 0.0100 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 15 0.0150 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 20 0.0200 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 25 0.0250 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 30 0.0300 1.000E-03 1.198E-01 4.176E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 35 0.0350 1.000E-03 1.198E-01 4.177E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 40 0.0400 1.000E-03 1.198E-01 4.177E-06 1.129E-01 3.763E+00 7.902E+01 6.722E+01 4.371E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03 Simulation finished after 41 time-steps Writing final snapshot at time t = 4.1000000000000016E-002 Wall clock time [hours] = 7.878E-03 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 17.32053 Maximum used memory per cpu [MBytes] = 13.711 Maximum used memory [MBytes] = 109.031 real 31.01 user 46.68 sys 4.74