Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) no cvsid is given i SVN: ------- v. ( ) no cvsid is given i initialize_mpicomm: enabled MPI WARNING: run: quad precision suppressed! 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.4777491000000000E-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] = 2.947E-03 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 9.20 Maximum used memory per cpu [MBytes] = 416.152 Maximum used memory [MBytes] = 416.152 real 14.24 user 13.53 sys 0.45