Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2016/04/06 21:39:26) xiangyu SVN: run.in,v v. 1.1 (2016/04/06 21:39:26) xiangyu initialize_mpicomm: enabled MPI 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 3-D run nxgrid, nygrid, nzgrid= 16 16 16 Lx, Ly, Lz= 0.50000000000000000 0.50000000000000000 0.50000000000000000 Vbox= 0.12500000000000000 input_persist_forcing: 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 input_persist_forcing: -10.000000000000000 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.0000000000000000 units_eos: cp, lnTT0, cs0, pp0, Rgas= 9700218.1368201599 -15.171368199678572 1.0000000000000000 21.428571428571445 83144726.887029916 initialize_gravity: constant x-grav= -9.8100000000000005 initialize_gravity: no y-gravity initialize_gravity: no z-gravity diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using lnrho and ss viscous force: nu*(del2u+graddivu/3+2S.glnrho) initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 1.2499999999999999E-007 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 initialize_particles: mass density per superparticle rhop_swarm= 4.0959999999999998E-003 initialize particles: set minimum collisional time-scale equal to two times the Courant time-step. initialize_particles_map: selected nearest-grid-point for particle-mesh method. initialize_particles_radius: mass per dust grain mpmat= 0.0000000000000000E+000 At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_energy: max(advec_cs2) = 3072.0000000000000 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> duu_dt: max(advec_uu) = 0.0000000000000000E+000 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt: Add drag force; tausp= 0.0000000000000000E+000 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: Constant gravity field in x-direction particle_gravity: No gravity in z-direction. particle_gravity: No radial gravity ---it---------t---------dt--------dvpm-------dvpmax-----ncoagpm-------apm--------apmin-------apmax-------admom0------admom1------admom3------admom24--- 0 0.00000 8.935E-04 0.0000E+00 0.0000E+00 0.0000E+00 1.0216E-05 4.4413E-06 2.5018E-05 7.6335E+08 7.7987E+03 9.1811E-07 -1.055E+02 50 0.04468 8.935E-04 1.0693E-01 7.4295E-01 0.0000E+00 1.0218E-05 4.4413E-06 2.5018E-05 7.6313E+08 7.7966E+03 9.1811E-07 -1.055E+02 100 0.08935 8.935E-04 1.0709E-01 7.4296E-01 0.0000E+00 1.0220E-05 4.4413E-06 2.5018E-05 7.6297E+08 7.7955E+03 9.1811E-07 -1.055E+02 150 0.13403 8.935E-04 1.0741E-01 7.4296E-01 0.0000E+00 1.0220E-05 4.4413E-06 2.5018E-05 7.6290E+08 7.7948E+03 9.1811E-07 -1.055E+02 200 0.17871 8.935E-04 1.0809E-01 7.5270E-01 0.0000E+00 1.0223E-05 4.4413E-06 2.5018E-05 7.6259E+08 7.7928E+03 9.1811E-07 -1.055E+02 250 0.22338 8.935E-04 1.0864E-01 7.7759E-01 0.0000E+00 1.0224E-05 4.4413E-06 2.5018E-05 7.6244E+08 7.7914E+03 9.1811E-07 -1.055E+02 300 0.26806 8.935E-04 1.0861E-01 7.7759E-01 0.0000E+00 1.0225E-05 4.4413E-06 2.5018E-05 7.6236E+08 7.7908E+03 9.1811E-07 -1.055E+02 350 0.31274 8.935E-04 1.0901E-01 7.7759E-01 0.0000E+00 1.0227E-05 4.4413E-06 2.5018E-05 7.6206E+08 7.7885E+03 9.1811E-07 -1.055E+02 400 0.35741 8.935E-04 1.0919E-01 7.7759E-01 0.0000E+00 1.0228E-05 4.4413E-06 2.5018E-05 7.6190E+08 7.7875E+03 9.1811E-07 -1.055E+02 450 0.40209 8.935E-04 1.0948E-01 7.7759E-01 1.0000E-04 1.0230E-05 4.4413E-06 2.5018E-05 7.6168E+08 7.7857E+03 9.1811E-07 -1.055E+02 Simulation finished after 451 time-steps Writing final snapshot at time t = 0.40298338706714787 Wall clock time [hours] = 1.552E-03 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 0.879 Maximum used memory per cpu [MBytes] = 27.648 Maximum used memory [MBytes] = 110.055 real 6.07 user 22.46 sys 0.36