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 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: 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= 64 1 64 Lx, Ly, Lz= 1.0471975511965976E-002 1.0471975511965976E-002 1.0471975511965976E-002 Vbox= 1.0966227112321507E-004 rsnap: read snapshot var.dat in 2.7107500000000002E-004 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 3.6757541328186911 6.2831853071795862 39.478417604357439 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 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_shear: turned shear acceleration off and handed it over to Particles_drag. initialize_particles: turned off and hand over shear acceleration to Particles_drag. initialize_particles: dust-to-gas ratio eps_dtog= 3.0000000000000000 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000 initialize_particles: mass per superparticle mp_swarm = 8.0319046232823537E-008 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 3.0000000000000000 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. initialize_particles_drag: tdrag = 1.5915494309189534E-002 initialize_particles_drag: reset mp_swarm = 8.0319046232823537E-008 initialize_particles_drag: reset rhop_swarm = 3.0000000000000000 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 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(advec_cs2) = 2949120000.0000010 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < p> Bcs for uy, x: < she>, y: < p>, z: < p> Bcs for uz, x: < she>, y: < p>, z: < p> duu_dt: max(advec_uu) = 71.955352855863524 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 shearing: Sshear,Sshear1= -9.4247779607693793 -9.4247779607693793 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000 dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F --it-------t--------dt-----dtv-------uxmin------------uxmax------------uymin------------uymax---------uzmin------uzmax---------rhom-----------rhomin----------rhomax-----nparmax-nmigmmax------vpxmin----------vpxmax----------vpymin-----------vpymax--------vpzmin-----vpzmax---npmax------rhopm----------rhopmin---------rhopmax----- 0 0.0000000 1.66E-05 1.00 1.177355970E-02 1.177366818E-02 -7.868703237E-02 -7.868694066E-02 -5.424E-08 5.424E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458340E-03 -3.92449257E-03 -7.849080524E-02 -7.849071399E-02 -5.196E-08 5.196E-08 1.0 3.000000000E+00 2.999997016E+00 3.000002984E+00 100 0.0016551 1.66E-05 1.00 1.177355942E-02 1.177366846E-02 -7.868703259E-02 -7.868694044E-02 -5.450E-08 5.450E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458358E-03 -3.92449239E-03 -7.849080544E-02 -7.849071379E-02 -5.218E-08 5.218E-08 1.0 3.000000000E+00 2.999997002E+00 3.000002998E+00 200 0.0033102 1.66E-05 1.00 1.177355916E-02 1.177366872E-02 -7.868703281E-02 -7.868694022E-02 -5.477E-08 5.477E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458376E-03 -3.92449220E-03 -7.849080563E-02 -7.849071360E-02 -5.240E-08 5.240E-08 1.0 3.000000000E+00 2.999996988E+00 3.000003012E+00 300 0.0049653 1.66E-05 1.00 1.177355890E-02 1.177366898E-02 -7.868703302E-02 -7.868694001E-02 -5.503E-08 5.503E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458394E-03 -3.92449202E-03 -7.849080583E-02 -7.849071340E-02 -5.262E-08 5.262E-08 1.0 3.000000000E+00 2.999996975E+00 3.000003025E+00 400 0.0066204 1.66E-05 1.00 1.177355865E-02 1.177366924E-02 -7.868703322E-02 -7.868693981E-02 -5.528E-08 5.528E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458412E-03 -3.92449184E-03 -7.849080602E-02 -7.849071321E-02 -5.283E-08 5.283E-08 1.0 3.000000000E+00 2.999996961E+00 3.000003039E+00 500 0.0082754 1.66E-05 1.00 1.177355840E-02 1.177366948E-02 -7.868703342E-02 -7.868693960E-02 -5.554E-08 5.554E-08 1.000000000E+00 1.000000000E+00 1.000000000E+00 1024.0 0.0 -3.92458432E-03 -3.92449164E-03 -7.849080622E-02 -7.849071301E-02 -5.305E-08 5.305E-08 1.0 3.000000000E+00 2.999996947E+00 3.000003053E+00 Simulation finished after 501 time-steps Writing final snapshot at time t = 8.2919934886555198E-003 Wall clock time [hours] = 2.243E-03 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.97 Maximum used memory per cpu [MBytes] = 29.395 Maximum used memory [MBytes] = 116.680 real 8.70 user 32.85 sys 0.53