SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.1 (2018/06/08 07:14:05) xiangyu
SVN: run.in,v v. 1.1 (2018/06/08 07:14:05) 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: ------- 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$
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= 0.12500000000000000 0.12500000000000000 0.12500000000000000
Vbox= 1.9531250000000000E-003
input_persist_forcing: 0.0000000000000000 0.0000000000000000 0.0000000000000000
input_persist_forcing: -10.000000000000000
rsnap: read snapshot var.dat in 1.7754399999999999E-004 seconds
setup_slices: slice_position = p
setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4
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.0309046517255039E-003
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 4.1351665567423561 5.0000000000000000 567.85714285714312 8.5714285714285712
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
initialize_gravity: no x-gravity
initialize_gravity: no y-gravity
initialize_gravity: constant gravz= -9.8100000000000005
diffusion: nothing (i.e. no mass diffusion)
WARNING: initialize_density: are rho_const or lnrho_const ok?!
inconsistency between the density constants from eos
(rho0 or lnrho0) and the ones from the density module
(rho_const or lnrho_const). It may damage your
simulation if you are using them in different places.
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
select_eos_variable: Using lnrho and ss
Supersaturation routine
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002
initialize_particles: mass per constituent particle mpmat= 4.1887902047863918E-012
initialize_particles: mass per superparticle mp_swarm = 7.8332133582218749E-009
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 1.6427439060621721E-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= 4.1887902047863918E-012
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: Constant 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_hydro: upwinding advection term
calc_pencils_energy: max(advec_cs2) = 1200490.8564055739
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
SOLVE dacc_dt
Bcs for acc, x: < p>, y: < p>, z: < p>
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
dvvp_dt: Add drag force; tausp= 0.0000000000000000
dnpswarm_dt_pencil: Calculate dnpswarm_dt
---it---------t---------dt-------urms-------epsK-------admom0--------admom1---------admom2------admom3------admom4------admom6-----admom12-----admom24--------accrms------------accm------------ttcrms------------ttcm-------condensationRaterms-condensationRatem-waterMixingRatiorms-waterMixingRatiom------ssatrms-----------ssatm------buoyancyrms--buoyancym--
0 0.00000 4.482E-04 0.000E+00 0.000E+00 2.6208E+10 2.6208496094E+05 2.6208E+00 2.6208E-05 2.6208E-10 2.6208E-20 -4.9582E+01 -1.0958E+02 1.5700000000E-02 1.5700000000E-02 2.9200000000E+02 2.9200000000E+02 0.0000000000E+00 0.0000E+00 0.0000000000E+00 0.0000E+00 0.0000000000E+00 0.0000000000E+00 2.1319E+00 2.1319E+00
20 0.00896 4.482E-04 2.537E-02 8.159E-03 2.6208E+10 2.6211013342E+05 2.6214E+00 2.6216E-05 2.6219E-10 2.6224E-20 -4.9581E+01 -1.0958E+02 1.5650646729E-02 1.5650644535E-02 2.9212277563E+02 2.9212277479E+02 3.6822149852E-03 3.6572E-03 1.1413912625E-01 1.0981E-01 1.0184102201E-02 1.0009504318E-02 1.0907E+00 1.0467E+00
40 0.01793 4.482E-04 3.304E-02 1.359E-02 2.6208E+10 2.6212240193E+05 2.6216E+00 2.6220E-05 2.6223E-10 2.6231E-20 -4.9581E+01 -1.0958E+02 1.5627459035E-02 1.5627457489E-02 2.9218045473E+02 2.9218045401E+02 1.7575912814E-03 1.7470E-03 1.1411066073E-01 1.0983E-01 5.2574288901E-03 4.9775036009E-03 1.0865E+00 1.0429E+00
60 0.02689 4.482E-04 3.921E-02 1.897E-02 2.6208E+10 2.6212858179E+05 2.6217E+00 2.6222E-05 2.6226E-10 2.6235E-20 -4.9581E+01 -1.0958E+02 1.5616191871E-02 1.5616191224E-02 2.9220847993E+02 2.9220847954E+02 8.8406059899E-04 8.6598E-04 1.1402661431E-01 1.0984E-01 2.8110849291E-03 2.5419035243E-03 1.0836E+00 1.0410E+00
80 0.03585 4.482E-04 4.350E-02 2.299E-02 2.6208E+10 2.6213176946E+05 2.6218E+00 2.6223E-05 2.6227E-10 2.6237E-20 -4.9581E+01 -1.0958E+02 1.5610538728E-02 1.5610538508E-02 2.9222254121E+02 2.9222254103E+02 4.6045633095E-04 4.4061E-04 1.1392700018E-01 1.0984E-01 1.5392892108E-03 1.3223763231E-03 1.0816E+00 1.0400E+00
Simulation finished after 81 time-steps
Writing final snapshot at time t = 3.6301992990234636E-002
Wall clock time [hours] = 3.725E-04 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 0.306
Maximum used memory per cpu [MBytes] = 32.207
Maximum used memory [MBytes] = 254.016
real 1.82
user 11.83
sys 0.35