Input Parameters

The SpecFWAT employs a YAML file to define the parameters for forward simulation, adjoint simulation, post-processing, and optimization. This file is typically named fwat_params.yml and is located in the DATA directory of your SpecFWAT project. Below shows section of the fwat_params.yml file:

NOISE Section

This section defines the noise parameters for the forward simulation and measurements of adjoint sources.

NOISE:
  MESH_PAR_FILE: DATA/meshfem3D_files/Mesh_Par_file # Mesh file
  SCOMPS: ['Z'] # Components of the source
  RCOMPS: ['Z'] # Components of the receiver
  CH_CODE: BX # Channel code
  IMEAS: 7 # Measurement type
  USE_RHO_SCALING: True # Whether to scale the density gradient
  NSTEP: 4500 # Number of time steps
  DT: 0.06 # Time step for the noise data
  SHORT_P: [6, 10, 20] # Short period of filters
  LONG_P: [15, 25, 40] # Long period of filters
  GROUPVEL_MIN: [2.3, 2.3, 2.5] # Approximate minimum group velocity
  GROUPVEL_MAX: [3.2, 3.5, 4.0] # Approximate maximum group velocity
  ADJ_SRC_NORM: False # Set the following to .true. to normalize adjoint sources across different bands
  USE_NEAR_OFFSET: False # Set the following to false if use only data > 1 average wavelength
  SUPPRESS_EGF: False # Set to .false. when the data are cross-correlation functions
  PRECOND_TYPE: 1 # 1: inner product of acceleration
  SIGMA_H: 5000
  SIGMA_V: 5000

Mesh Parameters

• MESH_PAR_FILE: Path to the mesh parameter file with the same format as meshfem3D. SpecFWAT uses internal mesh generator of Specfem3D to create the mesh based on the model parameters, but allow users to specify the mesh file rathar than a fixed mesh file DATA/meshfem3D_files/Mesh_Par_file of Specfem3D.

Solver Parameters

• RCOMPS: Components of the receiver, For tomography of isotropic media, it is usually set to ['Z'] for vertical component.

• SCOMPS: Components of the source used in ambient noise correlation. For isotropic inversions it is typically set to ['Z'].

• CH_CODE: Channel code, which is used to identify the channel in the data. It is usually set to BX for broadband data.

• NSTEP: Number of time steps for the forward/adjoint simulation. It overrides the NSTEP in DATA/Par_file, which allow users to specify the number of time steps for different data types.

• DT: Time step for the noise data, It overrides the DT in DATA/Par_file, which allow users to specify the time step for different data types.

Adjoint source

• IMEAS: Option of adjoint source measurements. It is used to determine the objective function. See measure_adj manual for more details of adjoint source measurements. The default value is 7, which corresponds to the cross-correlation time-shift adjoint source.

• SHORT_P: Short cut-off period of filters for the adjoint source measurements. It is a list of values, which are used to filter the data in different frequency bands.

• LONG_P: Long cut-off period of filters for the adjoint source measurements. It is a list of values with the same length as SHORT_P.

• GROUPVEL_MIN: Approximate minimum group velocity to determine the time window for the adjoint source measurements. It is a list of values with the same length as SHORT_P.

• GROUPVEL_MAX: Approximate maximum group velocity to determine the time window for the adjoint source measurements. It is a list of values with the same length as SHORT_P.

• ADJ_SRC_NORM: Set to True to normalize the adjoint sources across different bands.

• USE_NEAR_OFFSET: Set to False if you want to use only data with distance greater than 0.5 * wavelength in each period band. It is usually set to False in practice.

• SUPPRESS_EGF: Whether calculate difference to obtain empirical Green’s functions (EGF). Set to False if the data are cross-correlation functions.

Post-processing parameters of noise data

• PRECOND_TYPE: Type of preconditioning for gradient. 0 ($P_0$) or 1 ($P_1$) are accepted value based on inner product of acceleration, The default value is 1.

• SIGMA_H: Horizontal smoothing length in meter of gradient.

• SIGMA_V: Vertical smoothing length in meter of gradient.

• USE_RHO_SCALING: Whether to scale the density ($\rho$) gradient using the velocity gradient during the model update. The default value is True.

TELE Section

This section defines the parameters for the teleseismic FWAT.

TELE:
  TELE_TYPE: 1 # 1: teleseismic data, 2: receiver function, 3: teleseismic Cross-convolution
  MESH_PAR_FILE: DATA/meshfem3D_files/Mesh_Par_file # Mesh file
  RCOMPS: ['Z', 'R'] # Components of the receiver
  CH_CODE: BX # Channel code
  INJECTION_TYPE: 1 # 1: FK injection; 2: AxiSEM injection
  SAVE_FK: True # Save the FK wavefield (only when INJECTION_TYPE=1)
  COMPRESS_LEVEL: 0 # Compression level of the saved FK wavefield (only when INJECTION_TYPE=1)
  SUPPRESS_STF: True # Whether to suppress source time function convolution
  USE_LOCAL_STF: False # Whether to use a local (per-event) source time function
  NSTEP: 3000 # Number of time steps
  DT: 0.025 # Time step for the teleseismic data
  SHORT_P: [2] # Short period of filters
  LONG_P: [40] # Long period of filters
  TIME_WIN: [-5, 45] # Time window for the teleseismic data
  PRECOND_TYPE: 2 # 2: abs(Z); 3: root squared z
  SIGMA_H: 5000
  SIGMA_V: 5000
  RF:
    F0: [1.0] # Gaussian width for the RF
    MAXIT: 200 # Maximum number of iterations
    MINDERR: 0.001 # Minimum residual error when the RF converges
    TSHIFT: 5.0 # Time shift before P

Adjoint Source

The parameters in the TELE section with the same name as in the NOISE section are similar, but with different values. The following parameters are specific to the TELE section:

• TELE_TYPE: Type of objective function for teleseismic data. It can be

  • 1 for teleseismic waveform difference (Wang et al., 2021 ).
  • 2 for receiver function difference (Xu et al., 2023 ).
  • 3 for teleseismic cross-convolution difference.

• INJECTION_TYPE: Type of wavefield injection method for teleseismic simulation. It can be:

  • 1 (default) for FK (frequency-wavenumber) injection.
  • 2 for AxiSEM injection.

• AXISEM_BASE_PATH: Base path to the AxiSEM wavefield database. Only required when INJECTION_TYPE is 2.

• SAVE_FK: Set to True to save the FK wavefield for teleseismic data. Only applicable when INJECTION_TYPE is 1 (FK injection). The FK wavefield will be checked in the LOCAL_PATH (Set in DATA/Par_file) before running the forward simulation. If the FK wavefield is not found, it will be generated and saved in the LOCAL_PATH/FK_{event_name} directory, otherwise, it will be loaded from the existing FK wavefield directory to reduce time consumption. It will be useful for GPU acceleration.

• COMPRESS_LEVEL: Compression level of the saved FK wavefield in hdf5 format. Only applicable when INJECTION_TYPE is 1 (FK injection). Same as compression_opts in h5py.File. It can be set to 0 (no compression), 1 (fastest), 2, 3, 4, 5, 6, 7, or 9 (best compression). The default value is 0.

• SUPPRESS_STF: Whether to suppress source time function convolution when forward simulating teleseismic data. If set to True (default), the source time function will not be convolved. If set to False, the source time function will be convolved, and STF files named STF_{event_name}.sac must be prepared in the src_rec directory.

• USE_LOCAL_STF: Whether to use a per-event (local) source time function file from the src_rec directory. When set to True, SpecFWAT will look for STF_{event_name}.sac for each event individually. The default is False.

💡

The time 0 in the STF should correspond to the onset time of the earthquake. An example is shown below:

Source: SCARDEC Source Time Functions Database

• TIME_WIN: Time window for the teleseismic data. It is a list of two values, which are the start and end time before and after direct P arrival.

• PRECOND_TYPE: Type of preconditioning for gradient. 2 ($P_2$) or 3 ($P_3$) are accepted value based on Z-precondition. The default value is 2.

Receiver function parameters

For Receiver function adjoint tomography, the synthetic receiver function is calculated based on the iterative deconvolution method. The following parameters are used to control the receiver function inversion:

• RF.F0: Gaussian width for the receiver function. It is a list of values, which are used to filter the data in different frequency bands.

• RF.MAXIT: Maximum number of iterations for the receiver function inversion.

• RF.MINDERR: Minimum residual error when the receiver function converges.

• RF.TSHIFT: Time shift before P arrival for the receiver function inversion. It is used to align the receiver function with the direct P arrival.

LEQ Section

This section defines the parameters for Local Earthquake adjoint tomography (LEQ).

LEQ:
  MESH_PAR_FILE: DATA/meshfem3D_files/Mesh_Par_file # Mesh file
  NSTEP: 3000 # Number of time steps
  DT: 0.05 # Time step
  PRECOND_TYPE: 1 # 1: inner product of acceleration
  IMEAS: 7 # Measurement type (default: IMEAS_CC_TT_MT)
  RCOMPS: ['Z', 'R', 'T'] # Components of the receiver
  CH_CODE: BX # Channel code
  SHORT_P: [2, 5] # Short period of filters
  LONG_P: [10, 20] # Long period of filters
  ADJ_SRC_NORM: False # Normalize adjoint sources across bands
  USE_RHO_SCALING: True # Scale density gradient using velocity gradient
  SIGMA_H: 5000
  SIGMA_V: 5000
  WINDOW:
    WIN_TYPE: 1 # Window selection type (1–4)
    WIN_TYPE_1:
      SLIDING_WIN_LEN_FAC: 2.0
      MIN_VELOCITY: 2.4
      THRESHOLD_CORR: 0.8
      THRESHOLD_SHIFT_FAC: 0.3
      JUMP_FAC: 0.1
      MIN_WIN_LEN_FAC: 1.5
      MIN_SNR_WINDOW: 5.0
      MIN_ENERGY_RATIO: 5.0
      MIN_PEAKS_TROUGHS: 3
      IS_SPLIT_PHASES: False
    WIN_TYPE_2:
      MIN_VELOCITY: 2.4
      MIN_WIN_LEN_FAC: 1.5
      THRESHOLD_CORR: 0.8
      THRESHOLD_SHIFT_FAC: 2.0
      MIN_PEAKS_TROUGHS: 3
      THRESHOLD_DLNA: 1.0
      MIN_SNR_WINDOW: 5.0
      IS_SPLIT_PHASES: False
      MAX_TIME_BEFORE_FIRST_ARRIVAL: 10.0
      STALTA_WATER_LEVEL: 0.08
      C_0: 1.0
      C_2: 0.0
      C_3A: 4.0
      C_3B: 2.5
      C_4A: 2.0
      C_4B: 6.0
    WIN_TYPE_3:
      GROUPVEL_MIN: [2.3, 2.3]
      GROUPVEL_MAX: [3.5, 4.0]
    WIN_TYPE_4:
      PHASE: P
      TIME_WIN: [-5, 45]

Solver Parameters

The MESH_PAR_FILE, NSTEP, DT, RCOMPS, CH_CODE, SHORT_P, LONG_P, ADJ_SRC_NORM, USE_RHO_SCALING, SIGMA_H, and SIGMA_V parameters have the same meaning as in the NOISE section.

• PRECOND_TYPE: Type of preconditioning for gradient. 0 ($P_0$) or 1 ($P_1$) are accepted.

• IMEAS: Measurement type for adjoint sources. The default is IMEAS_CC_TT_MT (cross-correlation time shift via multi-taper).

WINDOW Sub-section

The WINDOW sub-section controls how time windows are selected for LEQ adjoint source measurements. The WIN_TYPE key selects among four strategies:

• WIN_TYPE:
  • 1 (WIN_SELECTOR_TYPE): Automated FLEXWIN-style selector.
  • 2 (WIN_FLEX_TYPE): Flexible window with STALTA-based criteria.
  • 3 (WIN_GROUPVEL_TYPE): Group-velocity based window.
  • 4 (WIN_ARRIVAL_TYPE): Arrival-time based window with a fixed time window around a named phase.

WIN_TYPE = 1 — Automated Selector (WIN_TYPE_1)

WINDOW:
  WIN_TYPE: 1
  WIN_TYPE_1:
    SLIDING_WIN_LEN_FAC: 2.0   # Sliding window length factor (relative to dominant period)
    MIN_VELOCITY: 2.4          # Minimum group velocity (km/s) for initial window
    THRESHOLD_CORR: 0.8        # Minimum cross-correlation coefficient to accept a window
    THRESHOLD_SHIFT_FAC: 0.3   # Maximum allowable time shift factor
    JUMP_FAC: 0.1              # Jump factor for discontinuous windows
    MIN_WIN_LEN_FAC: 1.5       # Minimum window length factor
    MIN_SNR_WINDOW: 5.0        # Minimum signal-to-noise ratio within window
    MIN_ENERGY_RATIO: 5.0      # Minimum energy ratio between signal and noise
    MIN_PEAKS_TROUGHS: 3       # Minimum number of peaks/troughs in window
    IS_SPLIT_PHASES: False     # Whether to split windows at phase arrivals

• SLIDING_WIN_LEN_FAC: Length of the sliding window as a factor of the dominant period. Default 2.0.
• MIN_VELOCITY: Minimum group velocity in km/s used to determine the earliest onset of the time window. Default 2.4.
• THRESHOLD_CORR: Minimum cross-correlation coefficient between observed and synthetic waveforms to accept a window. Default 0.8.
• THRESHOLD_SHIFT_FAC: Maximum allowable time-shift as a fraction of the dominant period. Default 0.3.
• JUMP_FAC: Factor controlling merging of adjacent windows. Default 0.1.
• MIN_WIN_LEN_FAC: Minimum window length as a factor of the dominant period. Default 1.5.
• MIN_SNR_WINDOW: Minimum signal-to-noise ratio required within the window. Default 5.0.
• MIN_ENERGY_RATIO: Minimum ratio of signal energy to noise energy. Default 5.0.
• MIN_PEAKS_TROUGHS: Minimum number of peaks and troughs within the window. Default 3.
• IS_SPLIT_PHASES: Whether to split windows at predicted phase arrivals. Default False.

WIN_TYPE = 2 — Flexible Window (WIN_TYPE_2)

WINDOW:
  WIN_TYPE: 2
  WIN_TYPE_2:
    MIN_VELOCITY: 2.4
    MIN_WIN_LEN_FAC: 1.5
    THRESHOLD_CORR: 0.8
    THRESHOLD_SHIFT_FAC: 2.0
    MIN_PEAKS_TROUGHS: 3
    THRESHOLD_DLNA: 1.0
    MIN_SNR_WINDOW: 5.0
    IS_SPLIT_PHASES: False
    MAX_TIME_BEFORE_FIRST_ARRIVAL: 10.0
    STALTA_WATER_LEVEL: 0.08
    C_0: 1.0
    C_2: 0.0
    C_3A: 4.0
    C_3B: 2.5
    C_4A: 2.0
    C_4B: 6.0

• MIN_VELOCITY, MIN_WIN_LEN_FAC, THRESHOLD_CORR, MIN_PEAKS_TROUGHS, MIN_SNR_WINDOW, IS_SPLIT_PHASES: Same meaning as in WIN_TYPE_1.
• THRESHOLD_SHIFT_FAC: Maximum allowable time-shift factor. Default 2.0 (larger than in the selector type).
• THRESHOLD_DLNA: Maximum allowable logarithmic amplitude difference. Default 1.0.
• MAX_TIME_BEFORE_FIRST_ARRIVAL: Maximum time (seconds) before the first predicted arrival to start the window. Default 10.0.
• STALTA_WATER_LEVEL: Water level for the STA/LTA ratio. Default 0.08.
• C_0, C_2, C_3A, C_3B, C_4A, C_4B: FLEXWIN-style tuning coefficients for the water-level criteria. Defaults: 1.0, 0.0, 4.0, 2.5, 2.0, 6.0.

WIN_TYPE = 3 — Group Velocity Window (WIN_TYPE_3)

WINDOW:
  WIN_TYPE: 3
  WIN_TYPE_3:
    GROUPVEL_MIN: [2.3, 2.3] # Minimum group velocity per filter band
    GROUPVEL_MAX: [3.5, 4.0] # Maximum group velocity per filter band

• GROUPVEL_MIN: List of minimum group velocities (km/s), one per filter band defined by SHORT_P/LONG_P.
• GROUPVEL_MAX: List of maximum group velocities (km/s), one per filter band.

WIN_TYPE = 4 — Arrival-based Window (WIN_TYPE_4)

WINDOW:
  WIN_TYPE: 4
  WIN_TYPE_4:
    PHASE: P          # Phase name for arrival time calculation
    TIME_WIN: [-5, 45] # Time window relative to the phase arrival [start, end] in seconds

• PHASE: Name of the seismic phase to use as the reference arrival (e.g. P, S). Default P.
• TIME_WIN: Two-element list specifying the start and end of the window in seconds relative to the predicted phase arrival.

ADJOINT_SOURCE Section

SpecFWAT employs ForAdjoint to measure adjoint sources of ambient noise and local earthquakes, which provides built-in measurement methods including:

• cross-correlation time shift and amplitude ratio
• waveform difference
• multi-taper phase shift and amplitude ratio
• Cross-convolution waveform difference
• Receiver function difference
• Exponentiated phase misfit

ADJOINT_SOURCE:
  ITAPER_TYPE: 1 # 1: Hanning, 2: Hamming, 3: Cosine, 4: Cosine P10:
  TAPER_PERCENTAGE: 0.3
  CC:
    TSHIFT_LIM: [-5.0, 5.0]
    DLNA_LIM: [-1.5, 1.5]
    CC_MIN: 0.7
    DT_SIGMA_MIN: 1.0
    DLNA_SIGMA_MIN: 0.5
  MT:
    NUM_TAPER: 5
    MT_NW: 4.0
    PHASE_STEP: 1.5
    TRANSFUNC_WATERLEVEL: 1e-10
    WATER_THRESHOLD: 0.02
    DT_FAC: 2.0
    ERR_FAC: 2.5
    DT_MAX_SCALE: 3.5
    MIN_CYCLE_IN_WINDOW: 3
    USE_MT_ERROR: False
  ENV:
    WTR_ENV: 0.2

Adjoint Source Tapering

• ITAPER_TYPE: Type of tapering window for the adjoint source measurements. It can be
  • 1 for Hanning window.
  • 2 for Hamming window.
  • 3 for Cosine window.
  • 4 for Cosine P10 window.
• TAPER_PERCENTAGE: Percentage of the tapering window. It is used to determine the length of the tapering window based on the total length of the time window.

Cross-correlation (CC) Measurement Parameters

• CC.TSHIFT_LIM: Time shift limits in seconds for cross-correlation measurement.
• CC.DLNA_LIM: Logarithmic amplitude ratio limits for cross-correlation measurement.
• CC.CC_MIN: Minimum cross-correlation coefficient to accept a measurement.
• CC.DT_SIGMA_MIN: Minimum standard deviation of time shift measurement in seconds.
• CC.DLNA_SIGMA_MIN: Minimum standard deviation of logarithmic amplitude ratio.

Multi-taper (MT) Measurement Parameters

• MT.NUM_TAPER: Number of tapers to use in multi-taper measurement.
• MT.MT_NW: bin width of multitapers (nw*df is the half bandwidth of multitapers in frequency domain, typical values are 2.5, 3., 3.5, 4.0)
• MT.PHASE_STEP: maximum step for cycle skip correction .
• MT.TRANSFUNC_WATERLEVEL: Waterlevel for the transfer function in multi-taper measurement.
• MT.WATER_THRESHOLD: The triggering value to stop the search. If the spectra is larger than 10*water_threshold it will trigger the search again, works like the heating thermostat.
• MT.DT_FAC: percentage of wave period at which measurement range is too large and MTM reverts to CCTM misfit.
• MT.ERR_FAC: percentage of error at which error is too large.
• MT.DT_MAX_SCALE: Used to calculate maximum allowable time shift
• MT.MIN_CYCLE_IN_WINDOW: Minimum number of cycles in the time window for multi-taper measurement.
• MT.USE_MT_ERROR: Whether to use multi-taper error for normalization.

Envelope (ENV) Measurement Parameters

• ENV.WTR_ENV: Waterlevel for envelope measurement.

MODEL_GRID Section

SpecFWAT update model parameters on a regular grid, which is easily to take sum of gradients of different data sets with different mesh. The size of each regular grid should be small to ensure the accuracy of gradient information. We recommended set the size of each regular grid to be at least half of element size of the mesh.

MODEL_GRID:
  REGULAR_GRID_MIN_COORD: [833950, -44274.0, -80000]
  REGULAR_GRID_INTERVAL: [2000, 2000, 1000]
  REGULAR_GRID_SIZE: [104, 48, 81]

• REGULAR_GRID_MIN_COORD: Minimum coordinate of the regular grid in meter. It is a list of three values, which are the minimum x, y, and z coordinates of the regular grid.

• REGULAR_GRID_INTERVAL: Interval of the regular grid in meter. It is a list of three values, which are the interval in x, y, and z directions.

• REGULAR_GRID_SIZE: Number of regular grid. It is a list of three values, which are the number of grid points in x, y, and z directions.

POSTPROC Section

This section defines the common parameters for the post-processing of the gradient after the adjoint simulation.

POSTPROC:
  INV_TYPE: [False, True, False] # Inversion type of noise, teleseismic, and other data
  JOINT_WEIGHT: [0.5, 1.0, 0.5]
  NORM_TYPE: 1 # Normalization type of gradient (1: default)
  TAPER_H_SUPPRESS: 5000
  TAPER_H_BUFFER: 15000
  TAPER_V_SUPPRESS: 0
  TAPER_V_BUFFER: 0
  IS_PRECOND: False

• INV_TYPE: Inversion type for each data type. It is a list of three boolean values controlling inversion of noise, teleseismic, and a third data type respectively.

• JOINT_WEIGHT: Weight for each data type in a joint inversion. It is a list of three values used to scale the gradient contribution of each data type.

• NORM_TYPE: Normalization type applied to the gradient during post-processing. The default value is 1.

• TAPER_H_SUPPRESS: Horizontal tapering length in meter to suppress the gradient at the margin of the mesh.

• TAPER_H_BUFFER: Horizontal tapering length in meter to buffer the gradient at the margin of the mesh.

• TAPER_V_SUPPRESS: Vertical tapering length in meter to suppress the gradient at the margin of the mesh.

• TAPER_V_BUFFER: Vertical tapering length in meter to buffer the gradient at the margin of the mesh.

• IS_PRECOND: Whether directly apply preconditioning to the gradient. If set to False, the preconditioner will be saved and apply on L-BFGS as rescale vector (Modrak and Tromp 2016 ).

MODEL_UPDATE Section

This section defines the parameters for the model update based on the gradient and optimization method.

MODEL_UPDATE:
  INIT_MODEL_PATH: initial_model.h5
  MODEL_TYPE: 1 # 1: vp,vs,rho; 2: L,N,Gc,Gs
  MODEL_TYPE_2: # Only required when MODEL_TYPE=2
    ALT_INV: False # Alternating inversion flag
    CURRENT_INV_TYPE: 1 # Current inversion type (1 or 2)
  OPT_METHOD: 2 # Optimization method, 1: SD; 2: LBFGS 
  ITER_START: 0
  LBFGS_M_STORE: 5
  MAX_SLEN: 0.02 # Maximum step length
  MAX_SHRINK: 0.618 # Maximum shrink factor
  MAX_SUB_ITER: 10
  DO_LS: False # Do line search
  C1: 0.01
  VPVS_RATIO_RANGE: [1.3, 2.5] # Min and max limitation of Vp/Vs ratio

• INIT_MODEL_PATH: Path to the initial model file. It should be a HDF5 file with the same format as DATA/tomo_files/tomography_model.h5, but the size could be different from the MODEL_GRID section. An interpolation will be performed to interpolate INIT_MODEL_PATH to MODEL_GRID size and save to optimize/model_M00.h5.

• MODEL_TYPE: Type of model parameters to be updated. It can be

  • 1 for updating $Vp$, $Vs$, and $\rho$.
  • 2 for updating azimuthal anisotropic parameters ($Vp$, $Vs$, $\rho$, $L$, $N$, $Gc'$, and $Gs'$).

• MODEL_TYPE_2: Sub-section required only when MODEL_TYPE is 2. Contains:

  • ALT_INV: Whether to use alternating inversion (alternate between isotropic and anisotropic updates). Default False.
  • CURRENT_INV_TYPE: Current inversion type for alternating inversion. Accepts 1 or 2. Default 1.

• OPT_METHOD: Optimization method for model update. It can be

  • 1 for steepest descent (SD).
  • 2 for L-BFGS.
  • 3 for conjugate gradient (CG).

• ITER_START: Starting iteration number of L-BFGS optimization.

• LBFGS_M_STORE: Number of previous iterations to store in L-BFGS optimization.

• MAX_SLEN: Maximum step length for model update.

• MAX_SHRINK: Maximum shrink factor for model update. It is used to control the step length reduce when the model update does not decrease the objective function.

• MAX_SUB_ITER: Maximum number of sub-iterations for model update.

• DO_LS: Whether to perform line search for model update. If set to True, the line search will be performed to find the optimal step length for model update, otherwise, the step length will be fixed to MAX_SLEN. The default value is False.

• C1: Constant for the Armijo condition in line search. It is used to determine whether the step length is sufficient to decrease the objective function. The default value is 0.01.

OUTPUT Section

This section controls the verbose output of the SpecFWAT workflow

OUTPUT:
  IS_OUTPUT_PREPROC: True # Output preprocessed data
  IS_OUTPUT_ADJ_SRC: True # Output adjoint sources
  IS_OUTPUT_EVENT_KERNEL: True # Output kernels
  IS_OUTPUT_SUM_KERNEL: True # Output sum of kernels
  IS_OUTPUT_INV_GRID: False # Output inverse grid
  IS_OUTPUT_DIRECTION: True

• IS_OUTPUT_PREPROC: Whether to output the preprocessed data in SAC format. If set to True, the preprocessed data will be saved in the solver/{model}.{simu_type}/{event_name}/OUTPUT_FILES/ directory.

• IS_OUTPUT_ADJ_SRC: Whether to output the adjoint sources in SAC format. If set to True, the adjoint sources will be saved in the solver/{model}.{simu_type}/{event_name}/SEM/ directory.

• IS_OUTPUT_EVENT_KERNEL: Whether to output the kernels for each event. If set to True, the kernels will be kept in the solver/{model}.{simu_type}/{event_name}/EKERNEL/ directory, otherwise, the kernels will be deleted after the post-processing step.

• IS_OUTPUT_SUM_KERNEL: Whether to output the sum of kernels after post-processing. If set to True, the sum of kernels will be saved in the optimize/SUM_KERNEL_{model}/ directory.

• IS_OUTPUT_INV_GRID: Whether to output the inverse grid. If set to True, the inverse grid will be saved in the optimize/ directory.

• IS_OUTPUT_DIRECTION: Whether to output the final descent direction of each iteration. If set to True, the direction will be saved in the optimize/ directory.