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$ WARNING: find_namelist: namelist "special_init_pars" is missing!! 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: gravity_r.f90,v v. 1.1 (2018/08/24 15:48:10) wlyra 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 8 Lx, Ly, Lz= 5.2000000000000002 5.2000000000000002 0.52000000000000002 Vbox= 14.060800000000002 rsnap: read snapshot var.dat in 2.9525700000000000E-004 seconds setup_slices: slice_position = m setup_slices: ix_loc, (video files) = 19 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.2027220936797545E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -4.6051701859880909 0.10000000000000001 1.0000000000000004E-002 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: non-smoothed newtonian gravity diffusion: (Dx*d^6/dx^6 + Dy*d^6/dy^6 + Dz*d^6/dz^6)rho diffusion: shock diffusion WARNING: initialize_density: For diffusion energy/momentum correction should use lmassdiff_fix=T! initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. WARNING: initialize_energy: llocal_iso=T. Make sure you have the appropriate INITIAL_CONDITION in Makefile.local! select_eos_variable: Using rho and cs2 resistivity: hyper3_aniso resistivity: shock initialize_shock: prenormalised shock_factor sum= 1.0000000000000000 viscous force(i): (nu.del6)ui + ((nu.uij5).glnrho) viscous force: nu_shock*(XXXXXXXXXXX) rt_int,rt_ext 0.40000000000000002 2.5000000000000000 initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000 initialize_particles: mass per superparticle mp_swarm = 2.8121600000000009E-004 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 0.16384000000000001 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=flk dxxp_dt: Particles boundary condition bcpy=rmv 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: Newtonian gravity from a fixed central object 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 ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 16) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomx ( 94) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomy ( 95) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil aa ( 151) 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_energy: max(advec_cs2) = 0.89294479277389083 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) = 5.6982620982701206 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: diffrho_shock= 5.0000000000000000 dlnrho_dt: diffrho_hyper3=(Dx,Dy,Dz)= 4.9999999999999997E-012 4.9999999999999997E-012 2.9999999999999998E-013 dlnrho_dt: max(diffus_diffrho ) = 0.41521311515592457 dlnrho_dt: max(diffus_diffrho3) = 4.5208832584840177E-006 daa_dt: SOLVE Bcs for Ax, x: < p>, y: < p>, z: < p> Bcs for Ay, x: < p>, y: < p>, z: < p> Bcs for Az, x: < p>, y: < p>, z: < p> daa_dt: iresistivity=hyper3-aniso shock daa_dt: max(diffus_eta) = 0.41521311515592457 daa_dt: max(diffus_eta2) = 0.0000000000000000 daa_dt: max(diffus_eta3) = 4.5208832584840177E-006 dspecial_dt: SOLVE dSPECIAL_dt dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt: Add drag force; tausp= 1.0000000000000000 pde: freezing variables for r < 0.40000000000000002 : T T T T T T T F pde: freezing variables for r > 2.5000000000000000 : T T T T T T T F ----it----------t-----------dt---------dtu---------dtnu--------dtb--------dtd--------rhomax------rhomin---------bzm-----------urupm--------npmax------rhopmax----nparmax-- 0 0.00000E+00 1.83659E-02 1.07163E+00 9.05205E-01 2.2540E-01 9.05205E-01 1.00000E+00 1.00000E+00 -0.0000000000 0.0000000000 3.00000E+00 1.804E-01 1.3000E+02 1 1.83659E-02 1.83202E-02 1.06897E+00 9.10437E-01 2.2489E-01 9.10437E-01 1.00800E+00 9.91218E-01 0.0000000000 0.0000000000 3.00000E+00 1.842E-01 1.3000E+02 2 3.66862E-02 1.82814E-02 1.06671E+00 9.15041E-01 2.2447E-01 9.15041E-01 1.01354E+00 9.83617E-01 0.0000000000 0.0000000000 3.00000E+00 1.841E-01 1.2900E+02 3 5.49676E-02 1.82494E-02 1.06484E+00 9.18981E-01 2.2412E-01 9.18981E-01 1.01704E+00 9.77203E-01 0.0000000000 0.0000000000 3.00000E+00 1.800E-01 1.2900E+02 4 7.32170E-02 1.82235E-02 1.06332E+00 9.22240E-01 2.2384E-01 9.22240E-01 1.01887E+00 9.71922E-01 0.0000000000 0.0000000000 3.00000E+00 1.723E-01 1.2900E+02 Simulation finished after 5 time-steps Writing final snapshot at time t = 9.1440517006895899E-002 Wall clock time [hours] = 3.662E-05 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 3.03 Maximum used memory per cpu [MBytes] = 34.281 Maximum used memory [MBytes] = 136.164 real 7.26 user 26.95 sys 0.66