SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.1 (2020/02/23 20:50:00) xiangyu
SVN: run.in,v v. 1.1 (2020/02/23 20:50:00) xiangyu
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: particles_dust.f90 v. 1.1 (2018/08/24 15:48:10) wlyra
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
The verbose level is ip= 14 (ldebug= F )
This is a 3-D run
nxgrid, nygrid, nzgrid= 16 16 16
Lx, Ly, Lz= 5.0000000000000000 5.0000000000000000 5.0000000000000000
Vbox= 125.00000000000000
input_persist_forcing: 0.0000000000000000 0.0000000000000000 0.0000000000000000
input_persist_forcing: -10.000000000000000
rsnap: read snapshot var.dat in 1.6681800000000000E-004 seconds
setup_slices: slice_position = p
setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4
units_general: unit_velocity= 100000.00000000000
units_general: unit_density= 6.7710000000000037E-023
units_general: unit_length= 3.0859999999999995E+018
units_general: unit_magnetic= 3.5449077018110318
units_eos: unit_temperature= 240.54441873595090
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -2.1517456033972402 0.34100000000000003 11.472283460000003 2.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
diffusion: nothing (i.e. no mass diffusion)
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
select_eos_variable: Using lnrho and cs2
initialize_shock: prenormalised shock_factor sum= 1.0000000000000000
initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002
initialize_particles: mass per constituent particle mpmat= 1.4849261275967752E-046
initialize_particles: mass per superparticle mp_swarm = 2.5257024657983126E-003
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 8.2762218399279108E-002
initialize particles: set minimum collisional time-scale equal to two times the Courant time-step.
initialize_particles_map: selected nearest-grid-point for particle-mesh method.
initialize_particles_radius: mass per dust grain mpmat= 1.4849261275967752E-046
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'
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
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: 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 rho1 ( 14) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil TT1 ( 59) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil ppvap ( 74) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil np ( 173) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil rhop ( 174) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil cc ( 189) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil cc1 ( 190) 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_hydro: upwinding advection term
calc_pencils_energy: max(advec_cs2) = 3.5721523200000016
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) = 0.0000000000000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < p>, z: < p>
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
dvvp_dt: Add drag force; tausp= 0.0000000000000000
---it--------t---------dt-------urms-------umax-------epsK-------rhom-------rho2m-----rhomin-----rhomax----shockmax-----vprms------vpmax-------apm-------ap2m-------apmin------apmax---
0 0.000E+00 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 1.000E-23 1.000E-46 1.000E-23 1.000E-23
20 9.524E+00 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 3.020E-23 9.121E-46 3.020E-23 3.020E-23
40 1.905E+01 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 5.286E-23 2.794E-45 5.286E-23 5.286E-23
60 2.857E+01 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 7.551E-23 5.702E-45 7.551E-23 7.551E-23
80 3.809E+01 4.762E-01 0.000E+00 0.000E+00 0.000E+00 4.933E+01 2.433E+03 4.933E+01 4.933E+01 0.000E+00 0.000E+00 0.000E+00 9.817E-23 9.637E-45 9.817E-23 9.817E-23
Simulation finished after 81 time-steps
Writing final snapshot at time t = 38.571146099636820
Wall clock time [hours] = 1.312E-04 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 0.204
Maximum used memory per cpu [MBytes] = 31.859
Maximum used memory [MBytes] = 249.598
real 1.92
user 12.47
sys 0.45