SVN: ------- v. ( ) $Id$
SVN: start.in,v v. 1.2 (2016/10/14 23:26:55) wlyra
SVN: run.in,v v. 1.2 (2016/10/14 23:26:55) 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: pointmasses.f90,v v. 1.1 (2019/02/02 03:54:41) wlyra
The verbose level is ip= 14 (ldebug= F )
This is a 2-D run
nxgrid, nygrid, nzgrid= 256 768 1
Lx, Ly, Lz= 1.2000000000000002 6.2831853071795862 0.20000000000000001
Vbox= 7.5398223686155044
rsnap: read snapshot var.dat in 7.7303670000000001E-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 -5.5859994389998171 5.0000000000000003E-002 1.5000000000000011E-003 0.40000000000000008
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
diffusion: mesh hyperdiffusion
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 nothing nothing nothing
heat conduction: hyperdiffusivity of ss
viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq
viscous force: nu_shock*(XXXXXXXXXXX)
sh: 1: Syntax error: Bad fd number
pde: ENTER
SVN: ------- v. ( ) $Id$
calc_pencils_hydro: call gij_etc
calc_pencils_energy: max(advec_cs2) = 865.83658027986576
WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh!
duu_dt: SOLVE
Bcs for ux, x: < a>, y: < p>, z: < p>
Bcs for uy, x: < a2>, y: < p>, z: < p>
Bcs for uz, x: < s>, y: < p>, z: < p>
duu_dt: max(advec_uu) = 482.79792988503829
dlnrho_dt: SOLVE
Bcs for lnrho, x: < a2>, y: < p>, z: < p>
dlnrho_dt: diffrho_hyper3_mesh= 20.000000000000000
dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000
dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000
denergy_dt: SOLVE denergy_dt
Bcs for ss, x: < a2>, y: < p>, z: < p>
denergy_dt: lnTT,cs2,cp1= -4.6697087071256620 6.2500000000000021E-003 0.0000000000000000
calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000
calc_heatcond_hyper3: added thdiff
At line 1154 of file pointmasses.f90
Fortran runtime warning: An array temporary was created for argument 'xxq' of procedure 'get_evr'
----it--------t-----------dt-----------rhomax---------rhomin----------urms-----------ux2m-----------uy2m------------xq1------------yq1------------xq2------------yq2-----------vxq1-----------vyq1-----------vyq2---------torqext_2------torqint_2---
0 0.000 7.21558E-04 1.00000E+00 1.00000E+00 9.99250E-01 0.00000E+00 9.98500E-01 1.26000E-05 -3.14159E+00 9.99987E-01 0.00000E+00 0.00000E+00 1.26000E-05 9.99987E-01 8.13152E-20 1.28749E-19
5 0.004 7.21558E-04 1.00002E+00 9.99999E-01 9.99250E-01 5.81655E-13 9.98500E-01 1.26000E-05 -3.13798E+00 9.99987E-01 3.60779E-03 2.73035E-11 1.26000E-05 9.99987E-01 6.10713E-12 6.00590E-12
10 0.007 7.21558E-04 1.00007E+00 9.99997E-01 9.99250E-01 2.32465E-12 9.98500E-01 1.26000E-05 -3.13438E+00 9.99987E-01 7.21558E-03 5.46915E-11 1.26000E-05 9.99987E-01 -1.38165E-11 -1.02103E-11
15 0.011 7.21558E-04 1.00017E+00 9.99993E-01 9.99250E-01 5.22572E-12 9.98500E-01 1.26000E-05 -3.13077E+00 9.99987E-01 1.08234E-02 8.21629E-11 1.26000E-05 9.99987E-01 9.14435E-12 1.73769E-11
20 0.014 7.21557E-04 1.00030E+00 9.99987E-01 9.99250E-01 9.28119E-12 9.98500E-01 1.26000E-05 -3.12716E+00 9.99987E-01 1.44312E-02 1.09717E-10 1.26000E-05 9.99987E-01 -3.30354E-11 -9.68500E-12
25 0.018 7.21557E-04 1.00044E+00 9.99980E-01 9.99250E-01 1.44870E-11 9.98500E-01 1.26000E-05 -3.12355E+00 9.99987E-01 1.80389E-02 1.37352E-10 1.26000E-05 9.99987E-01 -1.32815E-11 2.86737E-11
30 0.022 7.21557E-04 1.00067E+00 9.99972E-01 9.99250E-01 2.08387E-11 9.98500E-01 1.26000E-05 -3.11995E+00 9.99987E-01 2.16467E-02 1.65068E-10 1.26000E-05 9.99987E-01 -6.56673E-11 9.75895E-12
35 0.025 7.21557E-04 1.00084E+00 9.99962E-01 9.99250E-01 2.83314E-11 9.98500E-01 1.26000E-05 -3.11634E+00 9.99987E-01 2.52545E-02 1.92864E-10 1.26000E-05 9.99987E-01 -7.24493E-11 4.49498E-11
40 0.029 7.21556E-04 1.00117E+00 9.99952E-01 9.99250E-01 3.69599E-11 9.98500E-01 1.26000E-05 -3.11273E+00 9.99987E-01 2.88623E-02 2.20739E-10 1.26000E-05 9.99987E-01 -1.23903E-10 5.22231E-11
45 0.032 7.21556E-04 1.00143E+00 9.99940E-01 9.99250E-01 4.67186E-11 9.98500E-01 1.26000E-05 -3.10912E+00 9.99987E-01 3.24701E-02 2.48692E-10 1.26000E-05 9.99987E-01 -1.74388E-10 7.59112E-11
Simulation finished after 46 time-steps
Writing final snapshot at time t = 3.3191633335019489E-002
Wall clock time [hours] = 2.037E-03 (+/- 2.7778E-13)
Wall clock time/timestep/meshpoint [microsec] = 0.8109456
Wall clock time/timestep/local meshpoint [microsec] = 3.243782
Rhs wall clock time/timestep/local meshpoint [microsec] = 2.428758
Maximum used memory per cpu [MBytes] = 74.055
Maximum used memory [MBytes] = 295.547
real 9.00
user 29.47
sys 0.48