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
rsnap: read snapshot var.dat in 1.0844800000000000E-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.155E-04 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 3.09
Maximum used memory per cpu [MBytes] = 27.387
Maximum used memory [MBytes] = 54.461
real 0.78
user 1.00
sys 0.06