parouby
/
Pamhyr2
mirror of https://gitlab.com/pamhyr/pamhyr2
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

debug read rubar3 results

scenarios
Theophile Terraz 2025-10-22 15:01:10 +02:00
parent 08c414a66e
commit 861b78b928
1 changed files with 13 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
src/Solver/RubarBE.py
View File

@ -583,27 +583,30 @@ class Rubar3(CommandLineSolver):
for reach, lm in reachs: for reach, lm in reachs:
# First profile # First profile
e, s, q, z = read_data_line(f) h, s, q, z = read_data_line(f)
set_and_compute_limites(reach, 0, z, q, s) set_and_compute_limites(reach, 0, z+h, q, s)
# For each profile # For each profile
ind = 1 ind = 1
pe, ps, pq, pz = read_data_line(f) ph, ps, pq, pz = read_data_line(f)
while ind < lm - 2: while ind < lm - 2:
e, s, q, z = read_data_line(f) h, s, q, z = read_data_line(f)
set_and_compute_limites( set_and_compute_limites(
reach, ind, reach, ind,
(pz + z) / 2, (ph + pz + ph + h) / 2,
(pq + q) / 2, (pz + s) / 2,
(ps + s) / 2 (pq + q) / 2
) )
pe, ps, pq, pz = z, s, q, e ph = h
ps = s
pq = q
pz = z
ind += 1 ind += 1
# Last profile # Last profile
e, s, q, z = read_data_line(f) h, s, q, z = read_data_line(f)
set_and_compute_limites(reach, ind, z, q, s) set_and_compute_limites(reach, ind, z+h, q, s)
@timer @timer
def results(self, study, repertory, qlog=None, name="0"): def results(self, study, repertory, qlog=None, name="0"):