The steady state flow can also be called up directly without the graphical user interface via the command line, provided that a batch file *.bsi for the SITRA module or a batch file *.bge for the GEONEU module is present in the directory. Navigate the command line to directory in which the calculation is to be performed. Enter sitra and confirm with the Enter key, this will start a calculation with the SITRA module. The default batch file name for the SITRA module 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).
Entering geoneu and confirming with the Enter key starts a calculation with the GEONEU module. The default batch file name for the GEONEU module is geoneu.bge. This file is searched for in the directory. However, a different batch file can also be specified with the call geoneu filename (the extension ".bge" is automatically appended if required).
The batch files can also be created or modified manually using any editor.
Example of a batch file for a steady-state flow computation with the module GEONEU
out.g # Output file
1 1.0e-8 1.0e-6 # Equation solver (1-iterative,0-direct), residual, iter. difference
1 0 0 # Control lines? Velocity (0 average, 1 right, -1 left), vertical flow rate? (1-yes,0-no)
0 # Calculation of velocities? (1-yes, 0-no)
0 # Import K-values from EICHEN? (1-yes, 0-no)
1 5 1 0.5 # sat/unsat? (1/-1/0)(yes/contin/no),numb.,init.potentials EICH, atten.
0 # Calculation of pathlines ? (1-yes, 0-no)
Example of a batch file for a steady-state flow computation with the module SITRA
out.s # Output file
1 # detailed protocol (1-yes/0-no)
1 0 # Equation solver P/U (1-iterative, 0-direct, -1=multigrid)
1 0 # Calculation of control lines (1-yes, 0-no), IMITTEL
0 0 # Peclet numbers, hydr. cond. from EICHEN (1-yes, 0-no)
STROEMUNG #
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 program ignores the characters behind the sign `#'. This space can be used for individual remarks.
The presented batch file for the steady state flow can be downloaded on our Homepage.
Here you find a detailed description of the batch file for a steady state flow computation with GEONEU
1st line: name of the output file (starting at the first position!)
2nd line: flag for the solver: 0=direct, 1=iterative, -1=multigrid
3rd line: three flags: 1st flag for calculation of control lines (yes=1), 2nd flag for velocities for calculation of control lines (0=mean velocity, 1=velocity on the right hand side of the control line, -1=velocity on the left hand side of the control line), 3rd flag for calculation of vertical mass flow rates (3D), (only for 3D models without ceased elements! )
4th line: flag for calculation of velocities (yes=1), (automatically in the case of calculating control lines, mass flow rates or path lines)
5th line: flag for using calibrated hydraulic conductivity (yes=1), (only if calibration was running before!)
6th line: flag and parameter for iteration of thickness / free surface:
Three-dimensional and vertical models:
flag for iteration of the free surface(1 or -1 = yes/0=no), number of iteration (if flaf = 1 or -1), flag for initial values from EICH (flag=1), damping factor for iterations
Horizontal models with UNTE:
Flag for iteration of thickness (1 or -1=yes, 0=no), number of iterations, flag for protocol of iter. Thickness (1=yes/0=no), flag for protocol of iter. Potentials (1=yes/0=no), attenuation factor for iteration
With flag=-1 the iteration is performed, too. If possible, a former iteration is continued! For this case the use of calibrated hydraulic conductivities from EICHEN is not possible anymore!)
7th line: flag for pathline calculation (1=yes/0=no)
from the 8th line: input parameter for the path line calculation:
8th line: maximum length (m) of a path line
9th line: maximum flowing time (unit of time) of a path line
10th line: path lines equal (1) or against (0) the flow direction
11th line: porosity
12th line: snapping radius near wells (between 0.0 and 1.0)
13th line: 3 flags for number of pathlines with starting points in: in the centre of an element, at nodes, at any location:
From 14th line:
Start in the centres of the elements: element numbers (maximum 10 per line)
Start in nodes: first line: number of starting points per node, radius (m), next lines: node numbers (maximum 10 per line)
Start at any location: line with x, y (and z) coordinates (3D model) for each starting point
Here you find a detailed description of the batch file for a flow computation with SITRA
1st line: name of the output file (starting at the 1st column!)
2nd line: flag for complete record (flag=1)
3rd line: flag for the solver: 0=direct, 1=iterative
4th line: flag for computing control lines (1=yes/0=no), if flag = 1:
Flag for computing the velocities for the control lines:
0 = average between left and right velocity
1 = velocity from the right of the control line
-1 = velocity from the left of the control line
5th line: flag for the computing of Peclet numbers and (transient calculation only) Courant numbers (1=yes/0=no), flag for using the calibrated hydraulic conductivity from calibration (1=yes/0=no)
6th line: 'STROEMUNG' (only flow modelling) or 'TRANSPORT' (flow and transport modelling)
7th line: 3D model or vertical model: flag for a saturated/unsaturated calculation (1=yes/0=no)
Horizontal model using UNTE: flag for iteration of thickness (1=yes/0=no), else irrelevant line
8th line: 1. flag for transport: considering equal concentration at equal potential head (1=yes/0=no), 2. flag only important in 3D model or vertical model: density dependent calculation (1=yes/0=no), 3. flag: considering pressure or potential equation (1=pressure/0=potential)
9th line: Saturated/unsaturated calculation or iteration of thickness and/or density dependent calculation: number of iterations, attenuation factor, else irrelevant
Remark: For transient calculations with an iteration of thickness the number of iterations per time step is always 1 (no iteration!). Then, the attenuation factor should be equal 1.0 (thickness is only updated).
10th line: flag for steady state or transient flow and flag for steady state or transient transport: 3D model or vertical model: 0= steady state, 1= transient,
horizontal model with groundwater filled thickness:
0=steady state, 1=transient (without consideration of confined/unconfined conditions), 2=transient (with consideration of confined/unconfined conditions)
11th line: 2 flags for initial values of potential head and concentration/temperature (potential head for sat./unsat. calculation or transient flow) and concentration/temperature for density dep. calculation or transient transport: 2= data from file null (if file null not found then initial value = EICH/AKON), 1= initial potential head EICH or concentration/temperature AKON (if no values for EICH or AKON found then initial value = 0), 0= start value = 0.0
from the 12th line: input parameter for transient calculation
12th line: name of the transient input file (starting in 1st column!)
13th line: flag for warm start using the data in the file null (1=yes/0=no), (if file null not found then input is ignored!)
14th line: flag for definition of transient time steps:
0= flag for using the time steps from the transient input file, further input: factor for scaling down the time steps (integer), add. factor for scaling down time steps in transport eq. (optional)
1= constant time step length, further input (separated by blank or tab): time step length (double), unit of time (integer), (0=sec. 1=min. 2=hour. 3=day 4=month 5=year), add. factor for scaling down time steps in transport eq. (optional)
15th line: interpolation flags for transient boundary conditions for POTE, KNOT, FLAE, KONZ, 1KON,VORF and EFLA
For POTE, KNOT, FLAE, KONZ, 1KON, VORF and EFLA: 0= implicit use without interpolation, 1= implicit use with interpolation
For FLAE, EFLA and KNOT: -1= explicit use without interpolation, -2= explicit use with interpolation, flag for calculation with transient terrain settings as a result of subsurface mining (1=yes/0=no)
16th line: number of time steps or input for flooding: GRUB POTE or MENGE
17th line: flag for a complete record of the final results (file out66) to continue the transient calculation (flag=1)
18th line: flag for saving intermediate results to binary files: 0= flag for saving the results at every n-th time step, further input: n, 1= flag for saving the results at n defined time steps, further input: n, 2= flag for saving complete time-series plots at single nodes, further input: n
if Flag=1 in the 18th line:
from the 19th line: time steps in increasing order, maximum 10 per line, without the starting (0) and the last time step!
if Flag=2 in the 18th line:
from the 19th line: node numbers, maximum 10 per line
Stand: 11.06.2025 | 