The density-dependent mass transport calculation can also be called up directly via the command line without the graphical user interface, provided that a suitable batch file *.bsi and sitr file for the SITRA module is available in the directory. When entering the command line, navigate to the directory in which the calculation is to be carried out. Entering sitra and confirming with the Enter key starts 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, a different batch file can also be specified with the sitra file name call (the extension ".bsi" is automatically appended if required). The batch file can be created or modified manually using any editor
Example of a batch file of the density-dependent mass transport calculation
out.s # Output file
0 # 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 # Courant numbers, hydr. cond. from EICHEN (1-yes, 0-no)
TRANSPORT #
1 # satur./unsat. computation (1-yes,0-no)
0 1 0 # equal conc., density-related calc., p or h (1-yes/0-no)
3 0.5 # Number of iteration steps, factor of attenuation
1 1 # Flow/transport steady st./0 trans./1 trans.+conf./unconf./2
1 0 # Initial values POTE/KONZ from ASAT(3) null(2) EICH/AKON(1) =0(0)
instat.bis # Transient input file
0 # Warm start with null-file (1-yes, 0-no)
1 1 3 5 # time step 0->from file+factor, 1->time step+time step width+time unit
1 -1 -1 0 0 1 -1 1 # transient consideration of POTE,KNOT,FLAE,KONZ,1KON,VORF,EFLA, mining subsidence?
GRUB POTE # Mine with fixed potential
2 # Number of time steps
1 # Saving for optional continuation of the transient computation? (1->yes, 0->no)
0 1 # Save intermediary results 0-> every n.th step, 1 -> for n follow. TS, 2-> hydrographs for n nodes
The program ignores the characters behind the sign `#'. This space can be used for individual remarks.
The presented batch file for the model checking can be downloaded on our Homepage.
Here you find a detailed description of the batch file for a density-dependent mass transport calculation
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, -1=multigrid
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.03.2025 | 