SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) no cvsid is given i SVN: ------- v. ( ) no cvsid is given i 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$ SVN: particles_density. v. 20849 (2013-08-06 18:45:43) anders SVN: ------- v. ( ) $Id$ SVN: particles_dust.f90 v. 19206 (2012-06-30 21:40:24) sven.bingert SVN: ------- v. ( ) $Id$ The verbose level is ip= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 8 8 8 Lx, Ly, Lz= 0.20000000000000001 0.20000000000000001 0.20000000000000001 Vbox= 8.0000000000000019E-003 rsnap: read snapshot var.dat in 5.3191184997558594E-004 seconds setup_slices: slice_position = m setup_slices: ix_loc, (video files) = 7 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 3.6757541328186911 6.2831853071795862 39.478417604357439 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 Set z stratification: g_z = -gz_coeff^2 * z initialize_selfgravity: epicycle frequency kappa = 6.2831853071795862 WARNING: initialize_density_methods: density methods not yet implemented for density_stratified! initialize_density: mesh hyper-diffusion; diffrho_hyper3_mesh = 1.0000000000000000 initialize_density: Global density gradient with beta_glnrho_global= -0.10000000000000001 0.0000000000000000 0.0000000000000000 select_eos_variable: Using rho and cs2 viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq initialize_particles: dust-to-gas ratio eps_dtog= 2.0000000000000000E-002 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000 initialize_particles: mass per superparticle mp_swarm = 3.9166066791109385E-006 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000 initialize_particles: mass density per superparticle rhop_swarm= 0.25066282746310004 initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method. initialize_particles_sink: sink_radius= 2.5000000000000001E-002 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' 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 pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_energy: max(advec_cs2) = 189496.40450091567 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < p> Bcs for uy, x: < she>, y: < p>, z: < p> Bcs for uz, x: < she>, y: < p>, z: < p> coriolis_cartesian: add Coriolis force; Omega= 6.2831853071795862 duu_dt: max(advec_uu) = 12.222780051067222 dlnrho_dt: SOLVE the continuity equation. Bcs for rho, x: < she>, y: < p>, z: < p> denergy_dt: adding global pressure gradient force shearing: Sshear,Sshear1= -9.4247779607693793 -9.4247779607693793 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000 drhopswarm_dt_pencil: Calculate drhopswarm_dt 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' 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' At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 314 of file particles_sink.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= 6.2831853071795862 At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' dvvp_dt: No gravity in x-direction. dvvp_dt: Linear gravity field in z-direction. dvvp_dt: No radial gravity At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 314 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' ---it------t-------dt-----dtv-dtdragp----mass------urms------umax-----rhomin----rhomax---drhorms---drhomax--nblockmin-nblockmax--nparmin---nparmax------xpm-------xp2m-------ypm-------yp2m-------zpm-------zp2m------vpxm-------vpym-------vpzm------vpx2m-----vpy2m-----vpz2m-----npmax-----rhopm----rhopmax-- 0 0.000 1.66E-03 1.00 0.2498 7.987E-03 2.583E-01 2.974E-01 9.962E-01 9.999E-01 0.000E+00 0.000E+00 5.300E+01 6.200E+01 1.230E+02 1.300E+02 -1.828E-03 3.411E-03 -7.093E-04 3.284E-03 3.302E-03 3.394E-03 -1.097E-01 -2.272E-01 0.000E+00 1.220E-02 5.178E-02 0.000E+00 4.000E+00 2.507E-01 5.446E-01 At line 646 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 679 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 684 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 651 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' 1 0.002 1.66E-03 1.00 0.2498 7.987E-03 2.582E-01 2.947E-01 9.931E-01 1.003E+00 1.152E-03 3.321E-03 5.400E+01 6.200E+01 1.230E+02 1.310E+02 -1.620E-03 3.412E-03 -1.064E-03 3.279E-03 3.302E-03 3.394E-03 -1.097E-01 -2.274E-01 -2.128E-04 1.222E-02 5.185E-02 1.413E-05 4.000E+00 2.507E-01 5.370E-01 2 0.003 1.66E-03 1.00 0.2496 7.987E-03 2.580E-01 2.920E-01 9.914E-01 1.004E+00 1.608E-03 4.837E-03 5.400E+01 6.300E+01 1.230E+02 1.330E+02 -1.802E-03 3.412E-03 -2.196E-03 3.274E-03 3.302E-03 3.392E-03 -1.098E-01 -2.276E-01 -4.169E-04 1.223E-02 5.194E-02 5.474E-05 4.000E+00 2.507E-01 5.355E-01 At line 1036 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1041 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 975 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 980 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1345 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1378 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1383 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1350 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 646 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 679 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 684 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 651 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' 3 0.005 1.66E-03 1.00 0.2495 7.987E-03 2.580E-01 2.933E-01 9.917E-01 1.003E+00 1.352E-03 4.541E-03 5.400E+01 6.400E+01 1.210E+02 1.320E+02 -1.985E-03 3.412E-03 -2.154E-03 3.269E-03 3.301E-03 3.391E-03 -1.099E-01 -2.278E-01 -6.110E-04 1.225E-02 5.204E-02 1.193E-04 4.000E+00 2.507E-01 5.338E-01 4 0.007 1.66E-03 1.00 0.2494 7.987E-03 2.580E-01 2.959E-01 9.926E-01 1.002E+00 1.016E-03 3.559E-03 5.500E+01 6.500E+01 1.210E+02 1.340E+02 -1.776E-03 3.413E-03 -2.526E-03 3.264E-03 3.300E-03 3.388E-03 -1.100E-01 -2.281E-01 -7.949E-04 1.227E-02 5.215E-02 2.054E-04 4.000E+00 2.507E-01 5.376E-01 5 0.008 1.66E-03 1.00 0.2493 7.987E-03 2.579E-01 2.978E-01 9.927E-01 1.002E+00 1.068E-03 3.488E-03 5.600E+01 6.500E+01 1.190E+02 1.360E+02 -1.959E-03 3.413E-03 -2.485E-03 3.262E-03 3.298E-03 3.385E-03 -1.102E-01 -2.283E-01 -9.695E-04 1.230E-02 5.226E-02 3.108E-04 4.000E+00 2.507E-01 5.443E-01 At line 1345 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' At line 1378 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1383 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2' At line 1350 of file particles_mpicomm_blocks.f90 Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2' 6 0.010 1.66E-03 1.00 0.2494 7.987E-03 2.578E-01 2.965E-01 9.931E-01 1.003E+00 1.132E-03 3.652E-03 5.600E+01 6.700E+01 1.180E+02 1.370E+02 -1.751E-03 3.414E-03 -2.478E-03 3.254E-03 3.296E-03 3.381E-03 -1.103E-01 -2.285E-01 -1.137E-03 1.232E-02 5.235E-02 4.335E-04 4.000E+00 2.507E-01 5.502E-01 7 0.012 1.66E-03 1.00 0.2913 7.987E-03 2.578E-01 2.930E-01 9.920E-01 1.002E+00 1.089E-03 4.226E-03 5.600E+01 6.700E+01 1.190E+02 1.350E+02 -1.934E-03 3.414E-03 -1.266E-03 3.249E-03 3.294E-03 3.377E-03 -1.104E-01 -2.287E-01 -1.300E-03 1.234E-02 5.244E-02 5.715E-04 4.000E+00 2.507E-01 5.554E-01 8 0.013 1.66E-03 1.00 0.2915 7.987E-03 2.578E-01 2.937E-01 9.924E-01 1.002E+00 1.060E-03 3.777E-03 5.700E+01 6.700E+01 1.220E+02 1.370E+02 -2.118E-03 3.415E-03 -1.614E-03 3.244E-03 3.292E-03 3.372E-03 -1.105E-01 -2.289E-01 -1.460E-03 1.236E-02 5.251E-02 7.230E-04 4.000E+00 2.507E-01 5.599E-01 9 0.015 1.66E-03 1.00 0.2914 7.987E-03 2.577E-01 2.908E-01 9.936E-01 1.002E+00 1.084E-03 2.779E-03 5.700E+01 6.700E+01 1.230E+02 1.350E+02 -2.301E-03 3.415E-03 -1.178E-03 3.240E-03 3.289E-03 3.367E-03 -1.106E-01 -2.291E-01 -1.619E-03 1.238E-02 5.258E-02 8.865E-04 4.000E+00 2.507E-01 5.638E-01 Simulation finished after 10 time-steps Writing final snapshot at time t = 1.6598774037744242E-002 Wall clock time [hours] = 2.022E-04 (+/- 2.7778E-10) Wall clock time/timestep/(meshpoint+particle) [microsec] = 71.1 Maximum used memory per cpu [MBytes] = 10.090 Maximum used memory [MBytes] = 39.652 real 3.23 user 1.40 sys 0.63