Pencil1: iapp, nprocs, ncpus= 0 1 1 SVN: ------- v. ( ) $Id$ 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$ nxgrid, nygrid, nzgrid= 240 480 1 T WARNING: initialize_boundcond: invalid value of gamma, set it to 5/3! 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 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.3089969389957471E-002 0.0000000000000000E+000 (x,y,z)_ogrid(point2)= 3.0409771332038779E-003 -1.6362461737446838 0.0000000000000000E+000 Non-linear grid in radial direction - dx_rcyl, dx_rogrid: 5.8734793710668110E-005 2.2077187365501056E-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.88937751715686852 dx/dr_rogrid 0.94365885420207862 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= 0.0000000000000000E+000 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_map: particle_mesh = cic initialize_particles_radius: mass per dust grain mpmat= 2.8316221784356009E-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! init_uu: nothing init_lnrho: initlnrho(1) = const_rho init_lnrho: cs2bot,cs2top= 1.0000000000000000 1.0000000000000000 init_energy: initss(1) = const_TT initchem: nothing gaunoise_vect: ampl=0 for i1,i2= 1 3 gaunoise_ogrid: ampl=0 for i1,i2= 1 3 set_particles_radius: constant radius init_particles: Random particle positions within a box init_particles: Constant particle velocity init_particles: vpx0, vpy0, vpz0= 0.0000000000000000E+000 3000.0000000000000 0.0000000000000000E+000 start.x has completed successfully real 8.79 user 4.10 sys 0.26