.. _test_globk2018:

Run GLOBK to get repeatabilities and a combined H-file for the span
===================================================================

From the :file:`2018/` directory:

.. code-block:: console

   $ sh_glred -cmd

to get the :file:`globk.cmd` and :file:`glorg.cmd` files copied from :file:`~/gg/tables/` to :file:`2018/gsoln/`, then

.. code-block:: console

   $ sh_glred -s 2018 095 2018 097 -expt eura -netext G R E C -opt R H G T >& sh_glred.log

The script as commanded will translate the GAMIT plain text h-files in each day directory to GLOBK binary h-files (:file:`.glx`) and put them into the :file:`glbf/` directory (:console:`H` option);
create a :file:`.gdl` file for each day listing the h-file, run GLOBK for each day (:console:`G` option) using :file:`globk.cmd` and :file:`glorg.cmd`;
and generate time series plots (:console:`T` option using :console:`.pos` files, program :program:`tssum` and :program:`sh_plot_pos`).

For this test the :program:`globk`/:program:`glred` and :program:`glorg` command files created in :file:`gsoln/` by :console:`sh_glred -cmd` will work without editing;
however, in general there are changes you may need to make, for example, to define the reference frame with translation-only rather than tranlsation and rotation, and to change the list of reference (stabilization) sites.

The daily plots combining the GNSS measurements for each day, may be found in :file:`gsoln/plots_2018_095-2018_097/`.
You can view the :file:`.org` files to see what :program:`globk` did, paying attention to the stabilization iteration (what sites were removed), the :math:`\chi^2` increments in the :console:`EXPERIMENT LIST`, and the statitics of the stabilization (:console:`POS STATISTICS`).
In your own runs processing many days, it is convenient to get a quick summary by :command:`grep`'ing on :file:`.org` with :console:`USED` (from the :math:`\chi^2` list) and :console:`POS STAT`.

Next combine the days to get a single H-file for the 3-day span to use for the multi-year repeatabilities (time series) and velocities.

.. code-block:: console

  $ sh_glred -s 2018 095 2018 097 -expt eura -netext G R E C -ncomb 3 -globk_cmd_prefix COMB -opt G >& sh_glred_comb.log