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 (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 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 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.0000000000000000 0.0000000000000000 0.0000000000000000 input_persist_forcing: -10.000000000000000 rsnap: read snapshot var.dat in 1.7359999999999999E-004 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.0000000000000000 units_eos: cp, lnTT0, cs0, pp0, Rgas= 9700218.1368201599 -15.171368199678572 1.0000000000000000 21.428571428571445 83144726.887029916 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 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.0000000000000000 initialize_particles: mass per superparticle mp_swarm = 1.2499999999999999E-007 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 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.0000000000000000 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' 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' calc_pencils_hydro: max(advec_uu) = 0.0000000000000000 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> dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 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: 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.0689E-01 7.4295E-01 0.0000E+00 1.0217E-05 4.4413E-06 2.5018E-05 7.6328E+08 7.7981E+03 9.1811E-07 -1.055E+02 100 0.08935 8.935E-04 1.0696E-01 7.4296E-01 0.0000E+00 1.0218E-05 4.4413E-06 2.5018E-05 7.6320E+08 7.7976E+03 9.1811E-07 -1.055E+02 150 0.13403 8.935E-04 1.0727E-01 7.4296E-01 0.0000E+00 1.0218E-05 4.4413E-06 2.5018E-05 7.6313E+08 7.7970E+03 9.1811E-07 -1.055E+02 200 0.17871 8.935E-04 1.0810E-01 7.5270E-01 0.0000E+00 1.0220E-05 4.4413E-06 2.5018E-05 7.6290E+08 7.7953E+03 9.1811E-07 -1.055E+02 250 0.22338 8.935E-04 1.0850E-01 7.7759E-01 0.0000E+00 1.0221E-05 4.4413E-06 2.5018E-05 7.6274E+08 7.7942E+03 9.1811E-07 -1.055E+02 300 0.26806 8.935E-04 1.0841E-01 7.7759E-01 0.0000E+00 1.0223E-05 4.4413E-06 2.5018E-05 7.6259E+08 7.7931E+03 9.1811E-07 -1.055E+02 350 0.31274 8.935E-04 1.0874E-01 7.7759E-01 0.0000E+00 1.0223E-05 4.4413E-06 2.5018E-05 7.6259E+08 7.7931E+03 9.1811E-07 -1.055E+02 400 0.35741 8.935E-04 1.0905E-01 7.7759E-01 0.0000E+00 1.0223E-05 4.4413E-06 2.5018E-05 7.6259E+08 7.7931E+03 9.1811E-07 -1.055E+02 450 0.40209 8.935E-04 1.0953E-01 7.7759E-01 0.0000E+00 1.0226E-05 4.4413E-06 2.5018E-05 7.6229E+08 7.7902E+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] = 8.177E-04 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 0.463 Maximum used memory per cpu [MBytes] = 36.293 Maximum used memory [MBytes] = 144.273 real 3.32 user 11.99 sys 0.15