MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.2.0 phantom_version not found; no git information to display Using NATIVE inverse sqrt ERROR! compile settings: -DDUST currently not compatible with magnetic fields (-DMHD) Running on Intel(R) Xeon(R) CPU E5-2680 v3 @ 2.50GHz 24 cpus @ 2.50 Ghz, cache size 16384 KB Running in openMP on 16 threads --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1000000) 30.518 MB real(8) xyzh_soa(1000000,4) 30.518 MB real(8) vxyzu(4,1000000) 30.518 MB real(4) alphaind(2,1000000) 7.629 MB real(4) divcurlv(4,1000000) 15.259 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,1000000) 30.518 MB real(8) Bxyz(3,1000000) 22.888 MB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1000000) 167.847 MB real(8) dustevol(11,1000000) 83.923 MB real(8) ddustevol(11,1000000) 83.923 MB real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,11,1000000) 251.770 MB real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1000000) 30.518 MB real(8) dBevol(4,1000000) 30.518 MB real(4) divBsumm(1000000) 3.815 MB real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(11,1000000) 83.923 MB real(8) Bpred(4,1000000) 30.518 MB real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1000000) 976.562 KB integer(1) iphase_soa(1000000) 976.562 KB real(4) gradh(1,1000000) 3.815 MB real(8) tstop(22,1000000) 167.847 MB integer(4) ll(1000000) 3.815 MB integer(4) ibelong(1000000) 3.815 MB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,1000001) 7.629 MB integer(4) inodeparts(1000000) 3.815 MB real(8) xyzh_swap(1000000,4) 30.518 MB integer(4) inodeparts_swap(1000000) 3.815 MB integer(1) iphase_swap(1000000) 976.562 KB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) cellatid(1000001) 3.815 MB integer(4) ifirstincell(1000001) 3.815 MB kdnode nodeglobal(1000001) 49.591 MB kdnode node(1000001) 49.591 MB integer(4) nodemap(1000001) 3.815 MB --------------------------------------------------------- Total memory allocated to arrays: 1.370 GB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) --> TESTING M_4 cubic kernel checking wab0 = wab(0).................OK [max err = 0.000E+00, tol = 2.225-308] checking gradh0 = -3.*wab(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking dphidh0 = dphidh(0)...........OK [max err = 0.000E+00, tol = 2.225-308] checking radkern2 = radkern*radkern....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == wkern...........OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel == grkern..........OK [max err = 0.000E+00, tol = 2.225-308] checking w gradient equal to gradw.....OK [max err = 1.499E-07, tol = 2.000E-05] checking potential gradient = force....OK [max err = 6.722E-08, tol = 2.000E-05] checking get_kernel_grav1 == wkern.....OK [max err = 0.000E+00, tol = 2.225-308] checking get_kernel_grav1 == grkern....OK [max err = 0.000E+00, tol = 2.225-308] checking dphi/dh = phi - q*dphi/dq.....OK [max err = 8.105E-15, tol = 2.000E-07] <-- KERNEL TEST COMPLETE --> TESTING LINKLIST / NEIGHBOUR FINDING ----------------------------------------------------------------------------------------- random seed = -43587 32768 particles set in uniform random distribution ----------------------------------------------------------------------------------------- Test 1 of 2: building linked list... checking ncells > 0....................OK [got T should be T] checking neighbour number.............. checking nneigh (cached)...............OK [checked 32768 values] checking nneigh (no cache).............OK [checked 32768 of 32768 values] max nneigh = 8945 min nneigh = 1 mean = 2206.93 Test 2 of 2: building linked list... checking ncells > 0....................OK [got T should be T] Test 3: building linked list... particles in a line in x direction checking ncells > 0....................OK [got T should be T] particles in a line in y direction checking ncells > 0....................OK [got T should be T] particles in a line in z direction checking ncells > 0....................OK [got T should be T] <-- LINKLIST TEST COMPLETE --> TESTING DERIVS MODULE ---------------- particles set on 100 x 100 x 100 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 1.000E-02 dy: 1.000E-02 dz: 1.000E-02 ----------------------------------------------------------------------------------------- thread 0 npart = 1000000 hfact = 1.2 --> testing Hydro derivatives trial neigh mean : 367.52, real neigh mean = 57.00 ratio try/act= 6.45 trial neigh max : 548, max real neigh = 57 n neighbour calls : 155520, mean per part = 1.156 n density calcs : 2000000, mean per part = 2.000 dtmin = 1.048319032629973E-003 dtmax = 1.955273019938114E-002 dtmax/dtmin = 18.6515074045035 dtcour/dtf = 0.164449106870774 (wall: total = 19.62s force = 6.88s dens = 10.00s link = 2.75s) (cpu : total = 305.11s force = 111.18s dens = 158.52s link = 35.41s) (cpu/wall: total = 15.55 force = 16.17 dens = 15.85 link = 12.87) completed in 19.62 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.962E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradh.........................OK [max err = 3.167E-06, tol = 1.000E-05] checking force(x)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(y)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking force(z)......................OK [max err = 5.476E-04, tol = 1.000E-03] checking du/dt.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking n density calcs...............OK [is 2000000 should be 2000000] checking total nneigh..................OK [is 57000000 should be 57000000] --> testing artificial viscosity terms (individual alpha) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 5.121481238535053E-004 dtmax = 5.121481238535085E-004 dtmax/dtmin = 0.999999999999876 dtcour/dtf = 1.681859564884428E-002 (wall: total = 14.75s force = 6.88s dens = 5.12s link = 2.75s) (cpu : total = 227.35s force = 110.15s dens = 81.84s link = 35.36s) (cpu/wall: total = 15.41 force = 16.02 dens = 15.97 link = 12.86) completed in 14.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.670E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking art. visc force(x)............OK [max err = 2.705E-02, tol = 3.200E-02] checking art. visc force(y)............OK [max err = 1.147E-02, tol = 2.400E-02] checking art. visc force(z)............OK [max err = 9.451E-03, tol = 2.400E-02] --> testing ddivv/dt in Cullen & Dehnen switch trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 completed in 7.88 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.730E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking alphaloc......................OK [max err = 3.421E-04, tol = 3.500E-04] --> testing physical viscosity terms (two first derivatives) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 0.300000000000000 dtmax = 0.300000000000000 dtmax/dtmin = 1.00000000000000 dtcour/dtf = 447.734858356076 WARNING! force: viscosity constraining Courant timestep (wall: total = 15.50s force = 7.50s dens = 5.12s link = 2.88s) (cpu : total = 239.17s force = 121.76s dens = 81.97s link = 35.45s) (cpu/wall: total = 15.43 force = 16.23 dens = 15.99 link = 12.33) completed in 15.50 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlv(x)......................OK [max err = 1.478E-03, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.621E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvxdx.........................OK [max err = 8.214E-04, tol = 1.700E-03] checking dvxdy.........................OK [max err = 1.426E-16, tol = 2.500E-15] checking dvxdz.........................OK [max err = 1.426E-16, tol = 2.500E-15] checking dvydx.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvydy.........................OK [max err = 2.342E-16, tol = 2.500E-15] checking dvydz.........................OK [max err = 8.214E-04, tol = 1.000E-03] checking dvzdx.........................OK [max err = 4.243E-17, tol = 2.500E-15] checking dvzdy.........................OK [max err = 6.568E-04, tol = 1.500E-03] checking dvzdz.........................OK [max err = 4.328E-17, tol = 2.500E-15] checking viscous force(x)..............OK [max err = 3.949E-02, tol = 4.000E-02] checking viscous force(y)..............OK [max err = 2.994E-02, tol = 3.000E-02] checking viscous force(z)..............OK [max err = 2.730E-02, tol = 3.100E-02] checking mean nneigh...................OK [max err = 0.000E+00, tol = 2.225-308] checking max nneigh....................OK [got 57 should be 57] checking n density calcs...............OK [is 1000000 should be 1000000] checking total nneigh..................OK [is 57000000 should be 57000000] checking \sum v.dv/dt + du/dt = 0......OK [max err = 4.500E-07, tol = 1.520E-06] checking du/dt >= 0 for all particles..OK [got T should be T] --> testing dust evolution terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 1.048319035706207E-003 dtmax = 1.955273025675754E-002 dtmax/dtmin = 18.6515074045035 dtcour/dtf = 0.157834702373746 (wall: total = 21.00s force = 13.12s dens = 5.12s link = 2.75s) (cpu : total = 326.93s force = 209.95s dens = 81.51s link = 35.47s) (cpu/wall: total = 15.57 force = 16.00 dens = 15.90 link = 12.90) completed in 21.00 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divv..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking deps/dt.......................OK [max err = 1.826E-05, tol = 4.000E-05] checking du/dt.........................OK [max err = 4.482E-04, tol = 5.000E-04] checking deltavx.......................OK [max err = 2.282E-05, tol = 2.300E-05] checking energy conservation (dE=0)....OK [max err = 8.231E-16, tol = 6.500E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] checking dust mass conservation........OK [max err = 1.394E-17, tol = 1.000E-15] --> testing MHD derivatives (using B/rho directly) trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 4.072471057823637E-004 dtmax = 3.627823185039153E-003 dtmax/dtmin = 8.90816198207017 dtcour/dtf = 0.269845464080934 (wall: total = 19.75s force = 11.88s dens = 5.00s link = 2.88s) (cpu : total = 305.70s force = 188.74s dens = 81.45s link = 35.51s) (cpu/wall: total = 15.48 force = 15.89 dens = 16.29 link = 12.35) completed in 19.75 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB (symm)...................OK [max err = 8.214E-04, tol = 2.000E-03] checking dBx/dt........................OK [max err = 2.952E-04, tol = 2.000E-03] checking dBy/dt........................OK [max err = 1.161E-03, tol = 2.000E-03] checking dBz/dt........................OK [max err = 2.482E-04, tol = 2.000E-02] checking mhd force(x)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(y)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking mhd force(z)..................OK [max err = 2.500E-02, tol = 2.500E-02] checking div B (diff)..................OK [max err = 8.214E-04, tol = 1.000E-03] checking curlB(x)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(y)......................OK [max err = 8.219E-04, tol = 1.000E-03] checking curlB(z)......................OK [max err = 8.219E-04, tol = 1.000E-03] --> testing artificial resistivity terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 4.534818241332849E-004 dtmax = 3.167200767364767E-003 dtmax/dtmin = 6.98418458869349 dtcour/dtf = 0.278098673087194 (wall: total = 19.50s force = 11.62s dens = 5.12s link = 2.75s) (cpu : total = 302.16s force = 185.14s dens = 81.54s link = 35.48s) (cpu/wall: total = 15.50 force = 15.93 dens = 15.91 link = 12.90) completed in 19.50 s checking dBx/dt (resist)...............OK [max err = 0.000E+00, tol = 3.700E-02] checking dBy/dt (resist)...............OK [max err = 0.000E+00, tol = 3.400E-02] checking dBz/dt (resist)...............OK [max err = 0.000E+00, tol = 2.200E-01] checking \sum du/dt + B.dB/dt = 0......OK [max err = 0.000E+00, tol = 2.700E-03] --> testing div B cleaning terms trial neigh mean : 363.88, real neigh mean = 57.00 ratio try/act= 6.38 trial neigh max : 548, max real neigh = 57 n neighbour calls : 0, mean per part = 1.000 n density calcs : 1000000, mean per part = 1.000 dtmin = 1.949447112251919E-004 dtmax = 3.812184925514154E-004 dtmax/dtmin = 1.95552108162051 dtcour/dtf = 0.119550250157675 (wall: total = 19.50s force = 11.62s dens = 5.00s link = 2.88s) (cpu : total = 302.46s force = 185.51s dens = 81.49s link = 35.45s) (cpu/wall: total = 15.51 force = 15.96 dens = 16.30 link = 12.33) completed in 19.50 s checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking divB..........................OK [max err = 8.214E-04, tol = 1.000E-03] checking gradpsi_x.....................OK [max err = 4.261E-04, tol = 8.500E-04] checking gradpsi_y.....................OK [max err = 8.500E-04, tol = 8.500E-04] checking gradpsi_z.....................OK [max err = 3.729E-04, tol = 2.000E-03] --> testing Hydro derivs in setup with density contrast ---------------- particles set on 500 x 500 x 500 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-03 dy: 2.000E-03 dz: 2.000E-03 ----------------------------------------------------------------------------------------- ---------------- particles set on 500 x 500 x 500 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-03 dy: 2.000E-03 dz: 2.000E-03 ----------------------------------------------------------------------------------------- ---------------- particles set on 50 x 50 x 50 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-02 dy: 2.000E-02 dz: 2.000E-02 ----------------------------------------------------------------------------------------- thread 0 npart = 648432 in blob = 523984 to test = 73824 WARNING! maketree: maxlevel > max_indexed: will run faster if recompiled with NCELLSMAX=2*maxp, trial neigh mean : 430.91, real neigh mean = 57.47 ratio try/act= 7.50 trial neigh max : 88189, max real neigh = 988 n neighbour calls : 124015, mean per part = 1.191 n density calcs : 1452028, mean per part = 2.239 dtmin = 2.881763406465614E-004 dtmax = 2.154077551929841E-002 dtmax/dtmin = 74.7485913346109 dtcour/dtf = 1.32743250797427 (wall: total = 19.88s force = 8.62s dens = 9.12s link = 2.12s) (cpu : total = 310.66s force = 137.89s dens = 146.00s link = 26.77s) (cpu/wall: total = 15.63 force = 15.99 dens = 16.00 link = 12.60) completed in 19.88 s checking h (density)...................OK [max err = 3.510E-04, tol = 4.000E-04] checking divv..........................OK [max err = 3.289E-05, tol = 1.000E-03] checking curlv(x)......................OK [max err = 4.862E-05, tol = 1.500E-03] checking curlv(y)......................OK [max err = 1.126E-16, tol = 1.000E-03] checking curlv(z)......................OK [max err = 3.289E-05, tol = 1.000E-03] checking force(x)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking force(y)......................OK [max err = 6.909E-06, tol = 1.000E-03] checking force(z)......................OK [max err = 2.193E-05, tol = 1.000E-03] checking du/dt.........................OK [max err = 3.290E-05, tol = 1.000E-03] checking mean nneigh...................OK [max err = 0.000E+00, tol = 1.000E-17] checking max nneigh....................OK [got 988 should be 988] checking total nneigh..................OK [is 37263216 should be 37263216] <-- DERIVS TEST COMPLETE --> TESTING EQUATION OF STATE MODULE --> testing equation of state initialisation checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] Reading output_DE_z0.00x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.00x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.02x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.00.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.20.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.40.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.60.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading output_DE_z0.04x0.80.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading opacs.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) Reading /opacs.bindata in /home/dprice/phantom-nightly/phantom/data/eos/mesa/ (from PHANTOM_DIR setting) checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] Reading helm_data.tab in /home/dprice/phantom-nightly/phantom/data/eos/helmholtz/ (from PHANTOM_DIR setting) checking eos initialisation............OK [got 0 should be 0] checking eos initialisation............OK [got 0 should be 0] --> testing barotropic equation of state Barotropic eq of state: cs_ld = 3.162E-01 code units = 3.643E+04 cm/s Barotropic eq of state: cs = 3.162E-01 code units = 3.643E+04 cm/s Barotropic eq of state: rhocritT == rhoT = 2.514E-04 code units = 5.000E-19 g/cm^3 Barotropic eq of state: rhocrit0 == rho0 = 5.027E-04 code units = 1.000E-18 g/cm^3 Barotropic eq of state: rhocrit1 == rho1 = 5.027E+00 code units = 1.000E-14 g/cm^3 Barotropic eq of state: rhocrit2 == rho2 = 5.027E+04 code units = 1.000E-10 g/cm^3 Barotropic eq of state: rhocrit3 == rho3 = 5.027E+11 code units = 1.000E-03 g/cm^3 Barotropic eq of state: Barotropic eq of state: P = cs_ld*rho for rho/(g/cm^3) < rhoT Barotropic eq of state: P = 10**(log10(cs_bg**2) + M*(log10(rhoT)-log10(rho))) for rhoT <= rho/(g/cm^3) < rho0 Barotropic eq of state: P = cs*rho for rho0 <= rho/(g/cm^3) < rho1 Barotropic eq of state: P = cs*rho1*(rho /rho1)^ 1.400 for rho1 <= rho/(g/cm^3) < rho2 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^ 1.400*(rho /rho2)^ 1.100 for rho2 <= rho/(g/cm^3) < rho3 Barotropic eq of state: P = cs*rho1*(rho2/rho1)^ 1.400*(rho3/rho2)^ 1.100*(rho /rho3)^ 1.667 for rho3 <= rho/(g/cm^3) checking equation of state is continuous......... checking p/rho is continuous...........OK [checked 4999 of 5000 values] <-- EQUATION OF STATE TEST COMPLETE --> TESTING COOLING MODULE --> testing find_in_table routine checking table(i) < val < table(i+1)... <-- COOLING TEST COMPLETE --> TESTING KDTREE <-- KDTREE TEST COMPLETE --> SKIPPING SELF-GRAVITY TESTS (need -DGRAVITY) --> TESTING DUST MODULE --> testing drag initialisation checking drag initialisation...........OK [got 0 should be 0] checking drag initialisation...........OK [got 0 should be 0] checking deltav=0 gives Epstein drag...OK [got 1 should be 1] --> testing Epstein/Stokes drag transition checking ts is continuous into Stokes regime: rho=1.00E-19...........OK [max err = 3.172E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-18...........OK [max err = 3.172E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-17...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-16...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-15...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-14...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-13...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-12...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-11...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-10...........OK [max err = 6.244E-02, tol = 6.300E-02] checking ts is continuous into Stokes regime: rho=1.00E-09...........OK [max err = 6.244E-02, tol = 6.300E-02] checking Epstein drag formula matches non-linear solution..OK [max err = 3.865E-03, tol = 3.900E-03] --> testing DUST DRAG trial neigh mean : 5438.64, real neigh mean = 50.29 ratio try/act= 108.15 trial neigh max : 12675, max real neigh = 123 n neighbour calls : 46204, mean per part = 1.258 n density calcs : 785537, mean per part = 4.387 dtmin = 3.590489471945241E-004 dtmax = 3.297731851752349E-003 dtmax/dtmin = 9.18463033388343 dtcour/dtf = 63904.8670213282 (wall: total = 36.00s force = 14.25s dens = 21.25s link = 0.50s) (cpu : total = 575.11s force = 227.90s dens = 340.06s link = 7.15s) (cpu/wall: total = 15.98 force = 15.99 dens = 16.00 link = 14.30) checking acceleration from drag conserves momentum(x)......OK [max err = 1.974E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(y)......OK [max err = 1.814E-10, tol = 7.000E-07] checking acceleration from drag conserves momentum(z)......OK [max err = 3.948E-10, tol = 7.000E-07] checking acceleration from drag conserves energy...........OK [max err = 1.006E-07, tol = 7.000E-07] --> testing DUSTYBOX checking dust velocities match exact solution....OK [max err = 9.090E-05, tol = 1.000E-04] checking dust accel matches exact solution.......OK [max err = 2.848E-03, tol = 3.000E-03] checking gas velocities match exact solution.....OK [max err = 9.090E-05, tol = 1.000E-04] checking gas accel matches exact solution........OK [max err = 2.848E-03, tol = 3.000E-03] checking kinetic energy decay matches exact......OK [max err = 1.302E-05, tol = 1.000E-04] --> testing DUSTYDIFFUSE checking sum(epsilon_k) = epsilon......OK [max err = 0.000E+00, tol = 1.000E-14] checking dust diffusion matches exact solution...OK [max err = 2.430E-03, tol = 2.600E-03] <-- DUST TEST COMPLETE --> SKIPPING DUSTGROWTH TEST (REQUIRES -DDUST -DDUSTGROWTH) --> SKIPPING NON-IDEAL MHD TEST (REQUIRES -DNONIDEALMHD) --> TESTING READ/WRITE from dump file --------> TIME = 20.00 : full dump written to file test.dump <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000048F0F2 dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000096B05B testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000049284F dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000096B05B testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000493AC9 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000096B05B testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000494D47 dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000096B05B testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000497828 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 000000000096B05B testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000049C933 dump_utils_mp_wri 1681 utils_dumpfiles.f90 phantom 0000000000724391 readwrite_dumps_m 539 readwrite_dumps.F90 phantom 000000000096B05B testrwdump_mp_tes 158 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown >>> reading setup from file: test.dump on unit 9 FT:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 00:34:53.1 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000476D07 dump_utils_mp_rea 1225 utils_dumpfiles.f90 phantom 000000000072D9A3 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000047C50B dump_utils_mp_rea 1246 utils_dumpfiles.f90 phantom 000000000072D9A3 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000047E276 dump_utils_mp_rea 1253 utils_dumpfiles.f90 phantom 000000000072D9A3 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000047FFE5 dump_utils_mp_rea 1260 utils_dumpfiles.f90 phantom 000000000072D9A3 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004845FF dump_utils_mp_rea 1281 utils_dumpfiles.f90 phantom 000000000072D9A3 readwrite_dumps_m 782 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.30000000000000 External field found, Bext = 124.000000000000 125.000000000000 126.000000000000 ERROR reading values to verify conservation laws. Resetting initial values. ERROR reading grain size/density from file header time = 20.0000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(2,1001) 7.820 KB real(4) divcurlv(4,1001) 15.641 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,1001) 31.281 KB real(8) Bxyz(3,1001) 23.461 KB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1001) 172.047 KB real(8) dustevol(11,1001) 86.023 KB real(8) ddustevol(11,1001) 86.023 KB real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,11,1001) 258.070 KB real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1001) 31.281 KB real(8) dBevol(4,1001) 31.281 KB real(4) divBsumm(1001) 3.910 KB real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(11,1001) 86.023 KB real(8) Bpred(4,1001) 31.281 KB real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1001) 1001.000 B integer(1) iphase_soa(1001) 1001.000 B real(4) gradh(1,1001) 3.910 KB real(8) tstop(22,1001) 172.047 KB integer(4) ll(1001) 3.910 KB integer(4) ibelong(1001) 3.910 KB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,2003) 15.648 KB integer(4) inodeparts(1001) 3.910 KB real(8) xyzh_swap(1001,4) 31.281 KB integer(4) inodeparts_swap(1001) 3.910 KB integer(1) iphase_swap(1001) 1001.000 B integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) cellatid(2003) 7.824 KB integer(4) ifirstincell(2003) 7.824 KB kdnode nodeglobal(2003) 101.715 KB kdnode node(2003) 101.715 KB integer(4) nodemap(2003) 7.824 KB --------------------------------------------------------- Total memory allocated to arrays: 1.532 MB --------------------------------------------------------- reading particles 1: 1001, from block 1 lims= 1- 1001 forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000735146 readwrite_dumps_m 1261 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A25B0 dump_utils_mp_rea 1884 utils_dumpfiles.f90 phantom 000000000073617C readwrite_dumps_m 1288 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #2 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004A635E dump_utils_mp_rea 1963 utils_dumpfiles.f90 phantom 0000000000736708 readwrite_dumps_m 1305 readwrite_dumps.F90 phantom 000000000072EFE3 readwrite_dumps_m 898 readwrite_dumps.F90 phantom 000000000096CE22 testrwdump_mp_tes 205 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown got 11 sink properties from 10 sink particles ----------------------------------------------- ID| Mass | Racc | Macc | hsoft | ----------------------------------------------- 1| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 2| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 3| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 4| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 5| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 6| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 7| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 8| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 9| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 10| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | ----------------------------------------------- <<< finished reading (MHD) file --> checking header checking read without errors...........OK [got 0 should be 0] checking time..........................OK [max err = 0.000E+00, tol = 2.225-308] checking hfact.........................OK [max err = 0.000E+00, tol = 2.225-308] checking gamma.........................OK [max err = 0.000E+00, tol = 2.225-308] checking polyk.........................OK [max err = 0.000E+00, tol = 2.225-308] checking m(gas)........................OK [max err = 0.000E+00, tol = 2.225-308] checking m(dust).......................OK [max err = 0.000E+00, tol = 2.225-308] checking npartoftype(1)................OK [got 991 should be 991] checking npartoftype(2)................OK [got 10 should be 10] checking umass.........................OK [max err = 0.000E+00, tol = 2.225-308] checking udist.........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bexty.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextz.........................OK [max err = 0.000E+00, tol = 2.225-308] --> checking particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking vx............................OK [max err = 0.000E+00, tol = 2.225-308] checking vy............................OK [max err = 0.000E+00, tol = 2.225-308] checking vz............................OK [max err = 0.000E+00, tol = 2.225-308] checking u.............................OK [max err = 0.000E+00, tol = 2.225-308] checking alpha.........................OK [max err = 0.000E+00, tol = 1.175E-38] checking Bx............................OK [max err = 0.000E+00, tol = 1.175E-38] checking By............................OK [max err = 0.000E+00, tol = 1.175E-38] checking Bz............................OK [max err = 0.000E+00, tol = 1.175E-38] checking psi...........................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking dustfrac......................OK [max err = 0.000E+00, tol = 1.175E-38] checking particle type 1...............OK [checked 991 values] checking particle type 2...............OK [checked 10 values] --> checking sink particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking m.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking hsoft.........................OK [max err = 0.000E+00, tol = 2.225-308] checking maccreted.....................OK [max err = 0.000E+00, tol = 2.225-308] checking spinx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spiny.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spinz.........................OK [max err = 0.000E+00, tol = 2.225-308] checking tlast.........................OK [max err = 0.000E+00, tol = 2.225-308] checking vx............................OK [max err = 0.000E+00, tol = 2.225-308] checking vy............................OK [max err = 0.000E+00, tol = 2.225-308] checking vz............................OK [max err = 0.000E+00, tol = 2.225-308] --> checking read of single array from file x forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004ABF68 dump_utils_mp_rea 2103 utils_dumpfiles.f90 phantom 0000000000971F70 testrwdump_mp_tes 309 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ->y z vx vy vz u dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac dustfrac tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop tstop checking error flag....................OK [got 0 should be 0] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] --------> TIME = 20.00 : small dump written to file test.dump <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000048F0F2 dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 000000000072A94A readwrite_dumps_m 631 readwrite_dumps.F90 phantom 000000000096AFFA testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000049284F dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 000000000072A94A readwrite_dumps_m 631 readwrite_dumps.F90 phantom 000000000096AFFA testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000493AC9 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 000000000072A94A readwrite_dumps_m 631 readwrite_dumps.F90 phantom 000000000096AFFA testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000494D47 dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 000000000072A94A readwrite_dumps_m 631 readwrite_dumps.F90 phantom 000000000096AFFA testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000497828 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 000000000072A94A readwrite_dumps_m 631 readwrite_dumps.F90 phantom 000000000096AFFA testrwdump_mp_tes 156 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown --> checking read_dump >>> reading setup from file: test.dump on unit 9 ERROR! read_dump: file is not a Phantom full dump checking read_dump returns is_small_dump error code........OK [got 1978 should be 1978] --> checking read_smalldump >>> reading small dump file: test.dump on unit 9 ST:Phantom:1.2.0: (mhd+clean+dust) : 21/02/2019 00:34:53.5 WARNING! read_smalldump: *** VELOCITY WILL BE MISSING FROM SMALL DUMP FILES *** forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000476D07 dump_utils_mp_rea 1225 utils_dumpfiles.f90 phantom 0000000000731540 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 000000000096CD38 testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000047C50B dump_utils_mp_rea 1246 utils_dumpfiles.f90 phantom 0000000000731540 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 000000000096CD38 testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000047E276 dump_utils_mp_rea 1253 utils_dumpfiles.f90 phantom 0000000000731540 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 000000000096CD38 testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000047FFE5 dump_utils_mp_rea 1260 utils_dumpfiles.f90 phantom 0000000000731540 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 000000000096CD38 testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Read routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 00000000004845FF dump_utils_mp_rea 1281 utils_dumpfiles.f90 phantom 0000000000731540 readwrite_dumps_m 1022 readwrite_dumps.F90 phantom 000000000096CD38 testrwdump_mp_tes 201 test_rwdump.F90 phantom 00000000009C29CE test_mp_testsuite 251 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown n(gas) = 991 n(dust) = 10 npart = 1001 npart(total) = 1001 adiabatic eos: gamma = 1.29999995231628 comment: hfact was 1.12300002574921 now 1.12300000000000 comment: tolh was 9.999999747378752E-005 now 1.000000000000000E-004 comment: C_cour was 0.300000011920929 now 0.300000000000000 External field found, Bext = 124.000000000000 125.000000000000 126.000000000000 ERROR reading values to verify conservation laws. Resetting initial values. ERROR reading grain size/density from file header time = 20.0000000000000 --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1001) 31.281 KB real(8) xyzh_soa(1001,4) 31.281 KB real(8) vxyzu(4,1001) 31.281 KB real(4) alphaind(2,1001) 7.820 KB real(4) divcurlv(4,1001) 15.641 KB real(4) divcurlB(4,1001) 15.641 KB real(8) Bevol(4,1001) 31.281 KB real(8) Bxyz(3,1001) 23.461 KB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1001) 172.047 KB real(8) dustevol(11,1001) 86.023 KB real(8) ddustevol(11,1001) 86.023 KB real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,11,1001) 258.070 KB real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1001) 31.281 KB real(8) dBevol(4,1001) 31.281 KB real(4) divBsumm(1001) 3.910 KB real(8) fext(3,1001) 23.461 KB real(8) vpred(4,1001) 31.281 KB real(8) dustpred(11,1001) 86.023 KB real(8) Bpred(4,1001) 31.281 KB real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1001) 1001.000 B integer(1) iphase_soa(1001) 1001.000 B real(4) gradh(1,1001) 3.910 KB real(8) tstop(22,1001) 172.047 KB integer(4) ll(1001) 3.910 KB integer(4) ibelong(1001) 3.910 KB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,2003) 15.648 KB integer(4) inodeparts(1001) 3.910 KB real(8) xyzh_swap(1001,4) 31.281 KB integer(4) inodeparts_swap(1001) 3.910 KB integer(1) iphase_swap(1001) 1001.000 B integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) list(1001) 3.910 KB integer(4) cellatid(2003) 7.824 KB integer(4) ifirstincell(2003) 7.824 KB kdnode nodeglobal(2003) 101.715 KB kdnode node(2003) 101.715 KB integer(4) nodemap(2003) 7.824 KB --------------------------------------------------------- Total memory allocated to arrays: 1.540 MB --------------------------------------------------------- reading particles 1: 1001, from block 1 lims= 1- 1001 ERROR: missing velocity information from file WARNING: u not in file but setting u = (K*rho**(gamma-1))/(gamma-1) WARNING: alpha not found in file WARNING! sink particle velocities not found got 11 sink properties from 10 sink particles ----------------------------------------------- ID| Mass | Racc | Macc | hsoft | ----------------------------------------------- 1| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 2| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 3| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 4| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 5| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 6| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 7| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 8| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 9| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | 10| 1.04E+02 | 1.05E+02 | 1.07E+02 | 1.06E+02 | ----------------------------------------------- WARNING! div B cleaning field (Psi) not found in Phantom dump file: assuming psi=0 <<< finished reading (MHD) file --> checking header checking read without errors...........OK [got 0 should be 0] checking time..........................OK [max err = 0.000E+00, tol = 1.192E-07] checking hfact.........................OK [max err = 2.293E-08, tol = 1.192E-07] checking gamma.........................OK [max err = 3.668E-08, tol = 1.192E-07] checking polyk.........................OK [max err = 1.445E-08, tol = 1.192E-07] checking m(gas)........................OK [max err = 0.000E+00, tol = 1.192E-07] checking m(dust).......................OK [max err = 0.000E+00, tol = 1.192E-07] checking npartoftype(1)................OK [got 991 should be 991] checking npartoftype(2)................OK [got 10 should be 10] checking umass.........................OK [max err = 0.000E+00, tol = 2.225-308] checking udist.........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking xmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking ymax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmin..........................OK [max err = 0.000E+00, tol = 2.225-308] checking zmax..........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bexty.........................OK [max err = 0.000E+00, tol = 2.225-308] checking Bextz.........................OK [max err = 0.000E+00, tol = 2.225-308] --> checking particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking particle type 1...............OK [checked 991 values] checking particle type 2...............OK [checked 10 values] --> checking sink particle arrays checking x.............................OK [max err = 0.000E+00, tol = 2.225-308] checking y.............................OK [max err = 0.000E+00, tol = 2.225-308] checking z.............................OK [max err = 0.000E+00, tol = 2.225-308] checking m.............................OK [max err = 0.000E+00, tol = 2.225-308] checking h.............................OK [max err = 0.000E+00, tol = 2.225-308] checking hsoft.........................OK [max err = 0.000E+00, tol = 2.225-308] checking maccreted.....................OK [max err = 0.000E+00, tol = 2.225-308] checking spinx.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spiny.........................OK [max err = 0.000E+00, tol = 2.225-308] checking spinz.........................OK [max err = 0.000E+00, tol = 2.225-308] checking tlast.........................OK [max err = 0.000E+00, tol = 2.225-308] <-- READ/WRITE TEST COMPLETE --> ALLOCATING ALL ARRAYS --------------------------------------------------------- real(8) xyzh(4,1000000) 30.518 MB real(8) xyzh_soa(1000000,4) 30.518 MB real(8) vxyzu(4,1000000) 30.518 MB real(4) alphaind(2,1000000) 7.629 MB real(4) divcurlv(4,1000000) 15.259 MB real(4) divcurlB(4,1000000) 15.259 MB real(8) Bevol(4,1000000) 30.518 MB real(8) Bxyz(3,1000000) 22.888 MB real(8) dustprop(4,0) 0.000 B real(4) St(0) 0.000 B real(4) csound(0) 0.000 B real(8) abundance(5,0) 0.000 B real(4) temperature(0) 0.000 B real(8) dustfrac(22,1000000) 167.847 MB real(8) dustevol(11,1000000) 83.923 MB real(8) ddustevol(11,1000000) 83.923 MB real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,11,1000000) 251.770 MB real(8) xyzmh_ptmass(11,100) 8.594 KB real(8) vxyz_ptmass(3,100) 2.344 KB real(8) fxyz_ptmass(4,100) 3.125 KB real(8) fxyz_ptmass_sinksink(4,100) 3.125 KB real(4) poten(0) 0.000 B real(8) n_R(4,0) 0.000 B real(4) n_electronT(0) 0.000 B real(8) eta_nimhd(4,0) 0.000 B real(4) luminosity(0) 0.000 B real(8) fxyzu(4,1000000) 30.518 MB real(8) dBevol(4,1000000) 30.518 MB real(4) divBsumm(1000000) 3.815 MB real(8) fext(3,1000000) 22.888 MB real(8) vpred(4,1000000) 30.518 MB real(8) dustpred(11,1000000) 83.923 MB real(8) Bpred(4,1000000) 30.518 MB real(8) dustproppred(4,0) 0.000 B integer(1) iphase(1000000) 976.562 KB integer(1) iphase_soa(1000000) 976.562 KB real(4) gradh(1,1000000) 3.815 MB real(8) tstop(22,1000000) 167.847 MB integer(4) ll(1000000) 3.815 MB integer(4) ibelong(1000000) 3.815 MB integer(1) istsactive(1) 1.000 B integer(1) ibin_sts(1) 1.000 B integer(4) inoderange(2,1000001) 7.629 MB integer(4) inodeparts(1000000) 3.815 MB real(8) xyzh_swap(1000000,4) 30.518 MB integer(4) inodeparts_swap(1000000) 3.815 MB integer(1) iphase_swap(1000000) 976.562 KB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) list(1000000) 3.815 MB integer(4) cellatid(1000001) 3.815 MB integer(4) ifirstincell(1000001) 3.815 MB kdnode nodeglobal(1000001) 49.591 MB kdnode node(1000001) 49.591 MB integer(4) nodemap(1000001) 3.815 MB --------------------------------------------------------- Total memory allocated to arrays: 1.370 GB --------------------------------------------------------- --> TESTING STEP MODULE / boundary crossing ---------------- particles set on 50 x 50 x 50 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 2.000E-02 dy: 2.000E-02 dz: 2.000E-02 ----------------------------------------------------------------------------------------- t = 0.200000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.400000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.600000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 0.800000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.00000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.20000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.40000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.60000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 1.80000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] t = 2.00000000000000 dt = 0.200000000000000 checking h (density)...................OK [max err = 2.750E-04, tol = 3.000E-04] checking fx............................OK [max err = 0.000E+00, tol = 2.225-308] checking fy............................OK [max err = 0.000E+00, tol = 2.225-308] checking fz............................OK [max err = 0.000E+00, tol = 2.225-308] checking du/dt.........................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolx/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevoly/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dBevolz/dt....................OK [max err = 0.000E+00, tol = 2.225-308] checking dpsi/dt.......................OK [max err = 0.000E+00, tol = 2.225-308] <-- STEP TEST COMPLETE --> SKIPPING TEST OF IND TIMESTEP UTILS (need -DIND_TIMESTEPS) --> TESTING EXTERNAL FORCES MODULE ----------------------------------------------------------------------------------------- random seed = -43587 1000 particles set in uniform random distribution ----------------------------------------------------------------------------------------- checking star external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking corotate external force initialisation..OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 3.537E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 3.552E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking binary external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 5.620E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 4.658E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 5.188E-11, tol = 1.500E-03] checking prdrag external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking torus external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 2.704E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 1.051E-09, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking toystar external force initialisation...OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking external B field external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] No spiral arm spheroids set. No potential to be read in. checking spiral external force initialisation....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.338E-05, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 8.375E-05, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 6.864E-05, tol = 1.500E-03] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] ERROR! extern_neutronstar: Error opening ns-rdensity.tab checking neutronstar external force initialisation.........OK [got 1 should be 1] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking Einstein-prec external force initialisation.......OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 1.619E-10, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 1.526E-10, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 9.963E-11, tol = 1.500E-03] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.031E-11, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 9.824E-11, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 4.342E-11, tol = 1.500E-03] checking static sinusoid external force initialisation.....OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 7.196E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] ERROR! externalforces: Require number of particles per star : iexternalforce = 14 checking grav. wave inspiral external force initialisation...........OK [got 2 should be 2] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking disc gravity external force initialisation........OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 8.134E-10, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 0.000E+00, tol = 1.500E-03] checking corotating binary external force initialisation...OK [got 0 should be 0] checking fextx = -grad phi.............OK [max err = 4.960E-06, tol = 1.500E-03] checking fexty = -grad phi.............OK [max err = 4.744E-06, tol = 1.500E-03] checking fextz = -grad phi.............OK [max err = 1.824E-06, tol = 1.500E-03] --> testing accrete_particles routine checking accreted=was_accreted for all externf...OK [checked 16 values] --> testing velocity-dependent external force solvers checking corotate external force initialisation..OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking prdrag external force initialisation....OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking Lense-Thirring external force initialisation......OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.417E-16, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.184E-16, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.766E-16, tol = 1.000E-14] checking Einstein-prec external force initialisation.......OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.417E-16, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.184E-16, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 1.766E-16, tol = 1.000E-14] checking generalised Newtonian external force initialisation.........OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking corotating binary external force initialisation...OK [got 0 should be 0] checking fx(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fy(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] checking fz(x1,v1)=v12 + dt/2*[f(x1) + f(x1,v1)]...........OK [max err = 0.000E+00, tol = 1.000E-14] <-- EXTERNAL FORCE TESTS COMPLETE --> TESTING PTMASS MODULE --> testing integration of binary orbit dt for sinks = 1.250000000000000E-002 checking potential energy..............OK [max err = 0.000E+00, tol = 2.220E-16] checking angular momentum..............OK [max err = 0.000E+00, tol = 2.220E-10] nsteps per orbit = 711 norbits = 100 checking angular momentum..............OK [max err = 6.594E-15, tol = 3.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] checking total energy..................OK [max err = 2.440E-08, tol = 3.000E-08] --> skipping circumbinary disc test (-DPERIODIC is set) --> skipping circumbinary disc test (-DPERIODIC is set) --> testing softening in sink particle binary checking potential energy..............OK [max err = 3.391E-16, tol = 4.441E-16] nsteps per orbit = 1099 norbits = 10 checking angular momentum..............OK [max err = 6.324E-15, tol = 2.000E-14] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.225-308] checking total energy..................OK [max err = 1.021E-10, tol = 2.000E-09] --> testing accretion onto sink particles checking accretion flag................OK [got T should be T] checking isdead_or_accreted flag.......OK [got T should be T] checking x(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking y(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking z(ptmass) after accretion.....OK [max err = 0.000E+00, tol = 2.225-308] checking vx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking vz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fx(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fy(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking fz(ptmass) after accretion....OK [max err = 0.000E+00, tol = 2.225-308] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (uniform density) forrtl: warning (406): fort: (1): In call to GETNEIGH_POS, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000775E3C ptmass_mp_ptmass_ 964 ptmass.F90 phantom 00000000009A374D testptmass_mp_tes 594 test_ptmass.F90 phantom 00000000009C2A9A test_mp_testsuite 283 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.500E-02 -2.500E-02 -2.500E-02 0.000E+00 ) by accreting 93 particles: M= 3.321E-01 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] --> testing sink particle creation (sin) >>>>>> s t r e t c h m a p p i n g <<<<<< stretching to match density profile in r direction density at r = 0.000000000000000E+000 is 1.00000000000000 total mass = 6.960406353085308E-004 >>>>>> done forrtl: warning (406): fort: (1): In call to GETNEIGH_POS, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000775E3C ptmass_mp_ptmass_ 964 ptmass.F90 phantom 00000000009A374D testptmass_mp_tes 594 test_ptmass.F90 phantom 00000000009C2A9A test_mp_testsuite 283 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown ptmass_create: Testing 110 on thread 0 for ptmass creation... ptmass_create: OK created ptmass # 1 at (x,y,z,t)=(-2.765E-18 6.370E-19 1.870E-18 0.000E+00 ) by accreting 280 particles: M= 1.000E+00 checking nptmass=1.....................OK [got 1 should be 1] checking angular momentum..............OK [max err = 0.000E+00, tol = 1.000E-10] checking linear momentum...............OK [max err = 0.000E+00, tol = 2.220E-16] <-- PTMASS TEST COMPLETE --> TESTING GNEWTON MODULE --> testing relativistic orbit ---------- orbital parameters ---------- apocenter : 180.000 pericenter : 20.000 semi-major axis : 100.000 eccentricity : 0.800 initial velocity : 36.724 gravitational radius : 2.122 central mass : 1.000E+06 ---------------------------------------- integrating the orbit for one azimuthal period... ( C_force = 1.000000000000000E-004 ) checking energy........................OK [max err = 9.701E-06, tol = 1.000E-03] checking angular momentum..............OK [max err = 5.863E-07, tol = 1.000E-03] checking precession angle..............OK [max err = 5.762E-07, tol = 1.000E-05] <-- GNEWTON TEST COMPLETE --> TESTING COROTATION MODULE --> testing centrifugal force with rotation about z axis checking fx=Omega^2 x..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=Omega^2 y..................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] checking phi=1/2 Omega^2 R^2...........OK [max err = 2.321E-16, tol = 4.441E-16] --> testing coriolis force with rotation about z axis checking fx=2 Omega vy.................OK [max err = 0.000E+00, tol = 2.220E-16] checking fy=-2 Omega vx................OK [max err = 0.000E+00, tol = 2.220E-16] checking fz=0..........................OK [max err = 0.000E+00, tol = 2.220E-16] <-- COROTATE TEST COMPLETE --> SKIPPING DISC SETUP TESTS --> TESTING GEOMETRY MODULE --> checking transforms for coords and vectors are reversible checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 8.882E-16, L2 err = 7.326E-17 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 8.882E-16, L2 err = 1.433E-16 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 6.410E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 2.538E-16, L2 err = 1.465E-16 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 8.882E-16, L2 err = 1.036E-16 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 8.882E-16, L2 err = 1.282E-16 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 8.882E-16, L2 err = 1.794E-16 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 4.441E-16, L2 err = 1.110E-16 tol = 1.000E-15] --> checking transforms for coords and vectors are reversible at r=0 checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 1.225E-16, L2 err = 1.500E-32 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 4.441E-16, L2 err = 9.065E-17 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 4.441E-16, L2 err = 7.166E-17 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 1.000E-15, L2 err = 1.000E-30 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] --> testing get_coord_limits matches coord_transform checking cartesian x,y,z.............OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking cylindrical R,phi,z...........OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking spherical r,phi,theta.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking toroidal r,theta,phi.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking rotated x_1,x_2,x_3.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking flared cyl R,phi,zdash.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] checking log flared logR,phi,zd.......OK [max err = 0.000E+00, L2 err = 0.000E+00 tol = 1.000E-15] <-- GEOMETRY TEST COMPLETE --> testing Sedov blast wave ---------------- particles set on 16 x 16 x 16 uniform cubic lattice -------------- x: -0.500 -> 0.500 y: -0.500 -> 0.500 z: -0.500 -> 0.500 dx: 6.250E-02 dy: 6.250E-02 dz: 6.250E-02 ----------------------------------------------------------------------------------------- npart = 4096 particle mass = 2.441406250000000E-004 t = 0.38059E-03 dt = 2.822E-04 (tolv), np = 4096 t = 0.66275E-03 dt = 3.979E-04 (tolv) t = 0.10606E-02 dt = 4.818E-04 (tolv) t = 0.15424E-02 dt = 5.881E-04 (tolv) t = 0.21305E-02 dt = 6.514E-04 (courant) t = 0.27819E-02 dt = 6.608E-04 (courant) t = 0.34427E-02 dt = 6.721E-04 (courant) t = 0.41147E-02 dt = 6.851E-04 (courant) t = 0.47999E-02 dt = 7.000E-04 (courant) t = 0.54998E-02 dt = 7.166E-04 (courant) t = 0.62164E-02 dt = 7.349E-04 (courant) t = 0.69513E-02 dt = 7.550E-04 (courant) t = 0.77063E-02 dt = 7.768E-04 (courant) t = 0.84831E-02 dt = 8.002E-04 (courant) t = 0.92833E-02 dt = 8.254E-04 (courant) t = 0.10109E-01 dt = 8.521E-04 (courant) t = 0.10961E-01 dt = 8.805E-04 (courant) t = 0.11841E-01 dt = 9.104E-04 (courant) t = 0.12752E-01 dt = 9.418E-04 (courant) t = 0.13693E-01 dt = 9.746E-04 (courant) t = 0.14668E-01 dt = 1.009E-03 (courant) t = 0.15677E-01 dt = 1.044E-03 (courant) t = 0.16721E-01 dt = 1.080E-03 (courant) t = 0.17801E-01 dt = 1.118E-03 (courant) t = 0.18919E-01 dt = 1.156E-03 (courant) t = 0.20074E-01 dt = 1.194E-03 (courant) t = 0.21268E-01 dt = 1.233E-03 (courant) t = 0.22501E-01 dt = 1.271E-03 (courant) t = 0.23773E-01 dt = 1.311E-03 (courant) t = 0.25084E-01 dt = 1.352E-03 (courant) t = 0.26436E-01 dt = 1.392E-03 (courant) t = 0.27828E-01 dt = 1.432E-03 (courant) t = 0.29260E-01 dt = 1.473E-03 (courant) t = 0.30732E-01 dt = 1.513E-03 (courant) t = 0.32245E-01 dt = 1.553E-03 (courant) t = 0.33799E-01 dt = 1.593E-03 (courant) t = 0.35392E-01 dt = 1.632E-03 (courant) t = 0.37023E-01 dt = 1.670E-03 (courant) t = 0.38693E-01 dt = 1.707E-03 (courant) t = 0.40400E-01 dt = 1.742E-03 (courant) t = 0.42142E-01 dt = 1.776E-03 (courant) t = 0.43918E-01 dt = 1.809E-03 (courant) t = 0.45726E-01 dt = 1.842E-03 (courant) t = 0.47569E-01 dt = 1.876E-03 (courant) t = 0.49445E-01 dt = 1.909E-03 (courant) t = 0.51353E-01 dt = 1.942E-03 (courant) t = 0.53295E-01 dt = 1.974E-03 (courant) t = 0.55269E-01 dt = 2.006E-03 (courant) t = 0.57275E-01 dt = 2.038E-03 (courant) t = 0.59313E-01 dt = 2.069E-03 (courant) t = 0.61382E-01 dt = 2.100E-03 (courant) t = 0.63482E-01 dt = 2.131E-03 (courant) t = 0.65613E-01 dt = 2.162E-03 (courant) t = 0.67774E-01 dt = 2.192E-03 (courant) t = 0.69966E-01 dt = 2.223E-03 (courant) t = 0.72189E-01 dt = 2.253E-03 (courant) t = 0.74442E-01 dt = 2.284E-03 (courant) t = 0.76726E-01 dt = 2.291E-03 (courant) t = 0.79016E-01 dt = 2.268E-03 (courant) t = 0.81284E-01 dt = 2.251E-03 (courant) t = 0.83535E-01 dt = 2.241E-03 (courant) t = 0.85776E-01 dt = 2.054E-03 (courant) t = 0.87830E-01 dt = 2.032E-03 (courant) t = 0.89863E-01 dt = 2.017E-03 (courant) t = 0.91880E-01 dt = 2.009E-03 (courant) t = 0.93890E-01 dt = 2.009E-03 (courant) t = 0.95898E-01 dt = 2.015E-03 (courant) t = 0.97913E-01 dt = 2.028E-03 (courant) t = 0.99942E-01 dt = 5.849E-05 (dtprint) t = 0.10000 dt = 1.823E-03 (courant) --------> TIME = 0.1000 : full dump written to file test002 <-------- forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000048F0F2 dump_utils_mp_wri 1470 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 00000000008721E1 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C079D testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C2B99 test_mp_testsuite 323 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000049284F dump_utils_mp_wri 1494 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 00000000008721E1 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C079D testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C2B99 test_mp_testsuite 323 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000493AC9 dump_utils_mp_wri 1502 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 00000000008721E1 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C079D testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C2B99 test_mp_testsuite 323 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000494D47 dump_utils_mp_wri 1510 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 00000000008721E1 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C079D testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C2B99 test_mp_testsuite 323 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 0000000000497828 dump_utils_mp_wri 1534 utils_dumpfiles.f90 phantom 0000000000721170 readwrite_dumps_m 412 readwrite_dumps.F90 phantom 00000000008721E1 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C079D testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C2B99 test_mp_testsuite 323 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown forrtl: warning (406): fort: (1): In call to I/O Write routine, an array temporary was created for argument #1 Image PC Routine Line Source libifcoremt.so.5 00002AE6D45385A3 for_emit_diagnost Unknown Unknown phantom 000000000049C933 dump_utils_mp_wri 1681 utils_dumpfiles.f90 phantom 0000000000724391 readwrite_dumps_m 539 readwrite_dumps.F90 phantom 00000000008721E1 evolve_mp_evol_ 542 evolve.F90 phantom 00000000009C079D testsedov_mp_test 164 test_sedov.F90 phantom 00000000009C2B99 test_mp_testsuite 323 testsuite.F90 phantom 00000000009E08DF MAIN__ 93 phantom.F90 phantom 00000000004033FE Unknown Unknown Unknown libc-2.17.so 00002AE6D74B6C05 __libc_start_main Unknown Unknown phantom 00000000004032F9 Unknown Unknown Unknown input file test.in written successfully. Since code start: 70 timesteps, wall: 5.9s cpu: 94s cpu/wall: 16 Since last dump : 70 timesteps, wall: 5.9s cpu: 94s cpu/wall: 16 wall cpu cpu/wall frac step : 5.88s 93.95s 15.99 97.92% step (force): 1.88s 29.75s 15.87 31.25% step (dens) : 3.12s 47.96s 15.35 52.08% step (link) : 0.88s 15.56s 17.78 14.58% write_dump : 0.12s 1.41s 11.28 2.08% Etot= 9.669E-01, Ekin= 1.982E-01, Etherm= 7.686E-01, Epot= 0.000E+00 Emag= 0.000E+00 Linm= 1.544E-17, Angm= 1.141E-18 Centre of Mass = -1.166E-18, -4.825E-18, 0.000E+00 density (max)= 1.602E+00 (mean)= 1.159E+00 (max)= 9.455E+00 g/cm^3 Mgas = 1.000E+00 Mdust1 = 0.000E+00 Mdust2 = 0.000E+00 Mdust3 = 0.000E+00 Mdust4 = 0.000E+00 Mdust5 = 0.000E+00 Mdust6 = 0.000E+00 Mdust7 = 0.000E+00 Mdust8 = 0.000E+00 Mdust9 = 0.000E+00 Mdust10 = 0.000E+00 Mdust11 = 0.000E+00 alpha(max)= 1.000E+00 RMS Mach #= 6.243E-01 div B (max)= 0.000E+00, div B (mean)= 0.000E+00 h|div B|/B (max)= 0.000E+00, h|div B|/B (mean)= 0.000E+00 checking total energy..................OK [max err = 2.332E-04, tol = 4.700E-04] checking linear momentum...............OK [max err = 1.544E-17, tol = 7.000E-15] <-- testing complete total wall time = 9 min, 26.50 s (= 5.6650E+02s) total cpu time = 1 hr, 47 min, 36.37 s (= 6.4564E+03s) SUMMARY OF ALL TESTS: PASSED: 137 of 137 100.0% FAILED: 0 of 137 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found