Pencil1: iapp, nprocs, ncpus= 0 8 8 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2018/06/15 13:24:13) xiangyu SVN: run.in,v v. 1.1 (2018/06/15 13:24:13) 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.12500000000000000 0.12500000000000000 0.12500000000000000 Vbox= 1.9531250000000000E-003 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.0309046517255039E-003 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 4.1351665567423561 5.0000000000000000 567.85714285714312 8.5714285714285712 initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: constant gravz= -9.8100000000000005 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. select_eos_variable: Using lnrho and ss Supersaturation routine 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= 4.1887902047863918E-012 initialize_particles: mass per superparticle mp_swarm = 3.9166066791109380E-008 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 initialize_particles: mass density per superparticle rhop_swarm= 8.2137195303108618E-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= 4.1887902047863918E-012 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: Constant gravity field 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: 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: 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) = 1200490.8564055739 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 SOLVE dacc_dt Bcs for acc, x: < p>, y: < p>, z: < p> dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt: Add drag force; tausp= 0.0000000000000000E+000 dnpswarm_dt_pencil: Calculate dnpswarm_dt ---it---------t---------dt-------urms-------epsK------ncoagpm------admom0---------admom1---------admom2------admom3------admom4------admom6-----admom12-----admom24---------accrms-------------accm-------------ttcrms-------------ttcm-------condensationRaterms-condensationRatem-waterMixingRatiorms-waterMixingRatiom------ssatrms------------ssatm----------buoyancyrms--------buoyancym----- 0 0.00000 4.482E-04 0.000E+00 0.000E+00 0.0000E+00 5.2417E+09 5.2416992188E+04 5.2417E-01 5.2417E-06 5.2417E-11 5.2417E-21 -5.0281E+01 -1.1028E+02 1.5700000000E-02 1.5700000000E-02 2.9200000000E+02 2.9200000000E+02 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00 2.1318522625E+00 2.1318522625E+00 20 0.00896 4.482E-04 2.632E-02 2.633E-03 0.0000E+00 5.2417E+09 5.2423659803E+04 5.2430E-01 5.2437E-06 5.2444E-11 5.2457E-21 -5.0280E+01 -1.1028E+02 1.5687173832E-02 1.5687172474E-02 2.9203190978E+02 2.9203190927E+02 1.4947184855E-03 1.3053381812E-03 2.6031782488E-02 2.1964661151E-02 1.7912291155E-02 1.7852253298E-02 1.9193181683E+00 1.9143473192E+00 40 0.01793 4.482E-04 3.556E-02 4.718E-03 0.0000E+00 5.2417E+09 5.2429451422E+04 5.2442E-01 5.2454E-06 5.2467E-11 5.2492E-21 -5.0279E+01 -1.1028E+02 1.5676472726E-02 1.5676469743E-02 2.9205853423E+02 2.9205853298E+02 1.2236237909E-03 1.0907012978E-03 2.6032209020E-02 2.1971955907E-02 1.5691913989E-02 1.5527180830E-02 1.9175721055E+00 1.9125763325E+00 60 0.02689 4.482E-04 4.361E-02 7.044E-03 0.0000E+00 5.2417E+09 5.2434505188E+04 5.2452E-01 5.2470E-06 5.2487E-11 5.2522E-21 -5.0279E+01 -1.1028E+02 1.5667459367E-02 1.5667455546E-02 2.9208095812E+02 2.9208095635E+02 1.0271079748E-03 9.2539364395E-04 2.6019962391E-02 2.1978323269E-02 1.3831027814E-02 1.3573011293E-02 1.9160628557E+00 1.9110799414E+00 80 0.03586 4.481E-04 4.989E-02 9.103E-03 0.0000E+00 5.2412E+09 5.2435048704E+04 5.2459E-01 5.2483E-06 5.2508E-11 5.2559E-21 -5.0278E+01 -1.1027E+02 1.5659770673E-02 1.5659766650E-02 2.9210008499E+02 2.9210008296E+02 8.7800730150E-04 7.9415605942E-04 2.6063732293E-02 2.1983896835E-02 1.2237216914E-02 1.1908912266E-02 1.9148334910E+00 1.9098003365E+00 Simulation finished after 81 time-steps Writing final snapshot at time t = 3.6303644278855345E-002 Wall clock time [hours] = 6.810E-04 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 2.15 Maximum used memory per cpu [MBytes] = 34.371 Maximum used memory [MBytes] = 273.246 real 3.50 user 10.16 sys 2.35