The inverse calibration can also be called up directly without the SPRING graphical user interface via the command line, provided that suitable batch files *.bsi and *.biv are available in the directory. To do this, navigate to the directory in which the calculation is to be carried out. Enter sitra into the command line and confirm with the Enter key to start an inverse calibration. The default batch file name is sitra.bsi. This file is searched for in the directory. However, another batch file can also be specified with the sitra file name call (the extension "*.bsi" is automatically appended if required). The batch files can be created or modified manually using any editor.
For inverse calibration, the sitra.bsi file is first required, in which the batch file for the inverse calculation run (in the example: inv.biv) is specified.
Example of a batch file sitra.bsi (steady state calculation):
out.s # Output file
1 # detailed protocol (1-yes/0-no)
1 0 # Equation solver P/U (1-iterative, 0-direct)
1 0 # Calculation of control lines (1-yes, 0-no), IMITTEL
0 0 # Peclet numbers, hydr. cond. from EICHEN (1-yes, 0-no)
STROEMUNG+INVERS inv.biv# inverse modelling without transport
1 # satur./unsat. computation (1-yes,0-no)
0 0 0 # equal conc., density-related calc., p or h (1-yes/0-no)
5 0.5 # Number of iteration steps, factor of attenuation
0 0 # Flow/transport steady st./0 trans./1 trans.+conf./unconf./2
0 0 # Initial values POTE/KONZ from ASAT(3) null(2) EICH/AKON(1) =0(0)
The batch file inv.biv contains the parameters for the inverse calibration.
Example of a batch file inv.biv (inverse calculation), called from the steady state calculation:
out.inv # output file for inverse modelling record
invers.obs # observation data file
invers.csv # parameter file
1 # optimization with Levenberg-Marquard / 1: Quasi-Newton
10 # maximum number of iterations
KWER 1 0 # interpolation flags(if necessary) of KWER parameters / parameter d for distance weighting
SPEI 10 3 # interpolation (if necessary) of KSPE/parameter d for distance weighting
output = SENS;KOVAR
The value for the interpolation flag has the following meaning:
0: Kriging
1: Gaussian
10: distance weighting – Sampson
11: distance weighting 1/d
12: distance weighting 1/d^2
13: distance weighting 1/d^4