TOC  |      | Gpscom | Pages | Next Section | INDEX 

snxdif.html - 98Jul06

Snxdif

A program for comparing the site coordinates and EOP parameters in two Sinex files. Pages

Overview

The program Snxdif will read two Sinex files and determine by least squares adjustment a seven parameter transformation between the two sets of site coordinates. This program was written to be run interactively and prompts the user for input file names and options to define exactly how the comparison is to be done. However in actual practice, when the program is being run frequently, it is often usefull to prepare an input file with the answers to the questions and to feed that file to the program to avoid the tedium of answering the prompts. If earth orientation (EOP) data is in the file the program will produce a table of differences. If velocity data is available in the reference file the program will use that data to shift the coordinates to a common epoch. The program will prompt the user for a time to which to shift the coordinates of the reference data set (the second one specified). If the responce to this prompt is blanks or a carriage return the default will be to for the program to pick the mean time of parameters in the first data set. This program will time shift the coordinates of the second data set to the target time. Currently it is assumed that the times for the parameters in the first data set span a time which is very small or zero as compared to the velocities and thus do not have to be time shifted. Should this change the actions taken for the second data set must then be applied to the first data set as well. If no velocity data is available the program will simply compare coordinates with out regard to any time difference between the data in the Sinex files.

The current rules in Snxdif for matching site names from one Sinex file to another is that the name and the point code (the letter designation) must match. The sequence number is included in the match as a secondary condition because it is an arbitrary number under Sinex rules which is defined internally to the particular Sinex file. So Sinex files for the ITRF or any other agency would never have matches to those of NGS if the sequence number was included in the match. In the case of NGS Sinex files different sequence codes do have meaning and we could have more than one set of sequence codes in the same adjustment if the adjustment included a large time span. Thus following rules have been adopted. When the program is matching site names for two different files and it finds that the site name and point code matches it then searches to find a site where the sequence number also matches if no sequence number matches it will accept the first occurring name match.

This program was written in FORTRAN. Predominantly ANSI standard FORTRAN, a few isolated routines query the operating system for time and date, user and system identification. These "system" routines are inherently non-portable but, should be easily modified or disabled.

Input Data

The program normally requires two Sinex files as input, an output file name, specification of certian options and a list of site nemes that are to be included, or excluded depending on the option selected, in the adjustment for the seven parameter transformation. The following description descripes the prompts and responces for a typical comparison. As of the current writeing the following is not exaustive of all possible ways to run the program but should satisfy most needs.
 Program prompt -> TYPE SINEX FILE NAME  _ 
       responce ->  	     (give 1st file name)

 Program prompt -> READ ANOTHER DATA SET  Y/N/A _
       responce -> y
                      (if above answer is 'n' a TEST run is assumed
		       and alternate questions will follow.

 Program prompt -> TYPE SINEX FILE NAME  _ 
       responce ->       (give reference file name)

 Program prompt -> TYPE PRINT FILE NAME  _ 
       responce ->       (give print file name)

 Program prompt -> TYPE (freefield) Time for matching data sets _
       responce ->       (carriage return, or date and time)
                     date/time forms = yy mm dd hh mm ss.ss
                                     = yy doy hh mm ss.ss
                                     = yy mm dd
                                     = yy doy
                                     = mjd.ddddd
		     If this line is commented, comments must follow
		     column 32.

 Program prompt -> Rotate, Translate & Scale coordinates to match ?
       responce -> y
		     (if above answer is 'n' addional prompts will be
		      given to query the user to individually specify
		      whether to solve for rotation, translation or
		      scale. Experiment!)




         Before the next prompt a list of site names will be printed.


 Program prompt -> The sites selected from this list may be either
 Program prompt -> excluded or included in the adjustment. 
 Program prompt -> Will you specify an exclusion list? 
       responce -> y
                      (answering 'y' as shown means all site names
		       specified will be excluded from the adjustment
		       for the seven parameter transformation. The
		       program now requires the user to list site
		       names one per line and a blank line will end the
		       list. If only a blank line is entered then all
		       sites will be included in the adjustment.
		       A responce of 'n' will reverse the sence of the
		       list causing only listed sites to be included
		       in the adjustment.)


 Program prompt -> TYPE THE  NAMES OF SITES TO BE EXCLUDED(Included) 
 Program prompt -> FROM THE ADJUSTMENT FOR ROTATION & TRANSLATION 
 Program prompt -> A BLANK OR (CR) WILL END THE LIST 
 Program prompt -> 
 Program prompt -> THE FOLLOWING SITES WILL BE Excluded
 Program prompt -> Exclude SITE NAMED   _
       responce -> ALGO A
       responce ->                    (blanks or carriage return)

		      (the above 2 responces would cause the site
		       ALGO to be excluded from the seven parameter
		       adjustment.)

 Program prompt -> THE CENTER OF Adjustment CAN BE ABOUT THE ORIGIN,  '
 Program prompt -> the CENTROID , or ABOUT A SITE SELECTED HERE. '
 Program prompt -> TYPE BLANK OR (CR) FOR THE ORIGIN, ELSE SITE = _' 
       responce ->                    (blanks or carriage return)

                     (the above prompt and responce allow the user to
		      shift the origin of the coordinate system to the
		      centoid of the sites of to some particular site
		      before the seven parameter adjustment is done.
		      This is not generally used.)

When several comparisons of Sinex files are being made together it gets tiresome to answer the prompts when running the program. In this case the user can make a file containing the answers to the prompts which will be given. Then if the input file were named "snxdif.inp" the program could be run the as in the following example.

snxdif < snxdif.inp > junk.purge

In this case a file such as the examples below would be prepared. The comments shown could be included without harm. However on lines where a file name was being given trailing comments would cause problems.

Example input file ,this example would cause all sites to be
included in the seven parameter adjustment.

--------------------------------------------------------
gpscom.snx
y                        read another data set ?  
/gyp5/ITRF96.SNX.short
test.output
                                 Time for Matching
y                        Rotate, translate & scale?       
y                        specify exclusion list ?        
 
   
 

-----------------------------------------------------

Example input file ,this example would cause only the named
sites to be included in the seven parameter adjustment.
--------------------------------------------------------
gpscom.snx
y                        read another data set ?  
/gyp5/ITRF96.SNX.short
test.output
                                 Time for Matching
y                        Rotate, translate & scale?       
n                        specify exclusion list ?        
KOKB A  
ALGO A  
MDO1 A  
TROM A  
KOSG A  
WTZR A  
MATE A  
YELL A  
FAIR A  
GOLD B  
YAR1 A  
HART A  
MADR A  
 
   
 

-----------------------------------------------------
 TOC  |      | Gpscom | Pages | Next Section | INDEX 


snxdif.html
March 4, 1999
Bill Dillinger