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: particles_dust.f90 v. 1.1 (2018/08/24 15:48:10) wlyra SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 14 (ldebug= F ) This is a 1-D run nxgrid, nygrid, nzgrid= 64 1 1 Lx, Ly, Lz= 6.2830000000000004 6.2830000000000004 6.2830000000000004 Vbox= 6.2830000000000004 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.0000000000000000E+000 1.0000000000000000 1.0000000000000000 1.0000000000000000 select_eos_variable: Using lnrho and cs2 initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 1.6361979166666669E-002 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 initialize_particles: mass density per superparticle rhop_swarm= 0.16666666666666666 initialize particles: set minimum collisional time-scale equal to two times the Courant time-step. initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method. 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 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' 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: performing full pencil check (takes a while) pencil_consistency_check: possible overcalculation... pencil gpotself ( 186) is requested, but does not appear to be required! 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: performing full pencil check (takes a while) 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$ duu_dt: diagnostics ... dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt: Add drag force; tausp= 0.10000000000000001 ---it------t-------vpx2m------rhopmax---- 0 0.00 4.903E-09 1.000999E+00 5 0.10 4.985E-09 1.001009E+00 10 0.20 5.079E-09 1.001019E+00 15 0.30 5.178E-09 1.001029E+00 20 0.40 5.281E-09 1.001039E+00 25 0.50 5.386E-09 1.001049E+00 30 0.60 5.493E-09 1.001060E+00 35 0.70 5.603E-09 1.001070E+00 40 0.80 5.715E-09 1.001081E+00 45 0.90 5.829E-09 1.001092E+00 50 1.00 5.945E-09 1.001102E+00 55 1.10 6.064E-09 1.001113E+00 60 1.20 6.185E-09 1.001124E+00 65 1.30 6.308E-09 1.001136E+00 70 1.40 6.434E-09 1.001147E+00 75 1.50 6.562E-09 1.001158E+00 80 1.60 6.693E-09 1.001170E+00 85 1.70 6.827E-09 1.001181E+00 90 1.80 6.963E-09 1.001193E+00 95 1.90 7.102E-09 1.001205E+00 Simulation finished after 96 time-steps Writing final snapshot at time t = 1.9199999999999910 Wall clock time [hours] = 3.556E-05 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 2.98 Maximum used memory per cpu [MBytes] = 4.391 Maximum used memory [MBytes] = 4.391 real 0.66 user 0.64 sys 0.00