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 0-D run
 nxgrid, nygrid, nzgrid=           1           1           1
 Lx, Ly, Lz=   6.2831853071795862        6.2831853071795862        6.2831853071795862     
       Vbox=   1.0000000000000000     
 rsnap: read snapshot var.dat in    0.0000000000000000       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=   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_density: no hydro, turned off continuity equation
 initialize_density: 0-D run, turned off continuity equation
 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 cs2
 initialize_particles: Global pressure gradient with beta_dPdr_dust= -0.10000000000000001     
 initialize_particles: dust-to-gas ratio eps_dtog=   1.0000000000000000     
 initialize_particles: mass per constituent particle mpmat=   4.1887902047863905     
 initialize_particles: mass per superparticle mp_swarm =   1.0000000000000000E-002
 initialize_particles: number density per superparticle np_swarm=   0.0000000000000000     
 initialize_particles: mass density per superparticle rhop_swarm=   1.0000000000000000E-002
 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=   4.1887902047863905     
 pde: ENTER
SVN: -------            v.         (                   ) $Id$
 duu_dt: diagnostics ...
 dlnrho_dt: SOLVE
Bcs for      lnrho,  x: <       p>, y: <       p>,  z: <       p>
 dlnrho_dt: max(diffus_diffrho ) =   0.0000000000000000     
 dlnrho_dt: max(diffus_diffrho3) =   0.0000000000000000     
 SOLVE dlncc_dt
Bcs for         cc,  x: <       p>, y: <       p>,  z: <       p>
 dvvp_dt_pencil: calculate dvvp_dt
 dvvp_dt_pencil: ldraglaw_purestokes= F
 dap_dt_pencil: Calculate dap/dt

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
#0  0x765eb423e59 in ???
#1  0x765eb422e75 in ???
#2  0x765eb04532f in ???
	at ./signal/../sysdeps/unix/sysv/linux/x86_64/libc_sigaction.c:0
#3  0x60415d73652d in dap_dt_sweepup_pencil
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/particles_radius.f90:758
#4  0x60415d736e0b in __particles_radius_MOD_dap_dt_pencil
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/particles_radius.f90:653
#5  0x60415d70eeba in __particles_main_MOD_particles_pde_pencil
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/particles_main.f90:873
#6  0x60415d86dc3e in calc_all_rhs
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/equ.f90:1399
#7  0x60415d86f5fa in rhs_cpu
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/equ.f90:1541
#8  0x60415d874633 in __equ_MOD_pde
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/equ.f90:332
#9  0x60415d868298 in advance_substeps
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/timestep.f90:136
#10  0x60415d86987f in __timestep_MOD_time_step
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/timestep.f90:273
#11  0x60415d886eac in __run_module_MOD_timeloop
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/run.f90:544
#12  0x60415d889857 in run_start
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/run.f90:1406
#13  0x60415d88a61c in run
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/run.f90:1520
#14  0x60415d88ac34 in main
	at /home/pencil/defective/samples/0d-tests/coagulation-fragmentation/src/run.f90:1518
real 0.40
user 0.19
sys 0.03