SVN: ------- v. ( ) $Id$
WARNING: find_namelist: namelist "special_init_pars" is missing!!
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: particles_density. v. 20849 (2013-08-06 18:45:43) anders
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= 6.2831853071795866E-003 1.0000000000000000 6.2831853071795866E-003
Vbox= 3.9478417604357436E-005
rsnap: read snapshot var.dat in 1.3002599999999999E-004 seconds
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= 4.8108883747190188E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 4.0812192409268553 6.2831853071795862 23.687050562614463 0.40000000000000008
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
diffusion: mesh hyperdiffusion
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
initialize_density: Global density gradient with beta_glnrho_global= -0.10000000000000001 0.0000000000000000 0.0000000000000000
select_eos_variable: Using rho and ss
viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq
initialize_special: ux0, uy0 = 9.2791620582046985E-003 -0.12583302733752555
initialize_special: vpx0 = -1.4140910635015236E-003 -1.5730141989406352E-002
initialize_special: vpy0 = -0.12580297790242614 -0.12504652023805524
initialize_particles: Number of particle species = 2
initialize_particles: tausp_species = 6.7640850814055528E-003 1.5915494309189534E-002
initialize_particles: dust-to-gas ratio eps_dtog= 1.5000000000000000
initialize_particles: mass per constituent particle mpmat= 0.0000000000000000
initialize_particles: mass per superparticle mp_swarm = 1.8071785402385300E-009
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 0.18750000000000000
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'
drhopswarm_dt_pencil: Calculate drhopswarm_dt
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
dvvp_dt: Add Coriolis force; Omega= 6.2831853071795862
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) = 8192000000.0376654
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>
coriolis_cartesian: add Coriolis force; Omega= 6.2831853071795862
duu_dt: max(advec_uu) = 94.517383759778269
dlnrho_dt: SOLVE
Bcs for lnrho, x: < she>, y: < p>, z: < p>
dlnrho_dt: diffrho_hyper3_mesh= 1.4577334289680862E-003
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 29.696701559112856
denergy_dt: adding global pressure gradient force
dspecial_dt: special diagnostics
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
dvvp_dt_pencil: adding feedback from dust to gas
dvvp_dt: Add drag force; tausp= 0.0000000000000000
----it-------t-------dt----dtv----rhom-----drho2m-----urms------umax-------uxm--------uym------rdux2m----rduy2m----ruz2m----rduxduym---nparmax---nmigmmax----vpxm-------vpym----rhopdvpx2m-rhopdvpy2m-rhopvpz2m---npmax----drhop2m---rhopmax--
0 0.00000 9.9E-06 1.00 1.000E+00 3.757E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.283E-16 8.028E-16 9.283E-16 -7.108E-16 8.192E+03 0.000E+00 -8.572E-03 -1.254E-01 1.395E-16 1.487E-16 1.639E-16 8.000E+00 4.967E-13 1.500E+00
100 0.00099 9.9E-06 1.00 1.000E+00 3.727E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.325E-16 8.064E-16 9.325E-16 -7.141E-16 8.192E+03 0.000E+00 -8.572E-03 -1.254E-01 1.399E-16 1.492E-16 1.644E-16 8.000E+00 4.986E-13 1.500E+00
200 0.00199 9.9E-06 1.00 1.000E+00 3.766E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.365E-16 8.099E-16 9.365E-16 -7.171E-16 8.192E+03 6.400E+01 -8.572E-03 -1.254E-01 1.404E-16 1.497E-16 1.650E-16 8.000E+00 5.005E-13 1.500E+00
300 0.00298 9.9E-06 1.00 1.000E+00 3.778E-23 1.262E-01 1.262E-01 9.279E-03 -1.258E-01 9.404E-16 8.133E-16 9.404E-16 -7.201E-16 8.192E+03 0.000E+00 -8.572E-03 -1.254E-01 1.409E-16 1.502E-16 1.656E-16 8.000E+00 5.024E-13 1.500E+00
Simulation finished after 301 time-steps
Writing final snapshot at time t = 2.9886280746550956E-003
Wall clock time [hours] = 3.914E-03 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.27
Maximum used memory per cpu [MBytes] = 29.324
Maximum used memory [MBytes] = 115.914
real 14.58
user 57.08
sys 0.16