SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.2 (2019/02/11 21:08:12) wlyra SVN: run.in,v v. 1.2 (2019/02/11 21:08:12) 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= 32 32 1 Lx, Ly, Lz= 2.1000000000000001 6.2831840000000003 9.6000000000000002E-002 Vbox= 19.132295280000001 rsnap: read snapshot var.dat in 2.8569698333740234E-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 -4.6051701859880909 0.10000000000000001 1.0000000000000004E-002 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_gravity: non-smoothed newtonian gravity 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.9132295279999999E-002 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 1.0240000000000000 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' 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) = 3.8002779348098388 duu_dt: SOLVE Bcs for ux, x: < a: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) = 5.7004169022147577 duu_dt: max(advec_uu) = 20.131689030073787 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < p>, z: < p> dlnrho_dt: diffrho_shock= 10.000000000000000 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 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 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' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1378 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1345 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1383 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1350 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' ---it-------t----------dt----------umax--------rhomin--------rhom--------rhomax-------nshift--------rho2m--------uxuym--------urms---------orms--------rhopmax-------npmax------nparmax----nparmin------dtc-------ux2m-------uy2m---- 0 0.0000 9.309E-03 1.581E+00 1.000E+00 1.000E+00 1.000E+00 0.000E+00 1.000E+00 0.000E+00 8.305E-01 4.184E-01 5.581E+00 7.000E+00 2.5600E+02 2.4400E+02 4.5369E-02 0.0000E+00 6.8966E-01 At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 679 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 684 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 646 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 651 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1345 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1378 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1383 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1350 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' 10 0.0836 6.941E-03 1.583E+00 9.693E-01 9.999E-01 1.024E+00 0.000E+00 9.999E-01 2.490E-04 8.305E-01 4.184E-01 5.628E+00 7.000E+00 2.5500E+02 2.4400E+02 3.3825E-02 3.9658E-05 6.8961E-01 At line 646 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 679 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 684 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 651 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' 20 0.1623 6.680E-03 1.585E+00 9.294E-01 9.996E-01 1.047E+00 0.000E+00 9.994E-01 4.604E-04 8.305E-01 4.188E-01 5.558E+00 7.000E+00 2.5500E+02 2.3700E+02 3.2555E-02 9.8405E-05 6.8960E-01 30 0.2407 8.499E-03 1.586E+00 8.988E-01 9.991E-01 1.056E+00 0.000E+00 9.987E-01 6.563E-04 8.305E-01 4.195E-01 5.433E+00 7.000E+00 2.6000E+02 2.3900E+02 4.1423E-02 1.4416E-04 6.8961E-01 40 0.3192 8.326E-03 1.586E+00 8.788E-01 9.983E-01 1.075E+00 0.000E+00 9.975E-01 8.428E-04 8.305E-01 4.206E-01 5.277E+00 6.000E+00 2.6100E+02 2.3900E+02 4.0578E-02 1.7759E-04 6.8962E-01 50 0.3957 8.149E-03 1.586E+00 8.589E-01 9.974E-01 1.090E+00 0.000E+00 9.962E-01 1.019E-03 8.306E-01 4.220E-01 5.184E+00 6.000E+00 2.5800E+02 2.4400E+02 3.9715E-02 2.0260E-04 6.8963E-01 60 0.4738 7.971E-03 1.586E+00 8.393E-01 9.964E-01 1.106E+00 0.000E+00 9.946E-01 1.198E-03 8.306E-01 4.235E-01 5.028E+00 7.000E+00 2.5400E+02 2.4500E+02 3.8846E-02 2.2299E-04 6.8965E-01 70 0.5502 7.833E-03 1.586E+00 8.248E-01 9.953E-01 1.113E+00 0.000E+00 9.929E-01 1.367E-03 8.306E-01 4.251E-01 4.932E+00 7.000E+00 2.5500E+02 2.4500E+02 3.8175E-02 2.4033E-04 6.8965E-01 80 0.6273 7.750E-03 1.586E+00 8.171E-01 9.941E-01 1.123E+00 0.000E+00 9.911E-01 1.532E-03 8.306E-01 4.266E-01 4.808E+00 7.000E+00 2.5300E+02 2.4800E+02 3.7768E-02 2.5491E-04 6.8965E-01 90 0.7055 7.692E-03 1.587E+00 8.030E-01 9.930E-01 1.131E+00 0.000E+00 9.892E-01 1.695E-03 8.306E-01 4.282E-01 4.631E+00 7.000E+00 2.5200E+02 2.4600E+02 3.7486E-02 2.6789E-04 6.8964E-01 Simulation finished after 91 time-steps Writing final snapshot at time t = 0.71316974477638306 Wall clock time [hours] = 8.223E-04 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 16.1 Maximum used memory per cpu [MBytes] = 9.340 Maximum used memory [MBytes] = 36.762 real 5.43 user 3.85 sys 1.60