gamit 6月1日更新日誌（incremental_updates.20200601.tar.gz）

12345678901234567890123456789012345678901234567890123456789012345678901234567890

                    Incremental Updates for GAMIT/GLOBK 10.71

                      Last updated: Sunday, May 31, 2020 
 
kf/track/taiutc.f: Cleaned up code and fixed unit conflict problem. TAH 200531 
 
kf/track/{theory.f, get_miss.f, track_com.h}:  Modified behavior of USE_BLQ 
   command to allow solid-Earth tide model to be specified (IERS2010 or default 
   analytic model) and to turn off model completely (NOETIDE).   NOBLQ needs to 
   be given as argument to turn off ocean-tide model  (structured this way for '
   backwards compatibility) . TAH 200531 

kf/track/track.f: Fixed rounding error in num_epochs calculation. TAH 200527.

Thursday, May 28, 2020 

tables/svs_exclude.dat.nanu: Updated following USCG NANU 2020025. MAF 20200528

Summary of features added Wednesday, May 27, 2020 
- Changes to stop fatals when data occurs in a gap in station.info entries 
  (often associated with firmware updates. 
- Changes made to increase Glonass data yields due to issues with RINEX 2 headers.   
- New cddismgex archive to get RINEX 3 data from CDDIS MGEX area (needed for data 
  collected before MGEX data was saved in the same directories are GPS data.  
- Fixed some issues with occasional failures of parallel runs.  
 
gamit/arc/earthradTUM.f: Fix calculation of mjd but has no impact due to antenna 
  thrust not being computed In this routine. TAH 200526. 
 
gamit/model/update_coords.f: Tested antenna changes to see if messages should be 
  printed about changes.   TAH 200527. 
 
gamit/makex/makex.f: Fixed format reporting frequencies. TAH 200526. 
 
gamit/makex/rhead.f: For Glonass processing: Removed code that set P2 to 
  x2 when C2 available.  Made P2 be secondary choice to C2 so that files 
  with C2 in header but no C2 data, will include the P2 data in the x-file. TAH 200526. 
 
gamit/makex/rbatch.f: Added back feature to have blank line act as EOF in 
  batch file. TAH 200526. 
 
gamit/lib/sel_obtyp.f: Fixed initialization loop to 6 (from 4) for iobtypx; 
  Changed secondary choice to P2 for Glonass when C2 in header.  TAH 200526. 
 
gamit/lib/rstnfo.f: Updated to allow a gap between end and start times of 
  up 6hrs.  This happens with firmware updates where the  gap can be 1-minute 
  to a few hours.  There is probably no data in gap but GAMIT is still generating 
  epochs during the gap (makex and maybe model fatals). Cleaned up  indenting around 
  modified code. TAH 200524 
 
com/sh_sp3fit: Incorporated GNSS into decoding svs_exclude.dat.  Line format is now 
  Gnn, Rnn etc.  When gnss==G, old nn only designation allowed.  Note: doy is 
  3-digit.  TAH 200524 
 
com/sh_gnss_sel: Fixed issue with hardwired date in find command.  Updated 
  help message.  TAH 200524 
 
com/sh_preproc: Made lfile test name more unique by adding expt and doy to 
  temporary name.  (conflicts in parallel runs could lead to zero length 
  lfiles.) TAH 200523. 

Thursday, May 21, 2020 
 
kf/ctogobs/ctogobs_dd.f: Fixed error in replaced integer do loop. TAH 200520. 
 
tables/otl_FES2014bM_MRG.list: Merged FES2014 list from otl_FES2014bM_IGS.list 
  and otl_FES2014bM.list generated at Chalmers University. TAH 200520 
 
kf/track/{theory.f, get_miss.f, track_com.h}:  Modified behavior of USE_BLQ 
  command to allow solid-Earth tide model to be specified (IERS2010 or default 
  analytic model) and to turn off model completely (NOETIDE).   NOBLQ needs to 
  be given as argument to turn off ocean-tide model  (structured this way for 
  backwards compatibility). TAH 200520 
 

Wednesday 20 May 2020

com/sh_get_grids: New script from S. McClusky to download grids from
  ftp://everest.mit.edu to the current working directory. MAF 20200520

kf/update_pmu/igs_to_pmu.f: Added feature to have UT1 in the PMU file be
  generated by integrating the LOD values rather than using the globk estimates.
  This can be useful when generating Earth fixed SP3 files from inertial t-
  files. TAH 20200518

com/sh_gnss_sel: New script to add GNSS networks to sites.default.yyyy.ddd
  files. Assumes that the rinex files have been downloaded so that data types
  available in the files can be determined with scan_rinex. TAH 20200518

com/sh_plotvel: Added coast line specification -clw option. Default remains
  unchanged. TAH 20200515

com/sh_network_sel: Added feature to update existing sites.defaults file with
  additional net names added. Removed default apr_file name. Useful for having
  GPS and GNSS networks separate. TAH 20200515

com/sh_rx2apr: Forced rinex header coordinate extraction  if RINEX 3 files which
  svpos and svdiff can't read.  Cleaned up site name extraction when directory
  included in site file name. Cleaned up indentation. Remove spherical l-file
  output. TAH 20200518

com/sh_rx2apr: Added RINEX 3 trap to use approximate coordinates in RINEX header
  if file given with "-site" option is RINEX 3 or results of SVPOS if file given
  with "-site" option is RINEX 2 but file given with "-ref" option is RINEX 3.
  MAF 20200517

tables/svs_exclude.dat.nanu: Updated following USCG NANU 2020023. MAF 20200515

tables/rcvant.dat: Added STONEX S980 (STXS98), STONEX S980A (STX98A) and STONEX
  S990A (STX99A) receivers, and STXS990A (STX99A) antenna. MAF 20200514

Summary of added features and updates. Tuesday, May 12, 2020
- The Gipson 2017 VLBI model updated. Default: added UT1 libration term and
  keeps the IERS2010 model.
- New scan_rinex utility added assess observables in RINEX 2 and RINEX 3 files.
- Added pre-fit millisecond clock jump scanning and fixing in autcln. Normally
  not needed and default is not to do it. New status report option to look at
  pre-scan clock residuals. Useful if lots of BAD PREFIT data delete messages.
- Added lower frequency selection for GNSS processing (makex and sh_gamit -lfreq
  option) so G5, C7, E6/7/8 are options of lower frequency instead of G2, C6 and
  E5.
- Updated autcln to avoid errors in autcln/solve when some stations are missing
  a frequency band completely.

libraries/includes/freq_def.h: Fixed small error in E8 frequency (0.001E6 Hz)
  (impact MW-WL looked like LG in cview. TAH 20200512

gamit/makex/settyp.f: Added E8 to list of valid frequency. TAH 20200512

gamit/includes/model.h: Increased atxfrq to 4-characters and additional entry
  for primary and secondary choice of ANTEX frequency. TAH 20200512

gamit/utils/convert_antpcv.f: Updated read_antex call for new atxfrq definition.
  TAH 20200512

gamit/model/setup.f: Added selection of alternative lower frequency.
  TAH 20200512

gamit/model/get_svantpcv.f: Updated to allow primary and secondary choice of
  frequency read from ANTEX file. TAH 20200512

gamit/model/get_antpcv.f: Updated to allow primary and secondary frequency
  choice of frequency read from ANTEX file. TAH 20200512

gamit/model/get_antinfo.f: Modified format for reporting ANTEX frequencies used.
  A * next entry shows if primary or secondary frequency used for PCO/PCV.
  TAH 20200512

gamit/lib/sel_obtyp.f: Added additional data selections based choice of lower
  frequency. Added entries for G5, C7, E6/7/8; defaults are G2,C6 and E5.
  TAH 20200512

gamit/lib/read_antex.f: Passed modified atxfrq array now with 4 characters and
  three entries: (1) L1, (2) primary choice for L2, (3) secondary choice for L2
  where L2 if the lower frequency selection. TAH 20200512

tables/ftp_info: Added GFZ Potsdam as source for MGEX orbits. Currently SP3
  files need to be downloaded before sh_get_orbits is run. TAH 20200511

kf/ctogobs/{ctog_utils,read_cfheads}.f: Updated save non-zero frequencies when
  c-files are read. Solve has problems if the first c-file is missing a
  frequency band and autcln could have problem if last c-file is missing
  frequencies. This can happen with GNSS data missing standard frequencies at
  some sites. TAH 20200511

kf/ctogobs/set_freqs.f: Moved code that saves frequencies to ctog_utils.f and
  read_cfheads.f. TAH 20200511

kf/ctogobs/upd_cf_sd.f: Added code to write both frequencies bands to all c-
  files even if one band is missing (all data would have been deleted).
  TAH 20200511

gamit/makex/makex.f: Modified gnss input so that a lower frequency selection can
  be passed. For example C7 for the gnss runstring parameter will select C7 as
  lower frequency instead of default C6.  G5, E6, E7 and E8 are also options.
  TAH 20200511

gamit/lib/sel_obtyp.f: Added option for lower frequency selection. TAH 20200511

gamit/makex/rbatch.f: Reading of GNSS as 2-character with lower frequency option
  added. TAH 20200511

gamit/includes/global.h: Added gnsslf which contain GNSS and lower frequency
  selection (e.g., G or G5 for L5 lower frequency. TAH 20200511

gamit/makexp/makexp.f: Added passing gnsslf string into makex batch file.
  TAH 20200511

com/sh_makexp: Added -lfreq option to select lower frequency for G,C and E, e.g.
  "-gnss E -lfreq 7" will select L7 (E5b) as lower frequency for Galileo instead
  of L5 (E5a). Options allowed are G5, C7, E6/7/8. TAH 20200511

com/sh_preproc: Added -lfreq option (see sh_makexp). Fixed bug with g-file being
  checked with "^G" rather than "^$gnss". TAH 20200511

com/sh_gamit: Added -lfreq option (see sh_makexp); -remakex Y will be needed if
  older x-files already exist in processing directory. TAH 20200511

kf/ctogobs/prc_ctog_cmd.f: Added new option to prefit_clk command to remove
  millisecond jumps in clocks. Not a common problem but could be useful for some
  data. TAH 20200510

kf/ctogobs/ctog_utils.f: Initialization of prescan_msj option. (Default no
  scan). TAH 20200510

kf/ctogobs/prescan_clk.f: Added some commented out write statements to print
  information on clock fits. TAH 20200510

kf/ctogobs/ctogobs.f: New call to prescan_msj. TAH 20200510

kf/ctogobs/prescan_msj.f: New subroutine to scan and remove millisecond jumps.
  TAH 20200510

kf/ctogobs/ctogobs_com.h: Now logical added to implement prescan_msj. New
  version number for autcln. TAH 20200510

kf/ctogobs/fit_igs_clk.f: Limited updating linear fits to clocks with < 100 ns
  RMS. Millisecond jumps in clocks can cause high RMS scatters. Corrected
  indexing on frequencies (not issue since all frequencies are mapped to common
  frequency at this point). TAH 20200510

kf/ctogobs/upd_cf_sd.f: Removed a debug statement. TAH 20200510

kf/ctogobs/Makefile.generic: Added prescan_msj subroutine. TAH 20200510

gamit/utils/scan_rinex.f: New routine ported from GG11 to scan RINEX files to
  determine the numbers of observations of different types. Modified to add
  leading zero to PRNs<10 in output and generated a sorted output.
  RWK/TAH 20200509

gamit/includes/dimpar.h: Added max_ob11, maxgnss, max_sch for GG11 for
  scan_rinex program. TAH 20200509

gamit/utils/rinex.h: New include file from GG11 for use in scan_rinex.
  TAH 20200509

gamit/utils/Makefile.generic: Added scan_rinex to Makefile. TAH 20200509

com/link_gamit_bin: Added ../utils/scan_rinex. TAH 20200509

gamit/lib/{srotat,rotsnp_sofa}.f: Added module for hpeop_gipson and call to
  calc_hf_gip_xyu. This code is implementation of Gipson 2017 VLBI model. Both
  of these routines have same block of code for HF EOP but some time argument
  variable names differ. TAH 20200509

gamit/lib/hfeop_gipson.f90: New module implementing Gipson VLBI model. Libration
  and solid Earth terms were in a priori so this model needs these terms added
  and is compatible with Desai and Sibios implementation. TAH 20200509

gamit/lib/hfeop_desai.f90: Modified to remove functions and subroutines that are
  common to hfeop_gipson.f90. TAH 20200509

gamit/lib/hfeop_util.f90: Functions and subroutines common to the two HF EOP
  modules. TAH 20200509

gamit/lib/Makefile.generic: Modified to add new hfeop_gipson.f90 and
  hfeop_util.f90 modules and subroutines. TAH 20200509

tables/svs_exclude.dat.nanu: Updated following USCG NANU 2020020. MAF 20200508

Summary features added Tuesday 5 May 2020:
- Desai and Sibios HF EOP Model added, UT1 libration set in defaults.
- Fixed -autooff in htoglb (-a was OK).  Updated documentation to reflect that
  this option should be used for subnetworks processed in baseline mode.
- Fixed reading of FES2014 list and ocean tide potential models.  Critical
  update if these models used.
- Fixed reading of dates near the ends of years in UNR time series files.
- New MJD->date option in cplotx (field type 3).
- Added secondary list selection option in sh_network_sel (allows a primary list
  that should all be selected if they are available and their total number is
  less than the total network size.

tables/{sestbl.10.71,sestbl.10.71.MIT}: Set Etide model to IERS10 and added UT1
  Libration to Earth Rotation model. TAH 20200505

kf/{ctogobs/ctogobs_dd,svpos/spcsnr}.f: Replaced floating point do-loops with
  integer code. TAH 20200505

kf/htoglb/cresnx_file.f: Updated model names output to SINEX file. TAH 20200505

kf/gen_util/gamit_mod_n.f: Updated model names for UT1-Libration (UT1-LIBR) and
  Desai and Sibois (D&S2016). TAH 20200505

gamit/model/setup.f: Updated output documentation to reflect Desai and Sibois
  model, and the string assigned to speopmod character string. TAH 20200505

gamit/fixdrv/mdmake.f: Updated decoding of "Earth Rotation" sestbl. entry.
  TAH 20200505

gamit/orbits/{get_svcoords,partl,sp3tot,trot}.f: Replaced arbitrary setting of
  iut1pol bit mapped diurnal/semidiurnal models with value read from sestbl.
  TAH 20200505

gamit/lib/rot_gfile.f: Replaced setting iut1pol to zero with value read from
  sestbl. TAH 20200505

gamit/lib/{rotsnp_sofa.f,srotat.f}: Implemented new Desai and Sibois diurnal/
  semidiurnal ocean tide EOP model. TAH 20200505

gamit/lib/Makefile.generic: Added hfeop_desai.f90 and get_iut1pol.f routines.
  TAH 20200505

gamit/lib/get_iut1pol.f: New routine to get iut1pol value from sestbl. Called
  rather than arbitrary set values in different routines. TAH 20200505

gamit/lib/hfeop_desai.f90: New module to implement the Desai and Sibois model.
  TAH 20200505

gamit/arc/{earthradTUM.f,sbfn.f}: Replaced arbitrary setting of iut1pol bit-
  mapped diurnal/semidiurnal models with value read from sestbl. TAH 20200505

gamit/arc/satprop.f: Modified to remove repeated Block IIIA not modeling message
  when NCLE1 albedo model used. TAH 20200504

kf/htoglb/htoglb.f: Fixed decoding -a/-autooff option so both forms work. Only
  -a worked in earlier version. This would cause problems when baseline mode
  networks are combined. TAH 20200504

gamit/makex/{makek,makex}.f: Fix misnamed variable causing MAKEK to fail; fix
  inconsistent dimensions in both routines, with no effect; fix printing of the
  observables to the infor file, also with no effect. RWK 20200504

tables/rcvant.dat: Removed comment ("?") from beginning of records for several
  receivers and antennas without IGS codes but previously added to rcvant.dat by
  MIT. MAF 20200504

com/{links.arc,links.day,links.com}: Added test to only link otides.dat if no
  valid local file. Allows use of different ocean tide models (FES2004 versus
  FES2014b). TAH 20200503

gamit/arc/read_otides.f: Fixed logic in reading files with more than 18
  frequencies (FES2014b) and added output of model name to GAMIT.status, cleaned
  up debug output. Critical update when using FES2014. TAH 20200503

gamit/grdtab/get_otl_list.f: Updated reading of TANG RADI lines for FES2014
  format to ensure sites found (no sites found with old code). Critical update
  when using FES2014. TAH 20200502

tables/otl_FES2014bM_IGS.list: Updated model name to FES2014b from FES2014M,
  which had been in the file. TAH 20200502

com/sh_xtosess: Updated definition of "prn_list" to accommodate constellation
  prefix (G, R, E, C, etc.) to PRN number. RWK/MAF 20200430

kf/glinit/get_aprioris.f: Fixed mapping of ECOM-1 radiation parameters for old
  (pre 1.07) binary h-files. Impact when orbits are tightly constrained in RELAX
  solutions. TAH 20200429

kf/glist/glist.f: Increased number of digits in reported h-file version.
  TAH 20200429

tables/frames.dat: Added Caribbean ("CARIB18") and Puerto Rico-Virgin Islands
  ("PRVI18") reference frames relative to IGS14 from Wang et al. (2019),
  doi:10.3390/rs11060680. MAF 20200428

com/{sh_gamit,sh_make_sky_pngs}: Added .png extension to GMT's
  ps2raster/psconvert -F option and rename if GMT creates PNG file with double-
  extension ".png.png"). Y.-H. Chang/J. A. Garmen/MAF 20200426

kf/globk/{add_GGV.f,read_glb_mar.f,Makefile.generic: New routine, add_GGV, to
  add GGVersion information to the GLOBK description line at the end of
  read_glb_mar to allow better assessment of version of GG installed.
  TAH 20200426

help/htoglb.hlp: Expanded the documentation about correct procedures for
  combining sub-networks when baseline mode has been used in GAMIT processing.
  TAH 20200425.

kf/blsum/tscon.f: Updated reading dates from UNR txyz2 files. Old algorithm
  could have errors of a day near the end of the year and assign values to the
  wrong year. TAH 20200424

tables/svs_exclude.dat.nanu: Updated following USCG NANU 2020017. MAF 20200424

com/sh_network_sel: Added a secondary list of stations so that the core list can
  contain critical stations (the number of stations can be less than the total
  network size). After selection of all core stations, the secondary list is
  used for the priority of these stations. Added a -clear_list option to force
  the re-download of archive data holdings. TAH 20200422

Kf/utils/global_sel.f: Updated to handle non-"_GPS" site names in a priori file,
  changed reporting of short lists. TAH 20200422

help/head.snx: Updated DOMES codes. TAH 20200420

tables/rcvant.dat: Added SEPVC6150L (SEPVC6) and SEPVP6050_CONE (SEPVP6)
  antennas, and CHC P5E (CHCP5E) receiver. MAF 20200420

kf/utils/make_rms.f: Changed F6.1 to F6.2 format when writing statistics.
  TAH 20200417

kf/htoglb/mstinf.f: Changed dimensioning so that max_site sites can be in a
  station.info file to be merged with another one. E. D'Anastasio/TAH 20200417

tables/: Added igs14_2101.atx, ngs14_2101.atx and ngs14_2101_plus.atx, and
  linked antmod.dat to ngs14_2101_plus.atx; added igs_metadata_2101.snx and
  igs_metadata_2101_plus.snx, and linked igs_metadata.snx to
  igs_metadata_2101_plus.snx; updated svnav.dat.allgnss following
  igs_metadata_2101.snx. MAF 20200416

gamit/orbits/yawtab.f: Added one day to date for call to svnav_read because the
  start of the t-file may occur before the start time of a PRN number record.
  A. Koulali/R. Pickle/MAF 20200414

com/{sh_dem,sh_plotvel}: Updated to write new .grd-files to current working
  directory, rather than ~/gg/maps/; corrected definition of maximum longitude
  in "minmax" variable when passing topo_file name from sh_dem to grdimage for
  plotting, and updated logic to derive an illumination .grd-file from topo_file
  by requiring "topo_file" to be passed as the argument to the "-illu_file"
  option. MAF 20200413

com/sh_dem: Removed extraneous "-fg" option from GMT's grdcut command for
  SRTM3, SRTM30 and SRTM30_PLUS. S. Nistor/MAF 20200403

tables/sites.defaults: Appended to updates/source/tables.10.71.tar.gz, because
  it was missing from the original release, and incremental_updates/tables/.
  R. Romero Andrade/RWK/MAF 20200403

com/sh_sopac: Removed leading slash ("/") in definition of "sdir" variable to
  avoid double-slash ("//") typo in host URL. R. Pickle/MAF 20200402

com/links.arc: Changed links so that local ut1./pole. file will be retained
  during linking. TAH 20200401

com/sh_sp3fit: Changed logic to all allow multi-day orbit arcs of more than 3
  days (orbits modified to allow 9-day arcs). TAH 20200401

gamit/arc/read_arc_batch.f: Updated logic for reading arc.out file name and the
  hidden debug featutes that are included on that line. TAH 20200401

gamit/includes/orbits.h: Increased dimensioning to allow 9-day arcs.
  TAH 20200401

gamit/orbits/{orbdif,orbfit,orbrms,plt_postfit}.f: Modified outputs to add MJD
  to ends of residual values so that true times can be plotted. TAH 20200401

kf/plot/{dxaxis_mn,dxaxis_mx,get_xy_scale,make_line,plot_spec,report_scale}.f:
  Added new data type field to have MJD interpreted as date rather than
  converted to decimal year (help/{cplotx.hlp,ts_pos.plt,ts_pos2.plt} also
  updated accordingly). TAH 20200401

tables/nbody: Updated to use DE405 (little-endian) file for 1900-2050 from JPL
  (ftp://ssd.jpl.nasa.gov/pub/eph/planets/Linux/de405/lnx1900.405). MAF 20200401

tables/frames.dat: Added rotation rate vector for North America relative to
  IGS08 from Kreemer et al. (2018), doi:10.1029/2017JB015257. MAF 20200331

tables/hi.dat: Added DHTCR (= DHARP) height code for SEN671575-series antennas.
  P. Fang/Y. Bock/MAF 20200326

tables/rcvant.dat: Added STONEX S700A (STXS7A), STONEX S850A (STXS85) receivers,
  and STXS700A (STXS7A) antenna, following 2020-03-25 release of IGS
  rcvr_ant.tab. MAF 20200325

test_install/: Added missing "-netext" option to first sh_glred commands for
  2016 and 2014 in README; changed tssum product ID to conventional form of
  "mit.final_igs14" from "mit.final.igs14" in README and runtest. MAF 20200323

tables/: Updated igs_metadata_2097.snx and igs_metadata_2097_plus.snx; updated
  svnav.dat.allgnss following igs_metadata_2097.snx; updated eq.def with
  preliminary definitions for July 2019 Ridgecrest, California, earthquakes.
  MAF 20200323

tables/dcb.dat: Corrected open-ended end date for G037/G07 to 1993-334T23:59
  and removed erroneous extra line for G063/G01 (2019-305T00:00 to
  2019-304T23:59). S. Hreinsdottir/MAF 20200323

test_install/README: Corrected erroneous "-gnss" option in one sh_glred command
  to "-netext". MAF 20200319

tables/: Added finite end date for G034/G18 (2020-070) and dummy first
  record for G075/G18 (from 2020-073), until first estimates of DCB are
  available, in dcb.dat; updated igs_metadata_2097.snx and
  igs_metadata_2097_plus.snx to versions with corrected start date of G075/G18;
  updated svs_exclude.dat.nanu with NANU 2020014. MAF 20200319

tables/: Added igs14_2097.atx, ngs14_2097.atx and ngs14_2097_plus.atx, and
  linked antmod.dat to ngs14_2097_plus.atx; added igs_metadata_2097.snx and
  igs_metadata_2097_plus.snx, and linked igs_metadata.snx to
  igs_metadata_2097_plus.snx; updated svnav.dat.allgnss following
  igs_metadata_2097.snx. MAF 20200317

gamit/lib/rstnfo.f, kf/htoglb/{hfupd_ut.f,mstinf.h}: Changed length of comment
  to 36 characters from 30 (rstnfo.f), 24 (hfupd_ut.f) and 132 (mstinf.h)
  characters. TAH 20200317

com/{sh_plot_rinex,sh_rename_rinex}: Corrected calls of "CRX2RNX" (the original
  program name) to "crx2rnx" (our distribution name, in gamit/bin/).
  MAF 20200311

tables/rcvant.dat: Updated "PCN"-code based on old rcvant.dat.200109 for
  ROGUE SNR-8000 (P), TPS ODYSSEY_I (N), LEICA GS15 (C) and SEPT POLARXS (N)
  receivers. MAF 20200310

12345678901234567890123456789012345678901234567890123456789012345678901234567890