SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.1 (2011-03-24 20:03:33) wlyra
SVN: run.in,v v. 1.1 (2011-03-24 20:03:33) wlyra
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: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
The verbose level is ip= 20 (ldebug= F )
This is a 1-D run
nxgrid, nygrid, nzgrid= 1 1 16
Lx, Ly, Lz= 0.40000000000000002 1.6000000000000001 0.80000000000000004
Vbox= 0.80000000000000004
rsnap: read snapshot var.dat in 7.1017000000000001E-005 seconds
setup_slices: slice_position = p
setup_slices: ix_loc,iy_loc, (video files) = 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= 3.4363488390850127E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -3.6888794541139358 0.10000000000000001 7.1428571428571461E-003 0.28571428571428570
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.!
WARNING: initialize_gravity: You do not need gravity_simple for zero gravity!
initialize_gravity: no x-gravity
initialize_gravity: no y-gravity
initialize_gravity: no z-gravity
diffusion: mesh hyperdiffusion
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
select_eos_variable: Using rho and ss
initialize_energy: nheatc_max,iheatcond= 4 hyper3-mesh nothing nothing nothing
heat conduction: hyperdiffusivity of ss
viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq
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 = 8.0000000000000004E-004
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 1.6000000000000000E-002
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
pencil_consistency_check: checking pencil case
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=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
dvvp_dt: Add Coriolis force; Omega= 1.0000000000000000
dvvp_dt: No gravity in x-direction.
dvvp_dt: Linear gravity field in z-direction.
dvvp_dt: No radial 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 rho1 ( 14) 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 ee ( 55) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil lnTT ( 57) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil TT1 ( 64) is requested, but does not appear to be required!
pencil_consistency_check: possible overcalculation... pencil gg ( 91) 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
pencil_consistency_check: all tests passed
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_energy: max(advec_cs2) = 0.57794041309928257
WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh!
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>
coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000
duu_dt: max(advec_uu) = 0.0000000000000000
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < p>, z: < p>
dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < p>, y: < p>, z: < p>
denergy_dt: lnTT,cs2,cp1= -3.6888794541139358 1.0000000000000002E-002 1.1728390000000000
calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000
calc_heatcond_hyper3: added thdiff
----it--------t---------dt-------dtv-------dtu-------dtc-----nparmax-----zpm-------zp2m---
0 0.000 4.16E-02 1.00E+00 0.00E+00 7.79E-01 5.00E+00 1.23E-02 3.94E-02
1 0.042 4.16E-02 1.00E+00 4.48E-06 7.79E-01 5.00E+00 1.23E-02 3.93E-02
2 0.083 4.16E-02 1.00E+00 1.76E-05 7.79E-01 5.00E+00 1.22E-02 3.91E-02
3 0.125 4.16E-02 1.00E+00 3.90E-05 7.79E-01 5.00E+00 1.22E-02 3.88E-02
4 0.166 4.16E-02 1.00E+00 6.80E-05 7.79E-01 5.00E+00 1.21E-02 3.83E-02
5 0.208 4.16E-02 1.00E+00 1.04E-04 7.79E-01 5.00E+00 1.20E-02 3.78E-02
6 0.249 4.16E-02 1.00E+00 1.47E-04 7.79E-01 5.00E+00 1.19E-02 3.72E-02
7 0.291 4.15E-02 1.00E+00 1.95E-04 7.79E-01 5.00E+00 1.18E-02 3.64E-02
8 0.332 4.15E-02 1.00E+00 2.51E-04 7.79E-01 5.00E+00 1.17E-02 3.56E-02
9 0.374 4.15E-02 1.00E+00 3.15E-04 7.79E-01 5.00E+00 1.15E-02 3.47E-02
Simulation finished after 10 time-steps
Writing final snapshot at time t = 0.41550737663486215
Wall clock time [hours] = 1.287E-06 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 17.8
Maximum used memory per cpu [MBytes] = 25.805
Maximum used memory [MBytes] = 51.594
real 0.38
user 0.19
sys 0.06