SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) no cvsid is given i
SVN: ------- v. ( ) no cvsid is given i
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$
The verbose level is ip= 14 (ldebug= F )
This is a 2-D run
nxgrid, nygrid, nzgrid= 64 1 64
Lx, Ly, Lz= 1.0471975511965976E-002 1.0471975511965976E-002 1.0471975511965976E-002
Vbox= 1.0966227112321507E-004
rsnap: read snapshot var.dat in 1.1712500000000000E-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.2027220936797545E-008
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 3.6757541328186911 6.2831853071795862 39.478417604357439 1.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 rho and cs2
initialize_shear: turned shear acceleration off and handed it over to Particles_drag.
initialize_particles: turned off and hand over shear acceleration to Particles_drag.
initialize_particles: dust-to-gas ratio eps_dtog= 3.0000000000000000
initialize_particles: mass per constituent particle mpmat= 0.0000000000000000
initialize_particles: mass per superparticle mp_swarm = 8.0319046232823537E-008
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 3.0000000000000000
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.
initialize_particles_drag: tdrag = 1.5915494309189534E-002
initialize_particles_drag: reset mp_swarm = 8.0319046232823537E-008
initialize_particles_drag: reset rhop_swarm = 3.0000000000000000
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
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: 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_energy: max(advec_cs2) = 2949120000.0000010
duu_dt: SOLVE
Bcs for ux, x: < she>, y: < p>, z: < p>
Bcs for uy, x: < she>, y: < p>, z: < p>
Bcs for uz, x: < she>, y: < p>, z: < p>
duu_dt: max(advec_uu) = 71.955352855863524
dlnrho_dt: SOLVE
Bcs for lnrho, x: < she>, y: < p>, z: < p>
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
shearing: Sshear,Sshear1= -9.4247779607693793 -9.4247779607693793
shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
--it-------t--------dt-----dtv-------uxmin------------uxmax------------uymin------------uymax---------uzmin------uzmax---------rhom-----------rhomin----------rhomax-----nparmax-nmigmmax------vpxmin----------vpxmax----------vpymin-----------vpymax--------vpzmin-----vpzmax---npmax------rhopm----------rhopmin---------rhopmax-----
0 0.0000000 1.66E-05 1.00 1.177355970E-02 1.177366818E-02 -7.868703237E-02 -7.868694066E-02 -5.424E-08 5.424E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458340E-03 -3.92449257E-03 -7.849080524E-02 -7.849071399E-02 -5.196E-08 5.196E-08 1.0 3.000000000E+00 2.999997016E+00 3.000002984E+00
100 0.0016551 1.66E-05 1.00 1.177355942E-02 1.177366846E-02 -7.868703259E-02 -7.868694044E-02 -5.450E-08 5.450E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458358E-03 -3.92449239E-03 -7.849080544E-02 -7.849071379E-02 -5.218E-08 5.218E-08 1.0 3.000000000E+00 2.999997002E+00 3.000002998E+00
200 0.0033102 1.66E-05 1.00 1.177355916E-02 1.177366872E-02 -7.868703281E-02 -7.868694022E-02 -5.477E-08 5.477E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458376E-03 -3.92449220E-03 -7.849080563E-02 -7.849071360E-02 -5.240E-08 5.240E-08 1.0 3.000000000E+00 2.999996988E+00 3.000003012E+00
300 0.0049653 1.66E-05 1.00 1.177355890E-02 1.177366898E-02 -7.868703302E-02 -7.868694001E-02 -5.503E-08 5.503E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458394E-03 -3.92449202E-03 -7.849080583E-02 -7.849071340E-02 -5.262E-08 5.262E-08 1.0 3.000000000E+00 2.999996975E+00 3.000003025E+00
400 0.0066204 1.66E-05 1.00 1.177355865E-02 1.177366924E-02 -7.868703322E-02 -7.868693981E-02 -5.528E-08 5.528E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458412E-03 -3.92449184E-03 -7.849080602E-02 -7.849071321E-02 -5.283E-08 5.283E-08 1.0 3.000000000E+00 2.999996961E+00 3.000003039E+00
500 0.0082754 1.66E-05 1.00 1.177355840E-02 1.177366948E-02 -7.868703342E-02 -7.868693960E-02 -5.554E-08 5.554E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458432E-03 -3.92449164E-03 -7.849080622E-02 -7.849071301E-02 -5.305E-08 5.305E-08 1.0 3.000000000E+00 2.999996947E+00 3.000003053E+00
Simulation finished after 501 time-steps
Writing final snapshot at time t = 8.2919934886555198E-003
Wall clock time [hours] = 1.270E-03 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.11
Maximum used memory per cpu [MBytes] = 29.027
Maximum used memory [MBytes] = 115.289
real 4.94
user 18.51
sys 0.17