Hello riadh, thank you for your suggestions.
I changed the initial water depth as you told and the program runs for a brief period, but after some iterations it gives me the previous error. Besides, it correctly runs when the number of time steps and the graphic printout period are low.
Anyway, here are the parameters of my cas.file:
BOTTOM SMOOTHINGS =1
FRICTION COEFFICIENT =0.06
LAW OF BOTTOM FRICTION =4
TURBULENCE MODEL =3
/ EQUATIONS, BOUNDARY CONDITIONS
VELOCITY PROFILES =1;1
OPTION FOR LIQUID BOUNDARIES =1;1
/ EQUATIONS, INITIAL CONDITIONS
INITIAL DEPTH =0.1
INITIAL CONDITIONS ='CONSTANT DEPTH'
/ INPUT-OUTPUT, FILES
SECTIONS INPUT FILE ='control_section.txt'
LIQUID BOUNDARIES FILE ='hydro.liq'
GEOMETRY FILE ='geometry.slf'
SECTIONS OUTPUT FILE ='output3.txt'
STEERING FILE ='bisagnoSS.cas'
BOUNDARY CONDITIONS FILE ='bc_bisagno.cli'
RESULTS FILE ='results3.slf'
/ INPUT-OUTPUT, GRAPHICS AND LISTING
VARIABLES FOR GRAPHIC PRINTOUTS =U,V,B,H,F,S
MASS-BALANCE =true
GRAPHIC PRINTOUT PERIOD =50
/ NUMERICAL PARAMETERS
CONTINUITY CORRECTION =true
DEBUGGER =1
NUMBER OF TIME STEPS =300
TREATMENT OF THE LINEAR SYSTEM =2
SUPG OPTION =1;1;1;1
COMPATIBLE COMPUTATION OF FLUXES =true
TIME STEP =5.
INITIAL TIME SET TO ZERO =true
/NUMERICAL PARAMETERS, SOLVER
SOLVER =1
SOLVER ACCURACY =1.E-3
And this is my liquid boundaries file:
#
T Q(2) SL(1)
s m3/s m
0 20.0 3.0
5000 20.0 3.0
10000 20.0 3.0
15000 100.0 3.5
20000 100.0 3.5
25000 20.0 3.0
30000 20.0 3.0
35000 50.0 3.3
40000 20.0 3.0
45000 20.0 3.0
50000 20.0 3.0
#
Finally, is the stage discharge curve file an alternative to the column SL?
Thank you again for you attention :)