Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified 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 WARNING: run: quad precision suppressed! 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 5.6373000000000003E-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] = 7.708E-04 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 1.694777 Wall clock time/timestep/local meshpoint [microsec] = 13.55822 Rhs wall clock time/timestep/local meshpoint [microsec] = 9.219623 Maximum used memory per cpu [MBytes] = 29.574 Maximum used memory [MBytes] = 234.227 real 3.39 user 22.97 sys 1.04