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: ------- 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= 6.28318548 6.28318548 6.28318548 Vbox= 248.050232 input_persist_forcing: 0.00000000E+00 0.00000000E+00 0.00000000E+00 input_persist_forcing: -10.0000000 setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 35 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449078083038330 units_eos: unit_temperature= 1.2027220775801925E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 0.00000000E+00 1.00000000 1.00000000 1.00000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using lnrho and cs2 resistivity: constant eta Number of passive scalars = 1 viscous force: nu*(del2u+graddivu/3+2S.glnrho) initialize_particles: dust-to-gas ratio eps_dtog= 9.99999978E-03 initialize_particles: mass per constituent particle mpmat= 0.00000000E+00 initialize_particles: mass per superparticle mp_swarm = 8.26834148E-05 initialize_particles: number density per superparticle np_swarm= 0.00000000E+00 initialize_particles: mass density per superparticle rhop_swarm= 1.09226666E-02 initialize_particles_map: particle_mesh = cic pencil-criteria_magnetic: box divided along z dirn pencil_consistency_check: checking pencil case At line 542 of file particles_tracers.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 542 of file particles_tracers.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range' dxxp_dt: No gravity in z-direction. 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 cc ( 245) 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 At line 7670 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' At line 7672 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' At line 7716 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' At line 7718 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' 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) = 77.8146667 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) = 0.00000000E+00 dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> dlnrho_dt: max(diffus_diffrho ) = 0.00000000E+00 dlnrho_dt: max(diffus_diffrho3) = 0.00000000E+00 daa_dt: SOLVE Bcs for Ax, x: < p>, y: < p>, z: < p> Bcs for Ay, x: < p>, y: < p>, z: < p> Bcs for Az, x: < p>, y: < p>, z: < p> daa_dt: iresistivity=eta-const daa_dt: max(diffus_eta) = 0.389073312 daa_dt: max(diffus_eta2) = 0.00000000E+00 daa_dt: max(diffus_eta3) = 0.00000000E+00 SOLVE dlncc_dt Bcs for cc1, x: < p>, y: < p>, z: < p> dlncc_dt: pscalar_diff= 1.00000005E-03 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 gas velocity. At line 480 of file particles_stalker.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 480 of file particles_stalker.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range' --it-----t--------dt-------umax------urms------brms------jrms------jmax------orms------omax-----vAmax------abm-------abmn------abms------jbm-------oum------mass------bmx-------bmy-------bmz------ccmax------xp2m------yp2m------zp2m--- 0 0.000 4.53E-02 0.000E+00 0.000E+00 1.346E-03 2.832E-02 7.54E-02 0.000E+00 0.000E+00 4.019E-03 3.19E-10 -1.23E-09 1.55E-09 2.64E-08 0.00E+00 248.0478 3.607E-05 2.088E-05 2.502E-05 0.995185 3.285E+00 3.335E+00 3.298E+00 2 0.090 4.37E-02 3.709E-02 2.622E-02 1.220E-03 2.503E-02 6.66E-02 8.083E-02 1.143E-01 3.611E-03 2.77E-10 -1.02E-09 1.30E-09 2.14E-08 2.12E-03 248.0478 3.481E-05 2.007E-05 2.401E-05 0.995176 3.285E+00 3.335E+00 3.298E+00 4 0.177 4.24E-02 6.985E-02 3.671E-02 1.112E-03 2.227E-02 5.99E-02 1.117E-01 2.125E-01 3.256E-03 2.43E-10 -8.58E-10 1.10E-09 1.82E-08 4.10E-03 248.0482 3.366E-05 1.935E-05 2.312E-05 0.995453 3.285E+00 3.335E+00 3.299E+00 6 0.261 4.13E-02 9.761E-02 4.453E-02 1.020E-03 1.993E-02 5.43E-02 1.384E-01 3.030E-01 2.954E-03 2.16E-10 -7.30E-10 9.46E-10 1.57E-08 6.16E-03 248.0501 3.262E-05 1.871E-05 2.234E-05 0.995845 3.285E+00 3.335E+00 3.299E+00 8 0.343 4.08E-02 1.109E-01 5.098E-02 9.391E-04 1.794E-02 4.93E-02 1.581E-01 3.449E-01 2.704E-03 1.93E-10 -6.27E-10 8.20E-10 1.41E-08 8.05E-03 248.0503 3.161E-05 1.814E-05 2.161E-05 0.996214 3.285E+00 3.335E+00 3.299E+00 Simulation finished after 10 time-steps Writing final snapshot at time t = 0.42443747073411942 Wall clock time [hours] = 1.011E-03 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 5.80 Maximum used memory per cpu [MBytes] = 25.367 Maximum used memory [MBytes] = 25.367 real 83.24 user 83.07 sys 0.07