def lp_reservoir3():
    """
    Lake number 4185105 from Pocono test dataset
    an lp compute kernel to put under test
    """
    #print("lp3 test")

    water_elevation = 496.25399574057275
    lake_area = 1.95502996444702
    weir_elevation = 496.19599609375
    weir_coefficient = 0.4
    weir_length = 10.0
    dam_length = 10.0
    orifice_elevation = 495.210001627604
    orifice_coefficient = 0.1
    orifice_area = 1.0
    max_depth = 496.369995117188
    lake_number = 4185105
    initial_fractional_depth = 0.9
    args = [lake_area, max_depth, orifice_area,
            orifice_coefficient, orifice_elevation,
            weir_coefficient, weir_elevation, weir_length,
            initial_fractional_depth, 0.0, water_elevation]

    upstream_ids = array('l')
    k = MC_Levelpool(0, lake_number, upstream_ids, args)
    yield k
def lp_reservoir2():
    """
    an lp compute kernel to put under test
    """
    #print("lp2 test")

    water_elevation = 9.70
    lake_area = 15.0949
    weir_elevation = 9.626
    weir_coefficient = 0.4
    weir_length = 10.0
    dam_length = 10.0
    orifice_elevation = 7.733
    orifice_coefficient = 0.1
    orifice_area = 1.0
    max_depth = 9.96
    lake_number = 16944277
    initial_fractional_depth = 0.9
    args = [lake_area, max_depth, orifice_area,
            orifice_coefficient, orifice_elevation,
            weir_coefficient, weir_elevation, weir_length,
            initial_fractional_depth, 0.0, water_elevation]

    upstream_ids = array('l')
    k = MC_Levelpool(0, lake_number, upstream_ids, args)
    yield k