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
rsnap: read snapshot var.dat in 2.6125600000000002E-004 seconds
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
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
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)
sh: 1: Syntax error: Bad fd number
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_hydro: call gij_etc
calc_pencils_energy: max(advec_cs2) = 93.204484805507164
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.284019983021686
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.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 4.9518350650971099E-002
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < a2>, y: < a2>, z: < p>
denergy_dt: lnTT,cs2,cp1= 6.0800383749066729 1.2746403179318891E-002 0.0000000000000000
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.370E+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.370E+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.370E+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.370E+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.370E+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.370E+02 8.71197E-01 8.06383E+00 7.69498E+01 1.00000E-03
35 0.0350 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
40 0.0400 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
Simulation finished after 41 time-steps
Writing final snapshot at time t = 4.1000000000000016E-002
Wall clock time [hours] = 5.142E-04 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 1.130599
Wall clock time/timestep/local meshpoint [microsec] = 9.044788
Rhs wall clock time/timestep/local meshpoint [microsec] = 6.207180
Maximum used memory per cpu [MBytes] = 29.078
Maximum used memory [MBytes] = 230.645
real 2.20
user 14.90
sys 0.48