.. _intro_proc_sestbl.: :file:`sestbl.` --------------- The processing within GAMIT itself (:program:`arc`, :program:`model`, :program:`autcln`, :program:`solve`) is controlled by :file:`sestbl.`, :file:`sittbl.`, and :file:`autcln.cmd`. The :file:`sestbl.` entries that you are likely to consider in setting up your processing are shown below. The values are generally appropriate for regional or global processing for most data acquired after 1995. .. literalinclude:: sestbl. :lines: 17,6-8 :language: text By setting :content:`Choice of experiment = BASELINE`, you are fixing the orbits and omitting the orbital parameters from your GAMIT processing and output h-files. If you intend to combine your h-files with those from global processing at, e.g., MIT or SOPAC, or if you are processing a GNSS with poor orbits, then you should set :content:`Choice of experiment = RELAX` and apply constraints to the orbital parameters. The values given for orbital parameters (1 part in :math:`10^9` or about 2 cm) are reasonable for current IGS orbits (see http://acc.igs.org/igsacc_final.html for orbital accuracy over time). .. literalinclude:: sestbl. :lines: 20-21,30-31,29 :language: text The next four entries control the observations used and ambiguity resolution. :content:`LC_AUTCLN` and :content:`LC_HELP` are the usual choices with dual-frequency receivers and all but the shortest baselines. With the :content:`LC_AUTCLN` option, the widelane ambiguities are assigned and resolved in :program:`autcln` using the pseudoranges. With :content:`LC_HELP`, the wide-lane ambiguities are resolved by applying an ionospheric constraint (see :cite:t:`Dong_et_al_1998`). For data acquired with codeless receivers (before ~ 1994), :content:`LC_HELP` should be used. The sensitivity parameters for both wide-lane (with :content:`LC_HELP`) and narrow-lane ambiguities are set conservatively but may be varied if you need to squeeze greater accuracy out of short observing sessions. See :cite:t:`Dong_et_al_1998` and Section 3.5 of the `GAMIT Reference Manual `_ for a detailed discussion. For baselines less than a few km, using L1 and L2 independently or only L1 may reduce noise compared to using the ionosphere-free combination. To check for the size of ionospheric errors, it’s wise to try :content:`LC_HELP`, :content:`L1_ONLY`, :content:`L2_ONLY`, and :content:`L1,L2_INDEP` and compare the results. If your data were acquired with a single-frequency receiver, you should set :content:`L1_RECEIVER` in :file:`sestbl.` and also add the :content:`L1only` command to :file:`autcln.cmd` since :program:`autcln` by default will attempt to use both L1 and L2 in cleaning (even if :content:`L1_ONLY` is requested for :program:`solve`). .. literalinclude:: sestbl. :lines: 50-51,53,45-46,33-34,36-40,109,47 :language: text The next block of entries controls the modeling and estimation of the atmospheric delay. The first command specifies the mapping function to be used to project the zenith delay to the elevation angle of the satellites. Different functions must be used for the hydrostatic ("dry") delay and the "wet" delay due to water vapor. The default, sufficient for all but the most accurate studies of heights and for meteorological studies, is the "global mapping function" (GMF) developed by :cite:t:`Boehm_et_al_2006b` from fitting numerical weather model (NWM) data over 20 years. A more accurate reconstruction of the NWM data can be obtained by interpolating hydrostatic and wet mapping function coefficients as a function of time and location from the (large) global grid files compiled by the Vienna group :cite:p:`Boehm_et_al_2006a`. You have access to these values by downloading the VMF1 grid files for each year from ftp://everest.mit.edu/pub/GRIDS/, and setting :content:`Use map.grid = Y`, :content:`DMap = VMF1` and :content:`WMap = VMF1`. The next command controls the source of pressure (most importantly, but also temperature) for the *a priori* zenith hydrostatic delay (ZHD). The most accurate values, if they are available, are from local measurements of surface pressure, which can be written into RINEX met files and stored in a :file:`met/` directory for use by :program:`sh_gamit`. Setting the first option of :content:`Met obs source = RNX` tells GAMIT to use these for any station for which they are available. The next most accurate source would be the ZHD values from the VMF1 grid files. Since GAMIT reads the grid file for all sites used for the day and writes the values into the U-file, you select this option by setting :content:`Met obs source = UFL` and :content:`map.grid = Y`. (If you use a VMF1 station-list file, :file:`map.list`, the second option, :content:`GPT` or :content:`UFL` can be used for stations missing from the list file. Since the only source of NWM ZHD data currently tabulated for GAMIT is the VMF1 mapping-function grids, the :file:`met.grid` and :file:`met.list` files are not yet supported.) Finally, as with the mapping functions, the Vienna group has constructed from the NWM data an analytical model, designated "global pressure and temperature" (GPT2), of :cite:t:`Lagler_et_al_2013`. It reads from table :file:`gpt.grid` in :file:`~/gg/tables/` values of zenith hydrotstatic delay (ZHD), temperature, lapse rate, and dry and wet mapping functions as a function of latitude, longitude, and day-of-year that are averages from a fit to 10-year monthly averages from a global numerical weather model. The option STP implies standard constants (1013.25 hPa, 20°C), used prior to Release 10.3 but now effectively obsolete. To keep errors in height estimates below 2 mm, you need an accuracy of about 10 hPa in *a priori* pressure [Tregoning and Herring, 2006]. The :content:`Interval zen` command controls the number of zenith delay parameters estimated using the session. For geodetic studies, estimating values at 2-hour intervals with the constraints given is more than adequate. For meteorological studies, you may want to estimate more parameters and/or alter the constraints (see Chapter 7 of the `GAMIT Reference Manual `_). For atmospheric gradients, we allow a linear change during the session (two parameters each for N/S and E/W gradients). .. literalinclude:: sestbl. :lines: 42 :language: text The command to update the L-file assures that the second :program:`solve` solution will have adjustments within a linear range and allows :program:`sh_gamit` to apply updated coordinates to the processing of successive days. The criterion for updating is set by default to be adjustments larger than 30 cm (though this can be changed with a :file:`sestbl.` command, not shown). If you know that all of your *a priori* coordinates are accurate to within a few cm and don't want bad data to corrupt the coordinates in the L-file, you can set :content:`Update T/L files = NONE`. .. literalinclude:: sestbl. :lines: 24,22 :language: text The next two commands control how the phase data are weighted in the preliminary and final solutions. The :content:`Station Error` entry shown tells :program:`solve` to use elevation-dependent weighting for the phase date. In the preliminary solution, the assigned error is 10 mm with negligible elevation dependence. With :content:`AUTCLN postfit = Y` (default with :content:`Type of analysis = 1-ITER`) the second :program:`solve` run will assign weights to the phase data based on the actual scatter computed by :program:`autcln` in its "postfit" edit, and will contain both a constant and an elevation-dependent term, as recorded in the output print file :file:`autcln.post.sum`. Elevation-dependent weighting is recommended for almost all analyses; an exception might be tests of models (e.g. mapping functions) sensitive to low-elevation observations. .. literalinclude:: sestbl. :lines: 25,27,26 :language: text Related to the error weighting is the sampling of data for the solution. For cleaning purposes (:program:`model`, :program:`autcln`), the full sampling of the x-file (usually 30 s or 15 s) is used. However, since the phase errors are correlated over many minutes, it is not necessary to use sampling this frequent in :program:`solve`. The default decimation factor (:content:`4`, resulting in 2-minute sampling for the usual case) provides formal uncertainties that are usually reasonable within a factor of two while reducing time significantly. (See the discussion of uncertainties in :numref:`intro_gnss_solve`.) For the preliminary solution used to get decimeter-level coordinates for editing and to avoid non-linearity in the final adjustments, 5-minute sampling is sufficient. .. literalinclude:: sestbl. :lines: 74,60-62,66-67 :language: text There are over a dozen :file:`sestbl.` entries controlling the models used in processing the data. These are all given at the bottom of the template file in :file:`~/gg/tables/` and discussed in Section 3.2 of the `GAMIT Reference Manual `_. The three shown here are those most likely to need editing for current applications. Changing :content:`Antenna Model` to :content:`NONE` allows bypassing of the phase-center-variation (PCV) file :file:`antmod.dat`, an option advisable only if you are using identical antennas within a small regional network. The :content:`Tides applied` option is binary coded, allowing any combination of the various tidal models to be applied in :program:`model`. For most applications you may keep the default, which is to apply all of the tidal effects currently available in GAMIT (though the 2-bit for K1 is not necessary with the current IERS 2003 model for the solid-body tides, and the pole tide can be applied later in GLOBK). If your data are not strongly affected by ocean tidal loading ("OTL") or you do not wish to copy over the large ocean loading grid files, you might omit this model, setting :content:`Tides applied = 23` instead of :content:`31`. The OTL components are read from a station list file or a global grid, with station-list values taking precedent if :content:`Use otl.list = Y` and the station is within 10 km of a station in the file. If you have local stations close to an IGS station, it is advisable to use the grid exclusively so that you don't introduce artifacts in the relative motions of two stations.