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 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 (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 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.12500000000000000 0.12500000000000000 0.12500000000000000 Vbox= 1.9531250000000000E-003 input_persist_forcing: 0.0000000000000000 0.0000000000000000 0.0000000000000000 input_persist_forcing: -10.000000000000000 rsnap: read snapshot var.dat in 2.1089899999999999E-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.0309046517255039E-003 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 4.1351665567423561 5.0000000000000000 567.85714285714312 8.5714285714285712 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 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.0000000000000000 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_hydro: max(advec_uu) = 0.0000000000000000 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> 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 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.0000000000000000 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.2417E+09 5.2438927975E+04 5.2461E-01 5.2483E-06 5.2505E-11 5.2549E-21 -5.0278E+01 -1.1028E+02 1.5659770557E-02 1.5659766534E-02 2.9210008528E+02 2.9210008325E+02 8.7807976524E-04 7.9422123239E-04 2.6063732292E-02 2.1983896834E-02 1.2237194374E-02 1.1908888823E-02 1.9148334730E+00 1.9098003180E+00 Simulation finished after 81 time-steps Writing final snapshot at time t = 3.6303644278855345E-002 Wall clock time [hours] = 1.608E-04 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 0.507 Maximum used memory per cpu [MBytes] = 39.816 Maximum used memory [MBytes] = 314.359 real 1.03 user 5.62 sys 0.28