Pencil1: iapp, nprocs, ncpus= 0 8 8 SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "special_init_pars" is missing!! SVN: start.in,v v. 1.1 (2021/11/22 22:50:45) wlyra SVN: run.in,v v. 1.1 (2021/11/22 22:50:45) wlyra 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: particles_dust.f90 v. 1.1 (2018/08/24 15:48:10) wlyra 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= 32 32 32 Lx, Ly, Lz= 0.40000000000000002 1.6000000000000001 0.20000000000000001 Vbox= 0.12800000000000003 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 initialize_selfgravity: epicycle frequency kappa = 1.0000000000000000 diffusion: mesh hyperdiffusion diffusion: shock diffusion 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 shock nothing nothing heat conduction: hyperdiffusivity of ss heat conduction: shock viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq viscous force: nu_shock*(XXXXXXXXXXX) initialize_particles: Global pressure gradient with beta_dPdr_dust= -0.10000000000000001 initialize_particles: dust-to-gas ratio eps_dtog= 2.9999999999999999E-002 initialize_particles: mass per constituent particle mpmat= 0.0000000000000000E+000 initialize_particles: mass per superparticle mp_swarm = 1.6042420957638403E-007 initialize_particles: number density per superparticle np_swarm= 0.0000000000000000E+000 initialize_particles: mass density per superparticle rhop_swarm= 4.1068597651554302E-002 initialize particles: set minimum collisional time-scale equal to two times the Courant time-step. initialize_particles_map: selected triangular-shaped-cloud for particle-mesh method. initialize_particles_sink: sink_radius= 1.2500000000000001E-002 pde: ENTER SVN: ------- v. ( ) $Id$ At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' At line 934 of file particles_map.f90 Fortran runtime warning: An array temporary was created for argument 'x' of procedure 'real_to_index' calc_pencils_energy: max(advec_cs2) = 377.60386608391195 WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh! 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= 1.0000000000000000 duu_dt: max(advec_uu) = 3.9675517980296074E-002 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < p> dlnrho_dt: diffrho_shock= 1.0000000000000000 dlnrho_dt: diffrho_hyper3_mesh= 20.000000000000000 dlnrho_dt: max(diffus_diffrho ) = 5.8188604640949812E-002 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < p> denergy_dt: lnTT,cs2,cp1= -3.5357772971176513 1.1654440311231853E-002 0.0000000000000000E+000 calc_heatcond_shock: use shock diffusion calc_heatcond_shock: chi_shock= 1.0000000000000000 calc_heatcond_shock: added thdiff calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000E+000 calc_heatcond_hyper3: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F F F F dspecial_dt: SOLVE dSPECIAL_dt shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000E+000 dvvp_dt_pencil: calculate dvvp_dt dvvp_dt_pencil: ldraglaw_purestokes= F dvvp_dt_pencil: adding feedback from dust to gas dvvp_dt: Add drag force; tausp= 0.29999999999999999 drhopswarm_dt_pencil: Calculate drhopswarm_dt At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.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 particle_gravity: Calculating gravity particle_gravity: No gravity in x-direction. particle_gravity: Linear gravity field in z-direction. particle_gravity: No radial gravity dvvp_dt_selfgrav: add self-gravity At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 282 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 335 of file particles_selfgravity.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_quadratic_spline' At line 388 of file particles_sink.f90 Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_quadratic_spline' ---it-------t--------dt--------uxm--------uym-------uz2m-------ux2m-------uy2m-------uxuym------urms-------umax-------rhom-------rho2m-----rhomin-----rhomax------orms--------o2m-------ox2m-------oy2m-------oz2m--------ppm-------ethm--------ssm--------csm--------TTm-------TTmax------TTmin------nparmax-------npmax--------rhopm-------rhopmax------gpotenp-----potselfpm------gxgym--------gxgzm--------gygzm-------qtoomre----qtoomremin----nparsink----rhopinterp-- 0 0.000 8.760E-03 3.133E-07 -1.122E-09 1.001E-08 9.980E-09 9.909E-09 -3.378E-11 1.729E-04 4.838E-04 9.965E-01 9.955E-01 8.291E-01 1.144E+00 2.745E-02 7.535E-04 3.031E-04 3.712E-04 7.918E-05 7.143E-03 1.786E-02 4.763E-03 1.003E-01 2.515E-02 3.015E-02 2.185E-02 3.750E+03 1.800E+01 3.760E-02 4.394E-01 -1.284E-06 -1.402E-06 -8.410E-13 -3.840E-12 3.791E-12 4.040E+00 3.269E+00 0.000E+00 3.739E-01 1 0.009 8.768E-03 1.614E-06 -5.641E-09 8.859E-09 9.332E-09 9.010E-09 8.353E-11 1.649E-04 4.612E-04 9.965E-01 9.955E-01 8.254E-01 1.143E+00 2.613E-02 6.829E-04 2.748E-04 3.360E-04 7.200E-05 7.143E-03 1.786E-02 4.732E-03 1.003E-01 2.515E-02 3.018E-02 2.192E-02 3.761E+03 1.900E+01 3.760E-02 4.401E-01 -1.284E-06 -1.402E-06 4.068E-14 -3.792E-12 2.900E-12 4.040E+00 3.278E+00 0.000E+00 3.739E-01 2 0.018 8.776E-03 2.812E-06 -1.534E-08 7.542E-09 9.403E-09 8.264E-09 1.800E-10 1.588E-04 4.594E-04 9.965E-01 9.955E-01 8.222E-01 1.142E+00 2.494E-02 6.220E-04 2.502E-04 3.056E-04 6.626E-05 7.142E-03 1.786E-02 4.701E-03 1.003E-01 2.515E-02 3.020E-02 2.197E-02 3.757E+03 1.800E+01 3.760E-02 4.419E-01 -1.284E-06 -1.402E-06 9.291E-13 -3.769E-12 1.987E-12 4.040E+00 3.284E+00 0.000E+00 3.739E-01 3 0.026 8.785E-03 3.907E-06 -2.931E-08 6.196E-09 1.022E-08 7.656E-09 2.556E-10 1.552E-04 5.368E-04 9.965E-01 9.955E-01 8.194E-01 1.142E+00 2.386E-02 5.694E-04 2.285E-04 2.791E-04 6.177E-05 7.142E-03 1.786E-02 4.670E-03 1.003E-01 2.515E-02 3.021E-02 2.201E-02 3.761E+03 2.000E+01 3.760E-02 4.433E-01 -1.284E-06 -1.402E-06 1.822E-12 -3.774E-12 1.057E-12 4.039E+00 3.288E+00 0.000E+00 3.739E-01 4 0.035 8.794E-03 4.900E-06 -4.669E-08 4.950E-09 1.187E-08 7.173E-09 3.093E-10 1.549E-04 6.741E-04 9.965E-01 9.954E-01 8.171E-01 1.142E+00 2.289E-02 5.239E-04 2.095E-04 2.560E-04 5.835E-05 7.142E-03 1.786E-02 4.638E-03 1.003E-01 2.514E-02 3.021E-02 2.205E-02 3.761E+03 2.000E+01 3.760E-02 4.442E-01 -1.284E-06 -1.402E-06 2.717E-12 -3.810E-12 1.060E-13 4.039E+00 3.291E+00 0.000E+00 3.740E-01 5 0.044 8.804E-03 5.790E-06 -6.665E-08 3.905E-09 1.445E-08 6.813E-09 3.396E-10 1.586E-04 8.180E-04 9.965E-01 9.954E-01 8.153E-01 1.142E+00 2.201E-02 4.844E-04 1.927E-04 2.359E-04 5.583E-05 7.142E-03 1.785E-02 4.606E-03 1.003E-01 2.514E-02 3.020E-02 2.209E-02 3.760E+03 1.900E+01 3.760E-02 4.445E-01 -1.284E-06 -1.403E-06 3.611E-12 -3.878E-12 -8.645E-13 4.039E+00 3.295E+00 0.000E+00 3.740E-01 6 0.053 8.815E-03 6.579E-06 -8.840E-08 3.120E-09 1.808E-08 6.572E-09 3.443E-10 1.666E-04 9.659E-04 9.965E-01 9.954E-01 8.139E-01 1.142E+00 2.121E-02 4.501E-04 1.778E-04 2.182E-04 5.407E-05 7.142E-03 1.785E-02 4.574E-03 1.003E-01 2.514E-02 3.017E-02 2.213E-02 3.759E+03 2.000E+01 3.760E-02 4.444E-01 -1.284E-06 -1.403E-06 4.498E-12 -3.982E-12 -1.847E-12 4.038E+00 3.299E+00 0.000E+00 3.741E-01 7 0.062 8.825E-03 7.267E-06 -1.112E-07 2.610E-09 2.287E-08 6.454E-09 3.205E-10 1.787E-04 1.115E-03 9.965E-01 9.954E-01 8.131E-01 1.142E+00 2.050E-02 4.202E-04 1.645E-04 2.027E-04 5.293E-05 7.142E-03 1.785E-02 4.541E-03 1.003E-01 2.514E-02 3.014E-02 2.216E-02 3.763E+03 2.000E+01 3.760E-02 4.438E-01 -1.285E-06 -1.403E-06 5.375E-12 -4.121E-12 -2.842E-12 4.038E+00 3.299E+00 0.000E+00 3.741E-01 8 0.070 8.834E-03 7.854E-06 -1.343E-07 2.353E-09 2.894E-08 6.462E-09 2.654E-10 1.943E-04 1.263E-03 9.965E-01 9.953E-01 8.127E-01 1.143E+00 1.985E-02 3.940E-04 1.527E-04 1.890E-04 5.230E-05 7.141E-03 1.785E-02 4.507E-03 1.003E-01 2.514E-02 3.010E-02 2.217E-02 3.758E+03 2.000E+01 3.760E-02 4.428E-01 -1.285E-06 -1.403E-06 6.238E-12 -4.289E-12 -3.843E-12 4.038E+00 3.298E+00 0.000E+00 3.742E-01 9 0.079 8.844E-03 8.343E-06 -1.570E-07 2.297E-09 3.638E-08 6.602E-09 1.759E-10 2.128E-04 1.408E-03 9.965E-01 9.953E-01 8.129E-01 1.143E+00 1.927E-02 3.712E-04 1.422E-04 1.769E-04 5.209E-05 7.141E-03 1.785E-02 4.473E-03 1.003E-01 2.513E-02 3.005E-02 2.220E-02 3.758E+03 2.000E+01 3.760E-02 4.427E-01 -1.285E-06 -1.404E-06 7.083E-12 -4.477E-12 -4.846E-12 4.037E+00 3.298E+00 0.000E+00 3.743E-01 Simulation finished after 10 time-steps Writing final snapshot at time t = 8.8005663742706033E-002 Wall clock time [hours] = 1.091E-03 (+/- 2.7778E-13) Wall clock time/timestep/(meshpoint+particle) [microsec] = 6.26 Maximum used memory per cpu [MBytes] = 38.188 Maximum used memory [MBytes] = 303.391 real 4.63 user 15.33 sys 2.03