SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! WARNING: find_namelist: namelist "particles_breakup_init_pars" is missing!! SVN: start.in v. 43633 (2026-04-17 07:49:31) AxelBrandenburg WARNING: find_namelist: namelist "_pars" is missing!! SVN: run.in v. 43631 (2026-04-17 07:01:21) nilshau WARNING: run: quad precision suppressed! 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: 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$ The verbose level is ip= 14 (ldebug= F ) This is a 0-D run nxgrid, nygrid, nzgrid= 1 1 1 Lx, Ly, Lz= 6.2831853071795862 6.2831853071795862 6.2831853071795862 Vbox= 1.0000000000000000 rsnap: read snapshot var.dat in 0.0000000000000000 seconds setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449077018110318 units_eos: unit_temperature= 1.2027220936797545E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_density: 0-D run, turned off continuity equation diffusion: nothing (i.e. no mass diffusion) WARNING: initialize_density: are rho_const or lnrho_const ok?! inconsistency between the density constants from eos (rho0 or lnrho0) and the ones from the density module (rho_const or lnrho_const). It may damage your simulation if you are using them in different places. initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_hydro: 0-D run, turned off advection of velocity. select_eos_variable: Using rho and cs2 viscous force: nu*(del2u+graddivu/3+2S.glnrho) initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000 initialize_particles: mass per constituent particle mpmat= 4188.7902047863909 initialize_particles: mass per superparticle mp_swarm = 1.0000000000000000E-004 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 1.0000000000000000E-004 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= 4188.7902047863909 pde: ENTER SVN: ------- v. ( ) $Id$ At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' duu_dt: SOLVE Bcs for ux, x: < p>, y: < p>, z: < p> Bcs for uy, x: < p>, y: < p>, z: < p> Bcs for uz, x: < p>, y: < p>, z: < p> dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> At line 2314 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 2314 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range' dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt: Add drag force; tausp= 0.0000000000000000 dap_dt_pencil: Calculate dap/dt dnpswarm_dt_pencil: Calculate dnpswarm_dt 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: Constant gravity field in z-direction. particle_gravity: No radial gravity At line 221 of file particles_breakup.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 221 of file particles_breakup.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range' At line 222 of file particles_breakup.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'get_local_gas_density' At line 222 of file particles_breakup.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'get_local_gas_density' ---it--------t---------dt-------vpzmin-----vpzmax------vpzm--------npm-------np2m-------npmax-------apm-------ap2m-------apmin------apmax-----imskhm----npswarmm------zpm-----npar_found- 0 0.000E+00 5.000E-04 0.000E+00 0.000E+00 0.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 20 1.000E-02 5.000E-04 -9.810E+00 -9.810E+00 -9.810E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 40 2.000E-02 5.000E-04 -1.962E+01 -1.962E+01 -1.962E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 60 3.000E-02 5.000E-04 -2.943E+01 -2.943E+01 -2.943E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 80 4.000E-02 5.000E-04 -3.924E+01 -3.924E+01 -3.924E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 100 5.000E-02 5.000E-04 -4.905E+01 -4.905E+01 -4.905E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 120 6.000E-02 5.000E-04 -5.886E+01 -5.886E+01 -5.886E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 140 7.000E-02 5.000E-04 -6.867E+01 -6.867E+01 -6.867E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 160 8.000E-02 5.000E-04 -7.848E+01 -7.848E+01 -7.848E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 180 9.000E-02 5.000E-04 -8.829E+01 -8.829E+01 -8.829E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 200 1.000E-01 5.000E-04 -9.809E+01 -9.809E+01 -9.809E+01 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 220 1.100E-01 5.000E-04 -1.079E+02 -1.079E+02 -1.079E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 240 1.200E-01 5.000E-04 -1.177E+02 -1.177E+02 -1.177E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 260 1.300E-01 5.000E-04 -1.275E+02 -1.275E+02 -1.275E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 280 1.400E-01 5.000E-04 -1.373E+02 -1.373E+02 -1.373E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 300 1.500E-01 5.000E-04 -1.471E+02 -1.471E+02 -1.471E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 320 1.600E-01 5.000E-04 -1.569E+02 -1.569E+02 -1.569E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 340 1.700E-01 5.000E-04 -1.667E+02 -1.667E+02 -1.667E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 360 1.800E-01 5.000E-04 -1.765E+02 -1.765E+02 -1.765E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 380 1.900E-01 5.000E-04 -1.863E+02 -1.863E+02 -1.863E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 400 2.000E-01 5.000E-04 -1.961E+02 -1.961E+02 -1.961E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 420 2.100E-01 5.000E-04 -2.059E+02 -2.059E+02 -2.059E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 440 2.200E-01 5.000E-04 -2.158E+02 -2.158E+02 -2.158E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 460 2.300E-01 5.000E-04 -2.256E+02 -2.256E+02 -2.256E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 480 2.400E-01 5.000E-04 -2.354E+02 -2.354E+02 -2.354E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 500 2.500E-01 5.000E-04 -2.451E+02 -2.451E+02 -2.451E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 520 2.600E-01 5.000E-04 -2.549E+02 -2.549E+02 -2.549E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 540 2.700E-01 5.000E-04 -2.647E+02 -2.647E+02 -2.647E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 560 2.800E-01 5.000E-04 -2.745E+02 -2.745E+02 -2.745E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 580 2.900E-01 5.000E-04 -2.843E+02 -2.843E+02 -2.843E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+01 1.000E+02 1.000E+01 1.000E+01 0.000E+00 1.000E+00 0.000E+00 1 600 3.000E-01 5.000E-04 -2.941E+02 -2.941E+02 -2.941E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 620 3.100E-01 5.000E-04 -3.038E+02 -3.038E+02 -3.038E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 640 3.200E-01 5.000E-04 -3.135E+02 -3.135E+02 -3.135E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 660 3.300E-01 5.000E-04 -3.233E+02 -3.233E+02 -3.233E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 680 3.400E-01 5.000E-04 -3.330E+02 -3.330E+02 -3.330E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 700 3.500E-01 5.000E-04 -3.427E+02 -3.427E+02 -3.427E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 720 3.600E-01 5.000E-04 -3.524E+02 -3.524E+02 -3.524E+02 1.000E+00 1.000E+00 1.000E+00 2.240E+00 5.016E+00 2.240E+00 2.240E+00 0.000E+00 8.900E+01 0.000E+00 1 740 3.700E-01 5.000E-04 -3.619E+02 -3.619E+02 -3.619E+02 1.000E+00 1.000E+00 1.000E+00 6.002E-01 3.603E-01 6.002E-01 6.002E-01 0.000E+00 4.624E+03 0.000E+00 1 760 3.800E-01 5.000E-04 -3.713E+02 -3.713E+02 -3.713E+02 1.000E+00 1.000E+00 1.000E+00 6.002E-01 3.603E-01 6.002E-01 6.002E-01 0.000E+00 4.624E+03 0.000E+00 1 780 3.900E-01 5.000E-04 -3.804E+02 -3.804E+02 -3.804E+02 1.000E+00 1.000E+00 1.000E+00 1.963E-01 3.853E-02 1.963E-01 1.963E-01 0.000E+00 1.322E+05 0.000E+00 1 800 4.000E-01 5.000E-04 -3.881E+02 -3.881E+02 -3.881E+02 1.000E+00 1.000E+00 1.000E+00 7.679E-02 5.896E-03 7.679E-02 7.679E-02 0.000E+00 2.209E+06 0.000E+00 1 820 4.100E-01 5.000E-04 -1.664E+00 -1.664E+00 -1.664E+00 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 840 4.200E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 860 4.300E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 880 4.400E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 900 4.500E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 920 4.600E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 940 4.700E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 960 4.800E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 980 4.900E-01 5.000E-04 -4.556E-01 -4.556E-01 -4.556E-01 1.000E+00 1.000E+00 1.000E+00 6.197E-04 3.840E-07 6.197E-04 6.197E-04 0.000E+00 4.202E+12 0.000E+00 1 Simulation finished after 1000 time-steps Writing final snapshot at time t = 0.49999999999997269 Wall clock time [hours] = 3.139E-05 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.130E-02 Maximum used memory per cpu [MBytes] = 10.125 Maximum used memory [MBytes] = 10.125 real 0.14 user 0.14 sys 0.00