Exemple #1
0
advectiveFluxBoundaryConditions = {0: getAFBC_h, 1: getAFBC_u, 2: getAFBC_v}


def getDFBC_u(x, flag):
    if flag == 0:
        return lambda x, t: 0.0
    else:
        return None


def getDFBC_v(x, flag):
    if flag == 0:
        return lambda x, t: 0.0
    else:
        return None


diffusiveFluxBoundaryConditions = {0: {}, 1: {1: getDFBC_u}, 2: {2: getDFBC_v}}

#########################################
##### CREATE MODEL AND COEFFICIENTS #####
#########################################
bathymetry = {0: bathymetry_function}
LevelModelType = SW2DCV.LevelModel
coefficients = SW2DCV.Coefficients(g=g,
                                   bathymetry=bathymetry,
                                   cE=cE,
                                   LUMPED_MASS_MATRIX=LUMPED_MASS_MATRIX,
                                   LINEAR_FRICTION=LINEAR_FRICTION,
                                   mannings=mannings)
Exemple #2
0
advectiveFluxBoundaryConditions = {0: getAFBC_h, 1: getAFBC_u, 2: getAFBC_v}


def getDFBC_u(x, flag):
    if flag == 0:
        return lambda x, t: 0.0
    else:
        return None


def getDFBC_v(x, flag):
    if flag == 0:
        return lambda x, t: 0.0
    else:
        return None


diffusiveFluxBoundaryConditions = {0: {}, 1: {1: getDFBC_u}, 2: {2: getDFBC_v}}

#########################################
##### CREATE MODEL AND COEFFICIENTS #####
#########################################
bathymetry = {0: bathymetry_function}
LevelModelType = SW2DCV.LevelModel
coefficients = SW2DCV.Coefficients(g=g,
                                   bathymetry=bathymetry,
                                   cE=cE,
                                   LUMPED_MASS_MATRIX=LUMPED_MASS_MATRIX,
                                   mannings=mannings)
Exemple #3
0
g = physical_parameters['gravity']
LINEAR_FRICTION = physical_parameters['LINEAR_FRICTION']
mannings = physical_parameters['mannings']

# NUMERICAL PARAMETERS #
cE = numerical_parameters['cE']
LUMPED_MASS_MATRIX = numerical_parameters['LUMPED_MASS_MATRIX']

# ********************************** #
# ********** COEFFICIENTS ********** #
# ********************************** #
LevelModelType = SW2DCV.LevelModel
coefficients = SW2DCV.Coefficients(
    g=g,
    bathymetry={0: bathymetry} if bathymetry is not None else None,
    cE=cE,
    LUMPED_MASS_MATRIX=LUMPED_MASS_MATRIX,
    LINEAR_FRICTION=LINEAR_FRICTION,
    mannings=mannings)

# **************************************** #
# ********** INITIAL CONDITIONS ********** #
# **************************************** #
initialConditions = {
    0: initialConditions['water_height'],
    1: initialConditions['x_mom'],
    2: initialConditions['y_mom']
}

# ***************************************** #
# ********** BOUNDARY CONDITIONS ********** #