#!/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=100e-6 # 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 l-curve (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=1e-12 # 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