r'.!^
' ^y
NAV. -
.- SCHOOL
MOWTEP.tY,
, CALli-.
3J940
NAVAL POSTGRADUATE SCHOOL
Monterey, California
THESIS
APPLICATION OF ADDITIONAL SECONDARY
FACTORS TO LORAN-C POSITIONS
FOR
HYDROGRAPHIC OPERATIONS
by
Gerald Eo Wheat on
October 1982
Thesis
Advisor: Gers
.Id Bo Mills
Approved for public release; distribution unlimited
T2C8076
SCCUNtTV CUASSiriCATIOM OF THIS ^A<JC (Whtt Dmtm Btttmr*^)
REPORT DOCUMENTATION PAGE
3. OOVT ACCCSSIOM NO.
4. TITLE fanrf Su*(lffa)
Application of Additional Secondary Factors
Co LORAN-C Positions for Hydrographic
Operations
7. AuTMOni'«>
Gerald E. Wheaton
S. PCArOMMINC OnOANlZATION NAME ANO AOOKCtS
Naval Postgraduate School
Monterey, California 93940
n CONTWOLLIMC O^riCe NAMC AMO AOOKCSS
Naval Postgraduate School
Monterey, California
TT MONlToniNC ACCNCY NAMC * AOOnCSSff/ ^IH»rwnt Irom Caturotllng OlUem)
Naval Postgraduate School
Monterey, California 93940
READ INSTRUCTTDNS
BEFORE COMPLETTNO FORM
1 RECl^ltNT'S CATALOG NUMBEI
5. TYPE OF «CPO»T « PCmoO COVERCO
Master's Thesis
October 1982
«. PCWFOWMING ORG. HEPOWT Nu
Maen
•• CONTNACT OI» GHAnT NLMSERraj
"*• tSS^VA^^''*'*'^ '•'•OJtCT TASK
AHCA • WOKK UNIT NUMBEHS
12 (lEPOUT OATE
October 1982
*S NUMBER OF PAGES
129
IS. SeCUniTY CLASS, (ot thl» r*>artj
Unclassified
ISa. OeCLASSlFICATION/ OOWNSRAOINC
SCHEDULE
16. OISTRISUTION STATCMCNT (ol tMi R*p»H)
Approved for public release; distribution unlimited
17. OISTKISUTION STATEMENT (ol th» atattmct anlararf In Stack 30, II dlllmrmnt tram Rmport)
l«. SUPPLEMENTARY NOTES
19. KEY WORDS (Coftllnu* art rmirmf tld» II n*c««aarr an« l^mnltlt »V MocA nuai**r>
LORAN-C; Hydrographic Surveying; Additional Secondary Factors; Calibration;
20. ABSTRACT (Conllm— on ravwa* »ldm II naeaaaarr «•' ItfaMKfr *7 Moe* mmt^ar)
The application of LORAN-C in the hyperbolic mode as a positioning system
for hydrographic surveys was investigated. Observed LORAN-C time
differences from a field test conducted in Monterey Bay, California
were compared to calculated time differences determined from geographic
positions based on a microwave positioning system. Four methods
were used to determine the calculated time differences. The first
three methods were (1) applying only the seawater Secondary Factor,
DO I JAM 7j 1473 EDITION or 1 NOV •• IS OBSOLfTE
S/N 0 102-014- ««0I
1
IKCURITY CLASSIFICATION OF THIS RAOB (9hmn Oaf Mnfrad)
<*eij««»'
i'"^«" n«t« *■••»•«
(2) computing the time difference based on a Semi-Empirical TD Grid,
and (3) applying ASF Correctors from the DMAHTC LORAN-C Correction
Table. The final method applied multiple observed ASF Correctors
at five minute latitude and longitude intervals. By applying multiple
observed ASF Correctors, which was the most accurate method, a 38.3
meter 1 drms with a lane offset of 3 to 12 meters using the 9940 X-Y
LORAN-C combination was obtained. Based upon the results presented,
it may be possible to use LORAN-C for hydrographic surveys at scales
of 1:80,000,
DD Form 1473
1 Jan 73
S/N 0102-014-6601
•tcu«t»v CLAMiriCATioH or THU PAairmt*)' o««« t«««»»*»
Approved for public release: distribution unlimited
Application of Additional Secondary Factors to LORAN-C
Positions for Hydrographic Operations
by
Gerald E. WheatDn
Lieutenant , NOAA
B.S., California State University, Humboldt, 1975
Submitted in partial fulfillment of the
requirements for the degree of
MASTER OF SCIENCE IN OCEANOGRAPHY (HYDROGRAPHY)
from the
NAVAL POSTGRADUATE SCHOOL
October 1982
y^ I U II
a, I
MONTEREY, CALIF. 93S40
ABSTRACT
The application of LORAN-C in the hyperbolic mode as a
positioning system for hydrcgraphic surveys was investi-
gated. Observed LORAN-C time differences from a field test
conducted in Monterey Bay, California were compared to
calculated time differences determined from geographic posi-
tions based on a microwave positioning system. Four methods
were used to determine th? calculated time differences. The
first three methods were (1) applying only che seawater
Secondary Factor, (2) computing the time difference based on
a Semi-Empirical TD Grid, and (3) applying ASF Correctors
from the DMAHTC LORAN-C Correction Table. The final method
applied multiple observed ASF Correctors at five minute
latitude and longitude intervals. By applying multiple
observed ASF Correctors, which was the most accurate method,
a 38.3 me-cer 1 drms with a lane offset of 3 to 12 me-ers
using the 9940 X-Y LORAN-C combination was obtained. Based
upon rhe results presented, it may be possible to use
LORAN-C for hydrographic surveys at scales of 1:80,000.
TABLE OF CONTENTS
I. INTRODUCTION 9
A. USE OF LORAN-C 9
B. APPLICATION OP ASF CORRECTORS FOR NAVIGATION . 11
C. APPLICATION OF ASF TO HYDROGRAPHIC POSITIONING 13
D. OBJECTIVES 14
II. NATURE OF THE PROBLEM 18
A. THE PRINCIPLES OF LORAN-C 16
B. PHASE LAG 19
C. TD MODEL 22
D. SEMI-EMPIRICAL TD GRID CALIBRATION MODEL ... 24
1. Technical Approach 25
2. Generalized Range/Bearing Model 27
E. DMAHTC MODEL 31
1. Sea SF Model 31
2. Land SF Model . . . .- 31
3. Table Descri p-^.ion 34
F. ATTEMPTED DETERMINATION OF ASF CORRECTORS 3Y
HYDRO FIELD PARTIES 38
III. EXPERIMENTAL PROCEDURE 42
A. FIELD PROCEDURES 42
B. MICROWAVE SYSTEM POSITIONING 45
C. LORAN-C COMPUTATIONS 46
1. Seawa-er Secondary Factors (SF) 47
2. Semi-Empirical TD Grid 49
3. Calculated Table ASF Correctors 51
4. Observed ASF Correctors 52
IV. RESULTS 54
A. SEAWATER SECONDARY FACTORS 55
B. SEMI -EMPIRICAL T D GR ID 57
C. TABLE ASF CORRECTORS 59
D. MULTIPLE OBSERVED CORRECTORS 60
V. CONCLUSIONS 63
APPENDIX A: PROGRAM GPBYLQ 65
APPENDIX B: PROGRAM LORAN 78
APPENDIX C: PROGRAM LOPLC 88
APPENDIX D: PROGRAM LORTAB 101
APPENDIX E: PROGRAM ASFSEL 112
BIBLIOGRAPHY 123
INITIAL DISTRIBUTION LIST 127
LIST OF TABLES
I. Phase Retardation or Lag of Radio Waves 10
II. Coding Delay 9940 LORAN-C Chain 24
III. TD Bias (b) - jisec 29
IV. Seawater Coefficisnts 30
V. Land Coefficients 31
VI. Sea SF Model Coefficients 32
VII. Geographic Names and Positions 42
VIII. Microwave Positioning Equipment 45
IX. LORAN-C 9940 Chain Data 47
X. Data Set Parameter 55
XI. Seawater Secondary Factor Errors 56
XII. Semi-Empirical TD Sr id Correction Errors 58
XIII. Multiple LORAN-C Correction Table Errors 59
XIV. Multiple Observed Correction Errors 61
LIST OF FIGURES
2.1 Hyperbolic Fix (From Coast Guard LORAN-C User
Handbook, 1974) 18
2.2 Transmitted Radiation 19
2.3 Phase Lag 21
2.4 Location of West Coast LORAN-C Stations .... 26
2.5 Mixed Path TD Geometry 28
2.6 Composite Land-Ssa Path 33
2.7 99aO-W ASF Correction Table 35
2-8 9940-X ASF Correction Table 36
2.9 ■ 9940-Y ASF Correction Table 37
3.1 Location of Test Areas in Monterey Bay 44
I. INTRODUCriON
A. USE OF LORAN-C
In recent years there has been an increasing usage of a
LORAN-C receiver and the LORAN-C network as rhe primary
horizontal control for such scientific studies as deep ocean
dumpsites, marine fisheries studies [Rulon, 1979], bathyme-
tric surveys, and recently, a recDnn aissance hydrographic
survey. Examples of bathymetric surveys conducted by the
National Oceanic and Atmospheric Administration (NOAA) ,
which have used Loran-C for positioning are:
1) Su-100-1-79 Gulf of Alaska [NOAA H-9822, 1979],
2) SU-100-2-79 Gulf of Alaska [NOAA H-9823, 1979],
3) S-D902-WH-82 U.S. West Coast [NOAA Ship Surveyor,
1982].
Recently an attempt was made ♦:o use LORAN-C as the
sounding position control for a reconnaissance hydrographic
survey S-K902-WH-82 [NOAA, 1982]. This was a special survey
conducted by the NOAA Ship Whiting in May 1982. Special
surveys are field examinations of very limited extent or
scope and frequently require unique survey or data collec-
tion procedures [Umbach, 1976]. The purpose of this project
was to verify the existence and extent of reported shoaling
in three safety fairways in the Gulf of Mexico [NOAA, 1982].
The use of Loran-C as a positioning system for basic
hydrographic surveys has been very limited due to the abso-
lute accuracy of the long range system. A basic hydro-
graphic survey is defined as a survey which is so complete
that it need not. be supplemented by other surveys. "It must
be adequate to supersede for charting purposes all prior
surveys" [Umbach, 1976]. Variables which affact the accu-
racy of LORAN-C are signal propagation variations, weather,
and sky waves. The affects of weather and sky waves on
LORAN-C propagation are best described by Samaddar [1980]
and the American Practical Navigator, [DMA, 1977] respec-
tively.
Signal propagation variations are due to the phase
retardation of the signal as it passes over an all seawater
path, over land paths, or partial seawater-land paths as
compared to free space. Table I sammarizes phase retarda-
tion changes [Mortimer, 1978]. Errors due to an all sea
TABLE I
Phase Retardation or Lag of Radio Waves
Propagation Path
ReDr=s entation
Propa gation
V9l ocity
(kffl/secf
Vacuum
299792.5
Direct wave through 299 691
earth's atmosphere
Ground wave over
sea water
Ground wave over
rugged mountains
299560
298 899
Difference in Phase
Lag at 500 km
Compared with
Wave in Line
Abcv5 (m)
170
220
1,300
water path are known as the Secondary Factor (SF) and errors
due to a land path or mixsd path are known as the Additional
Secondary Factor (ASF) [Speight, 1982].
10
ASF Corrections in ths LORAN-C system can b9 as large as
plus or minus four microseconds, which is 600 meters on t.he
baseline. In other areas with the same LORAN-C coverage,
these corrections may be much larger due to the expansion of
the distance between adjacent hyperbolic lines of position.
For example, at 32° N and 80° W, using lattice pair 9960-X,
a four microsecond (jisec) error will offset the 9960-X line
of position approximately 2438 meters [Speight, 1982].
B. APPLICATION OF ASF CORRECTORS FOR NAVIGATION
To compensate for the Loran-C positional errors caused
by the ASF Correctors, the Secretary of Transportation
tasked the Defense Mapping Agency aydrographic/Topographic
Center (DMAHTC) , the National Oceanic and Atmospheric
Administration (NOAA) , and the United States Coast Guard
(USCG) with the job of determining and applying the
Additional Secondary Correctors for each Loran-C chain.
This task was published in the Department of Transportation
(DOT) National Plan for Navigation in the July 19, 1974
Federal Register. These corrections should provide 95^
assurance that a vessel could fix its position to a pred-
icted accuracy of 1/4 nautical mile (NM) within the U. S.
Coastal Confluence Zone (CCZ) and the Great Lakes. The CCZ
is defined as:
"the seaward approaches to land, the inner boundary of
which IS the narbor entrance and the cuter boundary of
which is 50 nautical miles offshore or the edge of the
Continental Shelf (100 fathom contour) whichever is
greater. "
The 1/4-NM accuracy reguirement also affects the
nautical chart. The National Ocean Survey (NOS) , which
publishes charts for the CCZ, engaged in a program with the
USCG and DMAHTC to provide the coastal navigator with charts
overprinted with lattices which meet 1/4-nm accuracy. The
OSCG, as operator of the LORAN-C radionavigation system.
11
conducts surveys to ensure that LORAN-C coverage exists
within the CCZ and will be reponsible for the verification
of 1/U-NM accuracy for all coastal LORAN-C service. In
conjunction with NOS , it assists in surveys of coastal
waters of the United States to allow production of LORAN-C
charts based on observed field data to meet the standards
set forth above [Speight, 1982].
DMAHTC, for LORAN-C civil need, prepares grid predic-
tions from its data base. Based on analysis and verifica-
xion of the predicted grid from a USCG and/or NOS survey, it
produces revisions to the initial grid predictions [ Speight,
1982]. At present, DMAHTC has provided NOS with ASF
Corrected LORAN-C Lattices which are overprinted on the NOS
Charts. Each chart with ASF Correctors applied contains one
of the following notes:
"The LORAN-C lines of oosition overprinted en this chart
have been prepared for use with groudwave sianals and
are presently compensatsd only for t heoretical' prooaga-
tion delays, wnich have not yst been verified by
observed data. Marinsrs are cautioned not to rely
entirely on the lattices in inshore waters. Skywave
corrections are not provided".
or
"The LOEAN-C lines of position overprinted on this chart
have been prepared for use with groundwave signals and
are compensated with propagation delays cofflDu.ed from
observed data. Mariners are cautioned not to rely
entirely on the lattices in inshore waters. Skywave
corrections are not provided" [Speight, 1982].
Presently, all of the NOS Charts of 1:80,000 to
1:120,000 scale covering the east coast. Gulf coast, and
Great Lakes show LORAN-C lattices that have been compensated
for Additional Secondary Factors. Most of 'h= lattices on
these charts have been constructed from DMAHTC data tapes
that provide adjusted LORAN-C readings for each rate at
every five minutes of latitude and longitude. A few
lattices were constructed using a single ASF Correction for
the entire chart area. Five minute data tapes were not
12
furnished by DMAHTC for constructing lines of position for
LORAN-C rates on the West Coast Charts. On these charts a
single average ASF Correction was used to adjust each
lattice [NOAA, Marine Chart Division, 1982].
In addition to supplying corrected LORAN-C lattices for
nautical charts, DMAHTC prepares, distributes, and periodi-
cally updates unclassified ASF LORAN-C Correction Tables
[Speight, 1982]. The ASF Correction Tables are for preci-
sion navigation, utilizing digital computers to convert
LORAN-C time differences to geographic coordinates [Speight,
198 2]. Presently, the ASF correctors found in the LOEAN-C
Correction Tables were determined using theoretical propaga-
tion delays. ASF correctors listed in the tables are going
to be updated with observed data and reprinted the first
quarter of 1983 [Wallace, 1982].
C. APPLICATION OF ASF TO HYDROGRAPHIC POSITIONING
Schr.ebele [1979] investigated the possibility of using
Loran-C as an electronic positioning system for hydrographic
surveying. He concluded that in Monterey Bay, California a
single Additional Secondary Factor (ASF) applied to offshore
lines of position gave a root mean sguare error (drms) of 66
meters for the West Coast 99U0 Y-W pair and a predicted U2
meter drms error for 9943 X-Y rates.
The 42 meter predicted drms is larger than Nelson's
[General Electric Co., 1979] findings in San Francisco Bay.
He demonstrated, in a dynamic mode, that the precision of
LORAN-C was 60.8 meters 2 drms (30.4 meter 1 drms) with a
worst case of 71.2 meters 2 drms (35.6 meter 1 drms) . A
mean difference or offset between the measured time differ-
ence and the calculated time difference was 34 nanoseconds
for the 9940-X rate and one nanosecond for the 9940-Y rate.
13
He also obtained a precision of 38.0 meters 2 drms (19.0
meters 1 drms) in the static mode. Nelson also states, that
the above precision is only achievable if the user has a
LORAN-C receiver which has the performance capabilities of
those used in the experiment. rh9 LORAN-C receiver must
have "comparable signal averaging time, extra notch filters,
and attenua-ion of the signal" [General Electric Co., 1979].
D. OBJECTIVES
The National Ocean Survey requires that hyperbolic
control systems used for hydrographic surveying exhibit a 1
drms of less than 0-5 millimeter at the scale of the survey
[Umbach, 1976]. Although this requirement is generally for
2 mHz phase comparision systems, it can be inferred that it
also applies to other hyperblic systems such as LORAN-C.
The scale routinely used" for coastal surveys is between
1:40,000 and 1:80,000 [Umbach, 1976] yielding an allowable
error of 20 to 40 meters not including systematic errors.
Schnebele [1979] concluded that hyperbolic LORAN-C, after
applying a single ASF Corrector, is unsuitable for basic
hydrographic surveying.
Whether or not applying multiple Additional Secondary
Factors (AS?) to LORAN-C lines of position will reduce the
drms sufficiently to meet the accuracy standards set by the
National Ocean Survey Hydrographic llanual will be ascer-
tained in this study. The term multiple ASF Correctors
refers to the applicatioQ of more than one corrector to
LORAN-C lines of position over a given area. The variable
ASF Correctors result from varying delays of the electromag-
netic wave as it propagates over different land segments.
Three methods of applying multiple ASF Correctors were
tested. The first method was the application of a
14
Semi-Empirical Time Difference Grid Calibration Model devel-
oped by The Analytic Science Corporation [1979]. The spon-
soring agency was the Onited States Coast Guard. The
Semi- Empirical Model applies Secondary Factors and
Additional Secondary Factors for each geodetic position
based on the distance over land, the distance over water,
and the total distance using mean se^ water and land conduc-
tivities.
The second method which was investigated applies ASF
Correctors found in the DMAHTC LORAN-C Correction Tables
[DMAHTC, 1981] to LORAN-C lines of positions. These correc-
tors were derived from the ground conductivities which have
been determined in the field by a Coast Guard calibration
team [U.S. Naval Oceanograp hie Office, 1982].
Finally, a third method was pursued. ASF Correctors,
which were determined by field observations, were applied to
the LORAN-C lines of position. These were derer mined by
computing the difference between the observed LORAN-C rates
and the expected time difference which was calculated using
four lines of position from a very accurate microwave posi-
tioning system. These ASF Correctors were determined at
five minute latitude and longitude intervals.
15
II. MIHRJ Ql IHI EIQBLEM
A. THE PRINCIPLES OF LORAN-C
To understand the problems associated with LORAN-C when
used during hydrcgraphic operations, one must first under-
stand its principles of operation. LORAN-C is a low
frequency, pulsed signal, hyperbolic, radio navigation
system, employing time difference measurements of signals
received by the navigator from an laast three ground tran-
smitting stations [Speight, 1982]. The stations ar9
comprised of a master transmitting station, two or more
secondary transmitting stations which are strategically
spaced several hundred miles apart and, if necessary, a
System Area Moniter (SAM) Station [U.S. Coast Guard, 197U].
System Area Monitor (SAM) stations associated rfith each
LORAN-C chain apply differential-type corrections to the
rates in real-time. SAM stations continuously monitor the
signals from all transmitters in the chain. If the oberved
time difference deviates by- more than 0.05 jusec from the
expected value, then the appropriate secondary adjusts its
emission delay time in order to remove the error [Schnebele,
1 97 9 ] .
The master and at least two secondary stations are
located such that the signals from the transmitting stations
can be received throughout the desired coverage area. The
master station is designated by the letter "M" and the
secondary stations or slave stations are aesignatisd W, X,
Y, or Z [.U.S. Coast Guard, 1974].
All stations transmit on the common frequency of 100
kHz. Interference between transmitters is avoided through
the use of time separation [Poppe, 1982]. After the master
16
station transmits a pulse, each secondary station delays its
own transmission for a fixed time, called the secondary
coding delay. This ceding delay is synchronized through the
use of cesium frequency standards at each station. The high
stability and accuracy of these standards permit each
station to derive its own time of transmission without
reference to another station [DMA, 1977]. Secondary coding
delays are predetermined by system propagation times and
equipment characteristics [ Laurila, 1976].
The pulse from the master transmitter is distinguished
from those of the secondaries through phase coding of the
pulses. Phase coding refers to the inversion of the nega-
tive and positive peaks of the sine wave comprising the 100
kHz carrier portion of the pulse. The purpose of the phase
coding is twofold:
"Firs-, it permits automatic discrimination between the
master and the various secondary stations, thereby
permitting all stations to be identified by their rela-
tive timing with respect to the master"
"Second, the phase coding nrcvides protection against
excessively long skywave delays which would cause the
late arrival of the proceeding pulse to coincide with
the leading edge or groun dwave portion of a oulse being
tracked" [Poppe, 1982].
The signals are received by a mobile receiver where the
differences in time of arrival of the master signal and
various secondary signals are measurad and displayed en the
indicator portion of the LORAN-C set. The accuracy of this
time difference is increased by phase comparision "of the
synchronized 100 kHz carrier within the master and secondary
pulses" [Laurila, 1976]. This measured time difference (TD
- in microseconds) represents a hyperbolic Line of Position
(LOP) [U.S. Coast Guard, 197a]. The intersection of two or
17
SECONDARY
( X )
THE LOCUS OF ALL POSITIONS WHERE THE
OBSERVED TIME DIFFERENCE BETWEEN
THE TIMES OF ARRIVAL OF THE M & X
SIGNALS IS CONSTANT.
OBSERVER
( HYPERBOLIC FIX )
SECONDARY
( Y)
J
Figure 2.1
Hyperbolic Fix (From Coast Guard LORAN-C User
Handbook, 197U)
more LORAN-C LOP*s defines the posi*:ion of the observer
(Figure 2.1). When plotted en a chart, the intersection of
the resultant hyperbolic lines defines a geographical
position [Speight, 1982].
18
B. PHASE LAG
In a vacuum, the velocity of radiated energy from an antenna
for LORAN-C is 29979 2,458 Icm/sec. Since radiated energy
cannot be shown pic^orially, the phase of "che transmitted
radiations is used. The lines of constant phase of the
transmitted radiation are shown in Figure 2.2 by the curved
lines labeled aa' , bb' , and other similar designa-ions.
They define the wave front as it proceeds outward from zh=i
antenna in all directions. The distance between each line
of constant phase is one wavelength (x) [Admiralty Manual,
1 96 5 ] .
\ (meters) = 29 9792.4 58 km/sec t frequency in kHz
"1
Antenna
A
< X ^\< X >
< X >
I
J
Figure 2.2 Transmitted Radiation
The velocity of the radiated energy in air depends on
:empera-ure, pressure, humidity, and the narure of the
19
surface over which the transmissions pass. The retardation
of a -cransmitted wave is known as phase lag. When low
frequencies are employed, such as LOHAN-C at 100, kHz, the
effects of change of temperature, pressure, and humidity are
swamped by the effects caused by changes in the na-ure of
the surface over which the transmissions are traveling
[Admiralty Manual, 1965]. The Genaral Electric Company,
TEMPO division, conducted a LORAN-C Signal Analysis
Experiment under the direction of the U.S. Coasr Guard.
This experiment was conducted along the U.S. Wesx Coast.
The General Electric Company recorded an overall change of
108 nsec and 116 nsec time of arrival from -tha master and
X-secondary stations respectively after a storm [Samaddar,
1980]. If ASF Correctors are as large as two microseconds
on the West Coast [DMAHTC, 1981] then the ASF Corrections
are 20 times larger in magnitude than weather effects for
the 9940 West Coast LORAN-C chain.
The change in transmission ratas or phase lag are a
result of the amount of energy transferred from -^he tran-
smitted radiation. This energy transfer depends en the
absorption qualities (inversely related to conductivity) of
the surface over which they are passing and their wavelength
(or frequency) . The lower the conductivity and the longer
the wavelength (or lower the frequency) the greater the
transfer of energy, and vice versa. Seawater has a rela-
tively high conductivity (5.0 mhos/meter). Land has a much
lower conductivity, which varies from marsh (fairly high) to
dry sand and rock (very low) [Admiralty Manual, 1965]. Two
excellant papers that discuss the electrical properties of
soil are those of Smith-Rose [1934] and Pressey, Ashwell,
and Fowler [1956]. Smith-Rose [1934] found that the conduc-
tivities for soil ranged from 0.18 mhos/meter for a grey
clay with salt to .00001 mhos/meter for granite.
20
Antenna
A
^0 120 160
Distance from Antenna (km)
Sea Level
Figure 2.3 Phasa Lag
Phase lag is illustrated in Figure 2.3, The lines of
constant phase , aa', bb ' , and cc' become distorted as they
progress along the sea surface. The dotted lines represent
the position of the lines of constant phase in the absence
of the sea surface. The wavelength (X) in meters, measured
at heights of several wavelengths above the sea, remains
about the same as the direct wave through the earth's atmo-
sphere at 299691 km/sec i- freguency in kHz. Near sea level
the absorption of energy retards -he progress of -he wave-
front, and makes the sea level wavelength ( x' ) less than A.
As the lines of constant phase progress away from, the
antenna the phase lag increases with distance. This is
known as the Secondary Factor (SF) [Admiralty Manual, 1965].
The most intriguing feature of phase lag occurs at the
coastal interface where there is an extreme change in the
conductivites between land and sea. Visualizing the wave
front in three dimensions, the lower part of the wave,
slowed by the drag of the ground, lags further and further
behind the upper part as the wave crosses the land. At the
21
coastline it suddenly encounters the much lower impedance of
the sea, and in a very short distance the bottom of the wave
tries to catch up with the top, as though the whole wave
front were an elastic balloon. This is known as "phase
recovery" [Eaton, 1979]. Phase recovery was verified during
tests on Decca transmissions across the south coast of
England by Pressey, Ashwell , and Fowler [1956]-
The determination of the Secondary Factor for seawater
is fairly direct since the conductivity of seawater (5.0
mhos/meter) is fairly constant. But for land the conduc-
tivity can vary depending on the type of soil and its water
content [Smith-Rose, 1934].
Phase lag for radiated energy over land can be deter-
mined two ways:
1) Assign an average land conductivity to the ASF
Model, For example, the average conductivity for
the soil on the west coast is 0.003 mhos/meter.
The average land conductivity will determine the
average phase retardation of the path [The
Analytic Science Corporation ,1979].
2) Determine every conductivity for each portion
of a line segment from the transmitter to the
receiver. The total of these conductivity
segments constituting a land-water profile will
determine the total phase retardation of the path
[Speignt, 1982].
C. TD MODEL
Positional fix accuracy using LORAN-C is primarily
dependent on a chart makers ability to accurately compute
the expected difference in time-of-acri val (TOA) of received
groundwave signals from the transmitting s-^ations. Time
22
differences (TD) , are the differences between the TO As of
the secondary and master transmitters.
TD = TOA - TOA (2.1)
i 1 m
i = Secondary Station
m = Master Station
TOA computations are dependent, upon an accurate knowledge of
the signal phase delay.
The phase delay of a groundwave signal is generally
expressed as:
T + SF (2.2)
nR
C
SF
where n is the surface refractive index, C is the speed of
light in a vacuum, R is the range between the raceiver posi-
tion and the trajismi-t; ting station. The primary phase delay,
T, is the computed travel time of the LOEAN-C pulse over a
distance equal to the tran smit ter-tD-receiver greax circle
path length, taking into account the veiccitv of electromag-
netic waves and the index of refraction of the atmosphere.
The secondary factor (SF) is a correction to the primary
phase delay and accounts for the phase lag. The dominant
term in (2.2) is the primary phase dslay (T) . The SF is
usually an order of magnitude smaller [The Analytic Science
Corporation, 1979].
Thus, time-of-arr ivals can be expressed as:
TOA = T + SF + CD (2.3)
i i i i
TOA = T + SF {2,'i)
m mm
where CD is the true emission delay or coding delay for the
LORAN-C chain [The Analytic Science Corporation, 1979]- The
coding delay is equal to a time delay plus a computed one
23
way baseline time (Be) which includes the secondary phase
correction for an all seawater path. The oneway baseline
time (Be) is equal to the distance between the masx.er and
secondary transmitters in meters divided by the propagation
velocity of LORAN-C through the earth's atmosphere (299.691
meters per microsecond [Navigation Department DMA, 1982]).
See Table II for Coding Delay values for the 9940 chain
TABLE II
Coding Delay 9940 LORAN-C Chain
Pair 9940-W:
9940-X:
9940-Y:
CD
CD
CD
+ Be = 1 1000
+ Be = 2 7000
+ Be = 4 0000
■•• 2796.90 = 13796.90
+ 1094.49 = 28094.49
+ 1967.27 = 41967,27
|isec
/isec
;us6c
1
1
1
J
[Riordan, 1979]. Combining equations 2.1, 2.3, and 2.4, the
true TD is given by equation 2.5 [The Analytic Science
Cor po rar io n , 1979].
TD = (T - T ) + (SF - SF ) + CD (2.5)
i i m i m i
D. SEMI-SaPIRICAL TD GRID CALIBRATION MODEL
The Semi-Empirical TD Model was developed by The
Analytic Science Corporation [1979] in Reading,
Massachusetts for the »est Coast 9940 LORAN-C chain.
Similar "time difference (TD) grid calibration techniques
have been successfully employed to develop an accurate
(approximately 100 nsec drms) calibra-ed grid for St. Marys
River LORAN-C chain", [The Analytic Science Corporation,
197 9].
24
^ • l^chnical Approach
The Semi-Empirical Model is based on Millington* s
empirical approach for computing the secondary factor over a
mixed (multiple-homogeneous segment) path which combines
land and sea phase delays. The generalized semi-empirical
polynomial functional form for the SF of the LORAN-C station
is given by:
K
2 L
SF = SF(T ,6)= I AT +S(C 5inl6 ♦ D ccsl B ) (2.6)
j j j k=-K^ k j 1=1 jl j jl j
where
j = secondary (W,X, or Y) cr master (M) station,
nS ,
T = = jth station-to user primary phase delay,
j c
R = jth station-to-user great-circle path length,
j
8 = user path bearing angl=^ at the jth station,
j
K , K and L are oositive integers,
12
C and D are the station-dependent coefficients
jl jl
of harmonic terms in the model,
A is the range-da pandent coefficient of -he model
which may in general be station-dependent.
Data from 27 coastal sites distributed along the
West Coast and 122 land-sea sitas distributed in the
Southern California CCZ (between Point Arguallo and San
Diego - see Figure 2.4) were used in a Kalman estimation
algorithm to compute the uncertain coefficients of the land
and sea models of the TD grid calibration algorithm. (An
explanation of Kalman filtering for the layman is presented
by Roger M. du Plessis [1967].) The calibrated algorithm
25
(X)
Middletown V /\
California \ ^Tn
Monterey Bay
Point Arguello
Figaro 2.4 Loca-iOx' of West Coast LORAN-C Statio:
26
was used to compute TDs at each data site and the TD resi-
duals (difference between measured aad calibrated TDs) were
examined. Adjustments were then made to the TD model struc-
ture in an attempt to further reduce the residuals. This
process of adjusting the aodel structure is repeated until
the residuals agree with the expected theoretical covariance
associated with the TD model. The model which exhibited the
"best" performance was selected as the West Coast TD grid
model.
2. Generalized Ranqe/B earing Model
The Generalized Range/Bearing (GRB) Model was
selected as the "best" semi-empirical calibration model for
the West Coast chain. The semi-empirical function is:
L
S? = AO + A1T + E (C sinlS + D cosl 6 ) usee (2.7)
J j 1=1 jl j jl j "^
where AO , A1 , C and D are the model coefficients, is the
path bearing angle measurad positive clockwise from north at
the jth ( W,X,Y or M) station and I is nhe path range to -he
jxh station. The GRB model is rela-cively complex and is
expected to exhibit superior performance. The extensive
mcdel is based on knowing the distance overland (TL) , -he
distance over water (TS) , the total distance (T) , and the
path bearing angle 3 (Figtirre 2.5).
It was noted that the calibrated mcdel was expected
to be accurate and applicable only over the extent of ranges
and bearing angles embodied in the calibration data. Hence,
outside the region covered by the calibration data the model
may not be as accurate as within the data coverage region.
Osing the GRB model, a drms value of approximately 0.3 jisec
was expected in areas where land data alone was used to
calibrate the model. Inclusion of sea calibration data
produced a drms value of 9.35 to 0.50 jis<5c.
27
r"
USER :
MASTER STATION
D'.
,^ECONDARY
STATION
Figure 2.5 Jlixed Path ID Geomezzj
The time difference (TD) is expressed by -he
following equation:
TD = (T - T ) + (SF. - SF ) + CD. + b. /isec (2.8)
i i m 11 11
28
where
n R.
T = |isec,
i c
n R
m
T = )isec,
m c
R = ith secondary s-at ion-to-user great-circle path length,
i
c = -he speed of an al ectromagnetic wave in a vacuum
= 2.99792458 X 10^ m/sec,
n = surface refractive index
= 1.000338,
R = master station-to- user great-circle path length,
m
CD.= coding delay found in Table II,
b = TD bias associated with the ith secondary
i
station (psec) (Table III) .
TABLE III
TD Bias (b) - jasec
r-
"T
TDW
-0.
BbU
1
TDX
-1.
,173
1
c-
TDY
-0.
353
1
-J
SF,= 0.5 {-SI + S2 + S3 - 34 + S5 + S6) (2.9)
1
The term SI is -he SF of a land path of length T
(;isec) from the jth station:
SI = 0.795/T + 0.439 + (0.00245) T usee. (2.10)
sj sj
The terms S2, S3, and S4 combine to make up the
secondary factior for the seawater path Isngihs.
29
S2 is the SF for the total path. S3 is the SF
using the seawater coefficients for the portion
with seawater, and S4 is the SF using the seawater
coefficients for the land path distance. The
seawater coefficients are found in Table IV.
TABLE IV
Seawater Coefficients
a1
a2
a3
b1
b2
b3
128.8
0. 187
0.000652
3.188 I
-0.594 I
0.000329 I
SF (T)
SF (T)
s
bl/T + b2 + (b3) T ;isec,
if 10 < T < 540 ;isec.
al/T + a2 ^- (a3) T
if T > 540 jasec.
;isec.
(2.1 la)
(2. lib)
Term S5 and S6 are the SFs of land paths of length
T and T . SF5 is the Secondary Factor for the
J Ij
total length using the land coefficients whereas
SF6 is the Secondary Factor for the distance over
land using •^. he land coefficients. The land
coefficients are found in Table V.
S5 = SF (T , 8 ) ,
L J J
S6 = SF (T ,6 ) ,
L Lj y
w he r e
SF (T , 3 ) = \0 ^ (A1) r + I (C sinl 3 + D cosl 3 ) (2.12)
L j j j 1=1 jl j jl j
30
TABLE V
Land Coefficients
AO
=
1.428
D
x2
=
0.9U2
A1
=
0.00158
C
=
0.0
C
=
0.0
c
=
0.588
Wl
y2
C
=
-0.711
D
=
0.0
w2
yi
D
=
0. 323
D
=
0.0
wl
y2
D
t=
0.0
c
=
1.010
w2
ml
C
=
0.0
c
=
-0.196
xl
in2
C
=
0.0
D
=
-0.893
x2
ml
D
=
0,0
D
=
-0.355
Xl
m2
._ J
S. DMAHTC MODEL
1 . S6 a S F Model
The equations for the Sea SF Model is:
SF = (B1/T) + 32 + (B3 T) usec, if ID < T < 537 jis^c, (2.13a)
SF = (A1/T) + A2 + (A3 T) a sec, if T > 537 jtisec, (2.13b)
where T is -he primary phase delay (or range) in microse-
conds (usee); Ak and Bk (k = 1, 2, and 3) are the sea model
coefficients used by DMAHTC in program TDGRID [Funakoshi,
1982]. The coefficients are found in Table VI.
2. Land SF Model
The solution used to resolve DMAHTC Taoles ASF
Corrections is called Millington's Method [DMAHTC, 1981].
31
TABLE VI
Sea SF Model Coefficients
A1
A2
A3
B1
iB2
|B3
129.04323
-0.40758
0.00064576813
2.741282
-0.011402
0.00032774815
— I
This method is based on the premise that the phase distor-
tion due to a composite land-sea path is the arithmetic
average of the phase distortion foand in the forward and
reverse paths of the propagated signal [ DMAHTC, 1981]. For
example, in Figure 2.6 two azimuths have been drawn on the
map and are labeled as 210° and 235o. Also placed on the
map are the proper ground conductivities which have been
determined in the field by the Coast Guard calibration team.
A great circle drawn on the appropriate chart or charts from
the LORAN-C Station coordinates to the area undsr considera-
tion spans various lengths of land and seawater. Each
length or segment will have a specific conductivity and
distance. The total of these conductivity segments, consti-
tuting a land-water profile, will determine the total phase
retardation along that path [U.S. Naval Oceanographic
Office, 1982]. All azimuths and distances are computed
based on the World Geodetic Systsm (WGS) datum [DMAHTC,
198 1]. The values of phase retardation for a given ground
conductivity are tabulated in the National Bureau of
Standards (N3S) Circular 57 3 [Speight, 1982].
The formula used to derive the ASF Correction for
the time difference for a master-slave transmitting station
pair is:
ASF Correction = (-Slave Error) - (-Master Error)
32
TRANSMITTER
(,003 mhos/m)
SEAWATER
(5.0 mhos/m)
^
Figure 2.6 Composite Land-Sea Path
The mean values derived for one station from the
forward and reverse solution of Millington's Me-hod are
subtracted from the Sea SF Model. The differences are
33
presented in an azimuth array. This array is a series of
geodetic azimuths radiating from the transmitter with
corrections computed at incremented distances along each
azimuth. After the valaes are computed in the azimuthal
array a compu-cer program rearranges them into a matrix form.
The matrix form is the arrangement of corrections into rows
and columns covering a specified geographical area at a
constant spacing. The purpose of the matrix is to enable
the corrections from two LORAN-C transmitters to be added
algebraically, combined into a single matrix, and arranged
in the desired tabular form. This tabular form is the body
of the table [DMAHTC, 1981].
3 • Ta b le Description
Each table contains a complete chain. Figures 2.5,
2.6, and 2.7 depict LORAN-C ASF Correctors for chains
9940-W, X, and Y for Monterey Bay, California. A table
section is prepared for each station pair (master station
and one slave station) in a LORAN-:: chain. As a rule the
limits of the table coverage are determined by the range of
the groundwave transmissions for t-hs LORAN-C chain. Each
page of corrections in the table covers an area three
degrees in latitude by one degree of longitude, with correc-
tions printed in increments of five minutes of arc. Rate
designation and page numbers are printed at the top of each
correction page. Those pages where latitude and longitude
limits contain both land and sea are included but correc-
tions apply only for the area covered by the a.S. Coastal
Confluence Zone (CCZ) . Large land bodies and areas ouside
the CCZ are represented by blank spaces on the page. ASF
correcticn values can be either positive or negative (posi-
tive values are shown without sign). Areas requiring no
correction show a zero value which in some cases in preceded
34
122"
LONGITUDE WEST
1
121
0' 55 50
45
40 35 30 25 20 15
10 5 0 •
39* O'
55
50
45
40
35
30
Z'j
20
15
10
-
L 0 \
A 38 0
T
I 55
T 50
U 45
0 40
E 35
30
25
20
N 15
0 10
R 0 5,
T 37 0
H
55
-1 .6 -1 .6
50
-1.6 -1 .5 -1.4
1 ^5
-1 .5 -1 .4 -1.6
i 40
-1.4 -1 .3 -1.5
,
i 35
-1.3
■ j
1 30
-1 .4
■
1 25
-1 .4 -1.5
20
-1 .2 -1.1
15
-1 .4 -1 .3 -1 .6
10
-1.3 -1.3 -1.7
-1.7
0 5
-1.3 -1.4 -1.7
-1 .7
-1.8
36 0'
-1 .2 -1 .6 -1 .6
• 1
-1 .7
-1.7 -1.7 -1.8
III 1
i
Figure 2.7
9940-v* ASF Correction Taole
by a negative sign indicating that the zero results from the
rounding off of a value slightly less than zero (indicates
the trend of the correction).
35
1
1
1
122"
LONGITUDE WEST
121"
1
1
0*
55
59
45 40 35 30 25 20 15
10 5 0 '
39° O'
55
! 50
i 45
40
35
*
30
'
25
20
15
10
>- 0 5,
A 38 0
-
T
1 55
1
T 50
1
U 45
1
0 40
f
E 35
1
30
•
1
25
1
1
20
N 15
1
1
0 10
1
1
R 0 5
T 37 0'
H
55
0.9
0.9
50
0.9
0.9
0.8
*
45
1 .0
1 .0
1 .0
40
1 .0
1.2
1 . 1
35
1.1
30
1 .0
25
1 . 1
1 .0
1
20
1.3
1 .4
15
1 .2
1.3
1 .1
10
1 .3
1.3
1 .0
1 .0
1
, 5
1 .3
1.3
1 . 1
1.1 0.9
j
36 0'
1 .3
1.3
1.2
1.2 1.0 1.0 0.7
1
1
I
" " " ■
Figure 2.8 99aO-X ASF Correction Tacls
The table can be entered directly by using the
ship's position determined to the nearest five minutes of
arc in la-irude and longitude either by dead reckoning or
36
LONGITUDE WEST
1
122'
121 'i
0' 55 50 45
40 35 30 25 20 15
10 5 0 '
0 >
39 0
55
50
45
40
35
30
■
25
20
15
10
A 38 0
-
T
I 55
T 50
•^
U 45
0 40
E 35
30
25
.
20
H 15
0 10
R 0 5
T 37 0
H
5 5
-0.2 -0.2
50
-0.3 -0.3 -0.3
45
-0.3 -0.2 -0.4
40
-0.4 -0.3 -0.6
35
-0.5
30
-0.7
25
-0.6 -0.8
20
-0.5 -0.6
15
-0.4 -0.3 -0.3
10
-0.2 -0.2 -0.2 -0.2
, 5
-0.1 -0. 1 -0.1 -0.1
-0.2
36 0'
0.0 0.0 0.0 0.1
-0.1 0.0 0.0
^
,
_ _ , J
Figure 2.9 9940-Y ASF Correction Table
scr^e other means. To find the page with the appropriate
correction, the Page Indexes of the table should be
utilized. These indexes show the limi-s and page number of
37
all pages in the table. To locate the number of the page on
which the desired correction is to be found ihe Page Index
is entered with the ship's position. In some cases the
ship's position will fall on the page limit in either lati-
tude or longitude or both. These positions are repeated on
both pages and either page may be used.
The ASF Correction is added algebraically to the
time difference for the LOR AN-C pair. Interpolation of this
data will not necessarily improve the accuracy due to the
method used to determine ASF Correctors [DMAHTC, 1981].
Since the correctors are computed in the azimuthal array and
are based on the conductivity and distance over which the
LORAN-C electromagnetic wave travels, the ASF Corrector
between the published AS? Correctors in the tables may not
be the linear interpolated values. For example, the ASF
Corrector for a distance of 500 m with an azimuth of 180° is
equal to 1.5 jisec. The ASF corrector for a distance of 500
m with an azimuth of 181^ is equal to 1.6 psec. The inter-
polated value between 1.5 and 1.6 is 1.55. The true ASF
Corrector for the 180° 30* azimuth is 1.U since the land
distance for the same azimuth is less than the land distance
for the 180O and 181 o azimuth. The LORAN-C signal passed
over a harbor [Dansford, 1982]-
F. ATTEMPTED DETERMINATION OF ASF CORRECTORS BY HYDRO FIELD
PARTIES
One of the major problems encoantered by hydrographic
survey operating units when using LORAN-C for position
control is the determination of the AS? Correctors for the
survey area. The four surveys mentioned in the Introduction
all made attempts to determine the correctors by comparing
the LORAN-C rates to a second source.
38
Bathymetric Surveys H-9822 [NOAA H-9822, 1979] and
H-9823 [NOAA H-9823, 1979] Gulf of Alaska, compared the
rates from an Internav LC-20U LORAN receiver to computed
rates from a position obtained from a JMR-1 Satellite
Navigation Receiver when available. Shore ties using radar
ranges, visual bearings, and sextant angles in comparison to
LORAN-C rates were also made prior to and after each survey.
The calibrations of LORAN-C rates were based on the satel-
lite positions only since the positional computation of
LORAN-C and JMR Doppler Satellite were made on the WGS 1972
datum whereas the land ties were based on the NAD 1927
datum.
Bathymetric survey SU- 40-7-82 which extended along the
Washington, Oregon, and California Coasts used LORAN-C as
navigational control. LORAN-C time differences were
compared with SATNAV positions. The report did not indicate
whether any correctors were applied [NOAA Ship Surveyor,
198 2]-
Finally, Hydrographic Survey S-K902-Wh-82,
Reconnaissance Survey of Safety Fairways, Gulf of Mexico
used LORAN-C as a pcsizioning control system. The posi-
tioning unit was an LC-20'i receiver. LORAN-C rates were
input via the HIDROPLOT Controller, a special purpose
input-outpu"*- interface which is the nucleus of the computer
system hardware [dmbach, 19 76]. Positions were computed and
plotted by Program RK121, LORAN-C Real-Time HYDROPLOT
[Backus, 1980].
ASF Correctors for LORAN-C were achieved by visual cali-
bration using three poirn: sextant fixes using charted oil
rigs as control in the vicinity of the survey area. A
three-point sextant fix is a convenient and accurate method
for determining the position of a hydrographic survey
vessel. Sextants are used to measure two angles between
39
three objects of known geographic positon. The center
object is common to both angles. The position of the obser-
vers taking rhe angles is fixed by the intersection of three
circular lines of position [Umbach, 1976].
These sextant angles were recorded and later tranf erred
onto their respective charts using a plastic rhree-arm
protractor. A plastic three-arm protractor is transparent
and made up of one fixed arm and two movable arms which
contain an etched line that is radial with the center of the
protractor [Ombach, 1976]. Sextant angles observed in pairs
for a resection fix with a common center mark may be plotted
directly by this instrument. When the three arms are placed
at the angles observed and fitted sd as to pass through the
plotted positions of the observed stations on the field
sheet, the hole at the center of the three-arm protractor is
the fixed positon of the vessel [Ingham, 1975].
Partial correctors for each area surveyed were defined
by comparing the observed rates and the determined rates
plotted on the nautical chart. The partial correctors were
applied via the HYDROPLOT Controller. However, even after
applying these correctors, the plotted position still disa-
greed with the ship's determined postion with respect to the
oil rigs. Ship's personnel attributed the discrepancies to
one or more of the following:
1) Accuracy of the charted rigs,
2) Weather effect on LORAN-C,
3) Time of day,
4) Propagation of signal over land path,
5) Three-arm protractor accuracy, and
6) Error in the conversion by the software
of the LORAN-C rates to latitude and
longitude [NO\A, 1982].
UO
There is an apparent need for a LORAN-C calibration
routine aboard NOAA ships which provides the ASF Correctors
for program RK121, LORAN-C Real-Time HYD20PL0T. The routine
should use the same geodetic distance computation found in
RK121 and use the same datum as that of the nautical chart
of the survey area.
The above mentioned discrepancies illustrate the defi-
ciencies in applying a single ASF Corrector -o LORAN-C data.
The accuracies for hydrography cannot be met using single
correctors because the errors are non-linear and systematic.
They cannot be distributed like residuals in a traverse.
Schnebele [1979] has already proven that single ASF area
correctors to LORAN-C positions do not meet the accuracy
s-andards of the NOS Hydrographic Manual.
Based on visual inspection of the DMAHTC LORAN-C
Correction Tables, ASF Correctors should be updated every
five minutes of Latitude or Longitude change. In Monterey
Bay, California, there is approximately 0.1 to 0.2 usee
difference for every five minures of change, a poten-ial
error of 55 to 110 meters.
41
III. EXPERIMENIAL PROCEDURE
A. FIELD PROCEDURES
In order to compare the use of differential LOHAN-C with
ASF multiple correctors, typical survey operations were
planned for the southern portion of Monterey 3ay,
California. This survey was conducted in con juncxion with a
comparative evaluation of multiple lines of positon for
selective positioning methods [Anderson, 1982]. Four micro-
wave ranging systems were set on Known geographic positions
TABLE VII
Geographic Names and Positions
Microwave System
Stations
Seaside 4 (1964)
Use Hon Ecc.
Geoceiver Sec.
Mussel Sec.
Park (1931)
Mulligan RM1
Range 7 (197 2)
Mussel (1932)
Geographic Position
(NAD 19 27)
35 0 36
1210 51
360 36
1210 52
36 0 36
121 0 53
36 0 36
1210 54
360 53
121 0 49
360 44
121 0 47
360 39
1210 49
36 0 37
1210 54
23.44596"
38.83281"
04,73031"
35.98040"
32.49281"
25.21162"
18.25484"
11.49661"
13.30600"
46.743 00"
5 6.49531"
52.31090"
02.47787"
08.58202"
18.15100"
11.49661"
^§l£§ !^§^d.
June 3-5,1982
June 3-5,1982
June 3-5,1982
June 3-5,1982
June 6-7,1982
Jane 6-7,1982
June 6-7, 1982
June 6-7, 1982
42
listed in Table VII. A series of tracklines were run in two
separate areas as shown on ?igure 3-1. To ensure that the
microwave postioning system was working properly, the equip-
ment was calibrated over known baselines of 1497.47 meters
and 7877.3 1 meters at the beginning and end of the project.
Trackline observations were only made during the daytime in
fair weather conditions so as to eliminate sky waves and
weather changes that influence LORAN-C signal propagation
characteristics [Samaddar, 1982]. The vessel used was the
126 foot R/V Acania which is operated by the Naval
Postgraduate School.
The positioning equipment consisted of a Micrologic
ML- 1000 LORAN-C receiver (0.0 1 usee resolution) and a
Trisponder Microwave System provided by Racal-DECCA Survey,
Inc. The Trisponder Microwave System consisted of four DNT1
Model 217C transponders, four DNT1 Model 21017 HP sector
antennas with 87o by 5° beam widths, cne DNT1 DDMU (Digital
Distance Measuring Unit), two Omni DVTI Model 2 1019 HP
antennas, a Houston Instruments Model DP3-M2D/RC3 plotter,
and a Texas Instruments 743 terminal (Table 7III) . The
manufacturer's published accuracy for the positioning equip-
ment is ± 1 m for a single range [ Ra::al-DECCA Survey, 1981].
Anderson [1982] discusses the accuracy of four lines of
position. The four Decca Trisponder distances were recorded
via a Texas Instruments 743 data terminal while the LORAN-C
rates were manually logged. The data was acquired at one
minute intervals while the ship maintained constant course
and speed. The recorded LORAN-C rates were 9940-W, X, and Y
of the West Coast chain.
To test the potential for calibrating the LORAN-C System
using the Semi-Smpirical Model, the correction tables, and
multiple observed field correctors, the positions derived
from the microwave system n easurement s were used to compute
43
122*
122
Figure 3.1 Location of Test irsas in Mcntsrey Bay
44
TABLE VIII
Microwave Positioning Equipment
Eauipment
S/N
DNT1 Model 217C Transponders (Code
(Code
(Code
(Code
72R]
74R
76R
78R
1
1
3323
3320
3321
3322
DNT1 Model 2107 HP Sector Antenna
(870 by 50 Beam Widths)
185
186
187
191
DNT1 DDMV Model 540
426
Omni, DNT1 Antenna Model 2 1019 HP
194
200
Houston Instruments Plotter Model
jP3-
12D/RC3
10722-10
Texas Instruments 7U 3 Terminal
34418
— _ _ _j
expected LORAN-C time differences at each point. The
difference or offset between these expected time differences
and the observed values were computed for the three methods.
The mean offset, standard error, and drms values were also
computed and compared.
3. MICROWAVE SYSTEM POSITIONING
The geographic position of the ship based on four lines
of position was determined using a computer program callei
GP3YLQ (Geographic Position by Least Squares) written by the
author (See Appendix A) . GPBYLQ contains subroutine LSQR
(Least Squares), which is a least squares adjustment written
by Paul R. Wolf, Ph.D. [1974] and revised by LCDR D. Leath
[1981]. Geographic postions were converted to X,Y (meters)
which in turn were converted to geographic position via
subroutines GPTOXY and XYTOGP, respectively [Wallace, 1974],
45
Subroutine GPTOXY and XYTOGP are based on the Modified
Transverse Mercator Grid srfhich was centered in the survey
area.
The Modified Transverse Mercator (MTM) projection is
used by the National Ocsan Survey and is similar to the
projection used in the Universal Transverse Msrcaror (UTM)
system. The main difference is that in the MTM a Central
Meridian is picked that is near the survey area instead of
being fixed at a particular meridian [Wallace, 1971].
Central Meridian (CMER) , False Easting (FEST) , and
Controlling Latitude (CLAT) are the three parameters which
define the MTM projection. CMER is the mean longitude
computed using the maximum and minimum longitudes of the
survey limits, FEST is the X-Coordinate that is assigned to
the Central Meridian, and CLAT is the distance in meters
from the equator to seme reference latitude [Wallace, 1971],
The Central Meridian, False Easting, and Controlling
Latitude used for Monterey Bay, California referenced to NOS
Chart 18685 are:
CMER = 1210 56' 00 .0",
FEST = 20000 .0,
CLAT = 4050000.0,
To be consistent with the National Ocean Survey charts
of the area, all computations were done relative to North
American Datum (NAD) 1927 geographic positions. All
programs were executed on an IBM 30 33 computer located at
W.R. Church Computer Center, Naval Postgraduate School,
Monterey, California.
C. LORAN-C COMPUTATIONS
The differences or offset between the observed and
computed LORAN-C rates using the Semi-Empirical TD Model,
ASF LORAN-C Correction Tables, or the Multiple Observed
U6
TABLE IX
LORAN-C 9 940 Chain Data
r
Station
Geogrj
aphic Position
(NAD 1927)
Master - Fallon, Nevada
390
1180
33'
49'
07.03"N
52.23"W
Slave - George. Washington
9946-W
470
1190
03'
44"
48.82"N
34.78"W
Slave - Middletown,
California 99 40-X
380
1220
46'
29'
' 57.49"N
40.04"W
Slave - Searchlight, Nevada
9940-Y
350
1140
19
48'
• 1 8 . 3 2 •• N
' 13.95»W
t -
_ J
Correctors were compared to the offsets between the observed
TD rates and the computed rates for which only the seawater
Secondary Factors (SF) were applied. The comparison of the
offsets between the four methods illustrates the improvement
in positional accuracy after applying ASF Correctors.
^ • Sea water Secondary Factors (SF)
Time differences using only the seawa'rer Secondary
Factors for each of the geographic positions were computed
using program LORAN written by the author {Appendix C) .
Seawater Secondary Factors (SF) were computed using formula
2.13 and the coefficents found in Table VI. All TOA
distances in meters were determined using subroutine INVER1.
INVER1 is a geodetic inverse routine using T. Vincenty's
modified Rainsford's method with Helmert's elliptical terms,
programed by LCDR L. Pfeifer, NOAA [1975]. Subroutine
INVER1 is accurate to 0.0 001 m halfway around the world
[Pfeifer, 1982]. All distances were converted to microse-
conds using 299.792458 m/|iS€C. Time differences (TD) were
47
computed from equation 2.5. North American Datum 1927
geographic positions were used for all computations. See
Table IX for the positions of LORAN-C 9940 transmitters
[Riordan, 1979].
To ensure that subroutine INVER1 was functioning
properly the distances between ths master and secondary
stations were compared to the NOS published baseline
distance [Riordan, 1979]. The published distances and -he
results from routine INVER1 are lisred below.
Published Baseline Commuted Baseline
Dis-ance (m) Distance z INVERT (m)
99U0-W 837,777.0929 837,777.115
99U0-X 327,886.3720 327,886.316
99aO-Y 589,298.5712 589,298.589
The difference between the publishei and computed baselines
ranged from 0.02 to 0.06 m.
Differences or offsets (x . ) were obtained by
subtracting the observed LORAN-C rates from the computed
values from the various methods. The mean difference or
offset (X) and standard deviation (s) in microseconds for
each rate were determined using equations 3.1 and 3,2
[Wonnacott, 1935]:
n
1
Mean offset (/isec) : x = - Z x (3.1)
n i
i=1
where: x = original observation in usee,
i
n = number of observations;
n
Standard 1
error: s = Z (x - x)2. (3.2)
(;isec) n-1 i
48
The mean offset in microseconds can be converted to meters
using equation 3.3 [Heinzen, 1977]:
Mean offset
(meters)
where w
1 =
X w
sin 01
(3.3)
= distance corresponding to one microsecond
on the baseline
= 149.396229 meters [Bigelow, 1963].
= one-half the angle between the
radius vectors from the position to
the master and secondary stations.
Equation 3.4 was used to compute drms values in meters
[Heinzen, 1977]:
drms=
sin ( + )
1 2
\1
(S^W)2
2
sin a
1
(S^W)2
+ +
2
sin a
2
2 pcos ( a + a ) s w
1 2 1
\
s w
2
(sin 01 ) (sin a )
(3.4)
where:
= correlation coefficient = 0.33,
s and a are as above with the subscript
denoting the appropria-^e
secondary station.
The correlation coefficient (p) is a resuli of the secondary
station having a common line of position with the master
station. Although often ignored, various authors assign
values ranging from 0.33 to 0.40. Bigelow [1963] chooses
P = 0.33.
2 . Semi-Empirical TD Grid
To determine if the Semi-Empirical TD Grid would
reduce either the drms value or the offset between the
observed and expected TD rates, program LOPLC (Line of
Position - LORAN-C) was written by the author (Appendix B) .
Program LOPLC computes the distance over land, distance over
water, the total distance, and the azimuth from north for
the transmitting station using Subroutine INVER1. The
49
land-sea distances were computed by selecting a point which
was located along the coast of Monterey Bay, Subroutine
SHORPT (Shore Point) interpolated a geodetic position from
23 geodetic points which outlined Monterey Bay, California,
The land/sea Secondary Factor was determined in Subroutine
SECFAC which is found in program LOPLC using eguations 2,8
through 2.12, All corapu taxions ware based on NAD 1927
geodetic datum.
Tests were made on program LOPLC using data found in
the Semi-Empirical TD Grid article [The Analytic Science
Corporation, 1979] using the WGS 1972 datum. Station TASC
55 located at latitude 3ao 34* 18,3" N and longitude 120o
39* 40. 3" W, was selected from The Analytic Science
Corpora-ion article. It was one of the stations used to
calibrate the coefficients for the Semi-Empirical TD equa-
tions discussed earlier. The only distances listed for TASC
55 were the individual distances ov=r land and over water.
The total distance between the transmitters and TASC 55 was
computed by adding the land and sea distances. The
following station-tc-site pazh segment lengths are listed
for TASC 55:
TASC 55 OI2R2
Distance Dist ance Distance Distance
Station (km) (km) (km) (km)
Master 540,730 35.548 575.278 576.083
X-Secondary 369,248 126,096 495.344 495.178
Y-Seccndary 525,659 15,832 541.491 541.308
Using the land distances, sea distance, and total distances
to TASC 55, and the computed azimth from subroutine INVER1,
program LOPLC produced offsets of 3.86 jisec and 1.13 jisf^c
for 9940-X and 9940-Y respectively when compared to the
expected time differences at TASC 55.
50
An attempt was made to determine if this discrepency
was due to program LOPLC. A comparison of the total
distances in the report between the transmitters and TASC 55
were compared to the total distances computed by INVEE1.
The difference in distance between Subroutine INVER1 and the
total distances from the Semi-Empirical TD 3rid acticle
ranged from 160 - 200 m corresponding to 0.5 to 0.65 usee.
This results in time difference errors of 0.04 usee for the
9940-X rate and a 0.10 ;isec for the 9940-Y rate, implying
that the offsets, 0.86 psec and 1.13 fisec, are caused in
part by the method in which the azimuth from north is deter-
mined. Unfortunately, azimuth data from TASC 55 was not
presented in the article. The Ar;alytic Science Corporation
has been contacted on numerous occasions in an attempt to
ascertain their method of determining distance and azimuth.
As of this date there has been no response.
Nevertheless, ttie data from Monterey Bay was
utilized in Program LOPLC to obtain results that could be
compared to that obtained by the other methods. If this
method is accurate enough, AS? Correctors could be deter-
mined via computer for each individual position without
using tables or field determined correctors. The mean
offset (X) and standard arror (s) in microseconds between
the observed and calculated rates w=re computed using equa-
tions 3.1 and 3.2, respectively. The mean offset (1) and
drms in meters were computed using equation 3.3 and 3.4,
respectively.
3. Calculated Table ASF Correctors
The offset between the observed LORAN-C rates and
the expected values with applied ASF Correctors from the
LORAN-C Correction Tables and the seawater Secondary Factor
(equation 2.13) were determined using program LORTAB which
51
was written by the author (Appendix D) . The ASF Corrector
for each position was selected using subroutine TABLE.
Subroutine TABLE, which is found in program LORTAB deter-
mines an AS? Corrector for each data point based on its
geodetic position. The ASF Correctors used in subroutine
TABLE (see Figures 2.5, 2.6, and 2.7) are located between
latitude 36° 35' N and 36o 55'N and longitude 121° 50' W and
1220 00' W. The difference in sign between the ASF
Correctors in subroutine TABLE and those found in ths
LORAN-C Correction Table is due to the difference in their
application. ASF Correctors from the tables are applied to
observed rates while ASF Correctors from Subroutine TABLE
are applied to the calculated LORAN-C rates. Negative ASF
Correctors from the LORAN-C Correction Table were applied to
the calculated time differences ro be consistent with the
application of Secondary Factors to the computed primary
phase delay.
As before equations 3.1 through 3.4 were used to
compute "he the mean offset , standard error, and drms. If
this application of LORAN-C Correction Tables is accurate
enough, it precludes the need to determine ASF Correctors in
the field.
'* • Observed ASF Cor rec tors
Observed ASF Correctors were determined using
Program ASFSEL (ASF Selection) which was written by the
author (Appendix E) , Schnebele's prior data [1979], and the
June 1982 data. Program ASFSEL (ASF Selec-ion) was writ-en
by the author. This program computes the ASF Correctors by
subtracting the observed LORAN-C rates from the expected
values. Only the seawater Secondary Factors from equation
2.13 have been applied to the calculated time differences.
52
The mean ASF Correctors for the LORAN-C rates were deter-
mined at every minute of latitude and longitude between 36°
50* N and 36° 35' N and 122° 04»W and 122o 49' W. See
j Appendix E for mean ASF Correctors at one minute intervals.
The ASF Correctors were then selected and assigned to
subroutine TABLS in Program L0RTA3 at five minute latitude
and longitude intervals. All values were entered to the
I nearest hundredth of a microsecond. Equations 3.1 through
3.4 were used as before for computations. This determina-
tion of ASF Correctors in the field, if accurate enough, may
allow the use of LORAN-C as a positioning system for hydro-
graphic surveying in the future.
53
IV. RESULTS
A total of 620 time differences (TD) and geographic
positions based on four lines of position were recorded in
the southern portion of Montery Bay in order to compare the
use of differential LORAN-C with ASF Multiple Correctors.
The data was divided into four sets. The first data set is
Schnebele's [1979] prior data consisting of a tozal of 130
data points collected on two separate days, Juna 12 and July
25, 1979 between latitude 36° 38» N and 36° 47' N and
longitude 121o 49' W and 12 20 02' W. The recorded LORAN-C
time difference rates were 9940-Y and 9940-W.
The second data set consists of 193 time differences
collected en June 7, 1982. The recorded LORAN-C rates were
also 9940-Y and 9940-W. Data set Number 2 is located
between latitude 36o 40' N and 36° 45' N and longitude 1210
54' M and 122° OO'W. See Figure 3.1 for the location of
test areas in Monterey Bay, California.
Data set Number 3 is located in the same area as data
set Number 2, between latitude 36o 40' N and 36° U5' N and
longitude 1210 W 54» and 122° 00' W. This set, consisting
of 128 points with recorded rates 9940-X and 9940-Y, was
recorded on June 6, 1982.
The final set, data set Number 4 was recorded between
June 3 and iJune 5, 1982. It contains 169 points located
between latitude 36o 36 N and 36o 39' N and longitude 121o
53' W and 121° 58' W. The recorded rates were 9940 -X and
9940-Y which are the same as data set three. The dana for
rates 9940-X and 9940-Y was Icepr separate so as to determine
if there was a significant difference between the offshore
(data set three) and inshore (data set four) drms values due
to phase recovery (see Chapter Two) .
54
TABLE X
Dara Set Parameter
Set
Number
Collection
Date
TD
R ates
N
Da
amber of
ta Points
Area
Lit
Limits
Lon
1
June 12, 1979
July 25, 1979
9 940-Y
9 9U0-W
130
36/38
36/47
121/49
122/02
2
June 7, 1982
9 940-Y
9940-W
193
36/40
36/45
121/54
122/00
3
June 6, 1982
9 940-X
9940-Y
128
36/40
36/45
121/54
122/00
4
June 5, 1982
9 94 0-X
9940-Y
169
36/36
36/39
121/53
121/58
J
Table X provides a convenient breakdown of the parame-
ters for each of the data sets. The table consists of the
data set number, the date the data was collected, the
LORAN-C time difference rates, the number of data poinxs,
and the area limits in latitude and longitude. It defines
the parameters for the four data sets of Tables XI (Seawater
Secondary Factor Errors) , XII (Semi-Empirical TD Grid
Correction Errors), XIII (Multiple LOP.AN-C Correc-ion Table
Errors) , and XIV (Multiple Observed Correction Errors) . All
basic data and computarion s are kept on file with the NOAA
hydrography instructor at the Naval Pos-graduate School,
Monterey, California. (Individuals seeking this information
should contact the Oceanography Department.)
A. SEAWATER SECONDARY FACTORS
LORAN-C rates were computed using only the seawater
Secondary Factor (equation 2. 13) . The LORAN-C computed time
differences are basically uncorrected rates since no
Additional Secondary Factors were applied. Offsets between
55
TABLE XI
Seawater Secondary Factor Errors
Rates
Mean
Offset
Mean
Offset
Standard
Error
(s-^sec)
a)
b)
c)
d)
Set
Y
W .
Set
Y
W
Set
X
Y
Se-c
X
Y
(Schnebele's data. 130 data points)
-0.508 -148.636 0.088
-1.2U1 -683.693 0.134
(June 7- 1982, 193 data points)
-0.526 -154.295 0.077
-1.283
-701.947
0. 124
(June 6, 1982, 128 data points)
1.565 535.449 0.059
-0.550
-161.404
0.083
(June 3-5, 1982, 169 data points)
1.582 550.089 0.055
-0.817 -237.109 0.187
drms
(1)
110.4
101.3
38.6
68. 1
the observed and calculated rates using only the seawater
Secondary Factor were computed to illustrate -he improvement
in posi-^,ion after applying the Semi-Empirical TD Model, ASF
Loran-C Correction Tables, or nhe Multiple Observed ASF
Correctors. The mean offset in microseconds and meters,
standard error, and the drms are found in Table XI.
Examples of offsets for several data points are listed in
Appendix 3 after program LORAN.
It was stated earlier in Chapter One that Schnebele
obtained a 66 m 1 drms using Y and W rates. The drms of 66
m was obtained using 48 data points which were located 10 km
or more offshore. The 110.4 m 1 dris for the Y and W rates
in Table XI is a result of combining the 130 inshore and
offshore positions. The increase from 66.0 m 1 drms for the
offshore positions to 110.4 m for combined offshore and
inshore positions indicatas that the application of a single
56
ASF Corrector is dependant on the size of the area. The 48
data points were located between latitude 36° 41 » N and 360
46» N and longitude 121o 55 • W and 122° 02« W.
Schnebele also obtained the 66 m 1 drms by adjusting the
LORAN-C observed time differences which were skewed due to
the ship's motion and the five second averaging interval of
the LORAN-C receiver. These caused the observed TDs to be
several seconds old in comparison to the microwave system
measurements [Schnebele, 1979]. Due tc the large amount.s of
data from the June 1982 survey operations, no deskewing was
done.
The offset, standard error, and drms for data set Number
2 (June 7, 1982) compares well to data set Number 1
(Schnebele's data - all). Also, the 38.6 m 1 drms for the
X-Y rates is between Schnebele's U2.0 m 1 drms prediction
[Schnebele, 1979] and Nelson's findings of 30 m 1 drms in
San Francisco Bay [General Electric 3o. , 1979]. The large
drms of 68.1 m for data set Number 4 (X-Y rates) is probably
due to the phase recovery of the electromagnetic wave from
The 99aO-Y transmitter located in Searchlight, Nevada.
Finally, if ASF Correctors are not applied, drms values
ranged from 38.6 m for 9940 X-Y to 101.4 m for 9940 Y-W for
data sets 2 and 3. The large offsets for all data sets
indicate a systematic error, ranging from 150 m to 700 m,
which precludes the use of zhis method for hydrographic
surveying.
B. SEMI -EMPIRICAL TD GRID
To determine if the Semi-Empirical TD Grid would reduce the
offset and drms for hydrographic surveying, program LOPLC
(Appendix C) was applied to data sets 1 through 4. The
Semi- Empirical TD Grid applies a Secondary Factor and an
Additional Secondary Factor to the primary phase delay based
57
or. the distance overland, the distance over water, and the
total distance using mean land and water conductivities.
The mean offset in meters and microseconds, standard error,
and drms are listed in Table XII. Examples of offsets for
several data points can be found in Appendix C after program
TABLE XII
Semi-Empirical TD Grid Correction Errors
Rates
Mean
Offset
(xz^sec)
Wean
Offset
(111)
Standard
Error
a)
b)
c)
d)
Set
Y
W
Set
Y
H
Set
X
Y
Set
X
Y
(Schnebele's data, 130 data points)
1.131 330.672 0.086
1. 145
630.278
(June 7, 1982 data, 193 data
1.108 324.975
1.105 603,926
(June 6, 1982 data, 128 data
1.131 386.940
0.5 89
(June 3-5,
1.0 50
0.341
172. 923
0. 139
points)
0.077
0. 126
goints)
.061
0.084
1982 data, 169 data points)
365.113 0.052
99.084 0.189
drms
(iS)
113.5
112.7
39.4
63.2
LOPLC. Earlier, Program LOPLC had been tested with data
found in the article by The Analytic Science Corporation
[1979], Program LOPLC produced offsets of 0.86 ;isec and
1.13 }is(^c for the time differences from rates 9940-X and
9940-Y, respectively. As previously stated, the large
offsets may be due to the method by which the distance and
azimuth were computed. The drms, for the data from Monterey
Bay, obtained with the Semi -Empirical TD Grid was similar to
the drms errors for seawater Secondary Factor .lodel. This
might be an indication that Program LOPLC is correct but the
58
bias needs to be adjusted to reduce the large offset between
the observed and calculated TD rates. The mean offset
ranged from 99 to 630 m« Bias reduction could be achieved by
applying land-sea data for the entire West coast to the
Semi-Empirical model. Again, the existence of large offset
values precludes the use of this method for hydrographic
surveying.
C. TABLE ASF CORRECTORS
To determine if the multiple ASF corrections from the
LORAN-C Correction Table would reducs the offset and drms to
meex the NOS accuracy standards, program LORTAB (Appendix D)
was applied to data sets 1 through !4 . ASF Correctors from
TABLE XIII
Multiple LORAN-C Correction Table Errors
Rates
Mean
Offset
(iZ^sec)
Mean
Offset
(iz^)
standard
Error
(sz^sec)
a) Set 1 (Schnebele's data, 130 lata points)
b)
c)
d)
Y
W
Set
Y
Set
X
Y
Set
No
eas
-0.229
0. 169
-^7,068"
92.540
(June 7- 1982 data, 192 data
-0.309 -90.614
0. 123
67. 189
(June 6, 1982 data, 128 data
0.559 19^.134
-0.319
-93.486
0. 104
0. 160
oints)
.085
0, 144
points)
0.073
0. 101
arms
(m)
123.6
1 16.9
47.3
4 (June 3-5, 1982 data, 169 data points)
ASF Corrections listed in Table for the south-
end of Monterey Bay next to the shore line.
the LORAN-C Correction Tables are determined from field
observation of land conductivities by the U.S. Coast Guard
59
Calibration Team [Marine Science Department, 1982]. The
results are listed in Tabls XIII.
When compared to the Saawater Secondary Factor Error in
Table XI, ^he drms values using LORAN-C Corrsction Tables
were increased slightly while the offsets were reduced
substantially. The offsets ranged from 65 - 200 m. Since
the drms ranged from 47.3 to 116.9 m for rates 9940 X-Y and
9940 Y-W respectively, the application of ASF Corrector from
the tables does not meet the NOS accuracy standard. Again,
large offsets and the increase in drms precludes the use of
this method for hydrographic surveying.
D. MULTIPLE OBSERVED CORRECTORS
To determine if multiple observed correctors would
diminish the offset and drms values, mean ASF Correctors
were selected at one minute latitude and longitude in-ervals
using Program ASFSEL. One minute ASF Correctors are shown
at the end of Program ASFSEL in Appendix E. From the one
minute grid, mean ASF Correctors were selected and entered
into Subroutine TABLE at five minute latitude and longitude
intervals in Program L0RTA3 . The following is an example of
the 9940-X ASF Correctors at five mi.iute latitude and longi-
tude intervals for the program:
122/00/00.0 121/55/00.0 121/50/00.0
36/50/00.0 -1.52 -1.63
36/45/00.0 -1.52 -1.56 -1.61
36/40/00.0 -1.58 -1.60
This is the same format used in the LORAN-C Correction
Tables. Program LORTAB was applied to the four data s*=ts.
The error results are listed in Table XIV.
The drms value obtained with multiple observed correc-
tors were all reduced when compared to the irms for the
60
TABLE XIV
Multiple Observed Correct icn Errors
-
Mean Mean Standard
Offset Offset Error
_ _ — _ ^
drms
Rates
1
(X-jisec) l-ffl) (s-^sec)
(Schnebels's data, 130 data points)
0.033 9.585 5.076
0.052 28.604 0.113
(I)
87.6
a)
Set
Y
W
b)
Set
Y
2
(June 7, 1982, 193 data poinds
0.028 8.318 0.073
-0-016 -8. 944 0.116
89.3
c)
Set
X
Y
3
(June 6, 1982, 128 data points)
-0.034 -11.490 0.055
0.010 2.913 0.086
38.3
d)
Set
X
Y
a
(June 3.5, 1982, 169 data points),
0.001 0.383 0.052
-0.008 -2.225 0.187
67.5
seawater Secondary Corrector. The idst: impressive reduction
in drms was within Schnebele's data which covered an area of
seven minutes of latitude and 12 minutes of longi-ude. The
drms for seawater Secondary Correctors was 110.4 m whereas
the drms for the same data using multiple observed correc-
tors was 87.6 m. This is a smaller drms than that of the
June 7, 1932 data (data set Number 2) which was obtained
three years later. It appears that 87.6 m 1 drms is nearly
the minimum error that can be obtained for the 9940 Y-W
rates in Monterey Bay after applying multiple observed AS?
Correctors at five minute latitude and longitude intervals.
For the LORAN-C rates 9940 X-Y, a 38, 3 m 1 drms was obtained
for same five minute area covered by data set Number 2 (June
7, 1982). LORAN-C rates 9940 X-Y were not obtained for the
same size area covered by data set Number 1 (seven minutes
of latitude and twelve minutes of longitude - 9940 Y-W) due
61
to the restriction of ship time and the length of time the
Racal-DECCA Trisponder alrctronic eqaipinent had been loaned.
The drms for the June 3-5, 1982 inshore data was only
reduced to 67.5 m from 68.1 m for seawater Secondary
Correctors. The small change in error at the coast is
probably a result of the erratic behavior of phase recovery
discussed earlier in Chapter Two. The drms value could
presumably be reduced if the correctors were applied at one
minute intervals. This would be a very costly method of
calibrating Loran-C for hydrographic surveying.
62
V. CONCLUSIONS
It was noted in Chapter One that the smallest scale
routinely used for coastal surveys Ls 1 ;80,000. This yields
an allowable error of 49 tn 1 drms with no systematic errors.
This paper determined whether or not applying multiple
Additional Secondary Factors (ASF) Correctors to LORAN-C
lines of position would reduce the drms sufficiently to meet
the accuracy standards set by the National Ocean Survey.
Three methods of applying multiple ASF Correctors were
tested. The first approach computes the time difference
based on a Semi-Empirical ID Grid. The Semi-Empirical Xodel
produced large offsets in the 9940-W, 99aO-X, and 99aO-Y
time differences. The offsets ranged from 99 to 630 m. The
drms for 9940 X-Y combination was 39.4 m and the drms for
9940 W-Y combination was 102.7 m.
The second method applies AS? Correctors found in the
DMAHTC LORAN-C Correction Tables to LORAN-C lines of
position. The application of the tables reduced the offset
in the LCRAN-C time differences. The offsets were between
67 and 191 m. The drms was increased to 47.3 m for the 9940
X-Y combination and 116.9 m for the 9940 W-Y pair.
The final and most accurate method applies multiple
observed ASF Correctors at five minute latitude and
longitude intervals to LORAN-C lines of position. This
method again reduced the offset ia the time difference.
This offset was between 3 and 12 m for the 9940 X-Y
combination. Part of the offset may have been a result of
the microwave positioning system. Reference is made to
Anderson's [1982] paper (in preparation) on the evaluation of
multiple lines of position.
63
The drms values were also reduced to 38.3 m for the 99U0
X-Y rates and 89.3 m for the 9940 W-Y combination. The 38.3
m 1 drms can be decreased by improving the sampling time for
LORAN-C receivers. Nelson obtained 30 m 1 drms for the 9940
X-I rates with special LORAN-C squipment used in San
Francisco Bay, California [General Electric Co., 1979].
Improving the sampling time for LORAN-C receivers used as
positioning equipment for hydrographic surveys should be
investigated.
With drms values of 38.3 m with the possibility of
obtaining 30.0 m 1 drms and offsets ranging from 3 to 12 m,
it may be possible to use LORAN-C for hydrographic surveys
at scales of 1:80,000 or less using multiple observed ASF
Correctors. The use of DMAHTC LORAN-C Correction Tables
should not be ignored. After updating these ASF Correctors
with observed data, the LORAN-C Correction Tables may allow
LORAN-C to be used as a positioning system for hydrographic
surveys.
64
APPENDIX A
PROGRAM GPBYLQ
C PROGRAM GPBYLQ
C
C GENERAL PROGRAM FOR DETERMINING GP FROM KNOWN STATION
C POSITIONS AND THE DISTANCES FROM THEM USING LEAST
C SQUARES. GP TO XY AND XY TO GP ARE DETERMINED 3Y
C SUBROUTINES GPTOXY AND XYTOG? WHICH ARE BASED ON THE
C MODIFIED TRANSVERSE MSRCATOR PROJECTION (MTM) .
C
C PROGRAMMED BY GERALD E. WHEAION, LT. NOAA
C
C LEAST SQUARES ADJUSTMENT BY PAUL R. WOLF, PH.D. AND
C REVISED BY D. LEATH, LCDR
C
C PROGRAM INPUT VARIABLE NAMES
C TITLE = ANY JOB IDENTIFICATION NAMES OR NUMBERS
C M AND N = THE NUMBER OF EQUATIONS (M) AND UNKNOWNS (N)
C XO AND YO = BEST QUESTIMATE OF THE POSTION
C STA(I,1) = X COORDINATE OF KNOW STATION
C STA(I,2) = Y COORDINATE OF KNOW STATION
C STA(I,3) = DISTANCE FROM KNOW STATION
C A (I, J) = THE COEFFICIENT MATRIX
C EL (I, J) = THE CONSTANT MATRIX
C QLL{I,J) = THE WEIGHT MATRIX (WEIGHTS ARE ENTERED AS
C 1 • S IF THE SOLUTION IS EQUALLY WEIGHTED)
C
IMPLICIT REAL*8 (A-H,0-Z)
DIMENSION XCORD (30) ,Y CORD (30) ,WT (3 0) ,ISNO(30)
COMMON /ISTAT/ STA (30 , 3) , FXY (30)
COMMON /LSQX/ A (30 ,30 ) , EL ( 30 , 1) , QLL ( 30 , 30) , AT (30 , 30) ,
65
1AQ(30,30), QXX(30,30) , AQL(30,10), X(30,1), V(30,1),
2VAE(30) , TITLE(80)
C
C READ AND WRITE OaTPOT TITLE
C
WRITS (6,509)
READ (5,710) TITLE
WRITE (6,710) TITLE
C
C DEFINE NUMBER OF SI3NALS (NOT GREATER THAN 30) AND
C NUMBER OF DATA SETS.
C
READ (5,502) NSIG,NDATA
502 FORMAT(I5,I10)
C
C DEFINE THE CENTRAL MERIDIAN (CMSR) , FALSE EASTING
C (FEST), AND CENTRAL LATITUDE (CLAT) .
C CMER IS EXPRESSED IN DEGREES, MINUTES, AND SECONDS.
C' FEST IS THE X-COORDINATS THAI IS ASSIGNED TO THE
C CENTRAL MERIDIAN AND IS EXPRESSED IN METERS.
C CLAT IS DEFINED AS THE CONTROLLING LATITUDE.
C IT IS USED TO REFERENCE THE Y-COORDIN ATES AND
C IS EXPRESSED IN METERS.
C
READ (5,503) ILO NC, ILM INC, RLSECC , FEST , CLAT
503 F0eMAT(1X,I3, 1X,I2,1X,F8. 5, 1X,?7. 1,F10.1)
CMER = ((IA3S (ILONC) * 60 + ILMINC) * 60) + RLSECC
C
C DEFINE VARIABLE FOR:
C NUMBER OF EQUATIONS (M)
C NUMBER OF UNKNOWNS (N)
C IPAGE = NUMBER OF LINES PER PAGE.
C
IPAGE = 1
66
M = 4
N = 2
C
C NULL WEIGHTS
C
DO 3 I a= 1,30,1
DO 3 JM=1,30,1
3 QLL(IM,JM) =0.0
C
C R2AD STATION NUMBERS , POSITION, AND WEIGHTS.
C CONVERT POSITIONS (3P) TO SECONDS AND THEN TO XY.
C
DO 12 J=1,NSIG,1
READ (5,80 0) I SNO (J) ,IL AT ,1 MIN, RS EC, JLON, JMIN , SSEC, WT (J)
800 FORMAT(1X,I3,IU,I3,F9 .5,I5,I3,F9.5,F5. 1)
C
RMAST = ((lABS(ILAT) * 60 + IilIN) * 60) + RSEC
RMASTL = ((lABS(JLON) * 60 + JMIN) * 60) + SSEC
C
CALL GPTOXY (RMA ST, RMA STL, XMETER, YMETER, FEST , CLAT , CMER)
XCORD(J) = XMETER
12 YCORD(J) = YMETER
C
C READ DATA (STATION NUMBERS AND THE DISTANCES)
C
15 DO 40 JC0UNT=1 ,NDATA ,1
READ(5,805) I F, IS, IT , lU ,NREC
805 FORM AT (515)
DO 16 IC0UNT=1 ,NSI3, 1
IF (ISNO (ICOUNT) .EQ. IF) IF=ICOUNT
IF (ISNO (ICOONT) .EQ. IS) IS = ICOUNT
IF (ISNO (ICOUNT) ,EQ. IT) IT=ICOUNT
16 I? (ISNO (ICOUNT) .EQ. 14) I4 = IC0UNT
C
67
STA(1,1) = XCORD(IF)
STA(1,2) = YCORD(IF)
QLL(1,1) = WT(IF)
C
STA(2,1) = XCOFD(IS)
STA(2,2) = YCORD(IS)
QLL(2,2) = WT(IS)
C
STA(3,1) = XCORD(IT)
STA(3,2) = YCORD(IT)
QLL(3,3) = WT(IT)
C
STA(a,1) = XCORD(IU)
STA(4,2) = YCORD(IU)
QLL(4,4) = WT(I4)
C
C READ THE DISTANCE RECORD AND LORAN HATS.
C
DO 38 KC00NT=1, NREC,1
READ (5,507) STA (1, 3) , STA(2,3) ,STA (3,3) ,STA (4,3) ,
1 RATE1,RATE2
507 FORMAT(6F10, 1)
C
C DETERMINE BEST GUESS COORDINATES XO AND YO
C WITH SUBROUTINE GUESS.
C
19 CALL GUESS (XO,YC)
C
IJUMP = 0
C
C COMPUTE FXY, A AND L MATRIX
C
20 DO 25 1=1, M
FXY(I) = DSQRT (DA3S( ( XO-STA (I, 1) ) **2
68
1 (Y0-STA(I,2) )**2) )
A (1,1) = (X0-STA(I,1) ) / (FXY(I))
A(I,2) = (Y0-STA(I,2) ) / (FXY (I) )
25 EL (1,1) = STA(I,3) - FXY(I)
C
C CALL SUBROUTINE LSQR
C
CALL LSQR (K,N)
C
C COMPUTE THE NEW QUESSTIMATE FOR XO AND YO
C
XO = XO + X (1 ,1)
YO = YO + X (2,1)
C
C EXIT IF STANDARD ARE MET USING IJUMP OR
C XO AND YO CUT OFF
C
IF (DABS (X (1,1) ) .LE.1. 00 .AND. D AB3 (X (2, 1) ) . LE. 1 . 00)
1 GO TO 35
IF (IJUMP .EQ. 10) GO TO 35
30 GO TO 20
C
C COMPUTE ERROR ELLIPSE
C
35 CONTINUE
CALL ELIPSE (SU, SV)
C
C CONVERT XY TO GP
C
CALL XYTOGP (XO, YO, SEC lAT, SECLON , FEST , CLAT, CMER)
CALL TODMS (SECL AT, IDE GP,I MINP, RSECP)
CALL TODMS (SECLON, JDE GP ,JMINP , 3SSC?)
C
C PAGE AND CONTINUE WITH NEXT SET OF OBSERVATIONS
69
c
WRITE (6,505) IDEGP ,IM INP,RSEC?, JDEGP, JMINP , SSECP ,
1 EATE1,RATE2
C WRITE(6,505) IDEGP ,IMINP, RSECP, JDEGP, JMINP, SSECP,
1 sa,sv
505 FORMAT (14, 13, F6. 2, 15, I3,F6 . 2, 2F9 . 2)
C IF(IPAGE .EQ. 50) WRITE (6,509)
IF(IPAGE .EQ. 50) IPAGE = 0
38 IPAGE = IPAGE + 1
C
C FORMAT STATEMENTS
C
710 FORMAT(80A1)
509 FORMAT (1H1)
UO CONTINUE
STOP
END
C======== = ======= = == ====== ====== = == = = = == ====== = = == = ======:
SUBROUTINE LSQR (M,N)
Q
IMPLICIT REAL*8 (A-H,0-Z)
COMMON ASQX/ A (30,30) ,SL (30, 1) ,QLL (30,30) , AT(30,30) ,
1AQ(30,30), QXX(30,30) , AQL(30,10), X(30,1), 7(30,1),
2VAR(30) , TITLE (80)
C
C COMPUTE A TRANSPOSE 3Y TRANSPOSING THE A MATRIX (AT)
C
DO 6 1 I=1,M
DO 61 J=1,N
61 AT (J, I) =A (I, J)
C
C USING STEPS (1) , (2) , AND (3) COMPUTE THE INVERSE
C OF THE TRANSPOSE (AT) * WEIGHTED MATRIC (QLL) *
C MATRIX A = QXX.
70
c
C (1) COMPOTE AQ = AT * QLL
C
DO 7 1 1 = 1 ,N
DO 71 J=1 ,M
AQ (I,J)=0.
DO 71 K = 1 ,M
71 AQ (I,J)=AQ(I,J) +(AT(I,K)*QLL(K, J) )
C
C (2) COMPUTE QXX = AQ * A
C
DO 81 1 = 1, N
DO 81 J=1,N
QXX(I,J)=0.
DO 81 K=1,M
81 QXX(I,J)=QXX (I, J) +AQ (I,K) *A (K,J)
C
C (3) INVEST QXX MATRIX
C
DO 3 07 K=1 ,N
DO 302 J=1,N
IF (J-K) 304,302, 304
304 QXX(K,J)=QXX (K, J)/QXX (K,K)
302 CONTINUE
QXX(K,K)=1./QXX (K,K)
DO 307 1=1, N
IF (I-K) 305,307,305
3 05 DO 3 03 J=1,N
IF (J-K) 306,303, 306
306 QXX(I,J)=QXX (I, J) -QXX (I,K) *QXX (K, J)
303 CONTINUE
QXX(I,K)=-QXX (I,K) *QXX (K, K)
307 CONTINUE
71
C USING STEPS (4) AND (5), COMPOTE THE UNKNOWNS X
C BY MULT THE INVERSE Q XX AND AQL.
C
C (U) COMPUTE AQL = A Q * EL
C
DO 101 1=1, N
AQL(I,1)=0.
DO 101 K= 1 , M
101 AQL(I,1)=AQL (I, 1) +A3 (I,K) *EL(K, 1)
C
C (5) COMPUTE X = QXX * AQL
DO 201 1=1, N
C
X (I, 1)=0.
DO 201 K=1,N
201 X (I, 1) =X(I,1) +QXX(I,K) *AQL (K, 1)
C
C (6) COMPUTE THE RESIDUAL (V = A * X -EL)
C
DO 301 1=1, M
V{I,1)=0.
DO 3 01 K=1,N
301 7 (I, 1)=7(I, 1) +A (I,K) *X(K, 1)
DO 1 1=1, M
1 V(I,1)=V(I,1)-EL(I,1)
C
C COMPUTE THE STANDARD DEVIATION OF UNIT WEIGHT SIGMA
C DM - NUMBER OF OBSERVATIONS
C DN - NUMBER OF KNKNOWNS
C
SIGMA=0.
DM=M
DN=N
DO 332 1=1, M
72
382 SIGMA=SIGMA+V (1,1) **2*QLL (1,1)
SIGMA=DSQRT (SIGMA/(D M-DM) )
C
C COMPUTE THE STANDARD DEVIATION OF THE ADJUSTED UNKNOWNS
C QXX - ARE THE ELEMENTS OF THE COVARIANCE MATRIX.
C
DO 446 1=1, N
446 VAR(I) =DSQRT (QXX(I ,1) *SIGMA**2)
C
5 10 CONTINUE
RETURN
END
C==r== ===== = == ============== = = ======== = = = = == = = = === = == = = ====== =
SUBROUTINE ELIP SE (SU, SV)
C
c
C SOLVE FOR THE SEMI MAJOR AND SEMI MINOR AXIS OF
C THE ERROR ELLIPSE
C
IMPLICIT REAL*8 (A-H,0-Z)
COMMON /LSQX/ A (30 ,30 ) , EL ( 30, 1) , QLL (30 , 30) , AT (30 , 30) ,
1AQ(30,30), QXX(30,30) , AQL(30,10), X(30,1), V(30,1),
2VAR(30) , TITLE(80)
C
SUS = .5*(QXX(1,1) + QXX(2,2) + DSQRT (DABS (QXX ( 1 , 1) -
1 QXX (2,2) +4.0*QXX(1,2) *QXX (2, 1) ) ) )
SVS = .5*(QXX(1,1) + QXX(2,2) - DSQRT (DABS (QXX (1 , 1) -
1 QXX (2,2) +4.0*QXX(1,2)*QXX (2,1) ) ) )
C
SU = DSQRT (SUS)
SV = DSQRT (SVS)
C
RETURN
END
73
C==:==s = =: == = ===== = = = = ======= ===== = = = = = = = === = == = == = = == = = = = == =
SD3R0UTINE G0ESS(XO,YO)
c
C SUBROUTINE GUESS DETERMINES THE BEST GUESS COORDINATES
C TO BE USED IN SUBROUTINE LSQR. USE RIGHT SIDE RULE FOR
C STATION ORDER.
C
IMPLICIT REAL*8 (A-H,0-Z)
COMMON /ISTAT/ STA (30 , 3) , FXY (30)
C
C DETERMINE DISTANCE BETWEEN STATION 1 AND STATION 2
C
D = DSQRT ((STA (2, 1) -STA (1,1) ) **2+ (STA (2,2)-3TA (1,2) )**2)
C
C DETERMINE ANGLE ALPHA BETWEEN XO , Y0/STA2/STA1
C
ALPHA=DARC0S ( (STA(2,3) **2-STA (1,3) **2 + D**2) /
1 (2.0*STA (2,3)*D) )
C
C DETERMINE ANGLE BROVO BETWEEN X-AXIS AND STA2-STA1
C
BROVO=DARSIN ( (STA (1,2)-STA (2,2) ) /D)
C
C DETERMINE X AND Y LENGTH
C
X=STA (2,3) *DCOS (ALPHA +BROVO)
Y = STA (2,3) *DSIN (ALPHA +BROVO)
C
C DETERMINE XO AND YO
C
XO = STA(2, 1) +X
YO = STA(2,2) ■^Y
in
RETURN
END
SUBROUTINE GPTOXY(SEC LAT, SECLON,XCO, YCO^FEST ,CLAT,CMER)
IMPLICIT REAL*8 (A-H,0-Z)
DATA E2,RK0,A / .00676 8658 DO ,. 99998D0 , 6378206 .UDO/
DATA RKGE0,W1,W2 /O. 0 4848 1 368D0 , 0. 11 422D0, 2 1 . 73607D0/
DATA W3, W 4/5 104. 57 33 8 DO, 63 673 99. 6 3 9D0/
DATA RADSEC /.O 000048 48 13681 1 1D0/
C
C
RADLAT = SECLAT * RADSEC
SINLAT = DSIN (RADLAT)
SIN2LA = SINLAT * SINLAT
COSLAT = DCOS (RADLAT)
C0S2LA = COSLAT * COSLAT
P = (CMER - SECLON) / 10000. ODO
V = A / DSQRT(1.0D0 - E2*SIN2LA)
TANCON = 1.0D0 - SIN2 LA/C0S2LA
S = W4 * (RADLAT - SI NLAT*C0SLAr/1 0. 0DO**6 *
* (W3 -COS2LA*(W2- W1*C0S2LA) ) )
T1 = S * RKO - CLAT
T2 = RKGEO * COSLAT * RKO * V
T3 = T2 * RKGEO / 2. 0 DO
T4 = T3 * SINLAT
T5 = T3 * RKGEO * COS 2LA / 3. ODO
T6 = T5 * TANCON
T7 = (4. ODO + TANCON) *T5 *RKGEO * SINLAT / 4, ODO
XCO = (T2 + (T6*P**2))*P + ?EST
YCO = (T7*p**4) + (T4*P**2) > T1
RETURN
END
75
SUBROUTINE XYTOGP(XCO , YCO, SECLAT, SECLON,FEST,CLAT,CMER)
IMPLICIT REAL*8 (A-H,0-Z)
C
DATA E2,A,SR /O .00676 86 58D0 ,6378206. UDO ,
1 0.0000048U81368D0/
DATA W1,W2,W3 /O. 2468 2D0, 30. 0233 5D0, 5078 .64 97700/
C
D = CLAT ♦ YCO
WO = 0. 15704998 1D0/10.0D0**6 * D
SINWO = DSIN(WO)
COSWO = DCOS (WO)
C0S2W0 = COSWO * COSWO
PHI1 = WO + SINWO*COS WO/1 0.0D0**6 *
* (W3-^COS2WO*{W2+W 1*C0S2W0) )
PHI2 = PHI1 / 0.99998D0
PHI3 = PHI2 / SR
Q = (XCO - FEST) / 10.0D0**6
V = A / DSQRT{1.0D0 - E2*DSIN {PHI2) **2)
T = DC0S(PHI2) * SR
C = V * 0.99998D0
T1 = 10.0D0**6 / (T*C)
T2 = (T1*10. 0DO**6) / (2.0D0*C)
T3 = (T2*10.0D0**6) / (3.0D0*C)
T4 = (T3*10.0D0**6) / (4.0D0*C)
DELLON = (T1-Q**2*T3* (2.0D0*DTAN (PHI2) **2+1 .ODO) ) *Q
SECLAT = ( ( (3.0D0*DTAN(PHI2) **2 + 5.0D0) *T4) *Q**2-T2)
1 *Q**2*DSIN (PHI2) + PHI3
SECLON = CMER - DELLON
RETURN
END
SENTRY
76
DATA SET EXAMPLE - PROGRAM GPBYLQ.
SHIP'S POSITION OBSERVED LORAN RATES
LATITUDE
(D-M-S)
LONGITUDE
(D-M-S)
9940-X
(fisec)
9940-Y
3 6 36 40.20
121 52 48,62
27508.79
42742.71
36 36 42.41
121 52 48.16
27508.81
42742.93
3 6 36 44.56
121 52 47.52
27508.83
42743.15
3 6 36 46.84
121 52 46.98
27508.76
42743. 26
36 36 49.09
121 52 46.5 9
27508.76
42743.48
77
APPENDIX B
PROGRAM LOR AN
C PROGRAM LORAN
C
C PROGRAM COMPOTES LINE OF POSITIONS FOR LORAN-C USING
C THE SECONDARY FACTOR (SF) BASED ON SEAWATER EM MODEL.
C
IMPLICIT REAL*8 (A-H,0-Z)
DIMENSION AXIS( 13) ,RF (13)
DIMENSION XSLAV1 (100 0) ,XSLAV2 (1000)
DATA RHOSEC,PI,UNCOV, RN/2. 06264 8062U71D05,
1 3. 1 4 15 92 6535 89 8D0, 2 9 9.79 245 800,1 .000338DO/
DATA XMEAN1,XMEAN2,VAR1,7AR2/O.OODO,O.OODO,
1 0.0ODO,00.0D0/
DATA XMEAN3,XMEAN4/0. 0ODO,O.00DO/
C
CCCCCC* ******************** *************
DATA AXIS/6. 3782064D06, 6. 378388 D06 , 6. 377397 1 55D06 ,
1 6.37816 D06, 6. 37816D06, 6, 378249145006,6. 378165D06,
2 6.378166006,6.378165 D06 , 6 . 378 1 45D06 , 6. 3775634D06,
3 6.378245006,6.3781350006/
DATA RF/6. 3565838 006,2.97 002,2.9915 28 12 85002,
1 2.98 2500 2,2. 98 2 4716 7 4270 02,2.9 3 46 50 02,2.9 8 25002,
2 2.983 002,2.983002,2-9825 002,6.3562569006,
3 2.933 002,2.9826002/
C*****ELLIPSOID OPTION NUMBER
C 1. CLARKE 1866 8. MERCURY
C 2. INTERNATIONAL (HAYFORD) 9. MARSHALL ISLAND
C 3. BESSEL 1841 10. NAVY 8D
C 4. AND (AUSTRALIAN) 11. AIRY
C 5, 1967 REFERENCE 12. KRASSOWSKI 1940
78
C 6. CLARKE 1880 MOD 13. WGS 1972
C 7. SAO
C
cccccc* ************ ******** *************
C CC1-2 = ELLIPSOID NUMBER (K)
C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC)
C
READ (5, 100) K,IREC
100 F0RMAT(lX,l2,ia)
TW0PI=2.*PI
A = AXIS(K)
F=1./RF (K)
IF (F.LT.3.D-3) F=(A-1./F)/A
C
C READ MASTER AND SLAVE STATIONS POSITIONS
C THE FIRST RECORD IS THE NUMBER OF MASTER AND SLAVE
C STATIONS FOR THE PARTICULAR CHAIN.
C
WRITE (6,202)
202 FORMAT(IHI)
WRITE (6,201)
201 FORMAT (1H )
READ (5, 105) ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM
WRITE (6,1 05) ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM
105 FORMAT (IX, 13, IX, 12, IX, F5. 2, IX, 14, IX, 12, IX, F5. 2)
RMAST = ( (lABS(ILATM) * 60 + IMINM) * 60 > RSECM) /
1 RHOSEC
IF(ILATM .LT. 0) RMAST = -RMAST
RMASTL = ( (TABS (ILONM) * 60 + ILMINM) *60 +RLSECM) /
1 RHOSEC
IF (ILONM .GT. 0) RMASTL = TWOPI - RMASTL
C
C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
79
READ (5, 101) ILATS,IMINS,RSECS,rLONS,ILMINS,ELSSCS,
1 DELAY1
WRITE (6,101) ILATS,IMINS,RSECSrILONS,ILMINS,RLSBCS,
1 DELAY1
101 FORMAT (IX, 13, IX ,12 , 1 X ,F5. 2, IX, lU , IX, 12, IX, F5 . 2 , F9 . 2)
RSLAV1 = ((lABS (ILATS) * 60 4- IMINS) * 60 + RSECS) /
1 RHOSEC
IF (ILATS .LT. 0) RSLAV1 = -RSLAV1
RSLAL1 = ((lABS (ILONS) * 60 + ILMINS) * 60 ■»• RLSECS) /
1 RHOSEC
IF (ILONS .GT- 0) RSLAL1 = TWOPI - RSLAL1
C
C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
C
READ (5, 101) ILATS, IMINS, RSECS, ILONS, ILMINS, RLSECS,
1 DELAY2
WRITE (6,101) ILATS, IMINS, RSECS, ILONS, ILMINS, RLSECS,
1 DELAY2
RSLAV2 = ((TABS (ILATS) * 60 ^ IMINS) * 60 ■•• RSECS) /
1 RHOSEC
IF (ILATS .LT. 0) RSLA V2 = -RSLA72
RSLAL2 = ( (IA5S (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2
WRITE (6,201)
C
C READ THE RECORDS POSITIONS AND OBSERVED RATES, THEN
C COMPUTE THE TOTAL DISTANCE.
C
DO 550 I=1,IREC,1
READ (5, 104) ILAT,IMIN ,RSEC ,ILON , ILMIN, RLSEC,
1 RATE1,RATS2
30
104 FOEMAT(1X,I3,I3,F6.3, 15,13, F6 . 3 , 3X, 2F1 0 .2)
C
C CONVERT LAT AND LONG TO RADIANS.
C
RPOST = ((lABS(ILAT) * 60 + IMIN) * 60 + RSEC) /
1 RHOSEC
IF(ILAT .LT. 0) RPOST = -RPOST
RPOSTL = ((lABS(ILON) * 60 + ILMIN) * 60 + RLSEC) /
1 RHOSEC
IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL
C
C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT
C
C MASTER
C
CALL INVER1 (A, F, RPOST , RPOSTL, RM AS T, RMASTL, FAZM , BAZM ,
1 DISTM)
DDISTM = (RN * DISTM) / UNCOV
CALL SECFAC {UDISTM,SFM)
C
C SLAVE1
C
CALL INVER1 (A,F ,RPOST , RPO STL,RSLA V1 , RSLAL1 , F AZ 1 , SAZ 1 ,
1 DIST1)
UDIST1 = (RN * DIST1) / UNCOV
CALL SECFAC (UDI ST1, SF 1)
C
C SLAVE2
C
CALL INVERT (A, F ,R?OSr , RPO STL,RSLAV2, ESLAL2, F AZ2, BAZ 2,
1 DIST2)
UDIST2 = (RN * DIST2) / UNCOV
CALL SECFAC (UDIST2,SF2)
81
C COMPOTE THE RATES AND COMPARE TO THE OBSERVED RATES
C
TDM1 = UDIST1 - ODISTM + SF1 - SFM ■•■ DELAYl
DIPF1 = TDM1 - EAT El
C
TDM2 = UDIST2 - UDIST M + SF2 - SFM + DSLAY2
DIFF2 = TDM2 - RATE2
C
C COMPUTE THE LANE WIDTH IN METERS BASED ON EQUATION
C 4.20 IN ELECTRONIC SURVEYING AND NAVIGATION -
C LAURILA, PAGE 94.
C
BR1 = DABS (FAZ1 - FAZ M)
3R2 = DABS(FAZ2 - FAZ M)
WIDTH1 = (DIFF1 * UNCOV * 0.5) / DSIN (BR1 * 0.5)
WIDTH2 = {DIFF2 * UNCOV * 0.5) / DSIN (BR2 * 0.5)
C
C WRITE THE POSITION OF VESSEL, COMPUTED RATE, OBSERV^ID
C RATE, AND THE DIFF BETWEEN THEM.
C
WRITE (6^200) ILAT,IMIN,RSEC,IL0N,ILMIN,RL3EC,RATS1,
DIFF1,RATE2,DIFF2
200 F0RMAT(1X,I3,I3,1X,F6.3,I5,I3,1X^F6.3,F12.2,F8.2,
1 F12.2,F8.2)
WRITE(6,201)
C
C XMEAN1 AND XMEAN2 ARE THE MEAN DIFFS BETWEEN THE
C COMPUTED RATS AND THE OBERSEHVED. XSLAV1 AND XSLAV2
C ARE THE STORED DIFFS.
C
XMEAN1 = XMSAN1 + DIFF1
XMEAN2 = XMEAN2 ■•■ DIFF2
XMEAN3 = XMEAN3 * WIDTHi
XMEAN4 = XMEAN4 + WIDTH2
82
XSLAV1 (I) = DIFF1
XSLAV2(I) = DIFF2
550 CONTINUE
C
C COMPUTE THE MIAN AND STANDARD DEVIATION
C
XMEAN1 = XMEAN1 / IREC
XMEAN2 = XMEAN2 / IREC
XMEAN3 = XMEAN3 / IREC
XMEAN4 = XMEAN4 / IREC
C
DO 6 00 I = 1, IREC, 1
VAR1 = 7AR1 + ((XSLA7 1(I) - XMSAN1)**2)
600 7AR2 = VAR2 + ((XSLA72(I) - XMEAN2)**2)
C
VAR1 = VAR1 / (IREC - 1.0)
VAR2 = VAR2 / (IREC - 1.0)
C
SD1 = DSQRT(VARI)
SD2 = DSQRT(VAR2)
C
WRITE (6,201)
WRITE (6,210) XMEAN1,SD1,XMEAN3
WRITE(6,201)
WRITE(6,211) XMSAN2,S D2,XMEANa
WRITE (6,202)
210 F0RMAT(1X,' SLAVE #1, MEAN = ',?10.3,
1 • STANDARD DEVIATION = • ,
2 F10.3,' DISTANCE IN METERS = ',F10.3)
211 F0RMAT(1X,' SLAVE #2, MEAN = ',F10.3,
1 ' STANDARD DEVIATION = ',
2 F10.3,' DISTANCE IN METERS = •,F10.3)
STOP
END
83
C=====: = = == = == == = = = = ======== == = = = = = = ==== = == = = = = = = = = = = == = == =
SUBROUTINE SECFAC (UTDIST, SF)
C
c
C THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR
C OTDIST = TOTAL DISTANCE
C SF = SECONDARY FACTOR
C
IMPLICIT REAL*8 (A-H, 0-Z)
C
c
C COEFFICIENTS
C
AO = 129.04323
A1 = -0.40758
A2 = 0.00064576813
C
BO = 2.741282
B1 = -0.011402
B2 = .00032774815
C
IF(UTDIST .GT. 537) 30 TO 10
SF = ( BO / UTDIST) + B1 + ( B2 * UTDIST)
GO TO 20
10 SF = ( AO / UTDIST) •»• A1 + ( A2 * UTDIST)
C
20 CONTINUE
RETURN
END
C== ======================== ================================:
SUBROUTINE INVSR1 (A, ? INV, GLAT1 , GLON 1 , GLAT2 , GL0N2 ,FAZ ,
1 BAZ,S)
C
c
34
Q 4e«:(c SOLUTION OF THE GEODETIC INVERSE PROBLEM AFTER
C *** T.VINCENTY MODIFIED RAINS FORD'S METHOD WITH HELMERT'S
C *** ELLIPTICAL TERMS. EFFECTIVE IN ANY AZIMUTH AND AT
C *** ANY DISTANCE SHORT OF ANTIPODAL STANDPOINT/FOREPOINT
C *** MUST NOT BE THE GEOGRAPHIC POLE
C
C *** A IS THE SEMI-MAJOR AXIS OF THE REFERENCE ELLIPSOID
C *** FINV IS THE FLATTENING (NOT RECIPROCAL) OF THE
C *** REFERNECE ELLIPSOID LATITUDES AND LONGITUDES IN
C *** RADIANS POSITIVE NORTH AND EAST FORWARD AZIMUTHS AT
C *** BOTH POINTS RETURNED IN RADIANS FROM NORTH GEODESIC
C *** DISTANCE S RETURNED IN UNITS OF SEMI-MAJOR AXIS A
C
C *** PROGRAMMED FOR CDC-66 00 BY LCDR L.PFSIFER NGS
C *** ROCKVILLE MD 18FEB75. MODIFIED FOR IBM SYSTEM 360
C *** BY JOHN G 3ERGEN NGS ROCKVILLE MD 7507.
C
IMPLICIT REAL*8 (A-H,0-Z)
DATA EPS/0. 5D-1 3/, PI/ 3. 14 1 59265 35 89 8D0/
TW0PI=2.*PI
R=1.-FINV
TU1 = R*DSIN(GLAT1) /DCOS (GLATI)
TU2=R*DSIN (GLAT2) /DCOS (GLAT2)
CU1=1 ./D3QRT(TU1*TU1+ 1.)
SD1=CU1*Tai
CU2=1 ./DSQRT(TU2*TU2+ 1.)
S=CU1*CU2
3AZ=S*TU2
FAZ=BAZ*TU1
X=GL0N2-GL0N1
100 SX=DSIN (X)
CX=DC0S (X)
TU1=CU2*SX
TU2=BAZ-SU1*CU2*CX
85
SY=DSQRT(Tai*T0 1+TU2*TU2)
CY=S*CX+FAZ
Y = DATAN2(SY,CY)
SA=S*SX/SY
C2A=-SA*SA-H.
CZ=FAZ+FAZ
IF(C2A.GT.O.) CZ=-CZ/C2A+CY
E=CZ*CZ*2.-1.
C= ((-3.*C2A + 4 .) *FINV4-a.) *C2A*FINV/16 .
D = X
X= ((E*CY*C + CZ) *SY*C-«-Y) *SA
X= (1 .-C) *X*FINV+GL0N2-GL0N1
IF (DABS (D-X) .GT.EPS) GOTO 100
FAZ=DATAN2(T01,TU2)
IF (FAZ.GE.TWOPI) FAZ= FAZ-T WOPI
IF (FAZ.LT.O.DO) FA2=F AZ+TWOPI
BAZ=DATAN2(CU1*SX,BAZ*CX-S01*Ca2) +PI
IF (BAZ.GE.THOPI) B AZ= BAZ-T WOPI
IF (BAZ.LT.O.DO) BAZ=B AZ+IWOPI
X = DSQRT ( (1./R/R-1.) *C2A+1.) +1.
X= (X-2.)/X
C=1.-X
C= (X*X/U. + 1.) /c
D=(0.375*X*X-1. ) *X
X=E*CY
S=1.-E-E
S= (( ( (SY*SY*4.-3.) *S*CZ*D/6.-X) *D/4.+CZ) *SY*D + Y) *C*A*R
RETURN
END
SENTRY
86
DATA SET EXAMPLE - PEOGRAM LORAN
SHIP'S POSITION OBSERVED LORAN RATES & ERRORS (E)
LATITUDE LONGITUDE 9940-Y EY 9940-W EW
(D-M-S) (D-M-S) ()1S6C) (;isec) (;isec) (usee)
36 U3 45.800 121 55 27.160 42789.31 -0.49 16294.04 -1.06
36 44 3.400 121 55 32.340 42791.13 -0.38 16293.46 -1.10
36 44 21.180 121 55 37.390 42793.04 -0.38 16292.73 -0.99
36 44 37.490 121 55 46.950 42795.13 -0.58 16292.03 -1.06
36 44 53.260 121 55 57.710 42796.93 -0.51 16291.43 -1.27
87
APPENDIX C
PROGRAM LOPLC
C PROGRAM LOPLC
C
C PROGRAM COMPUTES LINE OF POSITIONS FOR LORAN-C USING
C THE SECONDARY FACTOR (SF) AND ADDITIONAL SECONDARY
C FACTOR (ASF) THE SF AND ASF ARE BASED ON
C SEMI-EMPIRICAL TD GRID.
C
IMPLICIT REAL*8 (A-H,0-Z)
DIMENSION PHI (100) ,ELON (100) ,DISTM (100) ,DIST1 (100) ,
1 DIST2(10 0) ,AXIS (13) , RF (13) ,AZ1 (100) , AZ2 (1 00) , AZM (1 00)
DATA RHOSEC,PI, UNCOV/2. 06 264806 2U7 1D05,
1 3.1U15926535898D0,299.693D0/
C
cccccc* ****** ******** ****** *************
DATA AXIS/6. 3782064D06, 6. 378388 D06, 6. 377 397 1 55D06 ,
1 6.37816 D06,6. 37816D06^6.3782U91U5D06,6.378165D06,
2 6.378166D06,6. 378165 D06 , 6. 378 1 a5D06 , 6 . 377563UD06,
3 6.378245D06,6. 3781350D06/
DATA RF/6. 3565838 D06,2.97 D 02, 2 . 99 15 28 1285D02 ,
1 2.98 25D0 2,2. 98 247 167 4 27D 02, 2 .93 465D02 , 2 .9825D02 ,
2 2.983 D02,2.983D02,2 .9825 D02, 6 , 3562569D06 ,
3 2.983 D02,2.98 26D02/
C*****ELLIPSOID OPTION NUMBER
8. MERCURY
9. MARSHALL ISLAND
10. NAVY 8D
11. AIRY
12. KRASSOWSKI 19U0
13. WGS 1972
88
C
1.
CLARKE 1866
c
2.
INTERNATIONAL (HAYFORD)
c
3.
BESSEL 1841
c
4.
AND (AUSTRALIAN)
c
5.
1967 REFERENCE
c
6.
CLARKE 1830 MOD
C 7. SAO
C
cccccc* ************ ******** *************
C CC1-2 = ELLIPSOID NUMBER (K)
C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC)
C
READ(5,100) K,IREC,ISTA1,ISTA2
100 FOEMAT(1X,I2,I4 ,213)
TW0PI=2.*PI
A = AXIS(K)
F=1./RF (K)
IF (F.LT-3.D-3) F=(A-1./F)/A
C
C READ MASTER STATION POSITIONS.
C
READ (5, 105) ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM
105 FORMAT (1X,I3, 1X,I2,1X,F5.2,1X,I4,1X,I2,1X,F5.2)
RMAST = ( (IA3S (ILATM) * 60 + IMINM) * 60 + RSSCM) /
1 RHOSEC
IF (ILATM .LT. 0) RMAST = -RMAST
RMASTL = ((TABS (ILONM) * 60 + ILMINM) *60 +RLSECM) /
1 RHOSEC
IF(ILONM .GT. 0) RMASTL = TWOPI - RMASTL
C
C READ THE FIRST SLAVE STATION AND CHANGS THE LATITUDE
C AND LONGITUDE INTO RA DI ANS . DELAY IS THE CODING DELAY
C AND BIAS IS THE OFFSET IN MICROSECONDS.
C
READ (5, 101) ILATS,IMINS,RSSCS,ILONS,ILMINS,RLSECS,
1 DELAY1,3IAS1
101 FCRMAT(1X,I3, 1X,I2,1X,F5.2, 1X,I4, 1X,I2,1X,F5.2,
1 ?9.2,F7.3)
RSLA71 = ( (lABS (ILATS) * 60 + IMINS) * 60 * RSSCS) /
1 RHOSEC
39
IF(ILATS .LT. 0) RSLAV1 = -RSLAV1
RSLAL1 = ( (lABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL1 = TWOPI - RSLAL1
C
C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
C
READ (5, 101) ILATS,IMINS,RSECS, ILONS, ILMINS, RLSECS,
1 DELAY2,BIAS2
RSLAV2 = ((TABS (ILATS) * 60 + IMINS) * 60 + RSECS) /
1 RHOSEC
IF(ILATS .LT. 0) RSLAV2 = -RSLAV2
RSLAL2 = ( (lABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2
C
C READ COAST POINT LAT AND LONG AND CONVERT TO RADIANS.
C STORE LAT IN PHI AND LONG IN ELON.
C
DO 5 00 I=1,IRZC,1
READ (5,102) INUM,ILAT,IMIN,RSSC,ILON,ILMIN, RLSSC
102 FORM AT ( 1 X, 14, 14, 13, F7. 3, 1 5, 13, F 7. 3)
PHI (I) = ((lABS(ILAr) * 60 -<■ I2IIN) * 60 + RSSC) /
1 RHOSEC
IF(ILAT .LT. 0) PHI (I) = -PHI (I)
ELON (I) = ((IABS(IL3N) * 60 + ILMIN) * 60 ^- RLSEC) /
1 RHOSEC
IF(ILON .GT. 0) ELON (I) = TWOPI - ELON (I)
C
C COMPUTE DISTANCE FROM COAST POINT TO MASTER AND
C SLAVE STATIONS. SET UP COMPUTATION AND CALL INVER1
C
?1 = ?HI(I)
90
El = ELON(I)
CALL INVER1 (A, F, Pi, SI ,RMAST,R»1ASTL , AZF , AZB , S)
DISTM(I) = S
AZM(I) = AZF
CALL INVER1 (A,F,P1,E1 ,RSLAV 1 , RSLAL1 , AZF, AZB , S)
DIST1 (I) = S
AZ1(I) = AZF
CALL INVER1 (A,F,P1,E1 ,RSLA72, RSLAL2 , AZF , AZB , S)
DIST2 (I) = S
AZ2(I) = AZF
WRITE (6,210) I,INUM,ILAT,IMIN,RSEC,IICN,ILi5IN, RLSEC,
1 AZM (I) ,AZ1 (I) , AZ2 (I)
219 FORM AT (IX, 14, 14, 14, 1 3 ,F7. 3 , 15 , 13 , F7. 3 , 3F10 . 3)
500 CONTINUE
C
C READ THE NUMBER OF RECORDS AND THEIR POSITIONS.
C COMPUTE THE TOTAL DISTANCE OVER LAND AND SEA.
C TDISTM = TOTAL DISTANCE FROJl FASTER TO DATA POINT.
C TDISTS = TOTAL DISTANCE FROM SLAVE TO DATA POINT.
C
C COMPUTE THE FORWARD AND BACK AZIMUTHS
C FAZM AND FAZS = FORWARD AZIMUTH TO THE MASTER AND
C SLAVE STATION.
C EAZM AND 3AZS = BACK AZIMUTH TO THE MASTER AND
C SLAVE STATION.
C
READ (5, 103) JREC
103 FORM AT (IX, 14)
C
DO 550 1=1, JREC, 1
READ (5, 104) ILAT ,IilIN, RS EC , HON, ILMIN , RLSEC , RATS 1 , HAT E2
104 FORM AT (IX, 13,13, F6. 3, 15,13, F6. 3, F 9.3, F9. 3)
C
C CONVERT LAT AND LONG TO RADIANS.
91
RPOST = ((lABS(ILAT) * 60 + IMIN) * 60 + RSEC) /
1 RHOSEC
IF(ILAT .LT. 0) RPOSr = -RPOST
BPOSTL = ((lABS(ILON) * 60 + ILMIN) * 60 + RLSEC) /
1 RHOSEC
IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL
C
C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED
C POINT TO THE INTERPOLATED SHORE POINT AND WRITE.
C
ISTA = 1
CALL SHORPT (I RE C^RPOS T, RPOSTL, R MAST, RMASTL , UTDISM,
1 USDIST,ULDIST,FAZM,B AZM)
CALL SECFAC(UTDISM,aSDIST,ULDI3T,BAZM,ISTA,SFM)
C
CALL SHORPT (I REC,RPOS T, RP0STL,RSLAV1 , RSLAL1 ,UTDIS1,
1 USDIST,ULDIST,FAZ1,B AZ1)
CALL SECFAC (UTDIS 1 ,as DIST , ULDIST, BAZ1 ,ISTA 1 ,SF1)
C
CALL SHORPT (I RE C,RPOS T, RPOSTL, RSL A72 , RSLAL2 , UTDIS2,
1 USDIST,ULDIST,FAZ2,B AZ2)
CALL SECFAC (UTDIS2 ,as DIST , ULDIST, BAZ2 ,ISTA2 , SF2)
C
C COMPUTE THE RATES AND COMPARE TO THE OBSERVED RATES
C
TDM1 = UTDISI - UTDISM + SF1 - SFM + DELAY1 + BIASI
DIFF1 = TDMl - RATE1
C
TDM2 = UTDIS2 - UTDISM + SF2 - SFM + DELAY2 + EIAS2
DIFF2 = TDM2 - RATE2
C
C COMPUTE THE LANE WIDTH IN METERS BASED ON EQUATION
C 4.20 IN ELECTRONIC SURVEYING AND NAVIGATION - LAUEILA,
92
C PAGE 94.
C
BR1 = DABS (FAZ1 - FAZ M)
IF(BR1 ,GT. PI) BRI = TWOPI - BR1
BR2 = DABS (FAZ2 - FAZ M)
IF(BR2 .GT. PI) BR2 = TWOPI - BR2
WIDTH1 = (DIFF1 * UNCOV * 0.5) / DSIN(BR1 * 0.5)
WIDTH2 = (DIFF2 * UNCOV * 0.5) / DSIN (BR2 * 0.5)
C
C WRITE THE POSITION OF VESSEL, OBSERVED RATES,
C AND THE DIFFERENCES BETWEEN THEM.
C
WRITE (6,200) ILAT,IMIN,RSEC,IL0N,ILMIN,P.LSEC,RATS1,
1 DIFF1,RATE2, DIFF2
200 F0RMAT(1X,I2,I3,F7.3, Iit,l3,F7.3,?10.2,F6.2,F10.2,F6.2)
C
C XMEAN1 AND XMEAN2 ARE THE MEAN DIFFS BETWEEN THE
C COMPUTED RATE AND THE OBERSERVED. XSLAV1 AND XSLAV2
C ARE THE STORED DIFFS,
C
XM2AN1 = XMEAN1 + DIFF1
XMEAN2 = XMEAN2 + DIFF2
XMEAN3 = XMEAN3 + WIDTH 1
XMEANU = XMEANU + WIDTH2
XSLAV1 (I) = DIFF1
XSLAV2(I) = DIFF2
550 CONTINUE
C
C COMPUTE THE MEAN AND STANDARD DEVIATION
C
XMEAN1 = XMEAN1 / JREC
XMSAN2 = XMEAN2 / JREC
XMEAN3 = XMEAN3 / JREC
XMEAN4 = XMEANU / JREC
93
c
DO 600 I = 1, JREC, 1
VAR1 = VAR1 + ((XSLA7 1(I) - XMEAN1)**2)
600 VAR2 = VAR2 + ((XSLAV2(I) - XMEAN2)**2)
C
7AR1 = VAR1 / (JREC - 1.0)
VAR2 = VAR2 / {JREC - 1.0)
C
SD1 = DSQRT (VAR 1)
SD2 = DSQRT (VAR2)
C
WRITE(6,201)
WRITE (6,210) XMEAN1,SD1,XMEAN3
WRITE (6,201)
WRITE (6,211) XMEAN2,SD2,XMEANa
WRITE (6,202)
210 F0RMAT(1X,' SLAVE #1, MEAN = •,F10.3,
1 • STANDARD DEVIATION = •,
2 F10.3,' DISTANCE IN METERS = »,F10.3)
211 F0RMAT(1X,' SLAVE #2, MEAN = «,F10.3,
1 ' STANDARD DEVIATION = •,
2 F10.3,' DISTANCE IN METERS = •,?10.3)
WRITE(6,201)
201 FORMAT (ia )
202 FORMAT(IHI)
STOP
END
C==== ====================== ======= === = ==== ===== = ==== = ========
SUBROUTINE SHORPT (IRE C, RPOST , RPOSTL, RCOUT, RCONTL , UDIST ,
1 USD,ULD,FA,3A)
Q
C
C SUBROUTINE SHORPT WILL SELECT A POINT ALONG THE SHORE
C WHICH IS OUTLINED FROM NORTH BY SELECTED POINTS WITH
9a
C KNOWN LATITUDES AND LONGITUDES. THE SHORE POINT IS
C INTERPOLATED BETWEEN TWO KNOWN POINTS USING THE TOTAL
C DISTANCE BETWEEN THE POSITION AND THE CONTROL STATION
C AND THE AZIMUTH BETWEEN THE SHORE POINTS AND THE
C RECEIVERS POSITION.
C
IMPLICIT REAL*8 (A-H,0-Z)
DIMENSION FAZM(IOO)
COMMON/SHORE/PHI (100) , ELON (100) , UNCOV,RN, A , F
DATA PI/3, 1U159 26535398D0/
C
TWOPI = 2.0 * PI
C
CALL INVER1 (A,F ,RPOST ,RP0STL,R30NT, RCONTL, FA , B A, RDISTT)
UDIST = (RN * RDISTT) / UNCOV
C
RCOMP = 99999.99
C
DO 10 J=1,IREC, 1
P1 = PHI (J)
El = ELON (J)
CALL INVER1 (A, F, PI, El ,RCONT, 3C3NTL , AZF, A23, RDISTL)
CALL INVER1 (A,F ,RPOST ,RP0STL,P1 ,E1 , FAZ, BAZ , RDISTS)
USD = (RN * RDISTS) / UNCOV
ULD = (RN * RDISTL) / UNCOV
FAZM (J) = AZF
C
C COMPUTE THE DIFFERENCE BETWEEN THE TOTAL DISTANCE
C (UDIST) AND THE SUMMATION OF THE DISTANCE OVER THE WATER
C (USD) AND THE DISTANCE OVER THE LAND (ULD) . IT THE
C DISTANCE IS LESS THAN RCOMP, UPDAT RCOMP AND JSTA.
C JSTA IS THE CLOSEST POINT ALONG THE SHORELINE WHICH IS
C NEAR THE EM PROPAGATION PATH.
C
95
FDIFP = DABS(UDIST - (USD + ULD) )
IF(FDIFF .GT. RCOMP) GO TO 10
RCOMP = FDIFF
JSTA = J
10 CONTINUE
C
C NOW DETERMINE THE INTERPOLATED LATITUDE AND LONGITUDE
C SHORE POINT USING AZIMUTH PERCENTAGE.
C
TOPPER = JSTA - 1
ILOWER = JSTA ■«■ 1
AZMU = FAZM(IUPPER)
AZML = FA ZM (ILOWER)
IF (AZMU .GT. FA -AND. FA .GE. FAZM(JSTA)) ICH = lUPPER
IF (FAZM (JSTA) . GE. FA .AND. FA .GT. AZML) ICH = ILOWER
RADJ = 1.00 - DABS ((FAZM(ICH) - FA) / (FAZM (ICH) -
1 FAZM (JSTA) ) )
RNWPHI = PHI (JSTA) ^ ( (PHI (ICH) - PHI (JSTA) ) * RADJ)
IF (ELON (JSTA) . LE. ELON(ICH)) RADJ = 1.0 - RADJ
RNWELN = ELON (JSTA) + ((ELON (ICH) - ELON (JSTA)) * RADJ)
CALL TODMS (RNWPHI, IDG ,MIN, SEC)
RHOLD = TWOPI - RNWELN
CALL TODMS (RHOLD, IDGL,MINL,SSCL)
CALL INVER1 (A,F ,RNWPH I, RNWELN, RCONT, RCONTL , AZ, BZ , RDISTL )
CALL INVER1 (A,F,RPOST , RPOSTL,RN WPHI, RNWELN , AZ, 3Z, RDISTS )
USD = (RN * RDISTS) / UNCOV
ULD = (RN * RDISTL) / UNCOV
RETURN
END
96
SUBROUTINE SECF AC(UTD 1ST, USDIST , ULDIST, AZI , I STA , SF)
C
c
C THIS ROUTINE WILL COMPUTE THE LAND/SEA SECONDARY FACTOR
C UTDIST = TOTAL DISTANCE
C USDIST = DISTANCE OVER THE SEA WATER PATH
C ULDIST = DISTANCE OVER THE LAND PATH
C AZI = AZIMUTH FROM NORTH.
C
C MASTER = 1 (ISTA)
C W =2
C X =3
C Y =4
C
IMPLICIT REAL*8 (A-H,0-Z)
51 = (.795 / USDIST) + 0.439 + (.00245 * USDIST)
C
IF (UTDIST .GT. 540) GO TO 10
52 = (3.188 / UTDIST) - 0.594 + (.000329 * UTDIST)
GO TO 20
10 52 = (128.8 / UTDIST) + 0.187 + (.000652 - UTDIST)
C
20 IF(USDIST .GT. 540) GO TO 30
53 = (3.188 / USDIST) - 0.594 + (.000329 * USDIST)
GO TO 40
30 S3 = (128.8 / USDIST) + 0.187 + (.000652 * USDIST)
C
40 IF (ULDIST .GT. 540) GO TO 50
54 = (3.188 / ULDIST) - 0.594 + (.000329 * ULDIST)
GO TO 60
50 S4 = (128.8 / ULDIST) + 0.187 *■ (.000652 * ULDIST)
C
60 S5 = 1.428 + (.00158 * UTDIST)
S6 = 1.428 -•• (.00158 * ULDIST)
97
TAZI = 2-0 * AZI
GOTO (70,80,90,100) , ISTA
70 RHOLD = (1.010*DSIN (AZI) ) - ( . 1 96*DCOS (AZI) )
1 - (.893*DSIN(TAZI) ) - ( . 355*DCOS (TAZI) )
GO TO 200
80 RHOLD = (.323*DCOS (AZI) ) - (. 711 *DSIN (TAZI) )
GO TO 200
90 RHOLD = (.9a2*DCOS (TAZI) )
GO TO 200
100 RHOLD = (.588*DSIN (TAZI) )
C
200 S5 = S5 + RHOLD
S6 = S6 + RHOLD
SF = 0.5 * (S5 ^ S6 - SI + S2 + S3 - SU)
WRITZ (6,5 00) S1 ,S2,S3 ,SU,S4,S5,S6 ,S?
500 FORMAT(1X,7F15. 5)
RETURN
END
C== ============================================== ==========^
SUBROUTINE INVERT ( A, ? INV, GLAT1 , GLON 1 ,GLAT2 , GL0N2 ,FAZ,
1 3AZ,S)
C
See Appendix B for subroutine INVSR1.
98
short points around Monterey Bay, California. The point
are used to interpolate geodetic point for computation
of the distance over land and the distance ov^r sea.
NO. LATITUDE LONGITUDE
1 36 57 18.606 122 05 37.525
2 36 56 59.264 122 03 01.817
3 36 57 05.076 122 01 31.701
U 36 57 49.538 122 01 07.857
5 36 57 17.949 121 58 19.830
6 36 58 08.589 121 57 07.288
7 36 58 32.140 121 55 10. 083
8 36 58 01.498 121 53 57.390
9 36 56 46.115 121 52 22.313
10 36 55 38.140 121 51 24.399
11 36 53 13.806 121 49 46.743
12 36 49 38.384 121 47 48.395
13 36 47 39.241 121 47 10.818
14 36 46 27.554 121 47 39.637
15 36 44 56.717 121 47 52.416
16 36 41 14.439 121 48 32.642"
17 36 39 17.211 121 49 28.533
18 36 37 31.128 121 50 31.723
19 36 36 23.446 121 51 34.833
20 36 36 03.628 121 52 50.879
21 36 36 24.782 121 53 48.453
22 36 37 18.151 121 54 11.628
23 36 38 00.300 121 55 57.538
99
DATA SET EXAMPLE - PROGRAM LOPLC
SHIP'S POSITION OBSERVED LORAN RATES & ERRORS (E)
LATITUDE LONGITUDE 9940-Y EY 9940-W EM
(D-M-S) (D-M-S) (;isec) (;isec) (jjusec) (jasec)
36 43 45-800 121 55 27.160 42789.34 1.14 16294.04 1.32
36 44 3.400 121 55 32.340 42791.13 1.25 16293.46 1.29
36 44 21.180 121 55 37.390 42793.04 1.25 16292.73 1.40
36 44 37.490 121 55 46.950 42795.13 1.05 16292.03 1.33
36 44 53.260 121 55 57.710 42796.93 1.12 16291.43 1.11
100
APPENDIX D
PROGRAM LORTAB
C PROGRAM LORTAB
C
C PROGRAM COMPUTES LINE OF POSITION FOR LORAN-C USING
C SF SALT WATER CORRECTION FACTOR AND DMAHTC CALCULATED
C OR FIELD OBSERVED ASF CORRECTIONS.
C
C
IMPLICIT REAL*8 (A-H, 0-Z)
DIMENSION AXIS (13) ,RF (13)
DIMENSION XSLAV1 (1000) ,XSLAV2 (1000)
DATA RHOSEC,PI, UNCOV, RN/2 . 06264 8062471 DOS ,
1 3. 14 15 926535 89 800,29 9. 79 24 58D0 ,1 .00 03 3 8 DO/
DATA XMEAN1,XMEAN2,VAR1//AR2/0. OODO, O.OODO,
1 O.OODO,OO.ODO/
DATA XMEAN3,XMEAN4/0. OODO,O.OODO/
C
CCCCCC* ************ ******** ♦***********:*
DATA AXIS/6. 3782064006,6. 373338 006,6.377397155006,
1 6.37816 006,6.37816006,6.378249145006,6.373165006,
2 6,378166006,6.378165 006,6.37814 5006,6.3775634006,
3 6.378245006,6,3781350006/
DATA RF/6. 3565838 006,2.97 002,2.99152312 35002,
1 2.9 8 2500 2,2.98 24716 7 4 270 02,2.9 3 46 5002,2.9 8 2500 2,
2 2.983 002,2.983002,2.9825 002,6.3562569006,
3 2,983 002,2.9826002/
C*****ELLI?SOID OPTION NUMBER
C 1. CLARKE 1866 8. MERCURY
C 2. INTERNATIONAL (HAIFORD) 9. MARSHALL ISLAND
C 3. BESSEL 1841 10. NAVY 80
101
c
4. AND (AUSTRALIAN) 11. AIRY
c
5. 1967 REFERENCE 12. KRASSOWSKI 1940
c
6. CLARKE 1880 MOD 13. WGS 1972
c
7. SAO
c
cccccc* ************ ******** *************
1 c
CC1-2 = ELLIPSOID NUMBER (K)
c
CC3-5 = NUMBER OF POINTS ALONG COAST (IREC)
c
READ(5,100) K,IREC,ID1,ID2
100 F0RMAT(1X,I2,I4,2I3)
TW0PI=2.*PI
A=AXIS(K)
F=1./RF (K)
IF (F.LT.3.D-3) F=(A-1./F)/A
C
C READ MASTER AND SLAVE STATIONS POSITIONS
C THE FIRST RECORD IS THE NUMBER OF MASTER AND SLAVS
C STATIONS FOR THE PARTICULAR CHAIN.
C
WRITS (6,202)
202 FORMAT(IHI)
WRITE(6,201).
201 F0RMAT(1H )
READ (5, 105) ILATM,IMINM,RSSCM,ILONM,ILMINM,RLSECM
WRITE (6,105) ILATM,IMINM,RSSCM,ILONM,ILMINM,HLSECM
105 F0RMAT(1X,I3, 1X,I2,1X,F5. 2, 1X,I4, 1X,I2,1X,F5.2)
RMAST = ( (lABS(ILATM) * 60 + IMINM) * 60 + RSECM) /
1 RHOSEC
IF(ILATM .LT. 0) RMAST = -RMAST
RMASTL = ((TABS (ILONM) * 60 + ILMINM) *60 +RLSSCM) /
1 RHOSEC
IF (ILONM .GT. 0) RMASTL = TWOPI - RMASTL
102
C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
C
READ (5, 101) IIATS,IMINS,ESECS,ILONS,ILMINS,RLSECS,
1 DELAY1
WRITE (6,1 01) IIATS,IMINS,RSECS,ILONS,ILMINS,RLSECS,
1 DELAY1
101 FORMAT (IX, 13, IX ,12 , 1 X ,F5. 2, IX, I 4 , 1 X,I2, IX, F5 . 2 , F9 . 2)
RSLA71 = ( (lABS (ILATS) * 60 + IHINS) * 60 + 5SECS) /
1 RHOSEC
IF(ILATS .LT. 0) RSLAV1 = -PSLAV1
RSLAL1 = ((lABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL1 = TWOPI - R3LAL1
C
C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
C -
READ (5, 101) ILATS,IMINS,RSECS, ILONS, ILMINS, RLSECS,
1 DELAY2
WRITE (6,101) ILATS, iaiNS,RSECS, ILONS, ILMINS, RLSECS,
1 DELAY2
RSLAV2 = ((lABS (ILATS) * 60 + IMINS) * 60 + RSECS) /
1 RHOSEC
IF (ILATS .LT. 0) RSLAV2 = -RSLAV2
RSLAL2 = ((IA3S (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2
WRITE (6,201)
C
C READ THE RECORDS POSITIONS AND OBSERVED RATES, THEN
C COMPUTE THE TOTAL DISTANCE.
C
DO 5 50 I=1,IREC,1
103
READ (5, 10U) ILATrlMIN ,RSEC, ILOM , ILMIN, RLSEC ,
1 RATE1,RATE2
10U FORM AT ( IX, 13, 13, F6. 3, 15 ,1 3 ,F6 . 3 , 3X, 2F10 . 2)
C
C CONVERT LAT AND LONG TO RADIANS.
C
RPOST = ((lABS(ILAT) * 60 + IMIN) * 60 + RSEC) /
1 RHOSEC
IF(ILAT .LT. 0) RPOST = -RPOST
RPOSTL = ((IA3S(IL0N) * 60 + ILMIN) * 60 + RLSSC) /
1 RHOSEC
IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL
C
C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT
C
C MASTER
C
CALL INVER1 (A ,F , RPOST , RPOSTL, RM AST, RMASTL, FAZM ,BAZM,
1 DISTM)
UDISTM = (RN * DISTM) / UNCOV
CALL SECF AC (UDISTM, SFM)
C
C SLAVE1
C
CALL INVER1 (A,r , RPOST , RPO STL,RSLA V1 , RSLALI , FAZ 1 , BAZ 1 ,
1 DIST1)
UDIST1 = (RN * DIST1) / UNCOV
CALL SECFAC (UDIST1,SF 1)
C
C SLAVE2
C
CALL INVER1 (A, F , RPOST , RPO STL, RS LA V2, RSLAL2 , F AZ2 , 3A22,
1 DIST2)
UDIST2 = (RN * DIST2) / iJNCOV
104
CALL SSCFAC{0DIST2,SF2)
C
C DETERMINE THE ADDITIONAL SECONDARY CORRECTORS FROM THE
C LORAN-C CORRECTION TABLE FOR THE WEST COAST CHAIN 99U0
C
CALL TABLE(RPOST,RPOSTL,ID1,ASF1)
CALL TABLE(RPOST,RP0STL,ID2,ASF2)
C
C COMPUTE THE HATES AND COMPARE TO THE OBSERVED RATES
C
TDM1 = UDIST1 - UDI3TM + SF1 - SFM + ASF1 + DELAY1
DIFF1 = TDM1 - RATE1
C
TDM2 = UDIST2 - UDISTM + SF2 - SFM + ASF2 -»• DELAY2
DIFF2 = TDM2 - RATE2
C
C
C COMPUTE THE LANE WIDTH IN METERS BASED ON EQUATION
C 4.20 IN ELECTRONIC SURVEYING AND NAVIGATION -
C LAURILA, PAGE 94.
C
BR1 = DABS (FAZ1 - FAZ M)
BR2 = DABS(FA22 - FAZ M)
WIDTH1 = (DIFF1 * UNCOV * 0.5) / DSIN(BR1 * 0.5)
WIDTH2 = (DIFF2 * UNCOV * 0.5) / DSIN(BR2 * 0.5)
C
C WRITE THE POSITION OF VESSEL, COMPUTED RATE, OBSERVED
C RATE, AND THE DIFF BETWEEN THEM.
C
WRITE (6,200) ILAT,IMIN,RSEC,IL0N,ILMIN,RLSEC,RATE1,
DIFF1,RATE2,DIFF2
200 F0RMAT(1X,I3,I3,1X,F6 . 3,1 5,13, 1 X, F6. 3 ,F 12. 2 , F8 . 2,
1 F12.2,F8.2)
WRITE (6,201)
105
c
c
c
c
c
550
c
c
c
XMEAN1 AND XMEAN2 ARE THE MEAN DIFFS BETWEEN THE
COMPUTED RATE AND THE OBEHSERVED. XSLAVI AND XSLAV2
ARE THE STORED DIFFS.
XMEAN1 = XMEAN1 + DIFF1
XMEAN2 = XMEAN2 ••- DIFF2
XMEAN3 = XMEAN3 + MIDTH1
XMEAN4 = XMEANU + WIDTH2
XSLAV1 (I) = DIFF1
XSLAV2(I) = DIFF2
CONTINUE
COMPUTE THE MEAN AND STANDARD DEVIATION
XMEAN1 = XMEAN1 / IREC
XMEAN2 =-XMEAN2 / IREC
XMEAN3 = XMEAN3 / IREC
XMEAN4 = XMEANa / IREC
600
C
DO 6 00 I = 1, IREC, 1
VAR1 = VAR1 + ((XSLA7 1(I) - XMEAN1)**2)
VAR2 = VAR2 ■•• ((XSLA72(I) - XMEAN2)**2)
VAR1 = VAR1 / (IREC - 1.0)
VAR2 = VAR2 / (IREC - 1.0)
SD1 = DSQRT (VAR 1)
SD2 = DSQRT (VAR2)
WRITE (6,201)
WRITE(6,210) XMEAN1,SD1,XMEAN3
WRITE (6,201)
WRITS(6,211) XMEAN2,SD2,XMEAN4
106
WRIT2(6,202)
210 F0RMAT(1X,» SLAVE #1, MEAN = SFIO.S,
1 • STANDARD DEVIATION = «,
2 P10.3,' DISTANCE IN METERS = »,F10.3)
211 F0RMAT(1X,' SLAVE #2, MEAN = ',F10.3,
1 • STANDARD DEVIATION = »,
2 F10.3,' DISTANCE IN METERS = •,F10.3)
STOP
C========= = ========== ====== ==== === === = = === ====== = = = = = = ===
SUBROUTINE SECF AC(UTD 1ST, SF)
C
C
C THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR
C OTDIST = TOTAL DISTANCE
C SF = SECONDARY FACTOR
C
IMPLICIT REAL*8 (A-H,0-Z)
C
C
C COEFFICIENTS
C
AO = 129.0U323
A1 = -0.U0758
A2 = 0.00064576813
C
BO = 2.741282
B1 = -0.011402
B2 = .000327748 15
C
IF(UTDIST .GT. 537) GO TO 10
SF = ( BO / UTDIST) + B1 + ( 32 * UTDIST)
GO TO 20
10 SF = ( AO / UTDIST) + A1 + ( A2 * UTDIST)
C
107
20 CONTINUE
BETURN
END
C=="= = = == == = ==== = = ==:====== ===== = = = == = = = =:==:== = = = = = = = = = = = ==== =
SUBROUTINE TABLE(RLAT ,RLON,ID, ASF)
Q
C
C SUBROUTINE TABLE SELECTS THE PROPER ASF CORRECTOR FROM
C THE LORAN-C CORRECTION TABLE PUBLISHED BY THE DEFENSE
C MAPPING AGENCY.
C RLAT = POSITION LATITUDE IN SECONDS
C RLON = POSITION LONGITUDE IN SECONDS
C ID = LORAN-C CHAIN IDENTIFIER
C W = 1
C X = 2
C Y = 3
C ASF = ADDITIONAL SECONDARY FACTORS
C
C THE FOLLOWING TABLES OF ASF CORRECTORS ARE FOR
C MONTEREY BAY, CALIFORNIA - 9940 -W, -X, -Y,
C
IilPLICIT REAL*8 (A-H, 0-Z)
DIMENSION TA3LEW(3,5) ,TAB LEX (3 , 5) ,TABLEY (3 , 5)
C
DATA TABLSW/ 1. 6D0, 1 - 6D0 , 0. ODO ,
1 1. 6D0,1. 5D0,1.4D0,
2 1. 5D0,1, 4D0,1 .6D0,
3 1. 4D0,1. 3D0,1 .5D0,
4 1. 3D0,0. ODO, O.ODO/
C
DATA TABLEX/ -0 . 9D0, - 0. 9D0 , O.ODO,
1 -0.9D0,- 0.9D0,-0.3D0,
2 -1 ,0D0,- 1.0DO,-1.0DO,
3 -1 .ODO,- 1.2D0,-1 .1D0,
108
a -1.1D0, O.ODO, O.ODO/
c
DATA TABLEY/ 0. 2D0 , 0. 2D0,0 - ODO,
1 0. 3D0,0. 3D0,0.3D0,
2 0. 3D0,0. 2D0,0.4D0,
3 0. 4D0,0. 3D0,0.6D0,
4 0. 5D0,0. ODO, O.ODO/
DATA RHOSEC,PI/2,0626 4306 2a71DD5,3.1U159 26 5 358 98D0/
C
TWOPI = PI * 2. 0
C
C CONVEE RLAT AND RLON TO SECONDS
C
HLAT = ELAT * RHOSEC
HLON = TWOPI - RLON
HLON = HLON * RHOSEC
C
C STARTING LAT AND LONG FOR SEARCH
C LAT = 37/00/00.0 LONG = 122/05/00.0
C
c
C DETERMINE THE AS? CORRECTOR FOR THE LORAN-C COMBINATION
C
C LATITUDE
C
SLAT = 133200.0
SLON = 439500.0
RMID = 300.0
RDIFF = 150.0
DO 10 J=1,5, 1
SLAT = SLAT - RMID
ULAT = SLAT + RDIFF
VLAT = SLAT - RDIF?
109
10 IF (hLAT .LT. DLAT .AND. hLAT .3E. VLAT) GO TO 15
15 CONTINUE
C
C LONGITUDE
C
DO 30 1=1,3,1
SLON = SLON - RMID
ULON = SLON + RDIFF
VLON = SLON - RDIFF
30 IF (hLON .LT. ULON .iND. h LON .3T. VLON) GO TO 35
35 CONTINUE
C
C DETERMINE ASF CORRECTOR
C
ASF =0.0
IF (ID .EQ. 1) ASF = rABLEW(I,J)
IF (ID .SQ. 2) ASF = TABLEX(I,J)
IF (ID .EQ. 3) ASF = TABLEY(I,J)
RETURN
END
SUBROUTINE INVER1 ( A, F INV, GLAT1 , GL0N1 ,GLAT2 , GL0N2 ,FAZ,
1 BAZ,S)
See Appendix B for subroutine IN7ER1
110
DATA SET EXAMPLE - PROGRAM LORTAB
SHIP'S POSITION OBSERVED LORAN RATES & ERRORS (E)
LATITUDE LONGITODE 9940-Y EY 9940-W EW
(D-M-S) (D-M-S) (psec) (fisec) {fisec) (usee)
36 U3 45.800 121 55 27. 160
36 44 3.400 121 55 32.340
36 44 21. 130 121 55 37.390
36 44 37.490 121 55 46. 950
36 44 53. 260 121 55 57.710
42789.
.34
-0.
29
16294,
.04
0.
34
42791.
13
-0.
18
16293
.46
0.
30
42793,
04
-0,
18
16292,
.73
0.
41
42795.
13
-0.
38
16292,
.03
0.
34
42796.
93
-0.
31
16291,
.43
0.
13
111
APPENDIX E
PROGRAM ASFSEL
C PROGRAM ASFSEL
C
C PROGRAM DETERMINES OBSERVED ASF CORRECTORS BY SCANNING
C DATA AT 1 DEGREE LATITODE AND LONGITODE INTERVALS. THE
C ASF CORRECTORS ARE DETERMINED 3Y SUBTRACTING THE
C CALCULATED TD USING THE SEAWATER SECONDARY FACTOR FROM
C THE OBSERVED TD RATES.
C
IMPLICIT REAL*8 (A-H,0-Z)
DIMENSION AXIS ( 13) ,RF (13)
DIMENSION ASFCR 1 (16,2 6) ,ASFCR2(16 ,2 6) ,
1 IN01 (16,26) ,IN02(16, 26)
DATA RHOSEC,PI, UNCOV, RN/2, 06264 3062471D05,
1 3. 1 4 15 9265 35 89 8D0, 29 9.79 245800,1 .00 03 38D0/
DATA XMEAN1,XMEAN2,VAR1 //AR2/0. 0 0D0,0.00D0,
1 0.00D0,00.0D0/
DATA XMEAN3,XMEAN4/0. 00D0,0.00D0/
C
cccccc* ******************** *************
DATA AXIS/6. 3782064DO 6,6. 378383 D06 , 6. 3773971 55D06 ,
1 6.37816 D06, 6. 37816D06, 6. 378249145006,6. 378165D06,
2 6.373166006,6.378165 D06 , 6 . 378 1 45D06 ,6 . 3775634006 ,
3 6.378245006,6.3781350006/
DATA RF/6. 3565838 006,2.97 D 02 , 2 .99 152312 85D02 ,
1 2. 9 8 2500 2, 2. 98 247 16 7 4270 02,2.9 3 465002,2.9 8 2500 2,
2 2.983 002,2.983002,2.9825 002,6.3562569006,
3 2.983 002,2.9326002/
C*****ELLI?SOID OPTION NUMBER
C 1. CLARKE 1866 8. MERCURY
112
C 2. INTERNATIONAL (HAYFORD) 9. MARSHALL ISLAND
C 3. BESSEL 1841 10. NAVY 8D
C 4. AND (AUSTRALIAN) 11. AIRY
C 5. 1967 REFERENCE 12. KRASSCWSKI 194 0
C 6. CLARKE 1880 MOD 13. WGS 1972
C 7. SAO
C
cccccc* ****** ******** ****** *************
C CCl-2 = ELLIPSOID NUMBER (K)
C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC)
C
READ(5,100) K,IREC,ID1,ID2
100 F0RMAT(1X,I2,I4,2I3)
TW0PI=2.*PI
A = AXIS(K)
F=1./RF (K)
IF (F.LT.3.D-3) F=(A-1./F)/A
C
C READ MASTER AND SLAVE STATIONS POSITIONS
C THE FIRST RECORD IS THE NUMBER OF MASTER AND SLAVE
C STATIONS FOR THE PARTICULAR CHAIN.
C
WRITE (6,202)
202 F0RMAT(1H1)
WRITE (6,201)
201 FORMAT (1H )
READ (5, 105) ILATM,IMINM,RSECM,ILONM, ILMINM, RLSSCM
WRITE (6,105) ILATM,IMINM,RSSCM,ILONM,ILMINM,RLSECM
105 FCRMAT(1X,I3, 1X,I2,1X ,F5. 2, 1X,I4, 1X,I2, 1X,?5.2)
RMAST = ( (lABS(ILATM) * 6 0 > IMINM) * 60 + RSECM) /
1 RHOSEC
IF (ILATM .LT. 0) RMAST = -RMAST
RMASTL = ( (TABS (ILONM) * 60 + ILMINM) *60 +RLSECM) /
1 RHOSEC
113
IF(ILONM .GT- 0) RMASTL = TWOPI - RMASTL
C
C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
C
READ (5, 101) ILATS,IMINS,RSECS,ILONS,ILMINS,RLSECS,
1 DELAY1
WRITE (6,101) ILATS,IMINS,RSECS,ILONS,ILMINS,RLSECS,
1 DELAY1
101 FORM AT (IX, 13, IX, 12, IX, F5. 2, IX, 14, IX, 12, IX, F5. 2, F9. 2)
RSLAV1 = ((lABS (ILATS) * 60 + IMINS) * 60 + RSECS) /
1 RHOSEC
IF(ILATS .LT. 0) RSLAV1 = -RSLAV1
RSLAL1 = ((lABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL1 = TWOPI - RSLAL1
C
C READ THE SECOND SLAVS STATION AND CHANGE THE LATITUDE
C AND LONGITUDE INTO RADIANS.
C
READ (5, 101) ILATS, IMINS, RSECS, ILONS, ILMINS, RLSECS,
1 DELAY2
WRITE (6,101) ILATS, I ^ INS, RSECS, ILONS, ILMINS, RLSECS,
1 DELAI2
RSLAV2 = ( (TABS (ILATS) * 60 + IMINS) * 60 + RSECS) /
1 RHOSEC
IF (ILATS .LT. 0) RSLAV2 = -RSLAV2
RSLAL2 = ( (lABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) /
1 RHOSEC
IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2
WRITE(6,201)
C
C READ THE RECORDS POSITIONS AND OBSERVED RATES, THEN
C COMPUTE THE TOTAL DISTANCE.
11U
DO 55 0 I=1,IREC,1
READ (5, 104) ILAT,IMIN ,RSEC,ILON , ILMIN , RLSEC,
1 RATE1,RATE2
104 FOEMAT(1X,I3,I3,F6.3, 15,13 ,F6. 3 , 3X, 2F1 0. 2)
C
C CONVERT LAT AND LONG TO RADIANS.
C
RPOST = ((IA3S(ILAT) * 60 + IMIN) * 60 + RSEC) /
1 RHOSEC
IF(ILAT .LT. 0) RPOST = -RPOST
RPOSTL = ((IABS{ILON) * 6 0 -»• ILMIN) * 60 + RLSEC) /
1 RHOSEC
IF(ILON .GT. 0) RPOSTL = TMOPI - RPOSTL
C
C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT
C
C MASTER
C
CALL INVER1 (A ,F , RPOST , RPOSTL, RM AST ,RMASTL, FAZM , BAZM ,
1 DISTM)
UDISTM = (RN * DISTM) / UNCOV
CALL SECFAC (UDISTM,SFM)
C
C SLA7E1
C
CALL INVER1 (A,F ,RPOST ,RPO STL,RSLA71 , RSLAL1 , FAZ1 , 3AZ 1 ,
1 DIST1)
UDIST1 = (RN * DIST1) / UNCOV
CALL SECFAC (UDIST1,SF 1)
C
C SLAVE2
C
CALL INVER1 ( A, F, RPOST ,RPO STL,RS LAV2, RSLAL2, F AZ2, 3AZ2,
115
1 DIST2)
0DIST2 = (RN * DIST2) / UNCOV
CALL SECFAC (UDIST2,SF2)
C
C DETERMINE WHICH LAT AND LONG THE ASF CORRECTOR
C IS ASSIGNED TO.
C
CALL ASSIGN (RPOST, RPO STL, JN1, JN2)
C
C COMPUTE THE RATES AND COMPARE TO THE OBSERVED RATES
C
TDM1 = DDIST1 - UDISTM + SFI - SFM + DELAYI
DIFF1 = TDM1 - RATS1
C
TDM2 = UDIST2 - UDISTM + SF2 - SFK + DELAY2
DIFF2 = TDM2 - RATS2
C
C SUM THE DIFFERENCES TO THE MATRIX AND COUNT THE NUMBER
C OF ASF CORRECTORS FOR EACH BLOCK TO LATTER DETERMINE
C THE MEAN.
C
ASFCR1 (JN1, JN2) = ASFCR 1 ( JN 1, JN 2) + DIFFI
ASFCR2(JN1,JN2) = ASFCR2 ( JN1, JN2) + DIFF2
IN01 (JN1,JN2) = IN01 ( JN1, JN2) + 1
550 IN02 (JN1,JN2) = IN02 ( JN1, JN2) + 1
C
C DETERMINE THE MEAN ASF CORRECTOR FOR EACH LAT AND LONG
C
DO 650 1=1, 16,1
WRITE(6,201)
DO 600 J=1,26,1
IF(IN01(I,J) .EQ. 0) GO TO 580
ASFCR1(I,J) = ASFCR1(I,J) / IN01(I,J)
580 IF(IN02(I,J) .EQ. 0) GO TO 590
116
aSFCR2(I,J) = ASFCR2(I,J) / IN02(I,J)
590 IF(ASFCR1 (I,J) .EQ. 0.0) ASFCR1(I,J) = 9.99
600 IF(ASFCR2 (I,J) .EQ. 0.0) ASFCR2(I,J) = 9.99
650 CONTINUE
C
C WRITE THE CORRECTORS IN MATRIX FORMAT
C
C
IF(ID1 .EQ. 1) WRITE (6,1000)
IF(ID1 .EQ. 2) WRITE(6,1001)
IF(ID1 .EQ. 3) WRITE(6, 1002)
WRITE(6,201)
C
WRITE (6,300) ( (ASFCR1 (I,J) ,1 = 1 , 1 6) , J=1 , 26)
DO 700 L=1,5, 1
WRITE(6,201)
201 FORMAT (1H )
700 CONTINUE
C
IF(ID2 .EQ. 1) WRITE (6,1000)
IF(ID2 .EQ. 2) WRITE (6,100 1)
IF(ID2 .EQ. 3) WRITE(6,1002)
WRITE (6,201)
WRITE (6,300) ( (ASFCR2 (I,J) ,1 = 1,16) ,J = 1,26)
300 F0RMAT(16F6.2)
1000 FORMATC TABLE FOR 99U0-W »)
1001 FORMATC TABLE FOR 99 UO-X ')
1002 FORMATC TABLE FOR 9940-1 •)
STOP
END
C"== ===== = ================ == ======= = = === ====== =
SUBROUTINE SECF AC( UTD 1ST, SF)
C
C
117
C THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR
C UTDIST = TOTAL DISTANCE
C SF = SECONDARY FACTOR
C
IMPLICIT REAL*8 (A-H, 0-Z)
C
C
C COEFFICIENTS
C
AO = 129.04323
A1 = -0,40758
A2 = 0.00064576813
C
BO = 2,741282
B1 = -0.011402
32 = .00032774815
C
IF (UTDIST .GT. 537) GO TO 10
SF = ( 30 / UTDIST) + 31 + ( 32 * UTDIST)
GO TO 20
10 SF = ( AO / UTDIST) + A1 + ( A2 * UTDIST)
C
20 CONTINUE
RETURN
END
C== =============== = ================== = = === === = === = === = = =
SUBROUTINE ASSIGN (RLA T, RLON, I, J)
C
c
C SUBROUTINE ASSIGN SELECTS THE COLUMN AND ROW FOR
C THE LATITUDE AND LONGITUDE OF THE RECORD.
C RLAT = POSITION LATITUDE IN SECONDS
C RLON = POSITION LONGITUDE IN SECONDS
C I = COLUMN
118
C J = ROW
C
IMPLICIT REAL*8 (A-H, 0-Z)
DATA RHOSEC,PI/2.0 62 6a806 2U71D0 5, 3 . 1 U 1 59265 35898D0/
C
TWOPI = PI * 2. 0
C CONVER RLAT AND RLON TO SECONDS
C
HLAT = RLAT * PHOSEC
HLON = TWOPI - RLON
HLON = HLON * RHOSEC
C
C STARTING LAT AND LONG FOR SEARCH
C LAT = 37/05/00.0 LONG = 122/05/00.0
C
SLAT = 133200,0
SLGN = 439500.0
c
c
LATITUDE
c
J = 0
DO 10 IC = 1,26,1
SLAT = SLAT - 6 0.0
RULAT = SLAT + 30.0
RLLAT = SLAT - 30.0
10
IF (HLAT .LT. RULAT -AND
15
J = IC
C
C
LONGITUDE
HLAT .GE. RLLAT) GO TO 15
1 = 0
DO 20 IC = 1,16,1
SLON = SLON - 60.0
RLLON = SLON + 30.0
119
RRLON = SLON - 30.0
20 IF(HLON .LT. RLLON .AND. HLON .GE. HRLON) GO TO 25
25 I = IC
RETURN
END
C== ======================== ================================:
SUBROUTINE INV2R1 ( A, F INV, GLAT1 , GLON 1 ,GLAT2,GLON2 ,FAZ,
1 BAZ,S)
C
See Appendix B for subrcu-ine INVER1.
120
DATA SET EXAMPLE - PROGRAM ASFSEL
TABLE FOR 9940-Y
50»
1220 / 00*
US*
0.56 -0.51
-0.39
-0.41
-0.53 -0.44
-0.46 -0.40
-0.54 -0.39 -0.55 -0.45
-0. 63 -0.57 -0.49
-0.54 -0.50
-0.62
40'
35'
360
121
TABLE FOR 9940-Y
1210 / 55*
1210 / 50'
50"
-0.45 -G.U7 -0.41 -0.37
45* -0.51 -0.50
-0.38 -0.49 -0.50 -0.46
-0.54 -0.51 -0.49
-0.54 -0.57 -0.49
-0.55 -0.55 -0.50
40* -0.67 -0.61 -0.50
-0.64 -0.65
-0.60
-0.43 -0. 44 -0.47 -0.44
-0.52 -0.46
35'
122
BIBLIOGRAPHY
Admiralty Manual of H ydrographic Surve^in^, First Edi-icn,
7 . T7 Hy^rograpEer oT ^he Navy, ronton, England, 1965.
!!as^ers TEesis JTz Frep araTionj". "
School, Monterey, California, 1982,
Backus, Alexandra, "HK121, LOEAN-C R^al-Time Hydroplot,"
Language PAL 8, NOAA, Rockville, Maryland, December 10,
1980.
Bigelow, Henry W. , "Elecrronic Surveying: Accuracy of
Electronic Positioning System," Journal of the
Siirve^in^ and Mapping Division, October 1'9E3,
Dansford, Ed, Electronics Division, DMAHTC, Washington,
D.C., private communications, 9 September 1982.
Defense Mapping Agency. M-^-can Practical N§.'^-3^l0£/
DMA Publication Number "'57 DE'K JTy^rograpHIc"* Cen-er,
Washington, D.C, v. 1, 1977.
Defense Mappina Agency Hydr ographic/Topographic Center,
Washington, 6.C., DMA Stock No. LCPU322 1 2100-C, LORAN-C
Correction Table, West Coast, USA 99a0, 1981.
du Plessis, Roaer M. , Poor Man's Explanation of Kalman
Eaton, R.M,, Mortimer, A.R., and Gray, D.H., "Accurate
Chart Latticing for LORAN-C," International
H ydrographic Review, Monaco, LV1 (1), January T^T^".
Funakoshi, Larry, Private Communication Forwarding Codv
of FORTRAN Program TD3HID, DMAHTC, Coda:' N^rlE,
Washington, D.C, 26 February 1982.
General Electric Co., TEMPO division, Santa Barbara,
California, Report Number CG-D-U-80, LORAN-C Signal
Analysis Final Report, by Nelson, Larry W. and iJam^iTTT
Hur^, TJecemBer TT75.
Heinzen, Morris Rudy, H^dro^raghic Surve_ys: Geodetic
Ccntrcl Criteria, Hasher or science Thesis, Cornell
[rniversiry,"TTEaca, New York, 1977.
Ingham, Alan, Sea Surveying , Volume 1, John Wiley
and Sons, 1 9TS7
123
Laurila/ Sinio H. , Electronic Surveying and Navigation,
Wiley-Interscience PuETIcation, T975'. ""
Leath, Dudley W., FORTRAN Program LSQR, by Paul R. Wolf,
Naval Postgraduate School, Monterey, California, 1981.
Mortimer, A. R. , Eaton, R. M. , and Gray, D. K.,
"CaliiDration of the West Canadian LORAN-C Chain."
Canadian Aeronautics and Space Journal, v. 24, No. 3,
Hay73une T^TBT
National Oceanic and Atmospheric Administration, NOAA,
Marine Chart Division, OA/C321, UNCLASSIFIED Letter to
Chief Radionavigat ion Division, U.S. Coas- Guard,
G-NRN-3. Washington, D. C, , Subject: ASF Corrections,
3 May 19 8 2.
National Oceanic and Atmospheric Administration, National
Ocean Survey, Descriptive 3§2^£i (Hydrographic) ,
Registry Number: H-?H22 , AlasRl, Guir of IlasTca,
S"ou^east of~Kcfiak Island, Chief of Party: Capt J. G.
Grunwell, Capt B. I. Williams, 1979.
National Oceanic and Atmosoheric Administration, National
Ocean Survey, Descriotive Report (Hydrographic) ,
Registry Numoer: TF^HII . AlasSa, GuTT o? liaslca,
" ■ -. • - £ p^^^y. c^pt J. G.
National Oceanic and Atmospheric Administration, National
Ocean Survey, Descriptive Report, Registry Number
SiJS902-WH-82, Texas'; nuII"of "Mexico, rTeaIa~and ^"Sine
FanTcs, Cn^sr of Party: Cdr. Roy K. Matsushige, 1982.
Nav
L
igation Department, Defense Mapping Agency, UNCLASSIFIED
etter to: LCDR Jerry Mills, Naval Postgraduate School,
Monterey, California, Subject: Definition of a Lane on
a Baseline of the LORAN-C Radionavigat ion System in tEe
gygergglic 32l£/ August 9, HE2, ~
NOAA Ship Surveyor SI 32, Cruise Report 3U-aO-7-82. Ridges
and Fractures Tectonics Experiment 2, by Alexander
!!aTahcrf, TTEz-et ^lentist leg T~aQ3 Dr. Rcber Embley,
Chief Scientist Leg II., 1982.
Pfaifer, L., LCDR., Subroutine INVER1 , National Geodetic
Survey, Rockville, Maryland, 18 February 1975.
Pfsifer, L., private communications. National Geodetic
Survey, RocKville, Maryland, 10 September 1982.
Poppe, Martin Carl, Jr., "The LORAN-C Receiver A Functional
Description," Navigation, v. 29, No.1, Spring 1982.
124
Frsssey, B- G.- Ashwee, G.E., and Fowler, C.5., "Char.ga
of Phase with Distance of a Low Frequency Sroundwavs
Propagated Across a Coastline," Proceedinqs I.E.S.,
Paper No. 2082R, July 195 6.
Racal-DECCA Survey, Inc., Brochure. Titled 5U0 Digital
Distance Measuring Unit, Racal-DECCA Survey, Inc.,
TU^TJT~¥est T5IJIcelirZT Houston, Texas, 1931.
Riordan, R., National Ocean Survey, Code C64X2, Rockville,
Maryland, private communication Program LORAN, 24 May
197 9.
1979.
Samaddar, S. N., "The Theory of LORAN-C Ground Wave
Propagation - A Review," Navigation, 7. 26, No. 3,
Fall, 1979.
Samaddar S. N., "Weather Effect on LORAN-C Propagation,"
Navigation, v. 27, No. 1, Spring 19 80.
Schne
t
?
1979;
Smith-Rose, R.L., D.Sc, Ph.D., "Electric Measurements
on Soil with Alternating Currents," Proceedings I.E.E.
V. 75, August 1934. ~
Speight, John J., "DMAHTC Support to National Ocean Survey
LORAN-C Charting," Navisation, Spring 1982.
The Analytic Science Corporation, Reading, Massachusetts,
Ap2lication of Semi-Empirical TD Grid Calibration to
the~¥es^ ^oas^ Lg!?ITP"C Cliai.n, by Ra3Ea R. ^uota, July
4, 1976.
U.S. Coast Guard, CF-462, LORAN-C User Handbook, Department
of Transportation Coasx Suar^, August T^T^T.
U.S. Naval Oceanographic Office^ Marine Science Department,
Washington, Q.C. , Unpublished Manuscript, The
Preparation of Secondary Phase Correction Charts, rcr
ID'RO-r. , Transmi^t e^~75' FeSf uary~T7B"77
125
Wallace, J. L., HYDRO PL0T/HYD30L0G System Manual, National
Ocean Survey Technical Manual No. 2, September 1971.
Wallace, J. L., "HYDROPLOT/HYDROLOG Systems Manual," U.S.
Dept. of Commerce, National Oceanic and Atmospheric
Administration, NOS Technical Manual No. 2, v. 1,
Section 3, Sept 1972 (Revised 21 May 1974).
Wallace, J.L. , National Oceanic and Atmospheric
Administration, private communications, July 15, 1982.
Wolf, Paul R. , Solution Manual to Accompany Elements of
P hotoqra metrvT Hc(?r aw-!IiIT~Soolc Company, 1977.
Wonnacott, Thomas H. and Wonnacott, Ronald J., Introductiory
Statistics, 3rd ed. , John Wiley and Sons., 1915.
125
INITIAL DISTRIBUTION LIST
4. Chairman (Code 63 RD)
Department of Mereoi
Naval Postqraduat € Schoo
Department of Meteorology
Naval Postgraduate School
Mcn-erev, CA 939a 0
5. Director
Naval Oceanography Division (OP952)
Naval Observatory
34th and Massachusetts Avenue NW
Washington, D.C. 20390
6. Commander
Naval Oceanography Command
NSTL Station
Bay St. Louis, MS 39522
7. Commanding Officer
Naval Ocaanoaraphic Office
NSTL Station"
Bay St. Louis, MS 39522
8. Commanding Officer
Fleet Numerical Oceanography Cen^iar
Monterey, CA 93940
9. Commanding Officer
Naval Ocean Research and Development
Ac-tiivity
NSTL S-ation
Bay St. Louis, MS 39522
10. Commanding Officer
Naval Environmental Prediction Research
Facility
Monterey, CA 939 40
11. Chairman, Oceanography Department
0.3. Naval Academy
Annapokis, MD 21402
12. Chief of Naval Research
800 N. Quincv Street
Arling-on, VA 22217
127
No. Copies
1. Defense Technical Information Center 2
Cameron Station
Alexandria, VA 22314
2. Library, Code 014 2 2
Naval postgraduate School
Monterey, CA 93940
3. Chairman (Code 68Mr) 1
Depar-^ment of Oceanography
Naval Postgraduate School
Monterey, CA 9394 0
13. Office of Naval Research (Coda U20)
Naval Ocean Research and Development
Activity
NSTL Station
Bay St. Loais, MS 39522
14. Director (Code PPH)
Defense Mapping Agency
Bldg. 56, U.S. Naval Observatory
Washington, D.C. 20305
15. Director (Code HO)
Defense Mapping Avency Hydcrgraphic
Topographic Canter
6500 Brookes Lane
Washing-con, D.C. 20315
16. Director (Co4e PSD-MC)
Defense Mapping School
Ft. Belvoir, VA 22060
17. Director
National Ocean Survey (OA/C)
National Oceanic and Atmosoheric
Administration
Roclcville, MD 20 852
18.
Chief, Program Planning and Liaison (NC2)
Na-ional Oceanic and Atmospheric
Administration
Rockville, MD 20852
19. Associate Director, Marine Surveys and
Maps (OA/C 3).
National Oceanic and Atmospheric
Administration
Rockville, MD 20 852
20. Chief, Hydrcgraphic Surveys Division (OA/C35)
National Ocenaic and Atmospheric
Administration
Rockville, MD 20 852
21. Director
Pacific Marine Center - NOAA
1801 Fairview Avenue East
Seattle, WA 9810 2
22. Director
Atlantic Marine Center - NOAA
439 W. York Street
Norfolk, 7A 23510
23. Commanding Officer
NOAA Ship RAINIER
Pacific Marine Center, NOAA
1801 Fairview Avenue 2 ast
Seattle, WA 98102
128
24. Ccmmanding Officer
NOAA Ship FAIRWEATHER
Pacific Marine Center, NOAA
1801 Fairview Avenue Sast
Seattle, WA 98102
25. Ccmmanding Officer
NOAA Ship DAVIDSON
Pacific Marine Cen-cer, NOAA
1801 Fairview Avenue East
Seattle, WA 98102
26. Commandina Officer
NOAA Ship^MT. MITCHELL
Atlantic Marine Center, NOAA
439 West York Street
Norfolk, Virginia 23510
27. Commanding Officer
NOAA Ship WHITING
Atlantic Marine Center , NOAA
439 Wesx York Street
Norfolk, Virginia 23510
28. Commanding Officer
NOAA Ship PSIRCE
Atlantic Marine Center, NOAA
439 West York Srreet
Norfolk, Virginia 23510
29. Chief, Hydrographic Surveys Branch
Arilantic Marine Center , NOAA
439 WesT York S-ree::
Norfolk, Virginia 23510
129
1 0 O r> r- p
factors 1 f^^-^Mary
aesls 19^^P58
W4844 Wheaton
c.l Application of
additional secondary
factors to LORAN-C
positions for hydro-
graphic operations.
Application of additional secondary fact
3 2768 001 95034 8
DUDLEY KNOX LIBRARY