Setting initial conditions

Initial conditions are imposed after model construction. This can be easily done using the the set! function, which allows the setting of initial conditions using constant values, arrays, or functions.

julia> grid = RegularCartesianGrid(size=(16, 16, 16), extent=(1, 1, 1));

julia> model = IncompressibleModel(grid=grid);

julia> set!(model, u=0.1, v=1.5)

julia> grid = RegularCartesianGrid(size=(16, 16, 16), extent=(1, 1, 1));

julia> model = IncompressibleModel(grid=grid);

julia> ∂T∂z = 0.01;

julia> ϵ(σ) = σ * randn();

julia> T₀(x, y, z) = ∂T∂z * z + ϵ(1e-8);

julia> set!(model, T=T₀)