Pencil1: iapp, nprocs, ncpus= 0 1 1 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 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 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. 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! 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.0000000000000000E+000 (x,y,z)_ogrid(point2)= 3.0409771332038783E-003 -1.6362461737446836 0.0000000000000000E+000 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 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.0000000000000000E+000 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! 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.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 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.0000000000000000E+000 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 8313 of file solid_cells_ogrid.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic' At line 8314 of file solid_cells_ogrid.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'cyclic' 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 8318 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 8339 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 2353 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 2353 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'gp' of procedure 'interpolate_linear_range' At line 2353 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range' At line 4829 of file particles_dust.f90 Fortran runtime warning: An array temporary was created for argument 'ineark' of procedure 'calc_brownian_force' 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] = 2.618E-03 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 8.18 Maximum used memory per cpu [MBytes] = 415.934 Maximum used memory [MBytes] = 415.934 real 12.80 user 12.25 sys 0.32