SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
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= 16 16 16
Lx, Ly, Lz= 0.20000000000000001 0.20000000000000001 0.20000000000000001
Vbox= 8.0000000000000019E-003
rsnap: read snapshot var.dat in 1.1292100000000000E-004 seconds
setup_slices: slice_position = m
setup_slices: (video files) =
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
initialize_gravity: no x-gravity
initialize_gravity: no y-gravity
initialize_gravity: linear z-grav, nu= 1.0000000000000000
initialize_selfgravity: epicycle frequency kappa = 1.0000000000000000
diffusion: mesh hyperdiffusion
initialize_density: no need to read initial stratification for lanti_shockdiffusion=F.
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: Number of particle species = 2
initialize_particles: tausp_species = 0.50000000000000000 5.0000000000000001E-003
initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000
initialize_particles: mass per constituent particle mpmat= 0.0000000000000000
initialize_particles: mass per superparticle mp_swarm = 1.0026513098524003E-004
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 51.335747064442884
initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method.
initialize_particles_sink: sink_radius= 1.2500000000000001E-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'
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'
pde: ENTER
SVN: ------- v. ( ) $Id$
sh: 1: Syntax error: Bad fd number
calc_pencils_energy: max(advec_cs2) = 18424.999999999996
duu_dt: SOLVE
Bcs for ux, x: < she>, y: < p>, z: < s>
Bcs for uy, x: < she>, y: < p>, z: < s>
Bcs for uz, x: < she>, y: < p>, z: < a>
coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000
duu_dt: max(advec_uu) = 3.9999999999999996
dlnrho_dt: SOLVE
Bcs for lnrho, x: < she>, y: < p>, z: < s>
dlnrho_dt: diffrho_hyper3_mesh= 20.000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: adding global pressure gradient force
shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000
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 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 '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 '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'
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'
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'
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.
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: Calculate dvvp_dt
dvvp_dt: Add Coriolis force; Omega= 1.0000000000000000
dvvp_dt: No gravity in x-direction.
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: 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 565 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 565 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 565 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 565 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 565 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 571 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 619 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpibcast_real_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 720 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 720 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2'
At line 720 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2'
At line 671 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum_tmp' of procedure 'mpireduce_sum_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 720 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2'
At line 672 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'fsum' of procedure 'mpireduce_sum_arr2'
At line 712 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpisend_real_arr2'
At line 720 of file particles_sink.f90
Fortran runtime warning: An array temporary was created for argument 'bcast_array' of procedure 'mpirecv_real_arr2'
------it---------------t-----------------dt-------------nparmax-------------ux2m--------------uy2m--------------uz2m-------------uxuym--------------rhom-------------rhomin------------rhomax-------------vpxm--------------xpm---------------xp2m--------------zpm---------------zp2m-------------npmax-------------rhopm------------rhopmax-----------nparsink---------rhopinterp----
0 0.0000000000E+00 9.7017887701E-06 1.2500000000E+02 0.0000000000E+00 1.6971676641E-03 0.0000000000E+00 0.0000000000E+00 9.9811430117E-01 9.9501247919E-01 9.9997777802E-01 0.0000000000E+00 -4.7901681516E-04 3.3669032108E-03 -2.6921769025E-05 1.0454051865E-04 7.0000000000E+00 1.2533141373E+01 1.7115899214E+02 1.3000000000E+01 1.0975032884E+02
1 9.7017887701E-06 1.2025728300E-05 1.2400000000E+02 2.4005549969E-13 1.6971446894E-03 6.5229730409E-15 -1.3423729407E-09 9.9811430197E-01 9.9501249194E-01 9.9997922369E-01 -1.4931860727E-08 1.7508625222E-03 3.2261225036E-03 -1.1636308946E-05 1.0977078615E-04 7.0000000000E+00 1.1229694670E+01 1.3079069866E+02 9.0000000000E+00 8.4607797207E+01
2 2.1727517070E-05 1.2025730905E-05 1.2400000000E+02 1.1495511718E-12 1.6971177935E-03 3.2784725656E-14 -2.9797012544E-09 9.9811430298E-01 9.9501250808E-01 9.9998102422E-01 -3.2236188132E-08 1.7508625219E-03 3.2261225036E-03 -1.1636308944E-05 1.0977078611E-04 7.0000000000E+00 1.1229694670E+01 1.3078603622E+02 9.0000000000E+00 8.4605935719E+01
3 3.3753247976E-05 1.2025733134E-05 1.2400000000E+02 2.6548996376E-12 1.6970923027E-03 7.9410245293E-14 -4.5908544970E-09 9.9811430401E-01 9.9501252462E-01 9.9998283271E-01 -4.8313360462E-08 1.7508625214E-03 3.2261225036E-03 -1.1636308940E-05 1.0977078604E-04 7.0000000000E+00 1.1229694670E+01 1.3078137384E+02 9.0000000000E+00 8.4604073966E+01
4 4.5778981110E-05 1.2025734985E-05 1.2400000000E+02 4.6845014593E-12 1.6970680889E-03 1.4683449267E-13 -6.1788683375E-09 9.9811430507E-01 9.9501254153E-01 9.9998464770E-01 -6.3265160597E-08 1.7508625208E-03 3.2261225036E-03 -1.1636308934E-05 1.0977078593E-04 7.0000000000E+00 1.1229694670E+01 1.3077671152E+02 9.0000000000E+00 8.4602211949E+01
5 5.7804716095E-05 1.2025736458E-05 1.2400000000E+02 7.1801853399E-12 1.6970450390E-03 2.3566791626E-13 -7.7463416981E-09 9.9811430615E-01 9.9501255882E-01 9.9998646786E-01 -7.7184079144E-08 1.7508625199E-03 3.2261225036E-03 -1.1636308926E-05 1.0977078580E-04 7.0000000000E+00 1.1229694670E+01 1.3077204927E+02 9.0000000000E+00 8.4600349669E+01
6 6.9830452553E-05 1.2025737552E-05 1.2400000000E+02 1.0094732639E-11 1.6970230530E-03 3.4668909398E-13 -9.2955086718E-09 9.9811430726E-01 9.9501257650E-01 9.9998829200E-01 -9.0154209848E-08 1.7508625189E-03 3.2261225035E-03 -1.1636308916E-05 1.0977078563E-04 7.0000000000E+00 1.1229694670E+01 1.3076738709E+02 9.0000000000E+00 8.4598487125E+01
7 8.1856190105E-05 1.2025738272E-05 1.2400000000E+02 1.3389867824E-11 1.6970020421E-03 4.8083842964E-13 -1.0828297619E-08 9.9811430839E-01 9.9501259455E-01 9.9999011901E-01 -1.0225205529E-07 1.7508625177E-03 3.2261225035E-03 -1.1636308905E-05 1.0977078543E-04 7.0000000000E+00 1.1229694670E+01 1.3076272498E+02 9.0000000000E+00 8.4596624320E+01
8 9.3881928377E-05 1.2025738619E-05 1.2400000000E+02 1.7034625624E-11 1.6969819268E-03 6.3921275277E-13 -1.2346379852E-08 9.9811430954E-01 9.9501261298E-01 9.9999194789E-01 -1.1354725215E-07 1.7508625164E-03 3.2261225035E-03 -1.1636308892E-05 1.0977078520E-04 7.0000000000E+00 1.1229694670E+01 1.3075806294E+02 9.0000000000E+00 8.4594761254E+01
9 1.0590766700E-04 1.2025738598E-05 1.2400000000E+02 2.1004023096E-11 1.6969626362E-03 8.2306069248E-13 -1.3851209837E-08 9.9811431072E-01 9.9501263178E-01 9.9999377776E-01 -1.2410322430E-07 1.7508625150E-03 3.2261225035E-03 -1.1636308877E-05 1.0977078493E-04 7.0000000000E+00 1.1229694670E+01 1.3075340097E+02 9.0000000000E+00 8.4592897927E+01
Simulation finished after 10 time-steps
Writing final snapshot at time t = 1.1793340559398254E-004
Wall clock time [hours] = 2.405E-05 (+/- 2.7778E-13)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.70
Maximum used memory per cpu [MBytes] = 28.137
Maximum used memory [MBytes] = 222.582
real 0.47
user 1.22
sys 0.21