...
| Code Block | ||||
|---|---|---|---|---|
| ||||
indexBestFit = -1
bcBestFit = 1e100
roughnessBestFit = 1e100
devBestFit = 1e100
for i in range(1,len(timeSeriesToCompare)-1): # case 0 is data so don't use it.
case = listCalibration[i-1]
if case[2] < devBestFit:
indexBestFit = i
bcBestFit = case[0]
roughnessBestFit = case[1]
devBestFit = case[2] |
...
Once the best fit has been found, you can print the corresponding values to the Messages console.
| Code Block | ||||
|---|---|---|---|---|
| ||||
print "selected bc: " + str(bcBestFit) + "m3/s"
print "selected roughness " + str(roughnessBestFit) + "m^(1/3) s^-1" |