Pencil1: iapp, nprocs, ncpus= 0 2 2 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 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 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 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' initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000E-002 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 8.0000000000000004E-004 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 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 pencil_consistency_check: checking pencil case At line 2014 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 2014 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 At line 2014 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 2014 of file particles_dust_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' 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 ( 12) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil rho1 ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 15) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 53) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil lnTT ( 55) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil TT1 ( 62) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv ( 70) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil gg ( 89) 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.0000000000000000E+000 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.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 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.0000000000000000E+000 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] = 6.233E-04 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 8.631E+03 Maximum used memory per cpu [MBytes] = 8.250 Maximum used memory [MBytes] = 16.152 real 17.43 user 5.66 sys 10.54