Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "special_init_pars" is missing!! SVN: ------- v. ( ) no cvsid is given i SVN: ------- v. ( ) no cvsid is given i 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: particles_density. v. 20849 (2013-08-06 18:45:43) anders SVN: ------- v. ( ) $Id$ The verbose level is ip= 14 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 64 1 64 Lx, Ly, Lz= 6.2831853071795866E-003 1.0000000000000000 6.2831853071795866E-003 Vbox= 3.9478417604357436E-005 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= 4.8108883747190188E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 4.0812192409268553 6.2831853071795862 23.687050562614463 0.40000000000000008 diffusion: mesh hyperdiffusion initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_density: Global density gradient with beta_glnrho_global= -0.10000000000000001 0.0000000000000000E+000 0.0000000000000000E+000 select_eos_variable: Using rho and ss viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq initialize_special: ux0, uy0 = 9.2791620582046985E-003 -0.12583302733752555 initialize_special: vpx0 = -1.4140910635015245E-003 -1.5730141989406348E-002 initialize_special: vpy0 = -0.12580297790242614 -0.12504652023805524 initialize_particles: Number of particle species = 2 initialize_particles: tausp_species = 6.7640850814055528E-003 1.5915494309189534E-002 initialize_particles: dust-to-gas ratio eps_dtog= 1.5000000000000000 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 1.8071785402385300E-009 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 initialize_particles: mass density per superparticle rhop_swarm= 0.18750000000000000 initialize particles: set minimum collisional time-scale equal to two times the Courant time-step. initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method. 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' drhopswarm_dt_pencil: Calculate drhopswarm_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 dvvp_dt: Add Coriolis force; Omega= 6.2831853071795862 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: 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_energy: max(advec_cs2) = 8192000000.0376654 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < p> Bcs for uy, x: < she>, y: < p>, z: < p> Bcs for uz, x: < she>, y: < p>, z: < p> coriolis_cartesian: add Coriolis force; Omega= 6.2831853071795862 duu_dt: max(advec_uu) = 94.517383759778284 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < p> dlnrho_dt: diffrho_hyper3_mesh= 1.4577334289680860E-003 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 29.696701559112853 denergy_dt: adding global pressure gradient force dspecial_dt: special diagnostics shearing: Sshear,Sshear1= -9.4247779607693793 -9.4247779607693793 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000E+000 dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt_pencil: adding feedback from dust to gas dvvp_dt: Add drag force; tausp= 0.0000000000000000E+000 ----it-------t-------dt----dtv----rhom-----drho2m-----urms------umax-------uxm--------uym------rdux2m----rduy2m----ruz2m----rduxduym---nparmax---nmigmmax----vpxm-------vpym----rhopdvpx2m-rhopdvpy2m-rhopvpz2m---npmax----drhop2m---rhopmax-- 0 0.00000 9.9E-06 1.00 1.000E+00 3.757E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.283E-16 8.028E-16 9.283E-16 -7.108E-16 8.192E+03 0.000E+00 -8.572E-03 -1.254E-01 1.395E-16 1.487E-16 1.639E-16 8.000E+00 4.967E-13 1.500E+00 100 0.00099 9.9E-06 1.00 1.000E+00 3.727E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.325E-16 8.064E-16 9.325E-16 -7.141E-16 8.192E+03 0.000E+00 -8.572E-03 -1.254E-01 1.399E-16 1.492E-16 1.644E-16 8.000E+00 4.986E-13 1.500E+00 200 0.00199 9.9E-06 1.00 1.000E+00 3.766E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.365E-16 8.099E-16 9.365E-16 -7.171E-16 8.192E+03 6.400E+01 -8.572E-03 -1.254E-01 1.404E-16 1.497E-16 1.650E-16 8.000E+00 5.005E-13 1.500E+00 300 0.00298 9.9E-06 1.00 1.000E+00 3.778E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.404E-16 8.133E-16 9.404E-16 -7.201E-16 8.192E+03 0.000E+00 -8.572E-03 -1.254E-01 1.409E-16 1.502E-16 1.656E-16 8.000E+00 5.024E-13 1.500E+00 Simulation finished after 301 time-steps Writing final snapshot at time t = 2.9886280746550952E-003 Wall clock time [hours] = 3.106E-02 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 10.1 Maximum used memory per cpu [MBytes] = 12.199 Maximum used memory [MBytes] = 47.922 real 115.25 user 161.64 sys 23.11