open TELEMAC-MASCARET The mathematically superior suite of solvers

I am conducting a dam break study wherein I have modelled a reservoir and would like to route the water in the reservoir at a particular elevation downstream similar to the Malpasset example provided. Also I am coupled this model with SISYPHE. My problem is I am doing a simulation where the domain has a dry bed at the downstream of the dam so when i assign the "INITIAL CONDITION = ZERO DEPTH" the entire domain including the reservoir area where I Have assigned a water depth via blue kenue is also computed with a zero depth.
So instead i tried using the "INITIAL CONDITION = PARTICULAR" where we need to modify the subroutine condin.f. My problem is I only know how to input either an x coordinate value or y coordinate value in the condin file but I actually need to input a number of x & y coordinate value in condin.f subroutine.
As of now I did a simulation where i created an initial condition file with the reservoir water level reaching the required depth and used it as a cotinuation file by enabling COMPUTATION CONTINUED = YES and using the initial simulation as the PREVIOUS COMPUTATION FILE but the more accurate simulation would be with INITIAL CONDITION = particular.
I have read a lot of posts but haven't seen any post which describes how to modify the subroutine condin.f with multiple x and y coordinates and varying H values.
I have attached by files along with a result file that i have run using the Particular case modifying only the x coordinate.
Please help.

Hi
using particular initial condition or a manually made previous continuation file is exactly the same.
The accuracy of the simulation doesn't comes from the way you define initial condition!
If you really want to use the fortran, there is a subroutine called inpoly.f (in utils/bief) which allow you to select the nodes inside a polygon...

Hope this helps
Regards

Thank you for the reply sir

The problem when i create an initial condition file is not being able to represent the reservoir storage and elevation accurately.
While trying to create an initial condition file with the required initial elevation I had to edit the mesh at the dam location to the required elevation and then run the simulation but the reservoir storage as well as the flowrate at the dam location were underestimated which effects the downstream flow in the main simulation.
While creating an initial condition satisfying the required reservoir volume I am getting an underestimated initial elevation.

Regards
Sruthi

Hi
Not sure to well understand your problem...
Initial selafin file with the right free surface elevation --> storage and flowrate underestimated
Initial selafin file with right volume --> initial elevation underestimated

How this could be possible?

Greetings Sir,
In the first case with initial selafin file with required elevation I edited the interpolated mesh so that the elevation along the dam line is at the required level, ran the simulation and used this .slf file as a previous computation file for the new geometry file where i redited the mesh values back to the original values. the main file thus has a new selafin file with the same mesh nodes but different bottom topography values. The simulation runs but the results are not correct.

In the second case I initiated the simulation with a high PRESCRIBED FLOWRATE and PRESCRIBED ELEVATION to attain the reservoir required volume along with taking care of the fact that the water stays within the reservoir. But here the elevation at the dam location is not the required one.

Is there any other way i can represent the reservoir volume as well as the reservoir elevation at the the required levels at the same time.

Regards
Sruthi

Hi

Ok, I think your approach is not the best one

In your mesh, you have the real bottom of your area.
As it's a dam break study, you should consider there is no dam existing (total collapse of the dam)

In my opinion, working in BK, you add a new variable in your mesh named FREE SURFACE.
You could also duplicate your mesh to create a new initial condition file.
Then in the free surface variable, in the area of the reservoir, you set the value of the free surface variable to the reservoir elevation (before break)
And that's all!

Regards