SVN: ------- v. ( ) $Id$
SVN: ------- v. ( ) no cvsid is given i
SVN: ------- v. ( ) no cvsid is given i
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: ------- v. ( ) $Id$
SVN: ------- v. ( ) $Id$
The verbose level is ip= 14 (ldebug= F )
This is a 2-D run
nxgrid, nygrid, nzgrid= 128 256 1
Lx, Ly, Lz= 0.20000000000000001 0.40000000000000002 2.0000000000000000
Vbox= 8.0000000000000016E-002
rsnap: read snapshot var.dat in 1.0000000474974513E-003 seconds
setup_slices: slice_position = p
setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4
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= 4.8108883747190188E-009
units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 6.3969296552161463 20.000000000000000 240.00000000000006 0.40000000000000008
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 rho and ss
viscous force: nu*(del2u+graddivu/3+2S.glnrho)
By using fbc[x,y,z] I found the flow direction to be in the 2.0000000000000000 direction.
initialize_particles: dust-to-gas ratio eps_dtog= 1.0000000000000000
initialize_particles: mass per constituent particle mpmat= 2.8730912014629860E-008
initialize_particles: mass per superparticle mp_swarm = 4.0000000000000010E-003
initialize_particles: number density per superparticle np_swarm= 0.0000000000000000
initialize_particles: mass density per superparticle rhop_swarm= 1638.4000000000001
initialize particles: set minimum collisional time-scale equal to two times the Courant time-step.
initialize_particles_map: selected cloud-in-cell for particle-mesh method.
initialize_particles_radius: mass per dust grain mpmat= 2.8730912014629860E-008
sh: 1: Syntax error: Bad fd number
pencil_consistency_check: checking pencil case
At line 2313 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'xxp' of procedure 'interpolate_linear_range'
At line 2313 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'inear' of procedure 'interpolate_linear_range'
At line 2313 of file particles_map.f90
Fortran runtime warning: An array temporary was created for argument 'gp' of procedure 'interpolate_linear_range'
dxxp_dt: Calculate dxxp_dt
dxxp_dt: Particles boundary condition bcpx=p
dxxp_dt: Particles boundary condition bcpy=rmv
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
particle_gravity: Calculating gravity
particle_gravity: No gravity in x-direction.
particle_gravity: No gravity in z-direction.
particle_gravity: No radial gravity
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
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$
calc_pencils_energy: max(advec_cs2) = 326402499.99999988
duu_dt: SOLVE
Bcs for ux, x: < p>, y: < e2>, z: < p>
Bcs for uy, x: < p>, y: < e2>, z: < p>
Bcs for uz, x: < p>, y: < e2>, z: < p>
duu_dt: max(advec_uu) = 3439.4903954510964
dlnrho_dt: SOLVE
Bcs for lnrho, x: < p>, y: < e2>, z: < p>
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
dvvp_dt_pencil: calculate dvvp_dt
dvvp_dt_pencil: ldraglaw_purestokes= F
dvvp_dt: Add drag force; tausp= 0.0000000000000000
inlet_profile: uniform
inlet_YY_profile: uniform,
At line 5531 of file particles_dust.f90
Fortran runtime warning: An array temporary was created for argument 'part_pos' of procedure 'in_solid_cell'
---it----------t---------------dt---------------urms-------------umax-------------rhom----------------ozm------------------oz2m----------------vpxmax--------------vpymax--------------vpzmax-------------c_dragx-------------c_dragy-------
0 0.000000 1.5959793897E-05 5.0724107723 9.3062488478 1.0000000000 -8.08093282E-03 1.0372069552E+05 1.0792421268E-01 5.1005706624E+00 1.7864723862E-01 -5.4527607994E-03 8.6859838834E-01
10 0.000157 1.5360722548E-05 5.0675209861 9.7737254937 0.9999738337 -1.82013455E-03 8.3351549507E+04 1.0799875229E-01 5.0999457161E+00 1.7751175660E-01 -6.9433866025E-02 2.4747227511E+00
20 0.000311 1.5561133351E-05 5.0646483705 9.8435208743 0.9999275608 1.70355839E-03 7.3233643547E+04 1.0692140498E-01 5.0989685502E+00 1.7541776670E-01 5.2280257525E-02 2.1502691308E+00
30 0.000468 1.5904334171E-05 5.0625654390 9.1630038818 0.9998866582 1.61394303E-03 6.7214216205E+04 1.0509602520E-01 5.0979298633E+00 1.7308877100E-01 3.8787003496E-02 2.1542941656E+00
40 0.000628 1.6003888115E-05 5.0609051508 8.9766850719 0.9998551785 -9.21934338E-04 6.3304604192E+04 1.0378993004E-01 5.0970533100E+00 1.7057608641E-01 3.7526160098E-02 2.3671957780E+00
50 0.000789 1.6196264694E-05 5.0594402009 8.8721945663 0.9998373262 -3.36716931E-03 6.0665829991E+04 1.0221959987E-01 5.0954227043E+00 1.6793690431E-01 2.5748899706E-02 2.5425171830E+00
60 0.000951 1.6179860353E-05 5.0580828520 8.7110633022 0.9998254465 -3.89412670E-03 5.8875864702E+04 1.0089326862E-01 5.0942858101E+00 1.6509047001E-01 -2.8353266960E-02 2.5154539235E+00
70 0.001113 1.6311146372E-05 5.0568363983 8.7485931535 0.9998175744 -2.16006826E-03 5.7500653348E+04 9.9720051036E-02 5.0929658459E+00 1.6220677971E-01 -5.8804014716E-03 2.5181525895E+00
80 0.001277 1.6376894860E-05 5.0556141018 8.7595984582 0.9998059871 -2.36452614E-05 5.6381377836E+04 9.8291032861E-02 5.0924226952E+00 1.5932095321E-01 2.1652223508E-02 2.5204865642E+00
90 0.001441 1.6372577963E-05 5.0545104459 8.7471531174 0.9997925026 5.47962900E-04 5.5491352458E+04 9.6657371243E-02 5.0910655090E+00 1.5641668455E-01 1.5340784611E-02 2.4177485102E+00
Simulation finished after 91 time-steps
Writing final snapshot at time t = 1.4569666143881704E-003
Wall clock time [hours] = 1.161E-03 (+/- 2.7778E-07)
Wall clock time/timestep/(meshpoint+particle) [microsec] = 1.40
Maximum used memory per cpu [MBytes] = 57.750
Maximum used memory [MBytes] = 57.750
real 4.47
user 4.45
sys 0.01