Mixed Boundary Conditions

[SamDeAbreu]

Hi there!

I was wondering if it is possible to impose two different kinds of boundary conditions on one boundary. For instance, at z = 0:
ValueBoundaryCondition() if (x,y) satisfy some condition
NoFluxBoundaryCondition() otherwise

Thanks!

[glwagner]

This isn't possible, but it's always possible to compute the flux associated with a ValueBoundaryCondition, so this is the method that should be used, eg:

function discrete_flux_bc(i, j, grid, clock, fields, parameters)
    # compute gradient and diffusivity
    non_zero_flux = - diffusivity * gradient
    return ifelse(condition_satisified, non_zero_flux, zero(grid))
end