From 8dbb4a54d947c4edfa1b74d8e0b7ea5d2181a441 Mon Sep 17 00:00:00 2001 From: John Denker Date: Wed, 20 Nov 2024 11:54:12 -0700 Subject: Improve standards compliance: integer(kind=4) and real(kind=8). Note that integer*4 and real*8 are not actually part of the language standard, and never have been. --- htdp.f | 160 ++++++++++++++++++++++++++++++++--------------------------------- 1 file 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 -- cgit v1.2.3