%\documentclass[12pt]{report}
%\documentclass[12pt]{article}
%\documentclass[11pt]{report}
\documentclass[11pt]{article}
\setlength{\textheight}{8.7 in}
\setlength{\textwidth}{6.5 in}
\setlength{\topmargin}{1 cm}
\setlength{\oddsidemargin}{0.1 in}
\setlength{\evensidemargin}{-0.3 in}
\setlength{\headheight}{0 cm}
\setlength{\headsep}{0 cm}
\setlength{\topskip}{0 cm}
\pagestyle{empty}
\begin{document}
\title{
``HARPO: A versatile three-dimensional Hamiltonian ray-tracing program for
acoustic waves in an Ocean with Irregular Bottom'', by R. Michael Jones, J. P.
Riley, and T. M. Georges\\ NOAA special report, October 1986}
\date{Errata, \today}
\maketitle
\hspace*{-.25 in}page 47, line 7: Add ``with a copy of the input data file''
after ``raysets''.\\
\\
page 49, Fig. 2.24:  The implied decimal point for latitude of transmitter
should be between card columns 16 and 17.  The implied decimal point for
longitude of transmitter should be between card columns 22 and 23.  The implied
decimal point for imaginary part of wave polarization at transmitter should be
between card columns 75 and 76.  (20 April 2007)\\
\\
page 50, Fig. 2.25:  The imaginary part of wave polarization should be in card
columns 74-77, and the implied decimal point should be between card columns 75
and 76.  (20 April 2007)\\
\\
page 158, Table 7.25: Change ``XMIN0,YB'' to ``YMIN0,YB''\\
\\
page 158, Table 7.25: Change ``XMAX0,YT'' to ``YMAX0,YT''\\
\\
page 214:  The sentence ``Superimpose these raypath plots on the graph of the
previous runset:'' should read ``Superimpose these raypath plots on the graph
of the next runset:'' \\
\\
page 223, 6th line from bottom should read:\\
\hspace*{.5 in}the fractional increase of C with depth, $\epsilon_1 =$
\_\_\_\_\_\_\_\_(W158)\\
\\
page 223, last line should read:\\
\hspace*{.5 in}the fractional increase of C with depth, $\epsilon_2 =$
\_\_\_\_\_\_\_\_(W163)\\
\\
page 223: The formula for $\eta$ should contain a minus sign.\\
\\
page 224, 6th line from bottom should read:\\
\hspace*{.5 in}the fractional increase of C with depth, $\epsilon_1 =$
\_\_\_\_\_\_\_\_(W158)\\
\\
page 224, last line should read:\\
\hspace*{.5 in}the fractional increase of C with depth, $\epsilon_2 =$
\_\_\_\_\_\_\_\_(W163)\\
\\
page 224: The formulas for $\eta_1$ and $\eta_2$ should contain minus signs.\\
\\
page 225: Change ``W200'' to ``W150'', ``W201'' to ``W151'', and ``W202'' to
``W152''.\\
\\
page 235: Change the model check number for subroutine RBOTM to ``2.0''.\\
\\
page 236: Change the model check number for subroutine RVERT to ``3.0''.\\
\\
page 240: 9th line from the bottom, change ``HARPA'' to ``HARPO''.\\
\\
page 248: After line RAQC1670 of PROGRAM RAYTRC, insert\\
\\
\hspace*{.6 in}LHDRPG = LINES/LINSPP*LINSPP \hfill RAQC1675\\
\\
page 261:  line REXW0290 in FUNCTION READW should read:\\
\\
\hspace*{.6 in}PARAMETER (MXCMTS=83,BIGVAL=1.E9) \hfill REXW0290\\
\\
\\
page 261: line REXW0360 in FUNCTION READW should read:\\
\\
\hspace*{.6 in}LOGICAL NWOK,AB,UCON,UUCON,SETUCON \hfill REXW0360\\
\\
page 264: line REXW1980 in FUNCTION READW should read:\\
\\
\hspace*{.6 in}UUCON=SETUCON(X.Y,Z) \hfill REXW1980\\
\\
page 271: line TRWE0730 in SUBROUTINE TRACE should read:\\
\\
\hspace*{.6 in}IF(THERE.AND.FDOT.EQ.0.) RSIGN = SIGN(1.0,D2Z) \hfill TRWE0730\\
\\
page 273: In line TRWE1870 of SUBROUTINE TRACE, change SSURF TO ASURF.\\
\\
page 295: Line COPK0610 in SUBROUTINE CONBLK should read:\\
\hspace*{.5 in}DATA CUEF/1.d0,8.31436d-3,1.4d0/ \hfill COPK0610\\
\\
page 311: line OCBD1370 in SUBROUTINE OCNHD should read:\\
\\
1000  FORMAT(A1/A80,21X,2A10,' PAGE',I4) \hfill OCBD1370\\
\\
page 311: in line OCBD1540 in SUBROUTINE OCNHD, change ``F12.6'' to ``F18.6''.\\
\\
page 330: Change lines LAOT0750 and LAOT0760 to\\
C     LOOP FOR 9 MODELS AND PERTURBATIONS\\
\hspace*{.6 in}DO 1700 K=1,18\\
\\
\hspace*{.6 in}    Insert the following 4 lines after line LAOT0880\\
\hspace*{.6 in}IF(I.EQ.7) WRITE(LABEL,1600) (MODT(J),J=J1,J2)\\
\hspace*{.6 in}IF(I.EQ.8 .AND. J1.EQ.1)\\
\hspace*{.5 in}1 WRITE(LABEL,1600) (MODM(J),J=J1,J2)\\
\hspace*{.6 in}IF(I.EQ.9) WRITE(LABEL,1600) (MODP(J),J=J1,J2)\\
\\
page 335: Replace lines DRYS014 through DRYS019 in SUBROUTINE DRAWTKS with\\
\\
\hspace*{.6 in}IF(DTICV.GT.0.) THEN\\
\hspace*{.6 in}   YBP=YMID-AINT((YMID-YB)/DTICV)*DTICY\\
\hspace*{.6 in}   NTICX=(XR-XL)/DTICH+1.5\\
\hspace*{.6 in}ELSE\\
\hspace*{.6 in}   YBP=YB\\
\hspace*{.6 in}   NTICX=2\\
\hspace*{.6 in}ENDIF\\
C\\
\hspace*{.6 in}nticy=2\\
\hspace*{.6 in}IF(DTICV.GT.0.d0) NTICY=(YT-YBP+DTICV)/DTICV+0.5 ! added .5\\
\\
page 337: Change line PLGB0360 in SUBROUTINE PLTLB to\\
\hspace*{.6 in}WRITE(ANNOT,50) \\
\hspace*{.5 in}\&     IDINT((V-EARTHR)*F*dsign(1.d0,hb) + \\
\hspace*{.5 in}\& dsign(.5, (V-EARTHR)*F*dsign(1.d0,hb) ))\\
\\
page 337: Change line PLGB0400 in SUBROUTINE PLTLB to\\
60    WRITE(ANNOT,80) SNGL((V-EARTHR)*F*dsign(1.d0,hb))\\
\\
page 338: Change lines ARPC0100 and ARPC0110 in SUBROUTINE ARCTIC to\\
\\
\hspace*{.6 in}IF(TIC.NE.0.) NTIC=1.5+(THMAX-THMIN)/TIC \hfill ARPC0100\\
\hspace*{.6 in}NLINE=MAX0(3,100/NTIC) \hfill ARPC0110\\
\\
page 364: Line WGZ20140 in SUBROUTINE WGAUSS2 should be\\
\\
\hspace*{.6 in}DATA RECOGU/8.0/ \hfill WGZ20140\\
\\
page 381: Change line CSS10140 in SUBROUTINE CSMUNK1 to\\
\\
C\hspace*{.5 in}ETA = -2(Z-ZA)/H \hfill CSS10140\\
\\
page 384: Change line CSS20140 in SUBROUTINE CSMUNK2 to\\
\\
C\hspace*{.5 in}ETA = -2(Z-ZA)/H \hfill CSS20140\\
\\
page 389: (SUBROUTINE CTABLE)\\
\\
\hspace*{.5 in}Line CTUE0180, change CMX to CPX.\\
\hspace*{.5 in}Line CTUE0190, change CNTBL to CQTBL.\\
\hspace*{.5 in}Line CTUE0200, change CITBL to CLTBL.\\
\hspace*{.5 in}Line CTUE0210, change CFRMTBL to CIRMTBL.\\
\\
\hspace*{.5 in}After line CTUE0210:\\
\\
\hspace*{.5 in}Insert a copy of line CSS10370 through CSS10460 (from SUBROUTINE
\\
\hspace*{.5 in}CSMUNK1) on page 383.\\
\\
page 393: In lines CB2 0020 and CB2 0050 in SUBROUTINE NPSPEED, change B2 to
B4.\\
\\
page 415: Change lines GTLH0490 and GLTH0500 in SUBROUTINE GTANH to comments.\\
\\
page 428: Change line RVRT0190 in SUBROUTINE RVERT to \\
\\
\hspace*{.6 in}DATA RECORR/3.0/ \hfill RVRT0190\\
\\
page 429: SUBROUTINE GAMANG should be inserted before SUBROUTINE SMPANN in the middle of page 429.  This is a subroutine to calculate the side of a spherical triangle by law of cosines.\\
\\
%\small
%\scriptsize
\footnotesize
\hspace*{.5 in}SUBROUTINE GAMANG(THA,PHA) \hfill GAMANG 2\\
%      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
C\hspace*{.4 in}COMMON DECK ``GAMANG'' INSERTED HERE \hfill CGAMANG2\\
\hspace*{.5 in}COMMON/SPHGAM/SINLM0,COSLM0,GPH0,COSPHD,SINTH,COSTH \hfill  CGAMANG4\\
\hspace*{.5 in}COMMON/SPHGAM/GAMFUN,PGMTH,PGMPH,PGMTHTH,PGMPHPH,PGMTHPH \hfill  CGAMANG5\\
C \hfill GAMANG 4\\
\hspace*{.5 in}COSPHD=dCOS(PHA-GPH0) \hfill GAMANG 5\\
\hspace*{.5 in}SINTH=dSIN(THA) \hfill GAMANG 6\\
\hspace*{.5 in}COSTH=dCOS(THA) \hfill GAMANG 7\\
C \hfill GAMANG 8\\
\hspace*{.5 in}GAMFUN=SINLM0*COSTH+COSLM0*SINTH*COSPHD \hfill GAMANG 9\\
C \hfill GAMANG10\\
\hspace*{.5 in}PGMTH=-SINLM0*SINTH+COSLM0*COSTH*COSPHD \hfill GAMANG11\\
C \hfill GAMANG12\\
\hspace*{.5 in}PGMPH=-COSLM0*SINTH*dSIN(PHA-GPH0) \hfill GAMANG13\\
C \hfill GAMANG14\\
\hspace*{.5 in}PGMTHTH=-(SINLM0*COSTH+COSLM0*SINTH*COSPHD) \hfill GAMANG15\\
C \hfill GAMANG16\\
\hspace*{.5 in}PGMPHPH=-COSLM0*SINTH*COSPHD \hfill GAMANG17\\
C \hfill GAMANG18\\
\hspace*{.5 in}PGMTHPH=-COSLM0*COSTH*dSIN(PHA-GPH0) \hfill GAMANG19\\
C \hfill GAMANG20\\
\hspace*{.5 in}RETURN \hfill GAMANG21\\
\hspace*{.5 in}END \hfill GAMANG22\\
\normalsize
\\
We thank Arthur Newhall of WHOI for pointing out many of these errors.\\
\end{document}
