sh_gamit

Script to automate GAMIT processing

manual group:

Scripts

Description

sh_gamit is a script to automate GAMIT processing:

• sets up a directory structure for GAMIT

• gets tables and data from various archives around the world

• converts raw to RINEX if necessary

• runs GAMIT and saves phase sky-plots to check daily quality

To run for a particular experiment, create a tables/ directory within your experiment directory, copy into tables/ the following files from ~/gg/templates/ and modify them appropriately: autlcn.cmd process.defaults sestbl. sites.defaults sittbl. station.info Check other files in the templates directory. Edit them to suit the type of analysis you are running, and the options you like.

Example

$ sh_gamit -d 98 240 -orbit SIOF -expt tubi -eops usno -nopngs

Usage

sh_gamit <options>

where the following are some of the options available: (other options are set in the process.defaults file located in the tables/ directory)

-dir

Experiment directory or root directory for experiment processing. This must be an absoute path. [Default pwd]

-d yr days

where yr and days are the data to be processed e.g. 1997 153 156 178

-s yr d1 d2

where d1 is start day and d2 is stop day to be processed e.g. 1997 153 178

-r day

where day is number of days before current date.

-gnss

GNSS ( G R C E J I) default G

-lfreq #

Lower frequency selection depending on GNSS e.g., # = 5 for G L5; 6, 7 or 8 for E L6, L7 or L8; 5 or 7 for C L5 or L7 Options couples to GNSS selection. (Default, not specified and GAMIT selection used

-expt

4 char name of experiment/s to run. (list of expts can be used). [Default expt]

-orbit

Type of orbit to use. (igsf igsr igsu codf code codm codr emrf esaf gfzf gfzm grgm jaxm jplf mitf ngsf siof sior siou tumm wuhm [Default igsf] See sh_get_orbits for description. Optional second argument gives the day-of-year for the orbit if not the same as processing day

-orbt

4 char ID for g-files. Specify only to use a non-conventional, locally generated g-file [Default is lowercase of downloaded orbit type]

-srp_param

Radiation-pressure parameters to estimate in the fit to the input orbit (default allpr) (See sh_sp3fit for other options and interaction with radiation model in sestbl.)

-eops

Name of EOP series to use. (bull_b, bull_a, etc.) [Default usno]

-localeop

Use a local copy of EOP tables (series used is defined by -eops). No ftp update will be attempted. [Default off]

-remakex Y/C/N

Y - Remake x-files without checking. C - Remake x-file if orbit and session.info are incompatable and all RINEX available [Default C] N - Don’t remake x-files under any circumstances

-renamek Y/C/N

Y - Remove existing k-files and remake C - Do not remake k-files unless missing or X-files remade (default) N - Do not remake k-files – stop if missing

-remakej Y/N

Y - Remake j-file (default) N - Do not remake j-file (used with epoch-by-epoch j-files for SA with mismatched time tags)

-remakeg Y/N

Y - Remake g-file (default, unless usr g-file) N - Do not remake if it exists. Showld be used when processing sub-networks

-netext char

Add network suffix <char> to day directories (default if -gnss not in the command line is none (day number only): otherwise use the GNSS code, e.g., 301G)

-yrext

Add year prefix to day directories. Ie 1997_053

-pres

Plot gamit residuals as skyplots NO / YES (skyplot only) / ELEV (skyplot + phase vs elevation [Default N]

-nopngs

Do not create skyplot pngs.

-sessinfo

Session - sampling interval, number of epochs, start time (default 30 2880 0 0)

-rx_doy_plus

Search RINEX file names this many days into the future for current-day data

-rx_doy_minus

Search RINEX file names this many days into the past for current-day data

-minspan

Minimum span of data in rinex files in hours [Default is 0.1 (6 min)]

-mailto

Name of user and machine to send results summary (e.g. example@example.com)

-snrres

Calculate theoretical phase residuals from SNR. (requires raw data or rinex W/SNR)

-noacc

Do not exclude satellites base on accuracy codes in the .sp3 file. [Default exclude]

-ofit val Maximum orbit misfit allowed when fitting GMAIT orbit model to .SP3 file. [Default 0.1 m] -xsite list List of sites to be excluded from the processing. [format site:yyyy_ddd-yyyy_ddd OR site] -copt opts Types of files to be compressed in the day directory. [Default: x k ao D] -aopt opts Types of files to be archived under -archivepth <dir>. [Default: null] -dopt opts Types of files to be deleted from the day directory. [Default: c]

Possible d-, c-, and a- opts: D, H, ao, ac, as, b, c, d, e, g, h, i, j, k, l, m, met, o, p, ps, q, t, x, z, all Note: all for -dopt will also remove the day directories, so pairs logically with -aopt H h q as

-topt

opts Defines the list of table files to be relinked if they already exist [Default = none] none all nutabl. soltab. luntab. gdetic.dat guess_rcvant.dat leap.sec svnav.dat svs_exclude.dat tform.dat antmod.dat rcvant.dat dcb.dat core pole ut1

-c

Do not compress RINEX files when done [Default is to compress]

-h

Use Hatanaka RINEX compress algorithm. [Default No]

-archivepth dir Directory where files are to be archived [Default $procdir/archive] -xver char 1-character x-file version (6th character of x-file). [Default is single-digit year] -noftp Use this option if no ftp download connections are available. -met Link in RINEX met files from ../met [Default off] -ion [opts] Download IONEX files to ../ionex and link [Default igsg if no argument given]

Automatically switced on if “Ion model = GMAP” given in tables/sestbl.

-ftp_prog

Name of ftp program: ftp or ncftp [default is ftp]

-rinex_ftpsites List of ftp sites to be searched for rinex files [default sopac cddis and unavco] -subdir creates a temporary subdirectory to speed up uncompress when there are many rinex files in directory -aprfile Name of xyz coordinate .apr file to be used. [Default is itrf2000.apr] -xcheck Y/N Do not ftp rinex file if the xfile exists already [Default N] -noprocess opt Use this option if gamit processing should be avoided. [Default ftp].

Possible: ftp - ftp rinex, eop, brdc, orbits.

preproc - ftp + preprocess sites updating station.info and lfile if necessary.

-update_l (Y/N) Use this option to stop the l-file being updated by the .apr file [Default = Y] -no_lf_upd Do NOT update the master tables l-file (lfile.) when large coordinate adjustment occur - EVER.

[Default - update tables/lfile. if Update tolerance specified in sestbl. is exceeded

-fixed_rx_list Link ONLY rinex files for sites listed in sites.defaults [Default = N] -sd file Name of sites.defaults file to use [Default = $procdir/tables/sites.defaults] -metutil N/Z Output Zenith wet delay and precipitable water using zenith hydrostatic delays [Default N]

Z - Use GAMIT z-files to provide zenith hydrostatic delays (o-file interval is assumed)

-jclock

Source for j-file clocks: sp3 or brdc; brdc the default for now for GPS, forced to sp3 for non-GPS processing