aboutsummaryrefslogtreecommitdiff
path: root/htdp.f
diff options
context:
space:
mode:
Diffstat (limited to 'htdp.f')
-rw-r--r--htdp.f160
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
generated by cgit v1.2.3 (git 2.25.1) at 2024-12-26 06:51:13 +0000