SVN: ------- v. ( ) $Id$
SVN: start.in v. 12078 (2009-11-14 16:52:27) ajohan
SVN: run.in v. 12078 (2009-11-14 16:52:27) ajohan
bcx1,bcx2= p p p p p p p p p : p p p p p p p p p
bcy1,bcy2= p p p p p p p p p : p p p p p p p p p
bcz1,bcz2= p p p p p p p p p : p p p p p p p p p
initialize_mpicomm: enabled MPI
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: pointmasses.f90,v v. 1.1 (2019/02/02 03:54:41) wlyra
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= 10 (ldebug= F )
This is a 2-D run
nxgrid, nygrid, nzgrid= 64 64 1
Lx, Ly, Lz= 4.2000000000000002 4.2000000000000002 0.52000000000000002
Vbox= 17.640000000000001
rsnap: read snapshot var.dat in 5.0554500000000004E-004 seconds
rprint_list: nname= 18
rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1
rprint_list: nnamev= 0
rprint_list: nnamez= -1
rprint_list: nnamez= -1
rprint_list: nnamey= -1
rprint_list: nnamex= -1
rprint_list: nnamer= 0
rprint_list: nnamexz= -1
rprint_list: nnamexy= -1
rprint_list: nnamerz= -1
rprint_general: run through parse list
rprint_hydro: run through parse list
rprint_hydro: run through parse list
rprint_density: run through parse list
rprint_shock: run through parse list
rprint_viscosity: run through parse list
rprint_pointmasses: run through parse list
rprint_particles: run through parse list
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 -5.9914645471079817 5.0000000000000003E-002 2.5000000000000005E-003 1.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
timestep based on CFL cond; cdt= 0.40000000000000002
dVol_rel1= 2.4414062500000000E-004
dA_xy_rel1= 2.4414062500000000E-004
dA_xz_rel1= 1.5625000000000000E-002
dA_yz_rel1= 1.5625000000000000E-002
dL_y_rel1= 1.5625000000000000E-002
xmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
zmask_den= 1.0000000000000000
diffusion: (d^6/dx^6+d^6/dy^6+d^6/dz^6)rho
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.
xmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
zmask_hyd= 1.0000000000000000
WARNING: initialize_energy: llocal_iso=T. Make sure you have the appropriate INITIAL_CONDITION in Makefile.local!
select_eos_variable: Using rho and cs2
viscous force: nu*(del6u+S.glnrho)
viscous force: nu_shock*(XXXXXXXXXXX)
xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
rprint_pointmasses: run through parse list
choose_pencils: finding out which pencils are needed for the pencil case.
write_pencil_info: pencil information written to the file pencils.list.
initialize_particles: Number of particle species = 4
initialize_particles: tausp_species = 1.0000000000000000E-002 0.10000000000000001 0.29999999999999999 1.0000000000000000
initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002
initialize_particles: mass per constituent particle mpmat= 0.0000000000000000
initialize_particles: mass per superparticle mp_swarm = 1.7226562500000000E-004
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 4.0000000000000001E-002
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.
Lz= 0.52000000000000002
z= 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
life_support_on: needed for pencil consistency check.
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
At line 934 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index'
Entered inverse_laplacian_expandgrid
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline'
At line 282 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline'
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'fp_pt' of procedure 'gravity_pointmasses'
At line 1112 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'dfp_pt' of procedure 'gravity_pointmasses'
dxxp_dt: Calculate dxxp_dt
dxxp_dt: Particles boundary condition bcpx=rmv
dxxp_dt: Particles boundary condition bcpy=rmv
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
dvvp_dt_selfgrav: add self-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: performing full pencil check (takes a while)
pencil_consistency_check: possible overcalculation... pencil rho ( 13) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil grho ( 16) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil phix ( 92) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil phiy ( 93) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil rcyl_mn1 ( 97) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil uij5 ( 114) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil np ( 173) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil rhop ( 174) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil epsp ( 183) is requested, but does not appear to be required!
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: performing full pencil check (takes a while)
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
life_support_off: end of pencil consistency check.
pencil_consistency_check: all tests passed
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_energy: max(advec_cs2) = 0.56156081721680717
duu_dt: SOLVE
Bcs for ux, x: < p>, y: < p>, z: < p>
Bcs for uy, x: < p>, y: < p>, z: < p>
Bcs for uz, x: < p>, y: < p>, z: < p>
duu_dt: max(advec_uu) = 13.532882206920425
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < p>, z: < p>
dlnrho_dt: diffrho_shock= 1.0000000000000000E-002
dlnrho_dt: diffrho_hyper3= 5.0000000000000001E-009
dlnrho_dt: max(diffus_diffrho ) = 1.3369150792887986E-003
dlnrho_dt: max(diffus_diffrho3) = 0.12519416891898374
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
dvvp_dt_pencil: adding feedback from dust to gas
dvvp_dt: Add drag force; tausp= 0.0000000000000000
pde: freezing variables for r < 0.40000000000000002 : T T T T F F F F F
pde: freezing variables for r > 2.0000000000000000 : T T T T F F F F F
----it--------t---------dt-----dtdragp------rhom--------rhomax-------rhomin--------ux2m--------rhopm--------npmax-------rhopmax------nparmax-------vpy2m-------gpotenp--------xq1----------yq1----------xq2----------yq2-----
0 0.000 9.62E-04 1.00E+00 1.000E+00 1.000E+00 1.000E+00 4.1835E-01 1.000E-02 1.600E+01 3.257E-01 2.400E+02 5.081E-01 -4.437E-05 0.999000000 0.000000000 -0.001000000 0.000000000
1 0.001 9.62E-04 1.00E+00 1.000E+00 1.000E+00 1.000E+00 4.1835E-01 1.000E-02 1.600E+01 3.265E-01 2.400E+02 5.080E-01 -4.437E-05 0.998999538 0.000960828 -0.001000000 -0.000000962
2 0.002 9.62E-04 1.00E+00 1.000E+00 1.000E+00 1.000E+00 4.1835E-01 1.000E-02 1.600E+01 3.270E-01 2.480E+02 5.079E-01 -4.437E-05 0.998998152 0.001921752 -0.000999998 -0.000001924
3 0.003 9.62E-04 1.00E+00 1.000E+00 1.000E+00 1.000E+00 4.1835E-01 1.000E-02 1.600E+01 3.271E-01 2.480E+02 5.078E-01 -4.437E-05 0.998995841 0.002882745 -0.000999996 -0.000002886
4 0.004 9.62E-04 1.00E+00 1.000E+00 1.000E+00 9.999E-01 4.1835E-01 1.000E-02 1.600E+01 3.270E-01 2.480E+02 5.077E-01 -4.437E-05 0.998992605 0.003843787 -0.000999993 -0.000003848
5 0.005 9.62E-04 1.00E+00 1.000E+00 1.000E+00 9.999E-01 4.1835E-01 1.000E-02 1.600E+01 3.266E-01 2.480E+02 5.076E-01 -4.437E-05 0.998988445 0.004804864 -0.000999988 -0.000004810
6 0.006 9.62E-04 1.00E+00 1.000E+00 1.000E+00 9.999E-01 4.1835E-01 1.000E-02 1.600E+01 3.258E-01 2.480E+02 5.075E-01 -4.437E-05 0.998983360 0.005765966 -0.000999983 -0.000005772
7 0.007 9.62E-04 1.00E+00 1.000E+00 1.000E+00 9.999E-01 4.1835E-01 1.000E-02 1.600E+01 3.248E-01 2.480E+02 5.074E-01 -4.437E-05 0.998977350 0.006727083 -0.000999977 -0.000006734
8 0.008 9.62E-04 1.00E+00 1.000E+00 1.000E+00 9.999E-01 4.1834E-01 1.000E-02 1.600E+01 3.235E-01 2.480E+02 5.074E-01 -4.436E-05 0.998970416 0.007688211 -0.000999970 -0.000007696
9 0.009 9.62E-04 1.00E+00 1.000E+00 1.000E+00 9.999E-01 4.1834E-01 1.000E-02 1.600E+01 3.219E-01 2.480E+02 5.073E-01 -4.436E-05 0.998962556 0.008649345 -0.000999963 -0.000008658
Simulation finished after 10 time-steps
Writing final snapshot at time t = 9.6202497378347022E-003
wsnap_particles: written snapshot pvar.dat
written pointmass snapshot qvar.dat
wsnap_particles: written snapshot dpvar.dat
Wall clock time [hours] = 1.736E-05 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.22
Maximum used memory per cpu [MBytes] = 28.637
Maximum used memory [MBytes] = 227.223
real 3.43
user 13.51
sys 0.44