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: start.in,v v. 1.1 (2013-06-17 12:04:56) wlyra SVN: run.in,v v. 1.1 (2013-06-17 12:04:56) wlyra 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: gravity_r.f90,v v. 1.1 (2018/08/24 15:48:10) wlyra 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$ The verbose level is ip= 14 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 96 1 8 Lx, Ly, Lz= 2.1000000000000001 6.2831840000000003 9.6000000000000002E-002 Vbox= 1.8367007289947368 rsnap: read snapshot var.dat in 2.2848500000000001E-004 seconds setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 11 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 -2.3025848780656686 0.31622780000000000 0.10000002149284001 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: non-smoothed newtonian gravity diffusion: mesh hyperdiffusion diffusion: shock diffusion WARNING: initialize_density: For diffusion energy/momentum correction should use lmassdiff_fix=T! initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. WARNING: initialize_energy: llocal_iso=T. Make sure you have the appropriate INITIAL_CONDITION in Makefile.local! select_eos_variable: Using rho and cs2 initialize_shock: prenormalised shock_factor sum= 1.0000000000000000 viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq viscous force: nu_shock*(XXXXXXXXXXX) viscous force: nu*del2v 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 = 1.8367007289947369E-003 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 0.76800000000000002 initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method. particles_initialize_modules: reblocking particles At line 143 of file particles_map_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'ineargrid' of procedure 'get_brick_index' At line 143 of file particles_map_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'ineargrid' of procedure 'get_brick_index' sh: 1: Syntax error: Bad fd number pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 899.09316376215543 duu_dt: SOLVE Bcs for ux, x: < out>, y: < p>, z: < p> Bcs for uy, x: < a2>, y: < p>, z: < p> Bcs for uz, x: < s>, y: < p>, z: < p> dss_dt: max(advec_cs2) = 1798.1863275243109 duu_dt: max(advec_uu) = 0.24609708214271162 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < p>, z: < p> dlnrho_dt: diffrho_shock= 10.000000000000000 dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000 dlnrho_dt: max(diffus_diffrho ) = 1.9012160932860203 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 dspecial_dt: SOLVE dspecial_dt At line 2018 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 2018 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 2018 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 2018 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' dxxp_dt: Calculate dxxp_dt dxxp_dt: Particles boundary condition bcpx=rmv 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: No gravity in x-direction. dvvp_dt: No gravity in z-direction. dvvp_dt: Newtonian gravity from a fixed central object ---it-------t---------dt------umax--rhomin---rhom--rhomax--rho2m-----uxuym------urms-------orms---rhopmax-npmax-nparmax-nparmin-----dtc--------ux2m--------uy2m------uz2m--- 0 0.0000 1.227E-02 1.534 0.253 0.6231 3.953 0.6901 1.538E-05 7.346E-01 3.67E-01 8.298 6 527 473 9.1971E-01 9.687E-07 5.3967E-01 0.00000 25 0.0838 3.063E-03 1.534 0.208 0.6226 3.959 0.6989 1.143E-02 7.462E-01 2.25E+00 7.822 6 527 473 2.2963E-01 4.144E-03 5.5156E-01 0.00119 50 0.1602 3.098E-03 1.534 0.195 0.6218 3.951 0.7238 2.010E-02 7.551E-01 2.84E+00 6.700 6 524 473 2.3225E-01 7.049E-03 5.6217E-01 0.00103 75 0.2410 3.415E-03 1.534 0.191 0.6202 4.257 0.7586 2.753E-02 7.632E-01 3.04E+00 5.391 6 518 473 2.5597E-01 9.129E-03 5.7247E-01 0.00094 Simulation finished after 76 time-steps Writing final snapshot at time t = 0.24439575184057613 Wall clock time [hours] = 1.824E-04 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 4.89 Maximum used memory per cpu [MBytes] = 27.234 Maximum used memory [MBytes] = 54.281 real 1.21 user 1.59 sys 0.18