Options for invert
#!/bin/bash
#
# 2d/3d inversion tool invert
# (description of options)
# 
# OPTIONNAME=value # description (possible values)
#
#
# Global settings
#
DATAFILE=projectfile # defines the project filename (required)
DIMENSION=3 # defines the dimension of the problem (2 for 2d or 3 for 3d)
TOPOGRAPHY=0 # defines if topography is present (0 or 1)
TOPOPOINTS= # defines file which has additional coordinates for topography (x y z)
TOPOPOLY= # defines file which has additional polygons for topography (x0 y0 z0)
TIMESTEPS= # defines file which has the names of additional datafiles used in timestep inversion
PARAGEOMETRY= # defines polygon file for parametric geometry with extern script or program
REFRAKTOR= # defines file with information about refractor
CYLINDER=0 # defines cylindric geometry (0 or 1)
ELECTRODENODES=1 # defines that electrodes are represented as nodes
SPACECONFIG=0 # defines geometric factor (0 for half space, 1 for full space, 2 for mirrored half space)
UNDERWATER=0 # defines underwater survey, sets SPACECONFIG=2
#
# Data settings
#
INPUTOHM=0 # input data is in Ohm
OVERRIDEERROR=0 # overrides given errors with INPUTERRLEVEL and INPUTERRVOLTAGE (0 or 1)
INPUTERRLEVEL=3 # sets input error level (in percent) if no error given
INPUTERRVOLTAGE=100e6 # sets input voltage error (V) if no error given
FILTERVALS=0 # defines filter values (experimental)
#
# Inversion settings
#
RHOSTART=0 # resistivity of start model (0 = median of the data)
# can also be model specification file
MAXITER=20 # maximum number of iteration steps
SINGVALUE=1 # potential value at electrodes, for sensitivity (internal)
LAMBDA=20 # regularization parameter
LAMBDAOPT=0 # optimze lambda by using lcurve (0 or 1)
LAMBDADECREASE=1 # decrease lambda with each iteration
CONSTRAINT=1 # order 0, 1, 2 (experimental)
ZPOWER=0 # weight for vertical gradients
BLOCKYMODEL=0 # iteratively reweighted model (L1)
ROBUSTDATA=0 # iteratively reweighted data (L1)
LOWERBOUND=0.0 # lower resistivity bound (logarithmic barrier)
UPPERBOUND=0.0 # upper resistivity bound (0.0 = deactivated)
SENSMATUPDATE=0 # update sensitivity matrix with each iteration step (internal)
SENSMATDROPTOL=0 # only for very large problems (internal)
SENSMATMAXMEM=800 # only for very large problems (internal)
#
# Mesh settings
#
PARAMAXCELLSIZE=0 # maximum cell size volume (m3) (DIMENSION=3); area (m2) (DIMENSION=2) for para mesh
PRIMMAXCELLSIZE=0 # maximum cell size volume (m3) (DIMENSION=3); area (m2) (DIMENSION=2) for prim mesh
PARADEPTH=0 # maximum depth of parameter domain in meter (0 = automatic)
PARABOUNDARY=5 # boundary around electrodes in parameter domain (percent)
SPLINEBOUNDARY=0 # spline circle boundary instead of piecewise linear interpolation (experimental)
EQUIDISTBOUNDARY=0 # equidistant refined space between electrodes
BOUNDARY=500 # size of boundary area around parameter domain
MESHGEN=tetgen # 3d mesh generator
TETGENTOLERANCE=1e12 # tetgen tolerance limit
TETGENPRESERVEBOUNDARY=0 # tetgen should suppress splitting of boundary facets or segments
PARADX=0.0 # refinement for para grid (values 0.5 will be forced to 0.5)
PRIMDX=0.1 # refinement for prim grid (for 2d relative electrode spacing)
PARA2DQUALITY=33.0 # parameter grid (from 20 (bad) to 33.4 (good))
PRIM2DQUALITY=33.4 # primary grid (from 20 (bad) to 33.4 (good))
PARA3DQUALITY=1.5 # parameter grid (from 1.1 (good) to 2 (bad))
PRIM3DQUALITY=1.2 # primary grid (from 1.1 (good) to 2 (bad))
SURFACEQUALITY=30 # quality of topographical surface grid (from 20 (bad) to 33.4 (good))
SURFACEMAXTRISIZE=0.0 # maximal triangle area of paramatric surface grid
SURFACESMOOTH=0 # improve quality of topographical surface grid
ICDROPTOL=0.0 # if number of nodes 200k drop tolerance is set for ICCG solver
LINSOLVER=1 # sets linear solver, 0 for homebrew PCG (very slow), 1 for TAUCS (direct, ICCG), 2 for LDL (direct with AMD)
SECMESHREFINE=1 # use with caution (experimental)
SECP2MESH=0 # quadratic shapefunction for secmesh (experimental)
PRIMP2MESH=0 # use primary p2 mesh (experimental)
#
# Directory settings
#
MESHBASENAME=mesh # basename for mesh files
DIRMESHS=mesh # directory name for mesh files
DIRPOT=primaryPot # directory name for primary and interpolated potentials
DIRPRIMPOT=potentials # subdirectory name for primary potentials
DIRINTERPOLPOT=interpolated # subdirectory name for interpolated potentials
DIRFEM=femM # directory name for finite element solutions
DIRSENS=sensM # directory name for sensitivity matrix
OLDPRIMMESHSTYLE=0 # for internal use only
