sphinx.addnodesdocument)}( rawsourcechildren](docutils.nodessubstitution_definition)}(h.. |rel| replace:: 10.71h]h Text10.71}parenth sba attributes}(ids]classes]names]reladupnames]backrefs]utagnameh source lineKhh _documenthubh )}(h.. |date| replace:: |today| h]h2024 136}(hh)h(hh%Nh'Nubah}(h]h]h]dateah ]h"]uh$h h%h&h'Khhh(hubh target)}(h.. _intro_proc_sestbl.:h]h}(h]h]h]h ]h"]refidintro-proc-sestbluh$h8h'Khhh(hh%M/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.rstubh section)}(hhh](h title)}(h:file:`sestbl.`h]h literal)}(hhPh]hsestbl.}(hhTh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhhNubah}(h]h]h]h ]h"]uh$hLhhIh(hh%hFh'Kubh paragraph)}(hXThe 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.h](h$The processing within GAMIT itself (}(hhlh(hh%Nh'Nubhliteral_strong)}(h:program:`arc`h]harc}(hhvh(hh%Nh'Nubah}(h]h]programah]h ]h"]uh$hthhlubh, }(hhlh(hh%Nh'Nubhu)}(h:program:`model`h]hmodel}(hhh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthhlubh, }hhlsbhu)}(h:program:`autcln`h]hautcln}(hhh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthhlubh, }hhlsbhu)}(h:program:`solve`h]hsolve}(hhh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthhlubh) is controlled by }(hhlh(hh%Nh'NubhS)}(h:file:`sestbl.`h]hsestbl.}(hhh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhhlubh, }hhlsbhS)}(h:file:`sittbl.`h]hsittbl.}(hhh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhhlubh, and }(hhlh(hh%Nh'NubhS)}(h:file:`autcln.cmd`h]h autcln.cmd}(hhh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhhlubh. The }(hhlh(hh%Nh'NubhS)}(h:file:`sestbl.`h]hsestbl.}(hhh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhhlubh 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.}(hhlh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'KhhIh(hubh literal_block)}(hXChoice of Experiment = BASELINE ; BASELINE/RELAX./ORBIT Satellite Constraint = Y ; Y/N (next two lines are read and free-format. only all needed on next line) all a e i n w M D Y B 1UDC 1UDS 1UYC 1UYC 1UBC 1UBS 2UC 2US 4UC 4US SX SY SZ 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.001 0.001 0.001 h]hXChoice of Experiment = BASELINE ; BASELINE/RELAX./ORBIT Satellite Constraint = Y ; Y/N (next two lines are read and free-format. only all needed on next line) all a e i n w M D Y B 1UDC 1UDS 1UYC 1UYC 1UBC 1UBS 2UC 2US 4UC 4US SX SY SZ 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.001 0.001 0.001 }hjsbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl. xml:spacepreserveforcelanguagetexthighlight_args} linenostartKsuh$jh%hFh'K hhIh(hubhk)}(hXbBy 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).h](h By setting }(hj3h(hh%Nh'NubhS)}(h*:content:`Choice of experiment = BASELINE`h]hChoice of experiment = BASELINE}(hj;h(hh%Nh'Nubah}(h]h](code highlightcontenttexteh]h ]h"]languagejIuh$hRhj3ubhX, 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 }(hj3h(hh%Nh'NubhS)}(h':content:`Choice of experiment = RELAX`h]hChoice of experiment = RELAX}(hjRh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhj3ubhe and apply constraints to the orbital parameters. The values given for orbital parameters (1 part in }(hj3h(hh%Nh'Nubh math)}(h :math:`10^9`h]h10^9}(hjgh(hh%Nh'Nubah}(h]h]h]h ]h"]uh$jehj3ubh; or about 2 cm) are reasonable for current IGS orbits (see }(hj3h(hh%Nh'Nubh reference)}(h$http://acc.igs.org/igsacc_final.htmlh]h$http://acc.igs.org/igsacc_final.html}(hj{h(hh%Nh'Nubah}(h]h]h]h ]h"]refurij}uh$jyhj3ubh! for orbital accuracy over time).}(hj3h(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'KhhIh(hubj)}(hXtChoice of Observable = LC_AUTCLN ; L1_SINGLE/L1&L2/L1_ONLY/L2_ONLY/LC_ONLY/ ; L1,L2_INDEPEND./LC_HELP/LC_AUTCLN Ambiguity resolution WL = 0.15 0.15 1000. 99. 15000. ; used for LC_HELP only Ambiguity resolution NL = 0.15 0.15 1000. 99. 15000. ; Allow long baselines with LC_AUTCLN Ionospheric Constraints = 0.0 mm + 8.00 ppm h]hXtChoice of Observable = LC_AUTCLN ; L1_SINGLE/L1&L2/L1_ONLY/L2_ONLY/LC_ONLY/ ; L1,L2_INDEPEND./LC_HELP/LC_AUTCLN Ambiguity resolution WL = 0.15 0.15 1000. 99. 15000. ; used for LC_HELP only Ambiguity resolution NL = 0.15 0.15 1000. 99. 15000. ; Allow long baselines with LC_AUTCLN Ionospheric Constraints = 0.0 mm + 8.00 ppm }hjsbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.j+j,j-j.textj0}j2Ksuh$jh%hFh'KhhIh(hubhk)}(hX2The 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`).h](hNThe next four entries control the observations used and ambiguity resolution. }(hjh(hh%Nh'NubhS)}(h:content:`LC_AUTCLN`h]h LC_AUTCLN}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh and }(hjh(hh%Nh'NubhS)}(h:content:`LC_HELP`h]hLC_HELP}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhb are the usual choices with dual-frequency receivers and all but the shortest baselines. With the }(hjh(hh%Nh'NubhS)}(h:content:`LC_AUTCLN`h]h LC_AUTCLN}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh? option, the widelane ambiguities are assigned and resolved in }(hjh(hh%Nh'Nubhu)}(h:program:`autcln`h]hautcln}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh using the pseudoranges. With }(hjh(hh%Nh'NubhS)}(h:content:`LC_HELP`h]hLC_HELP}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhT, the wide-lane ambiguities are resolved by applying an ionospheric constraint (see }(hjh(hh%Nh'Nubh pending_xref)}(h:cite:t:`Dong_et_al_1998`h]h inline)}(hjh]hDong_et_al_1998}(hjh(hh%Nh'Nubah}(h]h](xrefcitecite-teh]h ]h"]uh$jhjubah}(h]id1ah]h]h ]h"]refdocintro/proc/control/sestbl refdomainjreftypet refexplicitrefwarn reftargetDong_et_al_1998uh$j h%hFh'Khjubh>). For data acquired with codeless receivers (before ~ 1994), }(hjh(hh%Nh'NubhS)}(h:content:`LC_HELP`h]hLC_HELP}(hj8h(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhE should be used. The sensitivity parameters for both wide-lane (with }(hjh(hh%Nh'NubhS)}(h:content:`LC_HELP`h]hLC_HELP}(hjKh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh) and narrow-lane ambiguities are set conservatively but may be varied if you need to squeeze greater accuracy out of short observing sessions. See }(hjh(hh%Nh'Nubj )}(h:cite:t:`Dong_et_al_1998`h]j)}(hj`h]hDong_et_al_1998}(hjbh(hh%Nh'Nubah}(h]h](jcitecite-teh]h ]h"]uh$jhj^ubah}(h]id2ah]h]h ]h"]refdocj, refdomainjlreftypej/ refexplicitrefwarnj2Dong_et_al_1998uh$j h%hFh'Khjubh and Section 3.5 of the }(hjh(hh%Nh'Nubjz)}(hG`GAMIT Reference Manual `_h]hGAMIT Reference Manual}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]nameGAMIT Reference Manualrefuri+http://geoweb.mit.edu/gg/docs/GAMIT_Ref.pdfuh$jyhjubh9)}(h. h]h}(h]gamit-reference-manualah]h]gamit reference manualah ]h"]refurijuh$h8 referencedKhjubh 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 }(hjh(hh%Nh'NubhS)}(h:content:`LC_HELP`h]hLC_HELP}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh, }(hjh(hh%Nh'NubhS)}(h:content:`L1_ONLY`h]hL1_ONLY}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh, }hjsbhS)}(h:content:`L2_ONLY`h]hL2_ONLY}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh, and }(hjh(hh%Nh'NubhS)}(h:content:`L1,L2_INDEP`h]h L1,L2_INDEP}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhf and compare the results. If your data were acquired with a single-frequency receiver, you should set }(hjh(hh%Nh'NubhS)}(h:content:`L1_RECEIVER`h]h L1_RECEIVER}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh in }(hjh(hh%Nh'NubhS)}(h:file:`sestbl.`h]hsestbl.}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh and also add the }(hjh(hh%Nh'NubhS)}(h:content:`L1only`h]hL1only}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh command to }(hjh(hh%Nh'NubhS)}(h:file:`autcln.cmd`h]h autcln.cmd}(hj-h(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh since }(hjh(hh%Nh'Nubhu)}(h:program:`autcln`h]hautcln}(hjBh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubhD by default will attempt to use both L1 and L2 in cleaning (even if }(hjh(hh%Nh'NubhS)}(h:content:`L1_ONLY`h]hL1_ONLY}(hjTh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh is requested for }(hjh(hh%Nh'Nubhu)}(h:program:`solve`h]hsolve}(hjgh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh).}(hjh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'KhhIh(hubj)}(hXHDMap = GMF ; GMF(default)/NMFH/VMF1 WMap = GMF ; GMF(default)/NMFW/VMF1 Use map.grid = Y ; VMF1 grid file with mapping functions and ZHD Met obs source = UFL GPT 50 ; hierarchical list: RNX ufile GPT/STP [humid value]; default GTP 50 if [humid value] < 0, use RNX or UFL if available Zenith Delay Estimation = Y ; Yes/No (default No) Interval zen = 2 ; 2 hrs = 13 knots/day (default is 1 ZD per day) Zenith Variation = 0.02 100. ; zenith-delay variation, tau in meters/sqrt(hr), hrs (default .02 100.) Elevation Cutoff = 0 ; default 0 to use value in autcln.cmd Atmospheric gradients = Y ; Yes/Np (default No) Number gradients = 2 ; number of gradient parameters per day (NS or ES); default 1 Gradient Constraints = 0.03 ; gradient at 10 deg elevation in meters; default 0.03 m Gradient Variation = .01 100 ; gradient variation Output met = N ; write the a priori met values to a z-file (Y/N) h]hXHDMap = GMF ; GMF(default)/NMFH/VMF1 WMap = GMF ; GMF(default)/NMFW/VMF1 Use map.grid = Y ; VMF1 grid file with mapping functions and ZHD Met obs source = UFL GPT 50 ; hierarchical list: RNX ufile GPT/STP [humid value]; default GTP 50 if [humid value] < 0, use RNX or UFL if available Zenith Delay Estimation = Y ; Yes/No (default No) Interval zen = 2 ; 2 hrs = 13 knots/day (default is 1 ZD per day) Zenith Variation = 0.02 100. ; zenith-delay variation, tau in meters/sqrt(hr), hrs (default .02 100.) Elevation Cutoff = 0 ; default 0 to use value in autcln.cmd Atmospheric gradients = Y ; Yes/Np (default No) Number gradients = 2 ; number of gradient parameters per day (NS or ES); default 1 Gradient Constraints = 0.03 ; gradient at 10 deg elevation in meters; default 0.03 m Gradient Variation = .01 100 ; gradient variation Output met = N ; write the a priori met values to a z-file (Y/N) }hjsbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.j+j,j-j.textj0}j2Ksuh$jh%hFh'K!hhIh(hubhk)}(hXu 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).h](hXThe 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 }(hjh(hh%Nh'Nubj )}(h:cite:t:`Boehm_et_al_2006b`h]j)}(hjh]hBoehm_et_al_2006b}(hjh(hh%Nh'Nubah}(h]h](jcitecite-teh]h ]h"]uh$jhjubah}(h]id3ah]h]h ]h"]refdocj, refdomainjreftypej/ refexplicitrefwarnj2Boehm_et_al_2006buh$j h%hFh'K%hjubhX% 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 }(hjh(hh%Nh'Nubj )}(h:cite:p:`Boehm_et_al_2006a`h]j)}(hjh]hBoehm_et_al_2006a}(hjh(hh%Nh'Nubah}(h]h](jcitecite-peh]h ]h"]uh$jhjubah}(h]id4ah]h]h ]h"]refdocj, refdomainjreftypep refexplicitrefwarnj2Boehm_et_al_2006auh$j h%hFh'K%hjubhX. You have access to these values by downloading the VMF1 grid files for each year from }(hjh(hh%Nh'Nubjz)}(h ftp://everest.mit.edu/pub/GRIDS/h]h ftp://everest.mit.edu/pub/GRIDS/}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]refurijuh$jyhjubh, and setting }(hjh(hh%Nh'NubhS)}(h:content:`Use map.grid = Y`h]hUse map.grid = Y}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh, }(hjh(hh%Nh'NubhS)}(h:content:`DMap = VMF1`h]h DMap = VMF1}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh and }(hjh(hh%Nh'NubhS)}(h:content:`WMap = VMF1`h]h WMap = VMF1}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhd. The next command controls the source of pressure (most importantly, but also temperature) for the }(hjh(hh%Nh'Nubh emphasis)}(h *a priori*h]ha priori}(hj0h(hh%Nh'Nubah}(h]h]h]h ]h"]uh$j.hjubh 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 }(hjh(hh%Nh'NubhS)}(h :file:`met/`h]hmet/}(hjBh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh directory for use by }(hjh(hh%Nh'Nubhu)}(h:program:`sh_gamit`h]hsh_gamit}(hjWh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh. Setting the first option of }(hjh(hh%Nh'NubhS)}(h:content:`Met obs source = RNX`h]hMet obs source = RNX}(hjih(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhX  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 }(hjh(hh%Nh'NubhS)}(h:content:`Met obs source = UFL`h]hMet obs source = UFL}(hj|h(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh and }hjsbhS)}(h:content:`map.grid = Y`h]h map.grid = Y}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh(. (If you use a VMF1 station-list file, }(hjh(hh%Nh'NubhS)}(h:file:`map.list`h]hmap.list}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh, the second option, }(hjh(hh%Nh'NubhS)}(h:content:`GPT`h]hGPT}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh or }(hjh(hh%Nh'NubhS)}(h:content:`UFL`h]hUFL}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh 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 }(hjh(hh%Nh'NubhS)}(h:file:`met.grid`h]hmet.grid}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh and }hjsbhS)}(h:file:`met.list`h]hmet.list}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh 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 }(hjh(hh%Nh'Nubj )}(h:cite:t:`Lagler_et_al_2013`h]j)}(hj h]hLagler_et_al_2013}(hj h(hh%Nh'Nubah}(h]h](jcitecite-teh]h ]h"]uh$jhjubah}(h]id5ah]h]h ]h"]refdocj, refdomainjreftypej/ refexplicitrefwarnj2Lagler_et_al_2013uh$j h%hFh'K%hjubh. It reads from table }(hjh(hh%Nh'NubhS)}(h:file:`gpt.grid`h]hgpt.grid}(hj+h(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh in }(hjh(hh%Nh'NubhS)}(h:file:`~/gg/tables/`h]h ~/gg/tables/}(hj@h(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubhX 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 }(hjh(hh%Nh'Nubj/)}(h *a priori*h]ha priori}(hjUh(hh%Nh'Nubah}(h]h]h]h ]h"]uh$j.hjubh- pressure [Tregoning and Herring, 2006]. The }(hjh(hh%Nh'NubhS)}(h:content:`Interval zen`h]h Interval zen}(hjgh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubhX; 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 }(hjh(hh%Nh'Nubjz)}(hG`GAMIT Reference Manual `_h]hGAMIT Reference Manual}(hjzh(hh%Nh'Nubah}(h]h]h]h ]h"]nameGAMIT Reference Manualj+http://geoweb.mit.edu/gg/docs/GAMIT_Ref.pdfuh$jyhjubh9)}(h. h]h}(h]id6ah]h]h ]gamit reference manualah"]refurijuh$h8jKhjubhz). For atmospheric gradients, we allow a linear change during the session (two parameters each for N/S and E/W gradients).}(hjh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'K%hhIh(hubj)}(hKUpdate T/L files = L_ONLY ; T_AND_L (default), T_ONLY, L_ONLY, NONE h]hKUpdate T/L files = L_ONLY ; T_AND_L (default), T_ONLY, L_ONLY, NONE }hjsbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.j+j,j-j.textj0}j2Ksuh$jh%hFh'K;hhIh(hubhk)}(hX?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`.h](h9The command to update the L-file assures that the second }(hjh(hh%Nh'Nubhu)}(h:program:`solve`h]hsolve}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubhA solution will have adjustments within a linear range and allows }(hjh(hh%Nh'Nubhu)}(h:program:`sh_gamit`h]hsh_gamit}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh 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 }(hjh(hh%Nh'NubhS)}(h:file:`sestbl.`h]hsestbl.}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjubh3 command, not shown). If you know that all of your }(hjh(hh%Nh'Nubj/)}(h *a priori*h]ha priori}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]uh$j.hjubh} 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 }(hjh(hh%Nh'NubhS)}(h":content:`Update T/L files = NONE`h]hUpdate T/L files = NONE}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh.}(hjh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'K?hhIh(hubj)}(hAUTCLN Command File = autcln.cmd ; Filename; default none (use default options) Station Error = ELEVATION 10 5 ; 1-way L1, a**2 + (b**2)(L**2) in mm, ppm, default = 10. 0. h]hAUTCLN Command File = autcln.cmd ; Filename; default none (use default options) Station Error = ELEVATION 10 5 ; 1-way L1, a**2 + (b**2)(L**2) in mm, ppm, default = 10. 0. }hj sbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.j+j,j-j.textj0}j2Ksuh$jh%hFh'KChhIh(hubhk)}(hXYThe 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.h](hjThe next two commands control how the phase data are weighted in the preliminary and final solutions. The }(hj2h(hh%Nh'NubhS)}(h:content:`Station Error`h]h Station Error}(hj:h(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhj2ubh entry shown tells }(hj2h(hh%Nh'Nubhu)}(h:program:`solve`h]hsolve}(hjMh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthj2ubh to use elevation-dependent weighting for the phase date. In the preliminary solution, the assigned error is 10 mm with negligible elevation dependence. With }(hj2h(hh%Nh'NubhS)}(h:content:`AUTCLN postfit = Y`h]hAUTCLN postfit = Y}(hj_h(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhj2ubh (default with }(hj2h(hh%Nh'NubhS)}(h$:content:`Type of analysis = 1-ITER`h]hType of analysis = 1-ITER}(hjrh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhj2ubh ) the second }(hj2h(hh%Nh'Nubhu)}(h:program:`solve`h]hsolve}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthj2ubhS run will assign weights to the phase data based on the actual scatter computed by }(hj2h(hh%Nh'Nubhu)}(h:program:`autcln`h]hautcln}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthj2ubh in its “postfit” edit, and will contain both a constant and an elevation-dependent term, as recorded in the output print file }(hj2h(hh%Nh'NubhS)}(h:file:`autcln.post.sum`h]hautcln.post.sum}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhj2ubh. 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.}(hj2h(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'KGhhIh(hubj)}(hDecimation Factor = 4 ; FOR SOLVE, default = 1 Quick-pre observable = LC_ONLY ; for 1st soln, default same as Choice of observable Quick-pre decimation factor = 10 ; 1st iter or autcln pre, default same as Decimation Factor h]hDecimation Factor = 4 ; FOR SOLVE, default = 1 Quick-pre observable = LC_ONLY ; for 1st soln, default same as Choice of observable Quick-pre decimation factor = 10 ; 1st iter or autcln pre, default same as Decimation Factor }hjsbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.j+j,j-j.textj0}j2Ksuh$jh%hFh'KNhhIh(hubhk)}(hXRelated 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.h](h`Related to the error weighting is the sampling of data for the solution. For cleaning purposes (}(hjh(hh%Nh'Nubhu)}(h:program:`model`h]hmodel}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh, }(hjh(hh%Nh'Nubhu)}(h:program:`autcln`h]hautcln}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh), 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 }(hjh(hh%Nh'Nubhu)}(h:program:`solve`h]hsolve}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjubh!. The default decimation factor (}(hjh(hh%Nh'NubhS)}(h :content:`4`h]h4}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjubh, 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 }(hjh(hh%Nh'Nubj )}(h:numref:`intro_gnss_solve`h]hS)}(hj)h]hintro_gnss_solve}(hj+h(hh%Nh'Nubah}(h]h](jstd std-numrefeh]h ]h"]uh$hRhj'ubah}(h]h]h]h ]h"]refdocj, refdomainj5reftypenumref refexplicitrefwarnj2intro_gnss_solveuh$j h%hFh'KRhjubh.) 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.}(hjh(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'KRhhIh(hubj)}(hXAntenna Model = AZEL ; NONE/ELEV/AZEL default = ELEV Tides applied = 31 ; Binary coded: 1 earth 2 freq-dep 4 pole tide (zero mean pole) ; 8 ocean 16 pole tide (IERS2010 mean pole) 32 atmosphere S1/S2 ; 64 pole tide (IERS20 secular pole) (31 default ITRF2014, 79 default ITRF2020). Use otl.list = N Use otl.grid = Y h]hXAntenna Model = AZEL ; NONE/ELEV/AZEL default = ELEV Tides applied = 31 ; Binary coded: 1 earth 2 freq-dep 4 pole tide (zero mean pole) ; 8 ocean 16 pole tide (IERS2010 mean pole) 32 atmosphere S1/S2 ; 64 pole tide (IERS20 secular pole) (31 default ITRF2014, 79 default ITRF2020). Use otl.list = N Use otl.grid = Y }hjQsbah}(h]h]h]h ]h"]sourceJ/chandler/home/floyd/public_html/gg/docs/source/intro/proc/control/sestbl.j+j,j-j.textj0}j2Ksuh$jh%hFh'KYhhIh(hubhk)}(hXNThere 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.h](hThere are over a dozen }(hjch(hh%Nh'NubhS)}(h:file:`sestbl.`h]hsestbl.}(hjkh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjcubhw entries controlling the models used in processing the data. These are all given at the bottom of the template file in }(hjch(hh%Nh'NubhS)}(h:file:`~/gg/tables/`h]h ~/gg/tables/}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjcubh% and discussed in Section 3.2 of the }(hjch(hh%Nh'Nubjz)}(hG`GAMIT Reference Manual `_h]hGAMIT Reference Manual}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]nameGAMIT Reference Manualj+http://geoweb.mit.edu/gg/docs/GAMIT_Ref.pdfuh$jyhjcubh9)}(h. h]h}(h]id7ah]h]h ]gamit reference manualah"]refurijuh$h8jKhjcubh`. The three shown here are those most likely to need editing for current applications. Changing }(hjch(hh%Nh'NubhS)}(h:content:`Antenna Model`h]h Antenna Model}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjcubh to }(hjch(hh%Nh'NubhS)}(h:content:`NONE`h]hNONE}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjcubh; allows bypassing of the phase-center-variation (PCV) file }(hjch(hh%Nh'NubhS)}(h:file:`antmod.dat`h]h antmod.dat}(hjh(hh%Nh'Nubah}(h]h]fileah]h ]h"]rolefileuh$hRhjcubhd, an option advisable only if you are using identical antennas within a small regional network. The }(hjch(hh%Nh'NubhS)}(h:content:`Tides applied`h]h Tides applied}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjcubh_ option is binary coded, allowing any combination of the various tidal models to be applied in }(hjch(hh%Nh'Nubhu)}(h:program:`model`h]hmodel}(hjh(hh%Nh'Nubah}(h]h]hah]h ]h"]uh$hthjcubhX. 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 }(hjch(hh%Nh'NubhS)}(h:content:`Tides applied = 23`h]hTides applied = 23}(hjh(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjcubh instead of }(hjch(hh%Nh'NubhS)}(h :content:`31`h]h31}(hj*h(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjcubhv. The OTL components are read from a station list file or a global grid, with station-list values taking precedent if }(hjch(hh%Nh'NubhS)}(h:content:`Use otl.list = Y`h]hUse otl.list = Y}(hj=h(hh%Nh'Nubah}(h]h](jFjGjHjIeh]h ]h"]languagejIuh$hRhjcubh 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.}(hjch(hh%Nh'Nubeh}(h]h]h]h ]h"]uh$hjh%hFh'K]hhIh(hubeh}(h](sestblhEeh]h](sestbl.intro_proc_sestbl.eh ]h"]uh$hGhhh(hh%hFh'Kexpect_referenced_by_name}j\h:sexpect_referenced_by_id}hEh:subeh}(h]h]h]h ]h"]sourcehFtranslation_progress}(totalK translatedKuuh$hcurrent_sourceN current_lineNsettingsdocutils.frontendValues)}(outputNhLN generatorN datestampN source_linkN source_urlN toc_backlinksentryfootnote_backlinksK sectnum_xformKstrip_commentsNstrip_elements_with_classesN strip_classesN report_levelK halt_levelKexit_status_levelKdebugNwarning_streamN tracebackinput_encoding utf-8-siginput_encoding_error_handlerstrictoutput_encodingutf-8output_encoding_error_handlerjerror_encodingutf-8error_encoding_error_handlerbackslashreplace language_codeenrecord_dependenciesNconfigN id_prefixhauto_id_prefixid dump_settingsNdump_internalsNdump_transformsNdump_pseudo_xmlNexpose_internalsNstrict_visitorN_disable_configN_sourcehF _destinationN _config_files]file_insertion_enabled raw_enabledKline_length_limitM'pep_referencesN pep_base_urlhttps://peps.python.org/pep_file_url_templatepep-%04drfc_referencesN rfc_base_url&https://datatracker.ietf.org/doc/html/ tab_widthKtrim_footnote_reference_spacesyntax_highlightlong smart_quotessmartquotes_locales]character_level_inline_markupdoctitle_xform docinfo_xformKsectsubtitle_xform image_loadinglinkembed_stylesheetcloak_email_addressessection_self_linkenvNubreporterNindirect_targets]substitution_defs}(hh h5h)usubstitution_names}(relhdateh5urefnames}refids}hE]h:asnameids}(j\hEj[jXjju nametypes}(j\j[juh}(hEhIjXhIj&jjsj^jjjjjjjjjjjju footnote_refs} citation_refs} autofootnotes]autofootnote_refs]symbol_footnotes]symbol_footnote_refs] footnotes] citations]autofootnote_startKsymbol_footnote_startK id_counter collectionsCounter}jKsRparse_messages](h system_message)}(hhh]hk)}(h9Duplicate explicit target name: "gamit reference manual".h]h=Duplicate explicit target name: “gamit reference manual”.}(hjh(hh%Nh'Nubah}(h]h]h]h ]h"]uh$hjhjubah}(h]h]h]h ]h"]jalevelKtypeINFOsourcehFlineKuh$jhhIh(hh%hFh'K9ubj)}(hhh]hk)}(h9Duplicate explicit target name: "gamit reference manual".h]h=Duplicate explicit target name: “gamit reference manual”.}(hj h(hh%Nh'Nubah}(h]h]h]h ]h"]uh$hjhj ubah}(h]h]h]h ]h"]jalevelKtypej sourcehFlineKuh$jhhIh(hh%hFh'Keubetransform_messages]j)}(hhh]hk)}(hhh]h7Hyperlink target "intro-proc-sestbl" is not referenced.}hj- sbah}(h]h]h]h ]h"]uh$hjhj* ubah}(h]h]h]h ]h"]levelKtypej sourcehFlineKuh$juba transformerN include_log] decorationNh(hub.