Improve standards compliance: integer(kind=4) and real(kind=8).pedantic

Note that integer*4 and real*8 are not actually part of the
language standard, and never have been.

1 files changed, 80 insertions, 80 deletions

diff --git a/htdp.f b/htdp.f
index b27b64e..e9a5985 100644
--- a/htdp.f
+++ b/htdp.f
@@ -11,7 +11,7 @@
 *
 **************************************************************
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
       character    HTDP_version*8
       character    Version_date*20
@@ -138,7 +138,7 @@ C    1      FORM='UNFORMATTED')
 *** Obtain parameters defining crustal motion model
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
       COMMON /TIMREF/ ITREF
 
@@ -249,7 +249,7 @@ C*** Set default reference epoch to Jan. 1, 2010
 ***   Region 65: Yangtze (YA)
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NMREGN = 65)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
       COMMON /BNDRY/ X(12000), Y(12000), NPOINT(70)
@@ -267,7 +267,7 @@ C*** Set default reference epoch to Jan. 1, 2010
 *** Determine the region containing a given point.
 
       IMPLICIT DOUBLE PRECISION(A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NMREGN = 65)
       COMMON /BNDRY/ X(12000), Y(12000), NPOINT(70)
       COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
@@ -323,7 +323,7 @@ C     NO MATTER HOW COMPLEX.
 C     ACCURACY -
 C     MAINTAINS FULL ACCURACY OF INPUT COORDINATES.
 C
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
       COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
       DIMENSION X(N),Y(N)
@@ -410,7 +410,7 @@ C
 C Convert horizontal velocities from mm/yr to rad/yr
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
       CALL RADII (YLAT,RADMER,RADPAR)
 c     write (*,*) "From RADR8T",RADMER,RADPAR
@@ -427,7 +427,7 @@ C
 C Compute the ITRF2008 velocity at a point in mm/yr
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NUMGRD = 10)
       parameter (NMREGN = 65)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
@@ -501,7 +501,7 @@ C*** then set its vertical velocity to 0.0
 ***   Snay 2003. "Introducing two spatial reference frames for regions of the Pacific Ocean." Surv. Land. Info. Sci., 63(1):5-12.
 
       IMPLICIT DOUBLE PRECISION (A-H, O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       PARAMETER (NUMPLATE = 52)
       
       DIMENSION WX(NUMPLATE), WY(NUMPLATE), WZ(NUMPLATE)
@@ -666,7 +666,7 @@ C*** then set its vertical velocity to 0.0
 ***  Print out a point-velocity (PV) record
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /FILES/ LUIN,LUOUT,I1,I2,I3,I4,I5, I6
 
       LATS = IDINT(SLAT*100.D0 + 0.5D0)
@@ -689,7 +689,7 @@ C*** then set its vertical velocity to 0.0
 ** Compute change in distance and change in azimuth.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
       COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
 
@@ -749,7 +749,7 @@ C      BAZ = Back azimuth   (radians, clockwise from north)
 C      S   = Distance (meters)
 C                                               
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)                               
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON/CONST/A,F,EPS2,EPS,AF,PI,TWOPI,RHOSEC                      
       DATA TOL/0.5D-14/                                                 
       R = 1.0D0-F                                                       
@@ -805,7 +805,7 @@ C
 C Convert Bluebook date to time in minutes
 
       IMPLICIT DOUBLE PRECISION (A-H, O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       CHARACTER    CARD*80
       CHARACTER   DATE*6
       COMMON /FILES/ LUIN,LUOUT, I1, I2, I3, I4, I5, I6
@@ -869,7 +869,7 @@ C Convert Bluebook date to time in minutes
 C Convert Bluebook date to time in years
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       character DATE*8
       COMMON /FILES/ LUIN,LUOUT, I1, I2, I3, I4, I5, I6
       READ(DATE,10) IYEAR,MONTH,IDAY
@@ -986,7 +986,7 @@ C
 
 *** Interactively obtain specifications for a grid.
 
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /FILES/ LUIN,LUOUT, I1, I2, I3, I4, I5, I6
       CHARACTER   NAMEG*10
       WRITE(LUOUT,200)
@@ -1076,7 +1076,7 @@ C
 *** Interactively obtain the specifications for a line.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /FILES/ LUIN,LUOUT,I1,I2,I3,I4,I5,I6
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
       CHARACTER   NAMEG*10
@@ -1174,7 +1174,7 @@ C *** PROGRAMMED FOR CDC-6600 BY LCDR L.PFEIFER NGS ROCKVILLE MD 20FEB75
 C *** MODIFIED FOR SYSTEM 360 BY JOHN G GERGEN NGS ROCKVILLE MD 750608
 C
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /CONST/ A,FINV,E2,EPI,AF,PI,TWOPI,RHOSEC
       DATA EPS/0.5D-13/
       R=1.D0-FINV
@@ -1225,7 +1225,7 @@ C
 *** Convert double precision real number to character*14
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       CHARACTER   CHAR14*14
 
       ORIG = ORIG*10000.D0
@@ -1241,7 +1241,7 @@ C
 ** Compute change in DX,DY,DZ-vector from time MIN1 to MIN2.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       character      PIDs*6
       parameter (nbbdim = 10000)
       COMMON /ARRAYS/ HT(nbbdim), LOC(nbbdim),PIDs(nbbdim)
@@ -1321,7 +1321,7 @@ C
 *         DUP = upward displacement in meters
 ************** 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
 *** Compute radii of curvature at fault center
       CALL RADII (EQLAT, RMER, RPAR)
@@ -1399,7 +1399,7 @@ C
 *            of normal dip slip
 *******************************************************************
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       LOGICAL VERT   
       PI = 3.141593D0
       TWOPI = PI + PI
@@ -1464,7 +1464,7 @@ C********1*********2*********3*********4*********5*********6*********7**
      &                  U3SS,U1DS,U2DS,U3DS)
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       LOGICAL VERT   
 
       YBAR = ETA*CDIP + Q*SDIP
@@ -1585,7 +1585,7 @@ C********1*********2*********3*********4*********5*********6*********7**
 C**** COMPUTES THE WEIGHTS FOR AN ELEMENT IN A GRID
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NUMGRD = 10)
       DIMENSION WEI(2,2)
       COMMON /CDGRID/ GRDLX(NUMGRD), GRDUX(NUMGRD), 
@@ -1673,7 +1673,7 @@ C********1*********2*********3*********4*********5*********6*********7**
 
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       DIMENSION VEL(2,2,3), B(*)
 
       DO 30 II = 0,1
@@ -1693,7 +1693,7 @@ C***************************************************
       SUBROUTINE RDEG (INPUT,VAL,CNEG)
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
       CHARACTER  INPUT*10
       CHARACTER  CNEG*1
@@ -1725,7 +1725,7 @@ C***************************************************
       INTEGER FUNCTION IUNGRD(IREGN, I, J, IVEC)
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NUMGRD = 10)
       COMMON /CDGRID/ GRDLX(NUMGRD), GRDUX(NUMGRD),
      1          GRDLY(NUMGRD), GRDUY(NUMGRD),
@@ -1743,7 +1743,7 @@ C***************************************************
 *** Convert velocities from vx,vy,vz to vn,ve,vu
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
       SLAT = DSIN(GLAT)
       CLAT = DCOS(GLAT)
@@ -1761,7 +1761,7 @@ C***************************************************
       SUBROUTINE TOVXYZ(GLAT,GLON,VN,VE,VU,VX,VY,VZ)
 *** Convert velocities from vn,ve,vu to vx,vy,vz
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
       SLAT = DSIN(GLAT)
       CLAT = DCOS(GLAT)
@@ -1781,7 +1781,7 @@ C***************************************************
 *** Predict displacements between two dates.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (numref = 19)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
       COMMON /FILES/ LUIN,LUOUT, I1, I2, I3, I4, I5, I6
@@ -2172,7 +2172,7 @@ C         IF(JN.EQ.'S' .OR. JW.EQ.'E')GO TO 320
 *** Compute velocities at specified locations
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (numref = 19)
       parameter (nrsrch = 0)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
@@ -2666,7 +2666,7 @@ C 129    FORMAT(F6.2, 1X, F6.2, 1X, F5.1)
 *** and ellipsoid height EHT (meters) in this reference frame.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       logical  Is_iopt_NAD83
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
 
@@ -2723,7 +2723,7 @@ c         CALL VTRANF(X,Y,Z,VX,VY,VZ, 1, IOPT)  !No longer NAD83
 *** datum with longitude positive west.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       LOGICAL FRMXYZ
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
 
@@ -2758,7 +2758,7 @@ c         CALL VTRANF(X,Y,Z,VX,VY,VZ, 1, IOPT)  !No longer NAD83
 *** Transform position between reference frames
 
       implicit double precision (a-h,o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
       parameter (nbbdim = 10000)
       parameter (rad2deg = 180.d0/3.14159265358979d0)
@@ -3522,7 +3522,7 @@ C*******************************************8***********************************
 ** Print updated and/or transformed parameters               
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
       CHARACTER    NAME24*24
       CHARACTER    LATDIR*1, LONDIR*1, LATDR*1, LONDR*1
@@ -3586,7 +3586,7 @@ C  However, the NAD 83 Pacific and Mariana frames use the IERS convention ITRF96
 C  Removed equivalence between NAD83(2011)and WGS84 original (Transit) in v3.4.0.
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /const/ a, f, e2, eps, af, pi, twopi, rhosec
@@ -4436,7 +4436,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** jopt --> input specifier of output reference frame
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /tranpa/ tx(numref), ty(numref), tz(numref), 
@@ -4486,7 +4486,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** jopt --> input specifier of output reference frame
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /tranpa1/ tx1(numref), ty1(numref), tz1(numref), 
@@ -4536,7 +4536,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** jopt --> input specifier of input reference frame
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /tranpa/ tx(numref), ty(numref), tz(numref), 
@@ -4586,7 +4586,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** jopt --> input specifier of input reference frame
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /tranpa1/ tx1(numref), ty1(numref), tz1(numref), 
@@ -4623,7 +4623,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 
 ** Write out options for reference frames
 
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       character    mframe*24, nframe*24
       dimension    nframe(24)
       dimension    iframe(24)
@@ -4753,7 +4753,7 @@ c      nframe(11)= 'SIO/MIT_92              ' !Included with ITRF91 in v3.4.0
 ** Update positions and/or observations to a specified date.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (numref = 19)
       parameter (nbbdim = 10000)
       CHARACTER    OLDBB*80, NEWBB*80, NAMEIF*80
@@ -5171,7 +5171,7 @@ c        ENDIF
 *** Retrieve geodetic coordinates from the Bluebook
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (nbbdim = 10000)
       CHARACTER    JN*1,JW*1
       CHARACTER    TYPE*4
@@ -5263,7 +5263,7 @@ C        HT(ISN) = HT(ISN) + GH
 *** Retrieve geodetic position from Bluebook with 5-digit SSN
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (nbbdim = 10000)
       CHARACTER   JN*1,JW*1
       CHARACTER   TYPE*3
@@ -5353,7 +5353,7 @@ C        HT(ISN) = HT(ISN) + GH
 *** Update Bluebook
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 
       parameter (nbbdim = 10000)
 
@@ -5681,7 +5681,7 @@ C    1          YLATT,YLONT,EHTNEW,DN,DE,DU,IOPT)
 *** Update 5-digit Bluebook
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (nbbdim = 10000)
       CHARACTER     CARD*212
       CHARACTER     SUBCRD*80
@@ -5980,7 +5980,7 @@ C    1          YLATT,YLONT,EHTNEW,DN,DE,DU,IOPT)
 *** Check if stations have *80* records
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (nbbdim = 10000)
       CHARACTER    CARD*80
       CHARACTER    PIDs*6
@@ -6009,7 +6009,7 @@ C*******************************************************************
 *** Update G-file of Bluebook
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       CHARACTER     CARD*120
       CHARACTER     TYPE*1
       CHARACTER     NRF*2
@@ -6210,7 +6210,7 @@ C  Initialize NAD 83 input/poutput logical variables as false:
 *** Update G-file of 5-digit Bluebook
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       CHARACTER     CARD*180
       CHARACTER     SUBCRD*80
       CHARACTER     TYPE*1
@@ -6358,7 +6358,7 @@ C        DECYR2 = DBLE(IYEAR2) + DBLE(MINO2 - MIN00)/525600.D0
 *** range is [-twopi to +twopi]
  
       implicit double precision(a-h,o-z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       common/CONST/A,F,E2,EP2,AF,PI,TWOPI,RHOSEC
  
     1 if(val.gt.twopi) then
@@ -6404,7 +6404,7 @@ C        DECYR2 = DBLE(IYEAR2) + DBLE(MINO2 - MIN00)/525600.D0
 *** a reference frame identifier in the Bluebook
 *** to a reference frame identifier in HTDP and back
 
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter ( numref = 19 )
       COMMON /REFCON/ IRFCON(40), JRFCON(numref)
 
@@ -6604,7 +6604,7 @@ C coordinates as if they were NAD 83(2011)coordinates.
 *** system used in the Bluebook to the
 *** system used in HTDP
 
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter ( numref = 19 )
       COMMON /REFCON/ IRFCON(40), JRFCON(numref)
 
@@ -6627,7 +6627,7 @@ C coordinates as if they were NAD 83(2011)coordinates.
 *** system used in HTDP to the system
 *** used in the Bluebook
 
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter ( numref = 19 )
       COMMON /REFCON/ IRFCON(40), JRFCON(numref)
 
@@ -6676,7 +6676,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** jopt2 --> (input) specifier of reference frame for output vector
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /tranpa/ tx(numref), ty(numref), tz(numref), 
@@ -6743,7 +6743,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** jopt2 --> (input) specifier of reference frame for output vector
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
 
       common /tranpa1/ tx1(numref), ty1(numref), tz1(numref), 
@@ -6798,7 +6798,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** Output longitude will be positive west.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       CHARACTER    NAME*24
       CHARACTER    COPT*1,LATDIR*1,LONDIR*1
       LOGICAL      FRMXYZ
@@ -7106,7 +7106,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** Interactively obtain velocity for a point
 
       IMPLICIT DOUBLE PRECISION (A-H, O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       CHARACTER   VOPT*1
       COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
 
@@ -7219,7 +7219,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** Input VN, VE, and VU are in mm/yr
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NDLOC = 2740)
 
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
@@ -7287,7 +7287,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** Predict displacements from time MIN1 to time MIN2.
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
 
       HT = HTOLD
@@ -7326,7 +7326,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 ** vu         output - upward velocity in mm/yr
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       logical  Is_iopt_NAD83
       COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC
 
@@ -7381,7 +7381,7 @@ C  The following 2 lines were added on 07/22/2015 after Rich found this bug
 *** Transform velocities from one reference frame to another
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       parameter (numref = 19)
       character    nameif*80,name24*80
       character    NAMEF*80
@@ -7563,7 +7563,7 @@ C  The parameters in common block tranpa are computed using the IGS values of IT
 C  The parameters in common block tranpa1 are computed using the IERS values of ITRF96==>ITRF97
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (numref = 19)
       common /tranpa/ tx(numref), ty(numref), tz(numref),
      &                dtx(numref), dty(numref), dtz(numref),
@@ -7627,7 +7627,7 @@ C  The parameters in common block tranpa are computed using the IGS values of IT
 C  The parameters in common block tranpa1 are computed using the IERS values of ITRF96==>ITRF97
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (numref = 19)
       common /tranpa1/ tx1(numref), ty1(numref), tz1(numref), 
      &                dtx1(numref), dty1(numref), dtz1(numref),
@@ -7686,7 +7686,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 *** Print transformed velocities
 
       implicit double precision (a-h, o-z)
-      implicit integer*4 (i-n)
+      implicit integer(kind=4) (i-n)
       character    name24*80
       COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
 
@@ -7719,7 +7719,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
        SUBROUTINE HEADER
 
        IMPLICIT DOUBLE PRECISION (A-H, O-Z)
-       IMPLICIT INTEGER*4 (I-N)
+       IMPLICIT INTEGER(kind=4) (I-N)
        character  HTDP_version*8
        COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6
        COMMON /VERSION/ HTDP_version                        
@@ -7741,7 +7741,7 @@ C  The parameters in common block tranpa1 are computed using the IERS values of
 ***    TEST       output - if (true) then there is an error
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       DIMENSION M(12)
       LOGICAL TEST
       CHARACTER   TOPT*1
@@ -7939,9 +7939,9 @@ C::9004.17, MSS, CHANGE ORDER YY MM DD
 C********1*********2*********3*********4*********5*********6*********7**
 C
       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
 C
-      INTEGER*4     A, B, C, D
+      INTEGER(kind=4)     A, B, C, D
 
       IYRP = IYR
 C
@@ -7999,7 +7999,7 @@ C      WRITE(*,*)
 *******
 
       IMPLICIT DOUBLE PRECISION (A-H, O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       LOGICAL INSIDE
 
       parameter (NUMPSG = 1)
@@ -8086,7 +8086,7 @@ C ROUTINE CHECKS IF THE POINT HAVING COORDINATES (POSX, POSY)
 C IS WITHIN THE REGION SPANNED BY THE GRID
 C
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       LOGICAL INSIDE
 
       INSIDE = .TRUE.
@@ -8150,7 +8150,7 @@ C********1*********2*********3*********4*********5*********6*********7**
 C**** COMPUTES THE WEIGHTS FOR AN ELEMENT IN A GRID
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NUMPSG = 1)
       DIMENSION WEI(2,2)
       COMMON /PSGRID/ PSGLX(NUMPSG), PSGUX(NUMPSG), 
@@ -8229,7 +8229,7 @@ C********1*********2*********3*********4*********5*********6*********7**
 
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       DIMENSION AMP(2,2,3), PS(*)
 
       DO 30 II = 0,1
@@ -8248,7 +8248,7 @@ C***************************************************
       INTEGER FUNCTION IPSGRD(IGRID, I, J, IVEC)
 
       IMPLICIT DOUBLE PRECISION (A-H,O-Z)
-      IMPLICIT INTEGER*4 (I-N)
+      IMPLICIT INTEGER(kind=4) (I-N)
       parameter (NUMPSG = 1)
       COMMON /PSGRID/ PSGLX(NUMPSG), PSGUX(NUMPSG),
      1          PSGLY(NUMPSG), PSGUY(NUMPSG),
@@ -8266,7 +8266,7 @@ C  In f90, use trim(name) to truncate all spaces after the name.
 C  But "trim" were not available in f77.
 
       implicit none
-      integer*4 i
+      integer(kind=4) i
       character name*80,scratch*80
 
       do i=80,1,-1
@@ -8282,8 +8282,8 @@ C-------------------------------------------------------------------------------
 
          implicit none
 
-         integer*4   i,length
-         real*8      x,y,z
+         integer(kind=4)   i,length
+         real(kind=8)      x,y,z
          character   name*80,record*120,record1*120,chars*120
          character   xxxx*80,yyyy*80,zzzz*80
 
@@ -8397,8 +8397,8 @@ C-------------------------------------------------------------------------------
 
          implicit none
 
-         integer*4   i,length
-         real*8      x,y,z,Vx,Vy,Vz
+         integer(kind=4)   i,length
+         real(kind=8)      x,y,z,Vx,Vy,Vz
          character   name*80,record*120,record1*120,chars*120
          character   xxxx*80,yyyy*80,zzzz*80
 
@@ -8595,8 +8595,8 @@ C-------------------------------------------------------------------------------
 
          implicit none
 
-         integer*4   i,length
-         real*8      x,y
+         integer(kind=4)   i,length
+         real(kind=8)      x,y
          character   name*24,record*120,record1*120,chars*120
          character   xxxx*80,yyyy*80
 
@@ -8682,8 +8682,8 @@ C-------------------------------------------------------------------------------
 
          implicit none
 
-         integer*4   i,length
-         real*8      x,y,vn,ve,vu
+         integer(kind=4)   i,length
+         real(kind=8)      x,y,vn,ve,vu
          character   name*80,record*120,record1*120,chars*120
          character   xxxx*80,yyyy*80,vnnn*80,veee*80,vuuu*80