Pencil1: iapp, nprocs, ncpus= 0 8 8 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2020/02/23 20:50:00) xiangyu SVN: run.in,v v. 1.1 (2020/02/23 20:50:00) 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: 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 3-D run nxgrid, nygrid, nzgrid= 16 16 16 Lx, Ly, Lz= 5.0000000000000000 5.0000000000000000 5.0000000000000000 Vbox= 125.00000000000000 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= 100000.00000000000 units_general: unit_density= 6.7710000000000014E-023 units_general: unit_length= 3.0859999999999995E+018 units_general: unit_magnetic= 3.5449077018110318 units_eos: unit_temperature= 240.54441873595090 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -2.1517456033972402 0.34100000000000003 11.472283460000003 2.0000000000000000 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 cs2 initialize_shock: prenormalised shock_factor sum= 1.0000000000000000 initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002 initialize_particles: mass per constituent particle mpmat= 1.4849261275967754E-046 initialize_particles: mass per superparticle mp_swarm = 2.5257024657983126E-003 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 initialize_particles: mass density per superparticle rhop_swarm= 8.2762218399279094E-002 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= 1.4849261275967754E-046 pencil_consistency_check: checking pencil case 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' 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' 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 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 rho ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rho1 ( 8) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil TT1 ( 53) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ppvap ( 68) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil csvap2 ( 69) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil np ( 166) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rhop ( 167) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cc ( 181) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cc1 ( 182) 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$ calc_pencils_hydro: upwinding advection term calc_pencils_energy: max(advec_cs2) = 3.5721523200000012 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 ---it--------t---------dt-------urms-------umax-------epsK-------rhom-------rho2m-----rhomin-----rhomax----shockmax-----vprms------vpmax-------apm-------ap2m-------apmin------apmax--- 0 0.000E+00 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 1.000E-23 1.000E-46 1.000E-23 1.000E-23 20 9.524E+00 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 3.020E-23 9.121E-46 3.020E-23 3.020E-23 40 1.905E+01 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 5.286E-23 2.794E-45 5.286E-23 5.286E-23 60 2.857E+01 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 7.551E-23 5.702E-45 7.551E-23 7.551E-23 80 3.809E+01 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 9.817E-23 9.637E-45 9.817E-23 9.817E-23 Simulation finished after 81 time-steps Writing final snapshot at time t = 38.571146099636820 Wall clock time [hours] = 2.204E-03 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 3.44 Maximum used memory per cpu [MBytes] = 12.520 Maximum used memory [MBytes] = 98.707 real 26.44 user 37.03 sys 6.93