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 Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ 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$ The verbose level is ip= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 32 Lx, Ly, Lz= 0.20000000000000001 0.20000000000000001 0.20000000000000001 Vbox= 8.0000000000000019E-003 rsnap: read snapshot var.dat in 4.0685899999999998E-004 seconds setup_slices: slice_position = m setup_slices: iz_loc,iz2_loc, (video files) = 19 35 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.8108898179852413E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 0.40546460810828949 1.0000000000000000 0.59999988000002402 0.40000011999997598 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: linear z-grav, nu= 1.0000000000000000 diffusion: (d^6/dx^6+d^6/dy^6+d^6/dz^6)rho initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. initialize_density: Global density gradient with beta_glnrho_global= -0.10000000000000001 0.0000000000000000 0.0000000000000000 select_eos_variable: Using rho and ss initialize_energy: cs2top,cs2cool,cs2cool2= 1.0000000000000000 0.0000000000000000 0.0000000000000000 initialize_energy: now set cs2cool=cs2top initialize_energy: nheatc_max,iheatcond= 4 nothing nothing nothing nothing heat conduction: nothing viscous force(i): nu_hyper/rho*(del6ui+der5(divu,i)) initialize_particles: Number of particle species = 4 initialize_particles: tausp_species = 0.10000000000000001 0.20000000000000001 0.29999999999999999 0.40000000000000002 initialize_particles: dust-to-gas ratio eps_dtog= 2.0000000000000000E-002 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000 initialize_particles: mass per superparticle mp_swarm = 2.0053026197048009E-007 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 0.82137195303108623 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. sh: 1: Syntax error: Bad fd number 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' 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= 1.0000000000000000 particle_gravity: Calculating gravity particle_gravity: No gravity in x-direction. particle_gravity: Linear 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_energy: max(advec_cs2) = 75224.999999999985 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < s> Bcs for uy, x: < she>, y: < p>, z: < s> Bcs for uz, x: < she>, y: < p>, z: < a> coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000 duu_dt: max(advec_uu) = 7.9999999999999991 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < s> dlnrho_dt: diffrho_hyper3= 4.0000000000000003E-015 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.18968670806249990 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < s> denergy_dt: lnTT,cs2,cp1= 0.40546460810828949 1.0000000000000000 1.1728390000000000 denergy_dt: adding global pressure gradient force shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000 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.50000000000000000 ---it-------t--------dt------nmigmax----nparmax----dtdragp-----ux2m-------uy2m-------uxuym------rhom------rhomin-----rhomax-------TTm-------TTmax------vpxm--------xpm-------xp2m--------zpm-------zp2m-------npmax------rhopm-----rhopmax----dedragp----decollp-- 0 0.000 1.309E-03 0.000E+00 2.500E+03 3.969E-01 2.461E-05 2.116E-03 -6.798E-05 9.971E-01 9.917E-01 1.000E+00 1.500E+00 1.500E+00 -8.526E-03 -2.131E-04 3.346E-03 1.111E-05 9.975E-05 1.000E+01 2.507E-01 4.062E+00 5.367E-04 0.000E+00 5 0.007 1.305E-03 1.000E+00 2.503E+03 3.773E-01 2.281E-05 2.112E-03 -6.632E-05 9.971E-01 9.909E-01 1.010E+00 1.500E+00 1.510E+00 -8.146E-03 -2.076E-04 3.346E-03 1.111E-05 9.974E-05 9.000E+00 2.507E-01 4.105E+00 4.835E-04 0.000E+00 10 0.013 1.302E-03 2.000E+00 2.505E+03 3.776E-01 2.152E-05 2.111E-03 -6.434E-05 9.971E-01 9.919E-01 1.007E+00 1.500E+00 1.507E+00 -7.805E-03 -1.796E-04 3.346E-03 1.112E-05 9.973E-05 9.000E+00 2.507E-01 4.119E+00 4.400E-04 0.000E+00 15 0.020 1.305E-03 1.000E+00 2.505E+03 3.786E-01 1.973E-05 2.111E-03 -6.199E-05 9.971E-01 9.919E-01 1.007E+00 1.500E+00 1.507E+00 -7.497E-03 -1.694E-04 3.346E-03 1.112E-05 9.971E-05 9.000E+00 2.507E-01 4.108E+00 3.946E-04 0.000E+00 Simulation finished after 16 time-steps Writing final snapshot at time t = 2.0873942813446639E-002 Wall clock time [hours] = 4.708E-04 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 2.48 Maximum used memory per cpu [MBytes] = 35.562 Maximum used memory [MBytes] = 140.969 real 3.26 user 11.15 sys 0.46