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
rsnap: read snapshot var.dat in 5.6476039999999996E-003 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] = 2.165E-05 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.79
Maximum used memory per cpu [MBytes] = 33.652
Maximum used memory [MBytes] = 134.227
real 5.26
user 16.77
sys 0.27