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.2831853071795862 6.2831853071795862 6.2831853071795862 Vbox= 248.05021344239853 input_persist_forcing: 0.0000000000000000 0.0000000000000000 0.0000000000000000 input_persist_forcing: -10.000000000000000 rsnap: read snapshot var.dat in 3.0000000006111804E-003 seconds 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.5449077018110318 units_eos: unit_temperature= 1.2027220936797545E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 0.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 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) pencil-criteria_magnetic: box divided along z dirn 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.2683404480799516E-005 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 1.0922666666666667E-002 initialize_particles_map: particle_mesh = cic sh: 1: Syntax error: Bad fd number pencil_consistency_check: checking pencil case At line 544 of file particles_tracers.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range' At line 544 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: checking dependence on pencil initialization pencil_consistency_check: results are independent of pencil initialization At line 8234 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' At line 8236 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' At line 8280 of file magnetic.f90 Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2' At line 8282 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: 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$ 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> dlnrho_dt: SOLVE Bcs for lnrho, x: < p>, y: < p>, z: < p> 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 SOLVE dlncc_dt Bcs for cc1, x: < p>, y: < p>, z: < p> dlncc_dt: pscalar_diff= 1.0000000000000000E-003 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.00E-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.0502 3.607E-05 2.088E-05 2.502E-05 0.995185 3.285E+00 3.335E+00 3.298E+00 2 0.080 4.00E-02 3.502E-02 2.476E-02 1.233E-03 2.537E-02 6.74E-02 7.630E-02 1.079E-01 3.653E-03 2.81E-10 -1.04E-09 1.32E-09 2.18E-08 1.89E-03 248.0502 3.494E-05 2.016E-05 2.412E-05 0.995173 3.285E+00 3.335E+00 3.299E+00 4 0.160 4.00E-02 6.650E-02 3.495E-02 1.132E-03 2.277E-02 6.11E-02 1.063E-01 2.023E-01 3.320E-03 2.49E-10 -8.87E-10 1.14E-09 1.87E-08 3.71E-03 248.0502 3.388E-05 1.948E-05 2.329E-05 0.995412 3.285E+00 3.335E+00 3.299E+00 6 0.240 4.00E-02 9.367E-02 4.273E-02 1.042E-03 2.048E-02 5.56E-02 1.328E-01 2.909E-01 3.020E-03 2.22E-10 -7.59E-10 9.82E-10 1.63E-08 5.67E-03 248.0502 3.288E-05 1.887E-05 2.254E-05 0.995766 3.285E+00 3.335E+00 3.299E+00 8 0.320 4.00E-02 1.072E-01 4.924E-02 9.610E-04 1.847E-02 5.06E-02 1.528E-01 3.334E-01 2.769E-03 1.99E-10 -6.54E-10 8.54E-10 1.46E-08 7.51E-03 248.0502 3.189E-05 1.830E-05 2.181E-05 0.996109 3.285E+00 3.335E+00 3.299E+00 Simulation finished after 10 time-steps Writing final snapshot at time t = 0.39999999999999991 Wall clock time [hours] = 2.097E-03 (+/- 2.7778E-07) Wall clock time/timestep/(meshpoint+particle) [microsec] = 12.0 Maximum used memory per cpu [MBytes] = 29.008 Maximum used memory [MBytes] = 29.008 real 12.57 user 12.16 sys 0.09