MPI parallelisation is OFF _ \ | | ___| _ \ | | | | __ \ _` | __ \ __| _ \ __ `__ \ \___ \ | | | | ___/ | | | ( | | | | ( | | | | | ___/ ___ | _| _| |_|\__,_|_| _|\__|\___/ _| _| _|_____/ _| _| _| ___ , _ |) (| \ _, o _ |\ /|/ \ ,_ o _ _ |/\_| | | |/ | /|/| | |/ |/ |__// | | / |/ \/ \/|/ (\__/ \/|_/ | |_/|/|_/|_/ | |/|/\__/|_/ (| Version: 1.2.0 phantom_version not found; no git information to display Using NATIVE inverse sqrt 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(3,1000000) 22.888 MB real(4) alphaind(2,1000000) 7.629 MB real(4) divcurlv(1,1000000) 3.815 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(2,0) 0.000 B real(8) dustevol(1,0) 0.000 B real(8) ddustevol(1,0) 0.000 B real(8) ddustprop(4,0) 0.000 B real(8) deltav(3,1,0) 0.000 B 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,1000000) 30.518 MB real(4) n_electronT(1000000) 3.815 MB real(8) eta_nimhd(4,1000000) 30.518 MB real(4) luminosity(0) 0.000 B real(8) fxyzu(3,1000000) 22.888 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(3,1000000) 22.888 MB real(8) dustpred(1,0) 0.000 B 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(2,1000000) 15.259 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: 609.415 MB --------------------------------------------------------- --> RUNNING PHANTOM TEST SUITE "Nobody cares how fast you can calculate the wrong answer." -- Richard West (former UKAFF manager) --> TESTING NON-IDEAL MHD ALGORITHMS --> testing wave dampening test NICIL: Version 1.2.4: 11 April 2018 NICIL: Copyright (c) 2015-2017 James Wurster NICIL: See LICENCE file for usage and distribution conditions NICIL: Reference: Wurster (2016) PASA, 33:e041. NICIL: Non-ideal terms used: Ambipolar diffusion NICIL: All resistivity coefficients are constant. NICIL: eta_ambi = 1.000E-02*|B|^2/rho cm^2 s^{-1} NICIL: Mean molecular mass: 2.310E+00 m_proton checking L2 error on wave damp test....OK [max err = 5.621E-05, tol = 7.150E-05] checking dt to ensure above valid default........OK [got T should be T] dtcourant, dtdiff: 4.51935681570E-03 2.88824049903E-02 checking initial courant dt............OK [max err = 2.897E-05, tol = 4.000E-04] checking initial dissipation dt from sts.........OK [max err = 9.756E-13, tol = 4.000E-04] --> testing standing shock & boundary particles ---------------- particles set on 32 x 15 x 15 uniform closepacked lattice -------------- x: -0.750 -> 0.00 y: -0.148 -> 0.148 z: -0.140 -> 0.140 dx: 2.344E-02 dy: 2.030E-02 dz: 1.914E-02 dy/dx: 0.8660 dz/dx: 0.8165 enforcing periodicity: nx, ny, nz = 32 x 14 x 15, np = 6720 dx: 2.344E-02 dy: 2.114E-02 dz: 1.860E-02 dy/dx: 0.9020 dz/dx: 0.7937 ----------------------------------------------------------------------------------------- ---------------- particles set on 87 x 12 x 12 uniform closepacked lattice -------------- x: 0.00 -> 2.50 y: -0.148 -> 0.148 z: -0.140 -> 0.140 dx: 2.884E-02 dy: 2.497E-02 dz: 2.354E-02 dy/dx: 0.8660 dz/dx: 0.8165 enforcing periodicity: nx, ny, nz = 87 x 12 x 12, np = 12528 dx: 2.884E-02 dy: 2.466E-02 dz: 2.325E-02 dy/dx: 0.8553 dz/dx: 0.8064 ----------------------------------------------------------------------------------------- NICIL: Version 1.2.4: 11 April 2018 NICIL: Copyright (c) 2015-2017 James Wurster NICIL: See LICENCE file for usage and distribution conditions NICIL: Reference: Wurster (2016) PASA, 33:e041. NICIL: Non-ideal terms used: Ohmic Resistivity, Hall Effect, Ambipolar diffusion NICIL: All resistivity coefficients are constant. NICIL: eta_ohm = 1.120E-09 cm^2 s^{-1} NICIL: eta_Hall = -3.530E-02*|B| cm^2 s^{-1} NICIL: eta_ambi = 7.830E-03*|B|^2/rho cm^2 s^{-1} NICIL: Mean molecular mass: 2.310E+00 m_proton checking density error on standing shock, compared to analytics......OK [max err = 2.072E-02, tol = 2.100E-02] checking v_x error on standing shock, compared to analytics..........OK [max err = 3.012E-02, tol = 3.050E-02] checking B_y error on standing shock, compared to analytics..........OK [max err = 1.002E-01, tol = 1.100E-01] checking dt to ensure above valid default........OK [got T should be T] checking Boundary particles are correctly initialised......OK [got T should be T] --> testing calculation of non-constant eta NICIL: Version 1.2.4: 11 April 2018 NICIL: Copyright (c) 2015-2017 James Wurster NICIL: See LICENCE file for usage and distribution conditions NICIL: Reference: Wurster (2016) PASA, 33:e041. NICIL: Using constant grain size. NICIL: Species abundance mass fraction NICIL: H 9.215E-01 0.747 NICIL: He 7.843E-02 0.252 NICIL: Na 1.601E-06 2.955E-05 NICIL: Mg 3.669E-05 7.161E-04 NICIL: K 9.874E-08 3.100E-06 NICIL: Mean molecular mass of cold gas: 2.310E+00 m_proton NICIL: Light element ion mass for CR's: 2.310E+00 m_proton NICIL: Metallic ion mass for CR's: 2.430E+01 m_proton NICIL: Dust-to-gas fraction: 1.000E-02 NICIL: Grain mass: 7.513E+09 m_proton NICIL: Grain radius: 1.000E-05 cm NICIL: Cosmic ray ionisation rate: 1.000E-17 s^{-1} NICIL: Including ionisation from Cosmic rays NICIL: Including thermal ionisation NICIL: Non-ideal terms used: Ohmic Resistivity, Hall Effect, Ambipolar diffusion Used rho,B_z,temp (cgs): 7.43448567310E-18 8.14587176850E-05 1.38195173036E+01 eta_ohm, eta_hall, eta_ambi (cgs): 1.14793940113E+10 3.40040077209E+14 5.26247580402E+17 checking calculated non-constant eta_ohm.........OK [max err = 1.249E-12, tol = 6.300E-05] checking calculated non-constant eta_hall........OK [max err = 1.148E-12, tol = 6.300E-05] checking calculated non-constant eta_ambi........OK [max err = 4.585E-14, tol = 6.300E-05] Used rho,B_z,temp (cgs): 4.66147407857E-03 1.94565874584E+02 7.69440525557E+03 eta_ohm, eta_hall, eta_ambi (cgs): 5.93454638765E+08 1.08059808926E+04 4.17918319187E-03 checking calculated non-constant eta_ohm.........OK [max err = 6.850E-14, tol = 6.300E-05] checking calculated non-constant eta_hall........OK [max err = 3.651E-12, tol = 6.300E-05] checking calculated non-constant eta_ambi........OK [max err = 7.553E-13, tol = 6.300E-05] <-- NON-IDEALMHD TEST COMPLETE <-- testing complete total wall time = 6 min, 6.38 s (= 3.6638E+02s) total cpu time = 1 hr, 37 min, 41.32 s (= 5.8613E+03s) SUMMARY OF ALL TESTS: PASSED: 3 of 3 100.0% FAILED: 0 of 3 0.0% ____ _ ____ ____ | _ \ / \ / ___/ ___| | |_) / _ \ \___ \___ \ | __/ ___ \ ___) |__) | |_| /_/ \_\____/____/ TEST SUITE PASSED sh: say: command not found