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: ------- 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= 240 480 1
Lx, Ly, Lz= 5.0000000000000003E-002 0.10000000000000001 1.0000000000000000E-002
Vbox= 5.0000000000000010E-003
rsnap: read snapshot var.dat in 1.0153262000000000E-002 seconds
T
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
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
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 rho and TT
heat conduction: nothing
heat conduction: nothing
heat conduction: nothing
Number of reactions= 1
Number of species= 5
viscous force: nu is calculated for a mixture
nxgrid_ogrid, nygrid_ogrid, nzgrid_ogrid= 48 240 1
Cylidner radius= 2.5000000000000001E-003
Cylindrical grid radius= 7.4999999999999997E-003
By using fbc[x,y,z] I found the flow direction to be in the 2 direction.
(x,y,z)_ogrid(point)= 4.0752373928797754E-003 -1.3089969389957545E-002 0.0000000000000000
(x,y,z)_ogrid(point2)= 3.0409771332038783E-003 -1.6362461737446836 0.0000000000000000
Non-linear grid in radial direction - dx_rcyl, dx_rogrid: 5.8734793710668110E-005 2.2077187365501013E-004
Theta grid spacing - r_cyl*dy_ogrid,r_int_outer*dy_ogrid,r_ogrid*dy_ogrid 6.5449846949787731E-005 1.9629176034748199E-004 1.9634954084936317E-004
dtheta/dr_surf 1.1143283702024811
dtheta/dr_rogrid 0.88937751715687030
dx/dr_rogrid 0.94365885420208051
dx/dtheta_rogrid 1.0610329539459631
Cartesian grid spacing - dx, dy, dz: 2.0833333333333335E-004 2.0876826722338206E-004 1.0000000000000000E-002
Timestep factor: 4
Cylinder boundary condition: Third order SBP boundary closures
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry_simple.
The values provided are for one-atomic ideal gas. Use at own risk!
initialize_mpicomm_ogrid: enabled MPI on overlapping grid
Interpolation zone: r_ogrid, r_int_outer, r_int_inner 7.4999999999999997E-003 7.4977929474027072E-003 6.8714881457325608E-003
interpolation_method==1: Linear interpolation used
particle_interpolate==3: Quadratic particle interpolation
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry_simple.
The values provided are for one-atomic ideal gas. Use at own risk!
initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002
initialize_particles: mass per constituent particle mpmat= 2.8316221784356013E-015
initialize_particles: mass per superparticle mp_swarm = 5.0000000000000008E-007
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 11.520000000000001
initialize particles: set minimum collisional time-scale equal to two times the Courant time-step.
initialize_particles_map: selected cloud-in-cell for particle-mesh method.
initialize_particles_radius: mass per dust grain mpmat= 2.8316221784356013E-015
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry_simple.
The values provided are for one-atomic ideal gas. Use at own risk!
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=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: 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(p%advec_cs2) = 50721411632298776.
duu_dt: SOLVE
Bcs for ux, x: < p>, y: < set:e2>, z: < p>
Bcs for uy, x: < p>, y: < set:e2>, z: < p>
Bcs for uz, x: < p>, y: < set:e2>, z: < p>
duu_dt: max(advec_uu) = 14370000.000000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < set:e2>, z: < p>
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: SOLVE denergy_dt
Bcs for lnTT, x: < p>, y: < set:e2>, z: < p>
denergy_dt: cs2 = 1103020644.7945871
dchemistry_dt: SOLVE dchemistry_dt
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
dvvp_dt: Add drag force; tausp= 0.0000000000000000
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry_simple.
The values provided are for one-atomic ideal gas. Use at own risk!
inlet_profile: uniform
inlet_YY_profile: uniform,
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry_simple.
The values provided are for one-atomic ideal gas. Use at own risk!
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_chemistry_simple.
The values provided are for one-atomic ideal gas. Use at own risk!
At line 8315 of file solid_cells_ogrid.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic'
At line 8316 of file solid_cells_ogrid.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic'
At line 8317 of file solid_cells_ogrid.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic'
At line 8320 of file solid_cells_ogrid.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic'
At line 8322 of file solid_cells_ogrid.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic'
At line 8341 of file solid_cells_ogrid.f90
Fortran runtime warning: An array temporary was created for argument 'bx_bound' of procedure 'boundary_x_8_6th'
---it-----------t-----------dt------urms------uxm------uym-----umax-------rhom-------c_dragx-------c_dragy-----nparmax--nparmin--nparsum-----Mamax-----
0 0.00000E+00 4.061E-09 3008.43 -0.00 2999.85 3977.06 0.00128313 2.1112E-09 1.8942E-04 0. 0. 0. 1.1975E-01
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 2370 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range'
At line 2370 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'gp' of procedure 'interpolate_linear_range'
At line 2370 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range'
At line 5169 of file particles_dust.f90
Fortran runtime warning: An array temporary was created for argument 'ineark' of procedure 'calc_brownian_force_li_ahmadi'
1 4.06149E-09 4.100E-09 3000.06 0.00 2989.10 3589.92 0.00128313 2.1508E-08 8.5367E-04 24. 24. 24. 1.0809E-01
2 8.16100E-09 4.100E-09 3000.06 -0.00 2989.09 3589.92 0.00128313 6.9959E-08 1.7774E-03 48. 48. 48. 1.0809E-01
3 1.22605E-08 4.100E-09 3000.06 -0.00 2989.09 3589.92 0.00128313 1.5380E-07 2.8941E-03 73. 73. 73. 1.0809E-01
4 1.63600E-08 4.100E-09 3000.06 -0.00 2989.09 3589.92 0.00128313 2.7558E-07 4.1643E-03 98. 98. 98. 1.0809E-01
5 2.04595E-08 4.100E-09 3000.06 -0.00 2989.09 3589.92 0.00128313 4.3043E-07 5.4870E-03 98. 98. 98. 1.0809E-01
6 2.45590E-08 4.100E-09 3000.05 -0.00 2989.08 3589.92 0.00128313 6.1225E-07 6.8102E-03 98. 98. 98. 1.0809E-01
7 2.86585E-08 4.100E-09 3000.05 -0.00 2989.08 3589.92 0.00128313 8.1258E-07 8.0942E-03 98. 98. 98. 1.0809E-01
8 3.27580E-08 4.100E-09 3000.04 -0.00 2989.07 3589.92 0.00128313 1.0216E-06 9.2981E-03 98. 98. 98. 1.0809E-01
9 3.68575E-08 4.100E-09 3000.03 -0.00 2989.06 3589.92 0.00128313 1.2304E-06 1.0396E-02 98. 98. 98. 1.0809E-01
Simulation finished after 10 time-steps
Writing final snapshot at time t = 4.0957054895072293E-008
Wall clock time [hours] = 1.904E-03 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 5.95
Maximum used memory per cpu [MBytes] = 416.309
Maximum used memory [MBytes] = 416.309
real 10.80
user 8.64
sys 0.24