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: ------- 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= 20 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 64 64 1 Lx, Ly, Lz= 6.2831853071795862 6.2831853071795862 6.2831853071795862 Vbox= 39.478417604357432 rsnap: read snapshot var.dat in 1.0000000474974513E-003 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_selfgravity: epicycle frequency kappa = 1.0000000000000000 select_eos_variable: Using lnrho and cs2 initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000 initialize_particles: mass per superparticle mp_swarm = 3.8553142191755302E-004 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 4.0000000000000001E-002 initialize particles: set minimum collisional time-scale equal to two times the Courant time-step. initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method. 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=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 dvvp_dt: Add Coriolis force; Omega= 1.0000000000000000 particle_gravity: Calculating gravity particle_gravity: No gravity in x-direction. particle_gravity: No gravity in z-direction. particle_gravity: No radial gravity dvvp_dt_selfgrav: add self-gravity At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' 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$ duu_dt: diagnostics ... shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000 dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt: Add drag force; tausp= 1.0000000000000000 ---it-------t--------dt-------vpxmax------vpymax------vpzmax-----rhopm----rhopmin---rhopmax---gpotenp-- 0 0.000 4.179E-03 0.000E+00 0.000E+00 0.000E+00 1.000000 0.999003 1.000998 -0.000003 5 0.021 4.179E-03 1.267E-04 1.321E-04 0.000E+00 1.000000 0.999000 1.001001 -0.000003 10 0.042 4.179E-03 2.458E-04 2.678E-04 0.000E+00 1.000000 0.998992 1.001009 -0.000003 15 0.063 4.179E-03 3.578E-04 4.076E-04 0.000E+00 1.000000 0.998980 1.001022 -0.000004 20 0.084 4.179E-03 4.631E-04 5.523E-04 0.000E+00 1.000000 0.998962 1.001040 -0.000004 25 0.104 4.178E-03 5.621E-04 7.028E-04 0.000E+00 1.000000 0.998940 1.001062 -0.000004 30 0.125 4.178E-03 6.550E-04 8.599E-04 0.000E+00 1.000000 0.998913 1.001089 -0.000004 35 0.146 4.178E-03 7.419E-04 1.024E-03 0.000E+00 1.000000 0.998882 1.001120 -0.000005 40 0.167 4.178E-03 8.229E-04 1.197E-03 0.000E+00 1.000000 0.998846 1.001156 -0.000005 45 0.188 4.178E-03 8.979E-04 1.380E-03 0.000E+00 1.000000 0.998806 1.001196 -0.000006 50 0.209 4.178E-03 9.667E-04 1.572E-03 0.000E+00 1.000000 0.998762 1.001241 -0.000007 55 0.230 4.177E-03 1.029E-03 1.776E-03 0.000E+00 1.000000 0.998713 1.001290 -0.000007 60 0.251 4.177E-03 1.085E-03 1.992E-03 0.000E+00 1.000000 0.998660 1.001343 -0.000008 65 0.272 4.177E-03 1.133E-03 2.222E-03 0.000E+00 1.000000 0.998602 1.001402 -0.000009 70 0.292 4.177E-03 1.173E-03 2.466E-03 0.000E+00 1.000000 0.998539 1.001465 -0.000010 75 0.313 4.176E-03 1.205E-03 2.727E-03 0.000E+00 1.000000 0.998471 1.001532 -0.000011 80 0.334 4.176E-03 1.227E-03 3.004E-03 0.000E+00 1.000000 0.998399 1.001605 -0.000013 85 0.355 4.176E-03 1.238E-03 3.300E-03 0.000E+00 1.000000 0.998321 1.001684 -0.000015 90 0.376 4.175E-03 1.237E-03 3.616E-03 0.000E+00 1.000000 0.998237 1.001768 -0.000016 95 0.397 4.175E-03 1.223E-03 3.953E-03 0.000E+00 1.000000 0.998148 1.001858 -0.000019 Simulation finished after 100 time-steps Writing final snapshot at time t = 0.41773351067090014 Wall clock time [hours] = 8.151E-03 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 2.76 Maximum used memory per cpu [MBytes] = 29.680 Maximum used memory [MBytes] = 29.680 real 30.10 user 30.03 sys 0.06