User Forum of Software BASEMENT

nicolagroff

User

Registered: 2017-08-07

Posts: 10

"Exception" problem with sediment transport

Hi!

I'm running a simulation with sediment transport with a gate and a weir condition.
I got this error:

EXCEPTION -> layer top and layer bottom are identical (down_hyd: 1202.63; top_hyd:
   1202.63)...this is not allowed., in file: ..\src\BASEplane\BPSedInput.cpp on
   line 467, Version: 2.8 R5771, compiled at: May 31 2018, 15:28:00

What's wrong??

This is the morfology block:

MORPHOLOGY {
			PARAMETER {
				porosity                       = 35
				density                        = 2650
				control_volume_type            = constant
				control_volume_thickness       = (0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1)
				control_volume_thickness_index = (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 25 26 30 99)
			}
			BEDMATERIAL {
				GRAIN_CLASS {
					diameters = (61)
				}
				MIXTURE {
					name            = M1
					volume_fraction = (100)
				}
				SOIL_DEF {
					name = S0
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S1
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S2
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S3
					LAYER {
						mixture          = M1
						bottom_elevation = -0.5
					}
				}
				SOIL_DEF {
					name = S4
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S5
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S6
					LAYER {
						mixture          = M1
						bottom_elevation = -0.5
					}
				}
				SOIL_DEF {
					name = S7
				}
				SOIL_DEF {
					name = S8
					LAYER {
						mixture          = M1
						bottom_elevation = -0.5
					}
				}
				SOIL_DEF {
					name = S9
					LAYER {
						mixture          = M1
						bottom_elevation = 0
					}
				}
				SOIL_DEF {
					name = S10
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S11
				}
				SOIL_DEF {
					name = S12
				}
				SOIL_DEF {
					name = S13
					LAYER {
						mixture          = M1
						bottom_elevation = 0
					}
				}
				SOIL_DEF {
					name = S14
					LAYER {
						mixture          = M1
						bottom_elevation = 0
					}
				}
				SOIL_DEF {
					name = S15
				}
				SOIL_DEF {
					name = S16
					LAYER {
						mixture          = M1
						bottom_elevation = -1
					}
				}
				SOIL_DEF {
					name = S25
					LAYER {
						mixture          = M1
						bottom_elevation = 0
					}
				}
				SOIL_DEF {
					name = S26
					LAYER {
						mixture          = M1
						bottom_elevation = 0
					}
				}
				SOIL_DEF {
					name = S30
					LAYER {
						mixture          = M1
						bottom_elevation = 0
					}
				}
				SOIL_DEF {
					name = S99
				}
				SOIL_ASSIGNMENT {
					index = (0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 25 26 30 99)
					soil  = (S0 S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 S13 S14 S15 S16 S25 S26 S30 S99)
					type  = index_table
				}
			}
			INITIAL {
				type = initial_mesh
			}
			BEDLOAD {
				FORMULA {
					bedload_formula      = smartjaeggi
					d30                  = 22
					d90                  = 293
					bedload_factor       = 0.25
					theta_critical       = (0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05 0.05)
					theta_critical_index = (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 25 26 30 99)
				}
				BOUNDARY {
					name              = SEZ_IN
					type              = transport_capacity
					string_name       = SEZ_IN
					fraction_boundary = 1
					mixture           = M1
					factor            = 0.25
				}
				BOUNDARY {
					name              = SEZ_OUT
					type              = IODown
					string_name       = SEZ_OUT
					fraction_boundary = 0.9
				}
				PARAMETER {
				}
				DIRECTION {
					lateral_transport_type = lateral_bed_slope
					lateral_index          = (5 6 7 9)
				}
				INNER_BOUNDARY {
					type         = open
					string_name1 = PRESA_M
					string_name2 = PRESA_V
				}
			}
			GRAVITATIONAL_TRANSPORT {
				index                      = (0 1 2 3 4 5 6 7 8 12 13 14 15 16 24 25 26)
				angle_failure_dry          = (40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40 40)
				angle_failure_wetted       = (30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30 30)
				angle_failure_deposited    = (25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25)
				gravity_transport_on_cells = partially wetted
			}
		}
	}
}

Thank you!

Nicola

Lukas Vonwiller

Developer

Registered: 2014-09-04

Posts: 130

Re: "Exception" problem with sediment transport

Hi Nicola

Does it work if you remove the Layers where you have a bottom_elevation = 0?

e.g. instead of

                SOIL_DEF {
                    name = S9
                    LAYER {
                        mixture          = M1
                        bottom_elevation = 0
                    }
                }

try

                SOIL_DEF {
                    name = S9
                }

This results in a fixed bed for soil S9.
Hope this helps!
cheers, Lukas

nicolagroff

User

Registered: 2017-08-07

Posts: 10

Re: "Exception" problem with sediment transport

Hi Lukas,

It works!

Thank you very much!

Have a nice day

Nicola