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: ------- 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$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
The verbose level is ip= 14 (ldebug= F )
This is a 3-D run
nxgrid, nygrid, nzgrid= 32 32 32
Lx, Ly, Lz= 0.128000006 0.128000006 0.128000006
Vbox= 2.09715241E-03
WARNING: input_persist_interstellar: t_next_SNI from run.in overwritten. Set l_persist_overwrite_tSNI=T to update!
WARNING: input_persist_interstellar: t_next_SNII from run.in overwritten. Set l_persist_overwrite_tSNII=T to update!
WARNING: input_persist_interstellar: x_cluster from run.in overwritten. Set l_persist_overwrite_xcluster=T to update!
WARNING: input_persist_interstellar: y_cluster from run.in overwritten. Set l_persist_overwrite_ycluster=T to update!
WARNING: input_persist_interstellar: z_cluster from run.in overwritten. Set l_persist_overwrite_zcluster=T to update!
WARNING: input_persist_interstellar: t_cluster from run.in overwritten. Set l_persist_overwrite_tcluster=T to update!
WARNING: input_persist_interstellar: lSNI from run.in overwritten. Set l_persist_overwrite_lSNI=T to update!
input_persist_interstellar: lSNI = F t_next_SNI = 0.00000000
WARNING: input_persist_interstellar: lSNII from run.in overwritten. Set l_persist_overwrite_lSNII=T to update!
input_persist_interstellar: lSNII = T t_next_SNII = 0.00000000
input_persist_forcing: 0.00000000 0.00000000 0.00000000
input_persist_forcing: -10.0000000
rsnap: read snapshot var.dat in 0.0000000000000000 seconds
units_general: unit_velocity= 97781.229948950713
units_general: unit_density= 1.6728000000000004E-024
units_general: unit_length= 3.0856775800000002E+021
units_general: unit_magnetic= 3.5449078083038330
unit temperature 1.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
initialize_eos: ENTER
initialize_gravity: no x-gravity
initialize_gravity: no y-gravity
initialize_gravity: constant gravz= 0.00000000
diffusion: nothing (i.e. no mass diffusion)
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
initialize_energy: cs2top,cs2cool,cs2cool2= 3.90849994E+37 0.00000000 0.00000000
initialize_energy: now set cs2cool=cs2top
initialize_energy: nheatc_max,iheatcond= 4 nothing nothing nothing nothing
heat conduction: nothing
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization.
The values provided are for one-atomic ideal gas. Use at own risk!
resistivity: constant eta
WARNING: initialize_magnetic: Resistivity coefficient eta is zero!
lalpha_aniso= F
lalpha_tensor= F
leta_tensor= F
initialize_magn_mf: lweyl_gauge= F
calc_angle_weights: correction_factor = 0.333333343
initialize_radiation: ndir = 2
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization.
The values provided are for one-atomic ideal gas. Use at own risk!
chiral_diffXX, chiral_diff, chiral_diffZZ = 0.00000000 0.00000000 0.00000000
initialize_dustvelocity: dust_chemistry = nothing
recalculated: md0= 1.00000000
initialize_dustvelocity: minmax(ad)= 0.620350480 0.620350480
initialize_dustvelocity: no betad calculation for draglaw=epstein_cst.
initialize_dustvelocity: dust geometry = sphere
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization.
The values provided are for one-atomic ideal gas. Use at own risk!
initialize_dustdensity: ldustcoagulation,ldustcondensation = F F
initialize_dustdensity: ldustcoagulation_simplified,ldustcondensation_simplified = F F
advec_ddensity: plain vanilla scheme
dust diffusion: nothing
gammacr1= 0.333333373
WARNING: initialize_cosmicrayflux: parameter tau was set to zero!
initialize_interstellar: t_next_SNI = 0.00000000E+00
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization.
The values provided are for one-atomic ideal gas. Use at own risk!
initialize_interstellar: mu not set by eos: 0.620000005
initialize_interstellar: eampl_SN, kampl_SN = 0.0000E+00 0.0000E+00
select_cooling: WSW cooling fct
initialize_interstellar: unit_Lambda 0.1811E+18
initialize_interstellar: unit_Gamma 0.3030E-06
WARNING: initialize_interstellar: using T/e instead of cv1 for diagnostics. Not yet implemented cv1 for ionization!
initialize_interstellar: t_next_SNI, t_next_SNII= 0.00000000E+00 0.00000000E+00
initialize_shock: prenormalised shock_factor sum= 1.00000000
WARNING: initialize_pointmasses: the masses do not sum up to one!
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization.
The values provided are for one-atomic ideal gas. Use at own risk!
initialize_particles: dust-to-gas ratio eps_dtog= 1.00000000
initialize_particles: mass per constituent particle mpmat= 0.00000000
initialize_particles: mass per superparticle mp_swarm = 4.01060570E-05
initialize_particles: number density per superparticle np_swarm= 0.00000000
initialize_particles: mass density per superparticle rhop_swarm= 626.657043
initialize particles: set minimum collisional time-scale equal to two times the Courant time-step.
initialize_particles_map: selected cloud-in-cell for particle-mesh method.
initialize_particles_radius: mass per dust grain mpmat= 0.00000000
WARNING: get_gamma_etc: gamma, cp, and cv are not constant in eos_ionization.
The values provided are for one-atomic ideal gas. Use at own risk!
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
dxxp_dt: Calculate dxxp_dt
dxxp_dt: Particles boundary condition bcpx=p
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
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
At line 273 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range'
At line 273 of file particles_selfgravity.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range'
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'
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: 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: 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
pencil_consistency_check: all tests passed
pde: ENTER
SVN: ------- v. ( ) $Id$
check_SN: ENTER
explode_SN: step, time = 1 0.00000E+00
calc_pencils_energy: max(advec_cs2) = 40506876.0
duu_dt: SOLVE
Bcs for ux, x: < she>, y: < p>, z: < p>
Bcs for uy, x: < she>, y: < p>, z: < p>
Bcs for uz, x: < she>, y: < p>, z: < p>
duu_dt: max(advec_uu) = 0.00000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < she>, y: < p>, z: < p>
dlnrho_dt: max(diffus_diffrho ) = 0.00000000
dlnrho_dt: max(diffus_diffrho3) = 0.00000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < she>, y: < p>, z: < p>
denergy_dt: lnTT,cs2,cp1= 9.21043491 216.036697 69.6056442
calc_heat_cool_interstellar: ENTER
daa_dt: SOLVE
Bcs for Ax, x: < she>, y: < p>, z: < p>
Bcs for Ay, x: < she>, y: < p>, z: < p>
Bcs for Az, x: < she>, y: < p>, z: < p>
daa_dt: iresistivity=eta-const
daa_dt_meanfield: SOLVE
daa_dt_meanfield: max(diffus_eta) = 0.00000000
daa_dt: max(diffus_eta) = 0.00000000
daa_dt: max(diffus_eta2) = 0.00000000
daa_dt: max(diffus_eta3) = 0.00000000
daatest_dt: SOLVE
Bcs for Axtest, x: < she>, y: < p>, z: < p>
Bcs for Aztest, x: < she>, y: < p>, z: < p>
SOLVE dlncc_dt
Bcs for cc, x: < she>, y: < p>, z: < p>
duud_dt: SOLVE duud_dt
Bcs for udx, x: < she>, y: < p>, z: < p>
Bcs for udy, x: < she>, y: < p>, z: < p>
Bcs for udz, x: < she>, y: < p>, z: < p>
dndmd_dt: SOLVE dnd_dt, dmd_dt, dmi_dt
Bcs for nd, x: < she>, y: < p>, z: < p>
SOLVE decr_dt
Bcs for ecr, x: < she>, y: < p>, z: < p>
decr_dt: no diffusion
decr_dt: max(diffus_cr) = 0.00000000
dfcr_dt: SOLVE
Bcs for Fecx, x: < she>, y: < p>, z: < p>
Bcs for Fecy, x: < she>, y: < p>, z: < p>
Bcs for Fecz, x: < she>, y: < p>, z: < p>
SOLVE dXY_dt
Bcs for XX_chiral, x: < she>, y: < p>, z: < p>
Bcs for YY_chiral, x: < she>, y: < p>, z: < p>
dXY_chiral_dt: max(diffus_chiral) = 0.00000000
shearing: Sshear,Sshear1= -0.00000000 -0.00000000
shearing: qshear,qshear0= 1.00000000 0.00000000
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
dnpswarm_dt_pencil: Calculate dnpswarm_dt
--it-------t--------dt-------urms----ekintot----ethtot-----rhom------ssm------TTmin-----TTmax------TTm------dtc-----dtu-----dtnu---dtchi--
0 0.0000E+00 1.41E-04 0.000E+00 0.00E+00 2.452E+00 1.000E+00 7.546E-01 1.000E+04 1.00E+04 1.000E+04 1.0000 0.0000 0.0000 0.0000
explode_SN: step, time = 2 0.14141E-03
1 1.4141E-04 1.43E-04 2.811E-02 8.28E-07 2.355E+00 1.000E+00 7.493E-01 7.628E+03 7.63E+03 7.628E+03 0.8819 0.1181 0.0000 0.0000
2 2.8421E-04 1.30E-04 5.656E-02 3.54E-06 2.319E+00 1.000E+00 7.469E-01 6.605E+03 7.18E+03 6.744E+03 0.7816 0.2184 0.0000 0.0000
3 4.1470E-04 1.22E-04 7.392E-02 6.63E-06 2.301E+00 1.000E+00 7.456E-01 5.977E+03 7.18E+03 6.288E+03 0.7312 0.2687 0.0000 0.0000
explode_SN: step, time = 5 0.53674E-03
4 5.3674E-04 1.21E-04 8.179E-02 8.55E-06 2.289E+00 1.000E+00 7.447E-01 5.576E+03 7.03E+03 6.004E+03 0.7196 0.2812 0.0000 0.0000
Simulation finished after 5 time-steps
Writing final snapshot at time t = 6.5811350214062259E-004
Wall clock time [hours] = 3.646E-04 (+/- 2.7778E-07)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 7.77
Maximum used memory per cpu [MBytes] = 77.844
Maximum used memory [MBytes] = 77.844
real 4.20
user 4.14
sys 0.04