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$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: pointmasses.f90,v v. 1.1 (2019/02/02 03:54:41) wlyra SVN: particles_dust.f90 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$ The verbose level is ip= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 32 Lx, Ly, Lz= 0.128000006 0.128000006 0.128000006 Vbox= 2.09715241E-03 WARNING: input_persist_interstellar: t_next_SNI from run.in overwritten. Set l_persist_overwrite_tSNI=T to update! WARNING: input_persist_interstellar: t_next_SNII from run.in overwritten. Set l_persist_overwrite_tSNII=T to update! WARNING: input_persist_interstellar: x_cluster from run.in overwritten. Set l_persist_overwrite_xcluster=T to update! WARNING: input_persist_interstellar: y_cluster from run.in overwritten. Set l_persist_overwrite_ycluster=T to update! WARNING: input_persist_interstellar: z_cluster from run.in overwritten. Set l_persist_overwrite_zcluster=T to update! WARNING: input_persist_interstellar: t_cluster from run.in overwritten. Set l_persist_overwrite_tcluster=T to update! WARNING: input_persist_interstellar: lSNI from run.in overwritten. Set l_persist_overwrite_lSNI=T to update! input_persist_interstellar: lSNI = F t_next_SNI = 0.00000000E+00 WARNING: input_persist_interstellar: lSNII from run.in overwritten. Set l_persist_overwrite_lSNII=T to update! input_persist_interstellar: lSNII = T t_next_SNII = 0.00000000E+00 input_persist_forcing: 0.00000000E+00 0.00000000E+00 0.00000000E+00 input_persist_forcing: -10.0000000 units_general: unit_velocity= 97781.229948950713 units_general: unit_density= 1.6728000000000004E-024 units_general: unit_length= 3.0856775800000002E+021 units_general: unit_magnetic= 3.5449078083038330 unit temperature 1.0000000000000000 initialize_eos: ENTER initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: constant gravz= 0.00000000E+00 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! WARNING: initialize_energy: cool+cool2==0 -> no meaningful cs2top_from_cool can be formed. Set to cs2top! initialize_energy: nheatc_max,iheatcond= 4 nothing nothing nothing nothing heat conduction: nothing WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! resistivity: constant eta WARNING: initialize_magnetic: Resistivity coefficient eta is zero! lalpha_aniso= F lalpha_tensor= F leta_tensor= F initialize_magn_mf: lweyl_gauge= F calc_angle_weights: correction_factor = 0.333333343 initialize_radiation: ndir = 2 WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! chiral_diffXX, chiral_diff, chiral_diffZZ = 0.00000000E+00 0.00000000E+00 0.00000000E+00 initialize_dustvelocity: dust_chemistry = nothing recalculated: md0= 1.00000000 initialize_dustvelocity: ad= 0.620350480 initialize_dustvelocity: No betad calculation for draglaw=epstein_cst initialize_dustvelocity: dust geometry = sphere WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! initialize_dustdensity: ldustcoagulation,ldustcondensation = F F initialize_dustdensity: ldustcoagulation_simplified,ldustcondensation_simplified = F F advec_ddensity: plain vanilla scheme initialize_dustdensity: diffnd_ndustspec= 0.00000000E+00 dust diffusion: nothing gammacr1= 0.333333373 WARNING: initialize_cosmicrayflux: parameter tau was set to zero! initialize_interstellar: t_next_SNI = 0.00000000E+00 WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! initialize_interstellar: mu not set by eos: 0.620000005 initialize_interstellar: eampl_SN, kampl_SN = 0.000E+00 0.000E+00 select_cooling: WSW cooling fct initialize_interstellar: unit_Lambda 0.181E+18 initialize_interstellar: unit_Gamma 0.303E-06 WARNING: initialize_interstellar: using T/e instead of cv1 for diagnostics. Not yet implemented cv1 for ionization! initialize_interstellar: t_next_SNI, t_next_SNII= 0.00000000E+00 0.00000000E+00 initialize_shock: prenormalised shock_factor sum= 1.00000000 WARNING: initialize_pointmasses: the masses do not sum up to one! WARNING: initialize_boundcond: invalid value of gamma, set it to 5/3! initialize_particles: dust-to-gas ratio eps_dtog= 1.00000000 initialize_particles: mass per constituent particle mpmat= 0.00000000E+00 initialize_particles: mass per superparticle mp_swarm = 4.01060570E-05 initialize_particles: number density per superparticle np_swarm= 0.00000000E+00 initialize_particles: mass density per superparticle rhop_swarm= 626.657043 initialize particles: set minimum collisional time-scale equal to two times the Courant time-step. initialize_particles_map: selected cloud-in-cell for particle-mesh method. initialize_particles_radius: mass per dust grain mpmat= 0.00000000E+00 WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization. The values provided are for one-atomic ideal gas. Use at own risk! pencil_consistency_check: checking pencil case dxxp_dt: Calculate dxxp_dt dxxp_dt: Particles boundary condition bcpx=p dxxp_dt: Particles boundary condition bcpy=p dxxp_dt: Particles boundary condition bcpz=p dxxp_dt: Set rate of change of particle position equal to particle velocity. dvvp_dt: Calculate dvvp_dt particle_gravity: Calculating gravity particle_gravity: No gravity in x-direction. particle_gravity: No gravity in z-direction. particle_gravity: No radial gravity dvvp_dt_selfgrav: add self-gravity At line 273 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 273 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range' At line 1154 of file pointmasses.f90 Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr' At line 1112 of file pointmasses.f90 Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses' At line 1112 of file pointmasses.f90 Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses' pencil_consistency_check: checking requested pencils pencil_consistency_check: performing small pencil check pencil_consistency_check: the small pencil check reported no problems pencil_consistency_check: checking dependence on pencil initialization pencil_consistency_check: results are independent of pencil initialization pencil_consistency_check: checking diagnostics pencils pencil_consistency_check: performing small pencil check pencil_consistency_check: the small pencil check reported no problems pencil_consistency_check: set lpencil_check_diagnos_opti=T in run.in to report pencil_consistency_check: optimisation possibilities for diagnostics pencils pencil_consistency_check: checking dependence on pencil initialization pencil_consistency_check: diagnostics are independent of pencil initialization pencil_consistency_check: all tests passed pde: ENTER SVN: ------- v. ( ) $Id$ check_SN: ENTER explode_SN: cvelocity_SN is uu_sedov, velocity profile = gaussian check_SNII: Old t_next_SNII = 0.00000000E+00 check_SNII: Next SNII at time = 0.19171357E-04 WARNING: sn_series.dat: new column SN_rate added 19.10.17 continuation of old data may need header and extra column appended! WARNING: sn_series.dat: new columns rhom, Ekin, Ecr 20.03.19 continuation of old data may need header and extra columns appended! WARNING: sn_series.dat: new columns added_Nsol 27.05.21 continuation of old data may need header and extra columns appended! explode_SN: step, time = 1 0.000E+00 explode_SN: dVol = 0.640E-07 explode_SN: SN type = 2 explode_SN: proc, l, m, n = 0 12 34 20 explode_SN: x, y, z = -0.030 0.058 0.002 explode_SN:remnant radius = 0.020 explode_SN: rho, TT = 0.100E+01 0.100E+05 explode_SN: maximum TT = 0.454E-04 explode_SN: Mean density = 0.955E+00 explode_SN: Total energy = 0.000E+00 explode_SN: CR energy = 0.000E+00 explode_SN: Added mass = 0.000E+00 explode_SN: Ambient Nsol = 0.791E+03 explode_SN: Sedov time = 0.502E-04 explode_SN: Shell speed = 0.159E+03 calc_pencils_energy: max(advec_cs2) = 40506876.0 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < p> Bcs for uy, x: < she>, y: < p>, z: < p> Bcs for uz, x: < she>, y: < p>, z: < p> duu_dt: max(advec_uu) = 0.00000000E+00 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.00000000E+00 dlnrho_dt: max(diffus_diffrho3) = 0.00000000E+00 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < p> denergy_dt: lnTT,cs2,cp1= 9.21043491 216.036697 69.6056442 calc_heat_cool_interstellar: ENTER daa_dt: SOLVE Bcs for Ax, x: < she>, y: < p>, z: < p> Bcs for Ay, x: < she>, y: < p>, z: < p> Bcs for Az, x: < she>, y: < p>, z: < p> daa_dt: iresistivity=eta-const daa_dt_meanfield: SOLVE daa_dt_meanfield: max(diffus_eta) = 0.00000000E+00 daa_dt: max(diffus_eta) = 0.00000000E+00 daa_dt: max(diffus_eta2) = 0.00000000E+00 daa_dt: max(diffus_eta3) = 0.00000000E+00 daatest_dt: SOLVE Bcs for Axtest, x: < she>, y: < p>, z: < p> Bcs for Aztest, x: < she>, y: < p>, z: < p> SOLVE dlncc_dt Bcs for cc, x: < she>, y: < p>, z: < p> duud_dt: SOLVE duud_dt Bcs for udx, x: < she>, y: < p>, z: < p> Bcs for udy, x: < she>, y: < p>, z: < p> Bcs for udz, x: < she>, y: < p>, z: < p> dndmd_dt: SOLVE dnd_dt, dmd_dt, dmi_dt Bcs for nd, x: < she>, y: < p>, z: < p> SOLVE decr_dt Bcs for ecr, x: < she>, y: < p>, z: < p> decr_dt: no diffusion decr_dt: max(diffus_cr) = 0.00000000E+00 dfcr_dt: SOLVE Bcs for Fecx, x: < she>, y: < p>, z: < p> Bcs for Fecy, x: < she>, y: < p>, z: < p> Bcs for Fecz, x: < she>, y: < p>, z: < p> SOLVE dXY_dt Bcs for XX_chiral, x: < she>, y: < p>, z: < p> Bcs for YY_chiral, x: < she>, y: < p>, z: < p> dXY_chiral_dt: max(diffus_chiral) = 0.00000000E+00 shearing: Sshear,Sshear1= -0.00000000E+00 -0.00000000E+00 shearing: qshear,qshear0= 1.00000000 0.00000000E+00 dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dnpswarm_dt_pencil: Calculate dnpswarm_dt --it-------t--------dt-------urms----ekintot----ethtot-----rhom------ssm------TTmin-----TTmax------TTm------dtc-----dtu-----dtnu---dtchi-- 0 0.0000E+00 1.41E-04 0.000E+00 0.00E+00 2.452E+00 1.000E+00 7.546E-01 1.000E+04 1.00E+04 1.000E+04 1.0000 0.0000 0.0000 0.0000 explode_SN: cvelocity_SN is uu_sedov, velocity profile = gaussian check_SNII: Old t_next_SNII = 0.19171357E-04 check_SNII: Next SNII at time = 0.45175306E-03 explode_SN: step, time = 2 0.141E-03 explode_SN: dVol = 0.640E-07 explode_SN: SN type = 2 explode_SN: proc, l, m, n = 0 16 32 17 explode_SN: x, y, z = -0.014 0.050 -0.010 explode_SN:remnant radius = 0.020 explode_SN: rho, TT = 0.100E+01 0.100E+05 explode_SN: maximum TT = 0.454E-04 explode_SN: Mean density = 0.974E+00 explode_SN: Total energy = 0.000E+00 explode_SN: CR energy = 0.000E+00 explode_SN: Added mass = 0.000E+00 explode_SN: Ambient Nsol = 0.807E+03 explode_SN: Sedov time = 0.507E-04 explode_SN: Shell speed = 0.158E+03 1 1.4141E-04 1.43E-04 2.811E-02 8.28E-07 2.355E+00 1.000E+00 7.493E-01 7.628E+03 7.63E+03 7.628E+03 0.8819 0.1181 0.0000 0.0000 2 2.8421E-04 1.30E-04 5.656E-02 3.54E-06 2.319E+00 1.000E+00 7.469E-01 6.605E+03 7.18E+03 6.744E+03 0.7816 0.2184 0.0000 0.0000 3 4.1470E-04 1.22E-04 7.392E-02 6.63E-06 2.301E+00 1.000E+00 7.456E-01 5.977E+03 7.18E+03 6.288E+03 0.7312 0.2687 0.0000 0.0000 explode_SN: cvelocity_SN is uu_sedov, velocity profile = gaussian check_SNII: Old t_next_SNII = 0.45175306E-03 check_SNII: Next SNII at time = 0.59069874E-03 explode_SN: step, time = 5 0.537E-03 explode_SN: dVol = 0.640E-07 explode_SN: SN type = 2 explode_SN: proc, l, m, n = 0 25 24 18 explode_SN: x, y, z = 0.022 0.018 -0.006 explode_SN:remnant radius = 0.020 explode_SN: rho, TT = 0.100E+01 0.629E+04 explode_SN: maximum TT = 0.454E-04 explode_SN: Mean density = 0.963E+00 explode_SN: Total energy = 0.000E+00 explode_SN: CR energy = 0.000E+00 explode_SN: Added mass = 0.000E+00 explode_SN: Ambient Nsol = 0.797E+03 explode_SN: Sedov time = 0.504E-04 explode_SN: Shell speed = 0.159E+03 4 5.3674E-04 1.21E-04 8.179E-02 8.55E-06 2.289E+00 1.000E+00 7.447E-01 5.576E+03 7.03E+03 6.004E+03 0.7196 0.2812 0.0000 0.0000 Simulation finished after 5 time-steps Writing final snapshot at time t = 6.5811350214062259E-004 Wall clock time [hours] = 5.382E-04 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 11.5 Maximum used memory per cpu [MBytes] = 77.656 Maximum used memory [MBytes] = 77.656 real 7.40 user 7.17 sys 0.13