Pencil1: iapp, nprocs, ncpus= 0 2 2 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 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 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 initialize_gravity: non-smoothed newtonian gravity diffusion: mesh hyperdiffusion diffusion: shock diffusion 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.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 1.8367007289947369E-003 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 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' pde: ENTER SVN: ------- v. ( ) $Id$ At line 2001 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 2001 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' calc_pencils_hydro: call gij_etc calc_pencils_energy: max(advec_cs2) = 899.09316376215554 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.1863275243111 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.9012160932860207 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 dspecial_dt: SOLVE dspecial_dt At line 2001 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 2001 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.0819 2.823E-03 2.701 0.210 0.6226 3.959 0.6981 1.462E-02 7.671E-01 2.37E+01 7.875 6 527 473 2.1164E-01 4.227E-03 5.8299E-01 0.00123 50 0.1489 2.672E-03 5.516 0.198 0.6220 3.953 0.7142 4.268E-02 8.422E-01 4.24E+01 7.111 6 524 473 2.0028E-01 1.035E-02 6.9754E-01 0.00144 75 0.2145 2.779E-03 5.269 0.196 0.6211 3.927 0.7251 7.352E-02 9.053E-01 5.24E+01 6.363 6 522 473 2.0834E-01 1.968E-02 7.9786E-01 0.00200 Simulation finished after 76 time-steps Writing final snapshot at time t = 0.21726431363288390 Wall clock time [hours] = 5.751E-04 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 15.4 Maximum used memory per cpu [MBytes] = 9.680 Maximum used memory [MBytes] = 18.980 real 3.81 user 2.97 sys 0.98