ali-ramadhan/gist:ec039bd71c21daff23c140b8aec180eb

## gistfile1.txt
ali-ramadhan@cyclops:~/Oceananigans.jl$ julia --project
               _
   _       _ _(_)_     |  Documentation: https://docs.julialang.org
  (_)     | (_) (_)    |
   _ _   _| |_  __ _   |  Type "?" for help, "]?" for Pkg help.
  | | | | | | |/ _` |  |
  | | |_| | | | (_| |  |  Version 1.1.0 (2019-01-21)
 _/ |\__'_|_|_|\__'_|  |  Official https://julialang.org/ release
|__/                   |

(Oceananigans) pkg> test
   Testing Oceananigans
 Resolving package versions...
Testing grids...
  Testing grid initialization...
  Testing grid constructor errors...
Testing fields...
  Testing field initialization...
  Testing field setting...
┌ Warning: Performing scalar operations on GPU arrays: This is very slow, consider disallowing these operations with `allowscalar(false)`
└ @ GPUArrays ~/.julia/packages/GPUArrays/0lvhc/src/indexing.jl:16
Testing halo regions...
  Testing initializing halo regions...
  Testing filling halo regions...
Testing operators...
  Testing function differentiation...
  Testing function interpolation...
  Testing 2D operators...
Testing Solvers...
  Testing FFTW planning...
  Testing divergence-free solution [CPU]...
  Testing divergence-free solution [GPU]...
  Testing analytic solution reconstruction...
[ Info: Error (ℓ²-norm), Float64, N=(32, 32, 32), m=(1, 1, 1): 0.0010429869370832128
[ Info: Error (ℓ²-norm), Float64, N=(32, 32, 32), m=(2, 2, 2): 0.004193186980784157
[ Info: Error (ℓ²-norm), Float64, N=(32, 32, 32), m=(3, 3, 3): 0.009515129976547675
[ Info: Error (ℓ²-norm), Float64, N=(48, 48, 48), m=(1, 1, 1): 0.0004631146916545021
[ Info: Error (ℓ²-norm), Float64, N=(48, 48, 48), m=(2, 2, 2): 0.001856639437966049
[ Info: Error (ℓ²-norm), Float64, N=(48, 48, 48), m=(3, 3, 3): 0.0041931869807848845
[ Info: Error (ℓ²-norm), Float64, N=(64, 64, 64), m=(1, 1, 1): 0.00026041644318885814
[ Info: Error (ℓ²-norm), Float64, N=(64, 64, 64), m=(2, 2, 2): 0.0010429869370831952
[ Info: Error (ℓ²-norm), Float64, N=(64, 64, 64), m=(3, 3, 3): 0.0023516875273211146
Testing Coriolis...
Testing surface waves...
Testing buoyancy...
Testing models...
  Testing doubly periodic model construction...
  Testing reentrant channel model construction...
  Testing non-dimensional model construction...
  Testing setting model fields...
Testing time stepping...
  Testing that time stepping works [CPU(), Float32, ConstantIsotropicDiffusivity]...
  Testing that time stepping works [CPU(), Float32, ConstantAnisotropicDiffusivity]...
  Testing that time stepping works [CPU(), Float32, SmagorinskyLilly]...
  Testing that time stepping works [CPU(), Float32, VerstappenAnisotropicMinimumDissipation]...
  Testing that time stepping works [CPU(), Float64, ConstantIsotropicDiffusivity]...
  Testing that time stepping works [CPU(), Float64, ConstantAnisotropicDiffusivity]...
  Testing that time stepping works [CPU(), Float64, SmagorinskyLilly]...
  Testing that time stepping works [CPU(), Float64, VerstappenAnisotropicMinimumDissipation]...
  Testing that time stepping works [GPU(), Float32, ConstantIsotropicDiffusivity]...
  Testing that time stepping works [GPU(), Float32, ConstantAnisotropicDiffusivity]...
  Testing that time stepping works [GPU(), Float32, SmagorinskyLilly]...
  Testing that time stepping works [GPU(), Float32, VerstappenAnisotropicMinimumDissipation]...
  Testing that time stepping works [GPU(), Float64, ConstantIsotropicDiffusivity]...
  Testing that time stepping works [GPU(), Float64, ConstantAnisotropicDiffusivity]...
  Testing that time stepping works [GPU(), Float64, SmagorinskyLilly]...
  Testing that time stepping works [GPU(), Float64, VerstappenAnisotropicMinimumDissipation]...
  Testing 2nd-order Adams-Bashforth...
  Testing recomputing w from continuity...
[ Info: Velocity divergence after recomputing w [CPU(), Float32]: min=-2.3841858e-7, max=2.3841858e-7, sum=4.7683716e-6, abs_sum=5.6505203e-5
[ Info: Velocity divergence after recomputing w [CPU(), Float64]: min=-4.440892098500626e-16, max=4.440892098500626e-16, sum=-8.881784197001252e-16, abs_sum=8.79296635503124e-14
[ Info: Velocity divergence after recomputing w [GPU(), Float32]: min=-2.3841858e-7, max=3.5762787e-7, sum=-1.1920929e-7, abs_sum=6.043911e-5
[ Info: Velocity divergence after recomputing w [GPU(), Float64]: min=-4.440892098500626e-16, max=4.440892098500626e-16, sum=-6.8833827526759706e-15, abs_sum=1.0502709812953981e-13
  Testing incompressibility...
[ Info: Velocity divergence after 1 time steps [CPU(), Float32]: min=-6.457305e-19, max=-6.457305e-19, sum=-2.4400536e-22, abs_sum=1.414861e-16
[ Info: Velocity divergence after 10 time steps [CPU(), Float32]: min=-6.826408e-19, max=-6.826408e-19, sum=-6.5854434e-21, abs_sum=1.646294e-16
[ Info: Velocity divergence after 100 time steps [CPU(), Float32]: min=-6.847821e-19, max=-6.847821e-19, sum=-2.1468129e-20, abs_sum=2.6427743e-16
[ Info: Velocity divergence after 1 time steps [CPU(), Float64]: min=-6.037244272213971e-19, max=-6.037244272213971e-19, sum=-1.9275293505266353e-22, abs_sum=1.284399769586337e-16
[ Info: Velocity divergence after 10 time steps [CPU(), Float64]: min=-7.735240399597831e-19, max=-7.735240399597831e-19, sum=1.666113807244092e-21, abs_sum=1.5625275138570133e-16
[ Info: Velocity divergence after 100 time steps [CPU(), Float64]: min=-6.005395833397209e-19, max=-6.005395833397209e-19, sum=8.58925289833789e-20, abs_sum=2.661629631159426e-16
[ Info: Velocity divergence after 1 time steps [GPU(), Float32]: min=-6.798761e-19, max=-6.798761e-19, sum=2.9982712e-22, abs_sum=1.3663385e-16
[ Info: Velocity divergence after 10 time steps [GPU(), Float32]: min=-8.801553e-19, max=-8.801553e-19, sum=4.8376236e-21, abs_sum=1.8132469e-16
[ Info: Velocity divergence after 100 time steps [GPU(), Float32]: min=-7.2684914e-19, max=-7.2684914e-19, sum=4.5685806e-20, abs_sum=2.7775364e-16
[ Info: Velocity divergence after 1 time steps [GPU(), Float64]: min=-6.28322264009662e-19, max=-6.28322264009662e-19, sum=-2.0630546222354245e-22, abs_sum=1.4153546122998988e-16
[ Info: Velocity divergence after 10 time steps [GPU(), Float64]: min=-5.88342308899259e-19, max=-5.88342308899259e-19, sum=6.221033481120677e-21, abs_sum=1.6029424585023134e-16
[ Info: Velocity divergence after 100 time steps [GPU(), Float64]: min=-7.438643342787266e-19, max=-7.438643342787266e-19, sum=-3.1245351358314914e-20, abs_sum=2.678438051320785e-16
  Testing tracer conservation in channel...
[ Info: Tracer conservation after 10 time steps [CPU(), Float32]: ⟨T⟩-T₀=-9.947598300641403e-14 °C
[ Info: Tracer conservation after 10 time steps [CPU(), Float64]: ⟨T⟩-T₀=7.815970093361102e-14 °C
[ Info: Tracer conservation after 10 time steps [GPU(), Float32]: ⟨T⟩-T₀=7.105427357601002e-15 °C
[ Info: Tracer conservation after 10 time steps [GPU(), Float64]: ⟨T⟩-T₀=1.4921397450962104e-13 °C
Testing boundary conditions...
Testing forcings...
  Testing forcing function initialization...
  Testing forcing function time stepping [CPU]...
  Testing forcing function time stepping [GPU]...
Testing turbulence closures...
  Testing closure instantiation...
  Testing calculation of nonlinear diffusivities...
    Calculating diffusivities for ConstantIsotropicDiffusivity [Float32, CPU()]
    Calculating diffusivities for ConstantAnisotropicDiffusivity [Float32, CPU()]
    Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float32, CPU()]
    Calculating diffusivities for TwoDimensionalLeith [Float32, CPU()]
    Calculating diffusivities for SmagorinskyLilly [Float32, CPU()]
    Calculating diffusivities for BlasiusSmagorinsky [Float32, CPU()]
    Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float32, CPU()]
    Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float32, CPU()]
    Calculating diffusivities for ConstantIsotropicDiffusivity [Float32, GPU()]
    Calculating diffusivities for ConstantAnisotropicDiffusivity [Float32, GPU()]
    Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float32, GPU()]
    Calculating diffusivities for TwoDimensionalLeith [Float32, GPU()]
    Calculating diffusivities for SmagorinskyLilly [Float32, GPU()]
    Calculating diffusivities for BlasiusSmagorinsky [Float32, GPU()]
    Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float32, GPU()]
    Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float32, GPU()]
    Calculating diffusivities for ConstantIsotropicDiffusivity [Float64, CPU()]
    Calculating diffusivities for ConstantAnisotropicDiffusivity [Float64, CPU()]
    Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float64, CPU()]
    Calculating diffusivities for TwoDimensionalLeith [Float64, CPU()]
    Calculating diffusivities for SmagorinskyLilly [Float64, CPU()]
    Calculating diffusivities for BlasiusSmagorinsky [Float64, CPU()]
    Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float64, CPU()]
    Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float64, CPU()]
    Calculating diffusivities for ConstantIsotropicDiffusivity [Float64, GPU()]
    Calculating diffusivities for ConstantAnisotropicDiffusivity [Float64, GPU()]
    Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float64, GPU()]
    Calculating diffusivities for TwoDimensionalLeith [Float64, GPU()]
    Calculating diffusivities for SmagorinskyLilly [Float64, GPU()]
    Calculating diffusivities for BlasiusSmagorinsky [Float64, GPU()]
    Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float64, GPU()]
    Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float64, GPU()]
  Testing constant isotropic diffusivity...
  Testing constant anisotropic diffusivity...
  Testing time-stepping with a tuple of closures...
Testing dynamics...
  Testing simple diffusion...
  Testing default model budgets with isotropic diffusion...
  Testing diffusion cosine...
  Testing passive tracer advection...
  Testing internal wave...
  Testing Taylor-Green vortex...
[ Info: Taylor-Green vortex test [CPU(), Float64, Nx=Ny=64, Nt=10]: Δu: (avg=2.33e-06, max=2.33e-06), Δv: (avg=2.33e-06, max=2.33e-06)
Testing diagnostics...
  Testing horizontal average [CPU]
  Testing horizontal average [GPU]
  Testing NaN Checker [CPU]
  Testing NaN Checker [GPU]
  Testing miscellaneous timeseries diagnostics [CPU]
  Testing miscellaneous timeseries diagnostics [GPU]
Testing output writers...
  Testing NetCDF output writer [CPU]...
  Testing JLD2 output writer [CPU]...
  Testing Checkpointer [CPU]...
  Testing NetCDF output writer [GPU]...
  Testing JLD2 output writer [GPU]...
  Testing Checkpointer [GPU]...
Running regression tests...
  Testing thermal bubble regression [CPU]
[ Info: Δu: min=-4.16063e-18, max=4.85011e-18, mean=1.46721e-22, absmean=2.63571e-19, std=4.85389e-19
[ Info: Δv: min=-4.5672e-18, max=5.20417e-18, mean=1.33432e-21, absmean=2.64667e-19, std=5.05737e-19
[ Info: Δw: min=-1.64257e-17, max=2.17383e-17, mean=3.80265e-21, absmean=3.16134e-18, std=4.4084e-18
[ Info: ΔT: min=-1.77636e-15, max=1.77636e-15, mean=1.73472e-18, absmean=1.47451e-17, std=1.61852e-16
[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
  Testing Rayleigh–Bénard tracer regression [CPU]
[ Info: Δu: min=-1.40998e-14, max=1.70697e-14, mean=8.91247e-18, absmean=1.06034e-15, std=1.68908e-15
[ Info: Δv: min=-1.01724e-14, max=1.04361e-14, mean=3.9698e-17, absmean=1.06877e-15, std=1.62495e-15
[ Info: Δw: min=-1.93179e-14, max=1.19904e-14, mean=-1.6183e-17, absmean=8.16755e-16, std=1.47607e-15
[ Info: Δb: min=-1.02141e-14, max=1.28786e-14, mean=3.16486e-19, absmean=7.39946e-16, std=1.2041e-15
[ Info: Δc: min=-9.76996e-15, max=1.03528e-14, mean=-9.8538e-18, absmean=3.56372e-16, std=6.30339e-16
  Testing oceanic large eddy simulation regression [VerstappenAnisotropicMinimumDissipation, CPU]
[ Info: Δu: min=0, max=0, mean=0, absmean=0, std=0
[ Info: Δv: min=0, max=0, mean=0, absmean=0, std=0
[ Info: Δw: min=0, max=0, mean=0, absmean=0, std=0
[ Info: ΔT: min=0, max=0, mean=0, absmean=0, std=0
[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
  Testing oceanic large eddy simulation regression [SmagorinskyLilly, CPU]
[ Info: Δu: min=0, max=0, mean=0, absmean=0, std=0
[ Info: Δv: min=0, max=0, mean=0, absmean=0, std=0
[ Info: Δw: min=0, max=0, mean=0, absmean=0, std=0
[ Info: ΔT: min=0, max=0, mean=0, absmean=0, std=0
[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
  Testing thermal bubble regression [GPU]
[ Info: Δu: min=-2.04934e-15, max=2.04991e-15, mean=-5.17854e-22, absmean=3.65673e-16, std=5.04904e-16
[ Info: Δv: min=-2.04934e-15, max=2.05037e-15, mean=2.25924e-22, absmean=3.63129e-16, std=5.02116e-16
[ Info: Δw: min=-5.55036e-15, max=5.99271e-15, mean=5.11368e-21, absmean=8.03186e-16, std=1.23162e-15
[ Info: ΔT: min=-1.77636e-15, max=1.77636e-15, mean=3.90313e-18, absmean=2.29851e-17, std=2.02051e-16
[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
  Testing Rayleigh–Bénard tracer regression [GPU]
[ Info: Δu: min=-2.53131e-14, max=3.56104e-14, mean=9.40878e-18, absmean=2.10036e-15, std=3.37324e-15
[ Info: Δv: min=-2.24959e-14, max=3.35287e-14, mean=3.68487e-17, absmean=2.15144e-15, std=3.30959e-15
[ Info: Δw: min=-3.69704e-14, max=2.26208e-14, mean=-2.77684e-17, absmean=1.65677e-15, std=2.97033e-15
[ Info: Δb: min=-4.09117e-14, max=1.73195e-14, mean=7.57425e-18, absmean=1.48783e-15, std=2.45749e-15
[ Info: Δc: min=-1.31561e-14, max=1.28786e-14, mean=-9.45242e-18, absmean=7.17898e-16, std=1.24366e-15
  Testing oceanic large eddy simulation regression [VerstappenAnisotropicMinimumDissipation, GPU]
[ Info: Δu: min=-7.66748e-16, max=8.53484e-16, mean=5.59111e-21, absmean=1.24728e-16, std=1.65753e-16
[ Info: Δv: min=-8.1532e-16, max=6.45317e-16, mean=-3.60471e-19, absmean=1.2172e-16, std=1.6096e-16
[ Info: Δw: min=-2.06953e-15, max=2.83541e-15, mean=-1.47763e-18, absmean=2.47498e-16, std=3.86671e-16
[ Info: ΔT: min=-1.06581e-14, max=1.06581e-14, mean=8.52846e-17, absmean=1.21957e-15, std=2.35339e-15
[ Info: ΔS: min=-2.13163e-14, max=2.84217e-14, mean=-1.82753e-17, absmean=2.23568e-15, std=4.52704e-15
  Testing oceanic large eddy simulation regression [SmagorinskyLilly, GPU]
[ Info: Δu: min=-9.15067e-16, max=7.89299e-16, mean=4.29501e-19, absmean=1.21116e-16, std=1.60125e-16
[ Info: Δv: min=-7.42462e-16, max=8.7777e-16, mean=-4.19436e-19, absmean=1.22525e-16, std=1.61836e-16
[ Info: Δw: min=-2.33841e-15, max=2.19703e-15, mean=-4.90932e-18, absmean=2.40275e-16, std=3.74179e-16
[ Info: ΔT: min=-1.42109e-14, max=1.06581e-14, mean=-1.21835e-18, absmean=1.11845e-15, std=2.27451e-15
[ Info: ΔS: min=-2.13163e-14, max=2.13163e-14, mean=9.13764e-17, absmean=2.20644e-15, std=4.45287e-15
Testing examples...
  Testing wind and convection-driving mixing example [CPU]
  Testing wind and convection-driving mixing example [GPU]
  Testing ocean convection with plankton example
  Testing simple diffusion example
  Testing internal wave example
  Testing two-dimensional turbulence example
  Testing Eady turbulence example
Testing abstract operations...
  Testing simple binary operations...
  Testing derivatives...
  Testing combined binary operations and derivatives...
  Testing combined binary operations and derivatives...
    Testing computation of abstract operations [Float32, CPU]...
      Testing compute! derivatives...
      Testing compute! unary operations...
      Testing compute! binary operations...
      Testing compute! multiary operations...
      Testing compute! kinetic energy...
      Testing horizontal averges...
    Testing computation of abstract operations [Float64, CPU]...
      Testing compute! derivatives...
      Testing compute! unary operations...
      Testing compute! binary operations...
      Testing compute! multiary operations...
      Testing compute! kinetic energy...
      Testing horizontal averges...
  Testing combined binary operations and derivatives...
    Testing computation of abstract operations [Float32, GPU]...
      Testing compute! derivatives...
      Testing compute! unary operations...
      Testing compute! binary operations...
      Testing compute! multiary operations...
      Testing compute! kinetic energy...
      Testing horizontal averges...
    Testing computation of abstract operations [Float64, GPU]...
      Testing compute! derivatives...
      Testing compute! unary operations...
      Testing compute! binary operations...
      Testing compute! multiary operations...
      Testing compute! kinetic energy...
      Testing horizontal averges...
Testing verification scripts...
  Testing stratified Couette flow verification [CPU]
Simulating stratified plane Couette flow

        N : 16, 16, 8
        L : 12.6, 6.28, 2
       Re : 4250.000
       Ri : 0.010
       Pr : 0.700
        ν : 0.000235
        κ : 0.000336
   U_wall : 1.000
   Θ_wall : 0.010

[ Info: Calculating JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u]...
[ Info: Calculation time: 174.274 ms
[ Info: Writing JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_fields.jld2...
[ Info: Writing done: time=20.204 ms, size=157.9 KiB, Δsize=153.2 KiB
[ Info: Calculating JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u]...
[ Info: Calculation time: 197.198 ms
[ Info: Writing JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_profiles.jld2...
[ Info: Writing done: time=3.872 ms, size=6.4 KiB, Δsize=1.8 KiB
[ Info: Calculating JLD2 output Symbol[:Nu, :Re_tau]...
[ Info: Calculation time: 426.566 ms
[ Info: Writing JLD2 output Symbol[:Nu, :Re_tau] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_statistics.jld2...
[ Info: Writing done: time=536.404 ms, size=9.7 KiB, Δsize=5.1 KiB
[10000000000000.00%] i: 1, t: 0.00, umax: (  1.57,  0.354,  0.439) m/s, CFL: 0.0002003, νκmax: (0.0963, 0.0912), νκCFL: (0.0001541, 0.0001459), next Δt: 5e-05, ⟨wall time⟩: 10.360 s
  Testing stratified Couette flow verification [GPU]
Simulating stratified plane Couette flow

        N : 16, 16, 8
        L : 12.6, 6.28, 2
       Re : 4250.000
       Ri : 0.010
       Pr : 0.700
        ν : 0.000235
        κ : 0.000336
   U_wall : 1.000
   Θ_wall : 0.010

[ Info: Calculating JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u]...
[ Info: Calculation time: 139.868 ms
[ Info: Writing JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_fields.jld2...
[ Info: Writing done: time=20.577 ms, size=157.9 KiB, Δsize=153.2 KiB
[ Info: Calculating JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u]...
[ Info: Calculation time: 215.113 ms
[ Info: Writing JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_profiles.jld2...
[ Info: Writing done: time=3.804 ms, size=6.4 KiB, Δsize=1.8 KiB
[ Info: Calculating JLD2 output Symbol[:Nu, :Re_tau]...
[ Info: Calculation time: 151.567 ms
[ Info: Writing JLD2 output Symbol[:Nu, :Re_tau] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_statistics.jld2...
[ Info: Writing done: time=19.365 ms, size=9.7 KiB, Δsize=5.1 KiB
[10000000000000.00%] i: 1, t: 0.00, umax: (  1.57,  0.354,  0.439) m/s, CFL: 0.0002003, νκmax: (0.0963, 0.0912), νκCFL: (0.0001541, 0.0001459), next Δt: 5e-05, ⟨wall time⟩: 27.217 s
Test Summary: | Pass  Broken  Total
Oceananigans  | 1941       7   1948
   Testing Oceananigans tests passed
	ali-ramadhan@cyclops:~/Oceananigans.jl$ julia --project
	_
	_ _ _(_)_ \| Documentation: https://docs.julialang.org
	(_) \| (_) (_) \|
	_ _ _\| \|_ __ _ \| Type "?" for help, "]?" for Pkg help.
	\| \| \| \| \| \| \|/ _` \| \|
	\| \| \|_\| \| \| \| (_\| \| \| Version 1.1.0 (2019-01-21)
	_/ \|\__'_\|_\|_\|\__'_\| \| Official https://julialang.org/ release
	\|__/ \|

	(Oceananigans) pkg> test
	Testing Oceananigans
	Resolving package versions...
	Testing grids...
	Testing grid initialization...
	Testing grid constructor errors...
	Testing fields...
	Testing field initialization...
	Testing field setting...
	┌ Warning: Performing scalar operations on GPU arrays: This is very slow, consider disallowing these operations with `allowscalar(false)`
	└ @ GPUArrays ~/.julia/packages/GPUArrays/0lvhc/src/indexing.jl:16
	Testing halo regions...
	Testing initializing halo regions...
	Testing filling halo regions...
	Testing operators...
	Testing function differentiation...
	Testing function interpolation...
	Testing 2D operators...
	Testing Solvers...
	Testing FFTW planning...
	Testing divergence-free solution [CPU]...
	Testing divergence-free solution [GPU]...
	Testing analytic solution reconstruction...
	[ Info: Error (ℓ²-norm), Float64, N=(32, 32, 32), m=(1, 1, 1): 0.0010429869370832128
	[ Info: Error (ℓ²-norm), Float64, N=(32, 32, 32), m=(2, 2, 2): 0.004193186980784157
	[ Info: Error (ℓ²-norm), Float64, N=(32, 32, 32), m=(3, 3, 3): 0.009515129976547675
	[ Info: Error (ℓ²-norm), Float64, N=(48, 48, 48), m=(1, 1, 1): 0.0004631146916545021
	[ Info: Error (ℓ²-norm), Float64, N=(48, 48, 48), m=(2, 2, 2): 0.001856639437966049
	[ Info: Error (ℓ²-norm), Float64, N=(48, 48, 48), m=(3, 3, 3): 0.0041931869807848845
	[ Info: Error (ℓ²-norm), Float64, N=(64, 64, 64), m=(1, 1, 1): 0.00026041644318885814
	[ Info: Error (ℓ²-norm), Float64, N=(64, 64, 64), m=(2, 2, 2): 0.0010429869370831952
	[ Info: Error (ℓ²-norm), Float64, N=(64, 64, 64), m=(3, 3, 3): 0.0023516875273211146
	Testing Coriolis...
	Testing surface waves...
	Testing buoyancy...
	Testing models...
	Testing doubly periodic model construction...
	Testing reentrant channel model construction...
	Testing non-dimensional model construction...
	Testing setting model fields...
	Testing time stepping...
	Testing that time stepping works [CPU(), Float32, ConstantIsotropicDiffusivity]...
	Testing that time stepping works [CPU(), Float32, ConstantAnisotropicDiffusivity]...
	Testing that time stepping works [CPU(), Float32, SmagorinskyLilly]...
	Testing that time stepping works [CPU(), Float32, VerstappenAnisotropicMinimumDissipation]...
	Testing that time stepping works [CPU(), Float64, ConstantIsotropicDiffusivity]...
	Testing that time stepping works [CPU(), Float64, ConstantAnisotropicDiffusivity]...
	Testing that time stepping works [CPU(), Float64, SmagorinskyLilly]...
	Testing that time stepping works [CPU(), Float64, VerstappenAnisotropicMinimumDissipation]...
	Testing that time stepping works [GPU(), Float32, ConstantIsotropicDiffusivity]...
	Testing that time stepping works [GPU(), Float32, ConstantAnisotropicDiffusivity]...
	Testing that time stepping works [GPU(), Float32, SmagorinskyLilly]...
	Testing that time stepping works [GPU(), Float32, VerstappenAnisotropicMinimumDissipation]...
	Testing that time stepping works [GPU(), Float64, ConstantIsotropicDiffusivity]...
	Testing that time stepping works [GPU(), Float64, ConstantAnisotropicDiffusivity]...
	Testing that time stepping works [GPU(), Float64, SmagorinskyLilly]...
	Testing that time stepping works [GPU(), Float64, VerstappenAnisotropicMinimumDissipation]...
	Testing 2nd-order Adams-Bashforth...
	Testing recomputing w from continuity...
	[ Info: Velocity divergence after recomputing w [CPU(), Float32]: min=-2.3841858e-7, max=2.3841858e-7, sum=4.7683716e-6, abs_sum=5.6505203e-5
	[ Info: Velocity divergence after recomputing w [CPU(), Float64]: min=-4.440892098500626e-16, max=4.440892098500626e-16, sum=-8.881784197001252e-16, abs_sum=8.79296635503124e-14
	[ Info: Velocity divergence after recomputing w [GPU(), Float32]: min=-2.3841858e-7, max=3.5762787e-7, sum=-1.1920929e-7, abs_sum=6.043911e-5
	[ Info: Velocity divergence after recomputing w [GPU(), Float64]: min=-4.440892098500626e-16, max=4.440892098500626e-16, sum=-6.8833827526759706e-15, abs_sum=1.0502709812953981e-13
	Testing incompressibility...
	[ Info: Velocity divergence after 1 time steps [CPU(), Float32]: min=-6.457305e-19, max=-6.457305e-19, sum=-2.4400536e-22, abs_sum=1.414861e-16
	[ Info: Velocity divergence after 10 time steps [CPU(), Float32]: min=-6.826408e-19, max=-6.826408e-19, sum=-6.5854434e-21, abs_sum=1.646294e-16
	[ Info: Velocity divergence after 100 time steps [CPU(), Float32]: min=-6.847821e-19, max=-6.847821e-19, sum=-2.1468129e-20, abs_sum=2.6427743e-16
	[ Info: Velocity divergence after 1 time steps [CPU(), Float64]: min=-6.037244272213971e-19, max=-6.037244272213971e-19, sum=-1.9275293505266353e-22, abs_sum=1.284399769586337e-16
	[ Info: Velocity divergence after 10 time steps [CPU(), Float64]: min=-7.735240399597831e-19, max=-7.735240399597831e-19, sum=1.666113807244092e-21, abs_sum=1.5625275138570133e-16
	[ Info: Velocity divergence after 100 time steps [CPU(), Float64]: min=-6.005395833397209e-19, max=-6.005395833397209e-19, sum=8.58925289833789e-20, abs_sum=2.661629631159426e-16
	[ Info: Velocity divergence after 1 time steps [GPU(), Float32]: min=-6.798761e-19, max=-6.798761e-19, sum=2.9982712e-22, abs_sum=1.3663385e-16
	[ Info: Velocity divergence after 10 time steps [GPU(), Float32]: min=-8.801553e-19, max=-8.801553e-19, sum=4.8376236e-21, abs_sum=1.8132469e-16
	[ Info: Velocity divergence after 100 time steps [GPU(), Float32]: min=-7.2684914e-19, max=-7.2684914e-19, sum=4.5685806e-20, abs_sum=2.7775364e-16
	[ Info: Velocity divergence after 1 time steps [GPU(), Float64]: min=-6.28322264009662e-19, max=-6.28322264009662e-19, sum=-2.0630546222354245e-22, abs_sum=1.4153546122998988e-16
	[ Info: Velocity divergence after 10 time steps [GPU(), Float64]: min=-5.88342308899259e-19, max=-5.88342308899259e-19, sum=6.221033481120677e-21, abs_sum=1.6029424585023134e-16
	[ Info: Velocity divergence after 100 time steps [GPU(), Float64]: min=-7.438643342787266e-19, max=-7.438643342787266e-19, sum=-3.1245351358314914e-20, abs_sum=2.678438051320785e-16
	Testing tracer conservation in channel...
	[ Info: Tracer conservation after 10 time steps [CPU(), Float32]: ⟨T⟩-T₀=-9.947598300641403e-14 °C
	[ Info: Tracer conservation after 10 time steps [CPU(), Float64]: ⟨T⟩-T₀=7.815970093361102e-14 °C
	[ Info: Tracer conservation after 10 time steps [GPU(), Float32]: ⟨T⟩-T₀=7.105427357601002e-15 °C
	[ Info: Tracer conservation after 10 time steps [GPU(), Float64]: ⟨T⟩-T₀=1.4921397450962104e-13 °C
	Testing boundary conditions...
	Testing forcings...
	Testing forcing function initialization...
	Testing forcing function time stepping [CPU]...
	Testing forcing function time stepping [GPU]...
	Testing turbulence closures...
	Testing closure instantiation...
	Testing calculation of nonlinear diffusivities...
	Calculating diffusivities for ConstantIsotropicDiffusivity [Float32, CPU()]
	Calculating diffusivities for ConstantAnisotropicDiffusivity [Float32, CPU()]
	Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float32, CPU()]
	Calculating diffusivities for TwoDimensionalLeith [Float32, CPU()]
	Calculating diffusivities for SmagorinskyLilly [Float32, CPU()]
	Calculating diffusivities for BlasiusSmagorinsky [Float32, CPU()]
	Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float32, CPU()]
	Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float32, CPU()]
	Calculating diffusivities for ConstantIsotropicDiffusivity [Float32, GPU()]
	Calculating diffusivities for ConstantAnisotropicDiffusivity [Float32, GPU()]
	Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float32, GPU()]
	Calculating diffusivities for TwoDimensionalLeith [Float32, GPU()]
	Calculating diffusivities for SmagorinskyLilly [Float32, GPU()]
	Calculating diffusivities for BlasiusSmagorinsky [Float32, GPU()]
	Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float32, GPU()]
	Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float32, GPU()]
	Calculating diffusivities for ConstantIsotropicDiffusivity [Float64, CPU()]
	Calculating diffusivities for ConstantAnisotropicDiffusivity [Float64, CPU()]
	Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float64, CPU()]
	Calculating diffusivities for TwoDimensionalLeith [Float64, CPU()]
	Calculating diffusivities for SmagorinskyLilly [Float64, CPU()]
	Calculating diffusivities for BlasiusSmagorinsky [Float64, CPU()]
	Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float64, CPU()]
	Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float64, CPU()]
	Calculating diffusivities for ConstantIsotropicDiffusivity [Float64, GPU()]
	Calculating diffusivities for ConstantAnisotropicDiffusivity [Float64, GPU()]
	Calculating diffusivities for AnisotropicBiharmonicDiffusivity [Float64, GPU()]
	Calculating diffusivities for TwoDimensionalLeith [Float64, GPU()]
	Calculating diffusivities for SmagorinskyLilly [Float64, GPU()]
	Calculating diffusivities for BlasiusSmagorinsky [Float64, GPU()]
	Calculating diffusivities for RozemaAnisotropicMinimumDissipation [Float64, GPU()]
	Calculating diffusivities for VerstappenAnisotropicMinimumDissipation [Float64, GPU()]
	Testing constant isotropic diffusivity...
	Testing constant anisotropic diffusivity...
	Testing time-stepping with a tuple of closures...
	Testing dynamics...
	Testing simple diffusion...
	Testing default model budgets with isotropic diffusion...
	Testing diffusion cosine...
	Testing passive tracer advection...
	Testing internal wave...
	Testing Taylor-Green vortex...
	[ Info: Taylor-Green vortex test [CPU(), Float64, Nx=Ny=64, Nt=10]: Δu: (avg=2.33e-06, max=2.33e-06), Δv: (avg=2.33e-06, max=2.33e-06)
	Testing diagnostics...
	Testing horizontal average [CPU]
	Testing horizontal average [GPU]
	Testing NaN Checker [CPU]
	Testing NaN Checker [GPU]
	Testing miscellaneous timeseries diagnostics [CPU]
	Testing miscellaneous timeseries diagnostics [GPU]
	Testing output writers...
	Testing NetCDF output writer [CPU]...
	Testing JLD2 output writer [CPU]...
	Testing Checkpointer [CPU]...
	Testing NetCDF output writer [GPU]...
	Testing JLD2 output writer [GPU]...
	Testing Checkpointer [GPU]...
	Running regression tests...
	Testing thermal bubble regression [CPU]
	[ Info: Δu: min=-4.16063e-18, max=4.85011e-18, mean=1.46721e-22, absmean=2.63571e-19, std=4.85389e-19
	[ Info: Δv: min=-4.5672e-18, max=5.20417e-18, mean=1.33432e-21, absmean=2.64667e-19, std=5.05737e-19
	[ Info: Δw: min=-1.64257e-17, max=2.17383e-17, mean=3.80265e-21, absmean=3.16134e-18, std=4.4084e-18
	[ Info: ΔT: min=-1.77636e-15, max=1.77636e-15, mean=1.73472e-18, absmean=1.47451e-17, std=1.61852e-16
	[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
	Testing Rayleigh–Bénard tracer regression [CPU]
	[ Info: Δu: min=-1.40998e-14, max=1.70697e-14, mean=8.91247e-18, absmean=1.06034e-15, std=1.68908e-15
	[ Info: Δv: min=-1.01724e-14, max=1.04361e-14, mean=3.9698e-17, absmean=1.06877e-15, std=1.62495e-15
	[ Info: Δw: min=-1.93179e-14, max=1.19904e-14, mean=-1.6183e-17, absmean=8.16755e-16, std=1.47607e-15
	[ Info: Δb: min=-1.02141e-14, max=1.28786e-14, mean=3.16486e-19, absmean=7.39946e-16, std=1.2041e-15
	[ Info: Δc: min=-9.76996e-15, max=1.03528e-14, mean=-9.8538e-18, absmean=3.56372e-16, std=6.30339e-16
	Testing oceanic large eddy simulation regression [VerstappenAnisotropicMinimumDissipation, CPU]
	[ Info: Δu: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: Δv: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: Δw: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: ΔT: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
	Testing oceanic large eddy simulation regression [SmagorinskyLilly, CPU]
	[ Info: Δu: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: Δv: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: Δw: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: ΔT: min=0, max=0, mean=0, absmean=0, std=0
	[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
	Testing thermal bubble regression [GPU]
	[ Info: Δu: min=-2.04934e-15, max=2.04991e-15, mean=-5.17854e-22, absmean=3.65673e-16, std=5.04904e-16
	[ Info: Δv: min=-2.04934e-15, max=2.05037e-15, mean=2.25924e-22, absmean=3.63129e-16, std=5.02116e-16
	[ Info: Δw: min=-5.55036e-15, max=5.99271e-15, mean=5.11368e-21, absmean=8.03186e-16, std=1.23162e-15
	[ Info: ΔT: min=-1.77636e-15, max=1.77636e-15, mean=3.90313e-18, absmean=2.29851e-17, std=2.02051e-16
	[ Info: ΔS: min=0, max=0, mean=0, absmean=0, std=0
	Testing Rayleigh–Bénard tracer regression [GPU]
	[ Info: Δu: min=-2.53131e-14, max=3.56104e-14, mean=9.40878e-18, absmean=2.10036e-15, std=3.37324e-15
	[ Info: Δv: min=-2.24959e-14, max=3.35287e-14, mean=3.68487e-17, absmean=2.15144e-15, std=3.30959e-15
	[ Info: Δw: min=-3.69704e-14, max=2.26208e-14, mean=-2.77684e-17, absmean=1.65677e-15, std=2.97033e-15
	[ Info: Δb: min=-4.09117e-14, max=1.73195e-14, mean=7.57425e-18, absmean=1.48783e-15, std=2.45749e-15
	[ Info: Δc: min=-1.31561e-14, max=1.28786e-14, mean=-9.45242e-18, absmean=7.17898e-16, std=1.24366e-15
	Testing oceanic large eddy simulation regression [VerstappenAnisotropicMinimumDissipation, GPU]
	[ Info: Δu: min=-7.66748e-16, max=8.53484e-16, mean=5.59111e-21, absmean=1.24728e-16, std=1.65753e-16
	[ Info: Δv: min=-8.1532e-16, max=6.45317e-16, mean=-3.60471e-19, absmean=1.2172e-16, std=1.6096e-16
	[ Info: Δw: min=-2.06953e-15, max=2.83541e-15, mean=-1.47763e-18, absmean=2.47498e-16, std=3.86671e-16
	[ Info: ΔT: min=-1.06581e-14, max=1.06581e-14, mean=8.52846e-17, absmean=1.21957e-15, std=2.35339e-15
	[ Info: ΔS: min=-2.13163e-14, max=2.84217e-14, mean=-1.82753e-17, absmean=2.23568e-15, std=4.52704e-15
	Testing oceanic large eddy simulation regression [SmagorinskyLilly, GPU]
	[ Info: Δu: min=-9.15067e-16, max=7.89299e-16, mean=4.29501e-19, absmean=1.21116e-16, std=1.60125e-16
	[ Info: Δv: min=-7.42462e-16, max=8.7777e-16, mean=-4.19436e-19, absmean=1.22525e-16, std=1.61836e-16
	[ Info: Δw: min=-2.33841e-15, max=2.19703e-15, mean=-4.90932e-18, absmean=2.40275e-16, std=3.74179e-16
	[ Info: ΔT: min=-1.42109e-14, max=1.06581e-14, mean=-1.21835e-18, absmean=1.11845e-15, std=2.27451e-15
	[ Info: ΔS: min=-2.13163e-14, max=2.13163e-14, mean=9.13764e-17, absmean=2.20644e-15, std=4.45287e-15
	Testing examples...
	Testing wind and convection-driving mixing example [CPU]
	Testing wind and convection-driving mixing example [GPU]
	Testing ocean convection with plankton example
	Testing simple diffusion example
	Testing internal wave example
	Testing two-dimensional turbulence example
	Testing Eady turbulence example
	Testing abstract operations...
	Testing simple binary operations...
	Testing derivatives...
	Testing combined binary operations and derivatives...
	Testing combined binary operations and derivatives...
	Testing computation of abstract operations [Float32, CPU]...
	Testing compute! derivatives...
	Testing compute! unary operations...
	Testing compute! binary operations...
	Testing compute! multiary operations...
	Testing compute! kinetic energy...
	Testing horizontal averges...
	Testing computation of abstract operations [Float64, CPU]...
	Testing compute! derivatives...
	Testing compute! unary operations...
	Testing compute! binary operations...
	Testing compute! multiary operations...
	Testing compute! kinetic energy...
	Testing horizontal averges...
	Testing combined binary operations and derivatives...
	Testing computation of abstract operations [Float32, GPU]...
	Testing compute! derivatives...
	Testing compute! unary operations...
	Testing compute! binary operations...
	Testing compute! multiary operations...
	Testing compute! kinetic energy...
	Testing horizontal averges...
	Testing computation of abstract operations [Float64, GPU]...
	Testing compute! derivatives...
	Testing compute! unary operations...
	Testing compute! binary operations...
	Testing compute! multiary operations...
	Testing compute! kinetic energy...
	Testing horizontal averges...
	Testing verification scripts...
	Testing stratified Couette flow verification [CPU]
	Simulating stratified plane Couette flow

	N : 16, 16, 8
	L : 12.6, 6.28, 2
	Re : 4250.000
	Ri : 0.010
	Pr : 0.700
	ν : 0.000235
	κ : 0.000336
	U_wall : 1.000
	Θ_wall : 0.010

	[ Info: Calculating JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u]...
	[ Info: Calculation time: 174.274 ms
	[ Info: Writing JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_fields.jld2...
	[ Info: Writing done: time=20.204 ms, size=157.9 KiB, Δsize=153.2 KiB
	[ Info: Calculating JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u]...
	[ Info: Calculation time: 197.198 ms
	[ Info: Writing JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_profiles.jld2...
	[ Info: Writing done: time=3.872 ms, size=6.4 KiB, Δsize=1.8 KiB
	[ Info: Calculating JLD2 output Symbol[:Nu, :Re_tau]...
	[ Info: Calculation time: 426.566 ms
	[ Info: Writing JLD2 output Symbol[:Nu, :Re_tau] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_statistics.jld2...
	[ Info: Writing done: time=536.404 ms, size=9.7 KiB, Δsize=5.1 KiB
	[10000000000000.00%] i: 1, t: 0.00, umax: ( 1.57, 0.354, 0.439) m/s, CFL: 0.0002003, νκmax: (0.0963, 0.0912), νκCFL: (0.0001541, 0.0001459), next Δt: 5e-05, ⟨wall time⟩: 10.360 s
	Testing stratified Couette flow verification [GPU]
	Simulating stratified plane Couette flow

	N : 16, 16, 8
	L : 12.6, 6.28, 2
	Re : 4250.000
	Ri : 0.010
	Pr : 0.700
	ν : 0.000235
	κ : 0.000336
	U_wall : 1.000
	Θ_wall : 0.010

	[ Info: Calculating JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u]...
	[ Info: Calculation time: 139.868 ms
	[ Info: Writing JLD2 output Symbol[:T, :w, :kappaT, :v, :nu, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_fields.jld2...
	[ Info: Writing done: time=20.577 ms, size=157.9 KiB, Δsize=153.2 KiB
	[ Info: Calculating JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u]...
	[ Info: Calculation time: 215.113 ms
	[ Info: Writing JLD2 output Symbol[:T, :w, :nu, :v, :kappaT, :u] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_profiles.jld2...
	[ Info: Writing done: time=3.804 ms, size=6.4 KiB, Δsize=1.8 KiB
	[ Info: Calculating JLD2 output Symbol[:Nu, :Re_tau]...
	[ Info: Calculation time: 151.567 ms
	[ Info: Writing JLD2 output Symbol[:Nu, :Re_tau] to stratified_couette_flow_data_Nxy16_Nz8_Ri0.01/stratified_couette_flow_Nxy16_Nz8_Ri0.01_statistics.jld2...
	[ Info: Writing done: time=19.365 ms, size=9.7 KiB, Δsize=5.1 KiB
	[10000000000000.00%] i: 1, t: 0.00, umax: ( 1.57, 0.354, 0.439) m/s, CFL: 0.0002003, νκmax: (0.0963, 0.0912), νκCFL: (0.0001541, 0.0001459), next Δt: 5e-05, ⟨wall time⟩: 27.217 s
	Test Summary: \| Pass Broken Total
	Oceananigans \| 1941 7 1948
	Testing Oceananigans tests passed