Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "special_init_pars" is missing!! SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ initialize_mpicomm: enabled MPI 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 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 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 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.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 2.8121600000000004E-004 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 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. pencil_consistency_check: checking pencil case At line 934 of file particles_map.f90At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' 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 ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 10) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomx ( 88) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil pomy ( 89) 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.89294479277389061 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.6982620982701233 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.41521311515592568 dlnrho_dt: max(diffus_diffrho3) = 4.5208832584840186E-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.41521311515592568 daa_dt: max(diffus_eta2) = 0.0000000000000000E+000 daa_dt: max(diffus_eta3) = 4.5208832584840186E-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.83267E-02 1.06935E+00 9.11056E-01 2.2497E-01 9.11056E-01 1.00795E+00 9.91183E-01 0.0000000000 0.0000000000 3.00000E+00 1.842E-01 1.3000E+02 2 3.66926E-02 1.83024E-02 1.06793E+00 9.17213E-01 2.2472E-01 9.17213E-01 1.01317E+00 9.83501E-01 0.0000000000 0.0000000000 3.00000E+00 1.841E-01 1.2900E+02 3 5.49950E-02 1.82861E-02 1.06698E+00 9.23186E-01 2.2457E-01 9.23186E-01 1.01585E+00 9.77067E-01 -0.0000000000 0.0000000000 3.00000E+00 1.800E-01 1.2900E+02 4 7.32811E-02 1.82723E-02 1.06617E+00 9.28527E-01 2.2444E-01 9.28527E-01 1.01625E+00 9.72017E-01 -0.0000000000 0.0000000000 3.00000E+00 1.722E-01 1.2900E+02 Simulation finished after 5 time-steps Writing final snapshot at time t = 9.1553375783758453E-002 Wall clock time [hours] = 1.896E-04 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 15.7 Maximum used memory per cpu [MBytes] = 16.957 Maximum used memory [MBytes] = 67.461 real 37.73 user 51.84 sys 7.48