Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified 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 WARNING: run: quad precision suppressed! 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 1.3813000000000001E-004 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] = 2.394E-06 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 33.1 Maximum used memory per cpu [MBytes] = 25.754 Maximum used memory [MBytes] = 51.270 real 0.64 user 0.42 sys 0.20