SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ 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: particles_dust.f90 v. 1.1 (2018/08/24 15:48:10) wlyra SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 14 (ldebug= F ) This is a 0-D run nxgrid, nygrid, nzgrid= 1 1 1 Lx, Ly, Lz= 6.2831853071795862 6.2831853071795862 6.2831853071795862 Vbox= 1.0000000000000000 rsnap: read snapshot var.dat in 0.0000000000000000 seconds 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 units_eos: unit_temperature= 1.2027220936797545E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_density: no hydro, turned off continuity equation initialize_density: 0-D run, turned off continuity equation diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using rho and cs2 initialize_particles: Global pressure gradient with beta_dPdr_dust= -0.10000000000000001 initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000 initialize_particles: mass per constituent particle mpmat= 4.1887902047863905 initialize_particles: mass per superparticle mp_swarm = 1.0000000000000000E-002 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 1.0000000000000000E-002 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= 4.1887902047863905 sh: 1: Syntax error: Bad fd number pde: ENTER SVN: ------- v. ( ) $Id$ duu_dt: diagnostics ... dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 SOLVE dlncc_dt Bcs for cc, x: < p>, y: < p>, z: < p> dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dap_dt_pencil: Calculate dap/dt dnpswarm_dt_pencil: Calculate dnpswarm_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 particle_gravity: Calculating gravity particle_gravity: No gravity in x-direction. particle_gravity: No gravity in z-direction. particle_gravity: No radial gravity ---it--------t---------dt-------dtfragp---dtsweepp-----rhom------rhomin-----rhomax------vpxm--------xpm-------xp2m--------zpm-------zp2m--------npm-------np2m-------npmax-------apm-------ap2m-------apmin------apmax----npswarmm------mpt-------mcct------dvp12m----dvp22mwnp- 0 0.000E+00 2.667E+01 0.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 1.000E+02 1.000E+04 1.000E+02 1.000E+00 1.000E+00 1.000E+00 1.000E+00 2.387E-03 1.000E+00 4.000E+00 1.000E-02 1.000E-02 10 3.249E+01 2.326E-06 0.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 1.000E+02 1.000E+04 1.000E+02 1.324E+00 1.754E+00 1.324E+00 1.324E+00 9.237E-04 8.988E-01 4.045E+00 1.000E-02 1.000E-02 Time step has become too short: dt = 7.4439867521649232E-010 Simulation finished after 16 time-steps Writing final snapshot at time t = 32.487410832055957 Wall clock time [hours] = 6.667E-06 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 14.9 Maximum used memory per cpu [MBytes] = 9.406 Maximum used memory [MBytes] = 9.406 real 0.08 user 0.06 sys 0.01