diff options
52 files changed, 1994 insertions, 887 deletions
diff --git a/HTDP-log.pdf b/HTDP-log.pdf Binary files differindex be8a8df..bd15a95 100755 --- a/HTDP-log.pdf +++ b/HTDP-log.pdf diff --git a/HTDP-user-guide.pdf b/HTDP-user-guide.pdf Binary files differindex 850c32f..a6fa0bc 100644 --- a/HTDP-user-guide.pdf +++ b/HTDP-user-guide.pdf @@ -28,5 +28,3 @@ initps.o: initps.f htdp: htdp.o initvl.o initeq.o initbd.o initps.o $(FC) $(CFLAGS) $(INCLUDES) -o $@ $? $(LIBS) - -#MADE diff --git a/htdp.exe b/htdp.exe Binary files differnew file mode 100644 index 0000000..9d7a1d8 --- /dev/null +++ b/htdp.exe @@ -3,16 +3,18 @@ ************************************************************** * NAME: HTDP (Horizontal Time-Dependent Positioning) * -* WRITTEN BY: Richard A. Snay & Chris Pearson +* WRITTEN BY: Richard Snay, Chris Pearson, Jarir Saleh, +* and Michael Dennis * * PURPOSE: Transform coordinates across time * and between reference frames * -**************************************************************** +************************************************************** IMPLICIT DOUBLE PRECISION (A-H,O-Z) IMPLICIT INTEGER*4 (I-N) - character HTDP_version*10 + character HTDP_version*8 + character Version_date*20 CHARACTER OPTION*1 character cont*5 @@ -20,9 +22,10 @@ COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6 COMMON /VERSION/ HTDP_version -C You must change HTDP version here if necessary +C You must change HTDP version and date here, if necessary - HTDP_version = 'v3.2.8' + HTDP_version = '3.2.9' + Version_date = 'May 17, 2020' *** Introduce variables for file id's @@ -61,17 +64,17 @@ C 1 FORM='UNFORMATTED') *** Initialize conversion table between reference frame identifiers CALL SETRF - WRITE(LUOUT,5) HTDP_version + WRITE(LUOUT,5) HTDP_version, Version_date 5 FORMAT( 1 '********************************************************'/ - 1 '* HTDP (Horizontal Time-Dependent Positioning) *'/ - 1 '* SOFTWARE VERSION ',a10 / - 1 '* *'/) + 1 ' HTDP (Horizontal Time-Dependent Positioning) '/ + 1 ' SOFTWARE VERSION: ', a8 / + 1 ' VERSION DATE: ', a20 /) WRITE(LUOUT,501) 501 FORMAT( - 1 '* AUTHORS: Richard Snay, Chris Pearson & Jarir Saleh *'/ - 1 '* Email: ngs.cors.htdp@noaa.gov *'/ - 1 '* *'/ + 1 ' AUTHORS: Richard Snay, Chris Pearson, Jarir Saleh, '/ + 1 ' and Michael Dennis '/ + 1 ' Email: ngs.cors.htdp@noaa.gov '/ 1 '********************************************************'/) WRITE(LUOUT,10) 10 FORMAT( @@ -90,11 +93,12 @@ C 1 FORM='UNFORMATTED') 26 FORMAT(' ***************************************'/ 1 ' MAIN MENU:',/ 6 ' 0... Exit software.',/ - 7 ' 1... Estimate displacements between two dates.'/ - 8 ' 2... Estimate velocities.'/ - 9 ' 3... Update positions and/or observations ' - & ,'to a specified date.'/ - & ' 4... Transform positions between reference frames. '/ + 7 ' 1... Estimate horizontal displacements between two dates.'/ + 8 ' 2... Estimate horizontal velocities.'/ + 9 ' 3... Transform observations to a specified reference ' + & ,'frame and/or date.'/ + & ' 4... Transform positions between reference frames ' + & ,'and/or dates.'/ & ' 5... Transform velocities between reference frames. ') 30 READ(LUIN,35,err=52,iostat=ios) OPTION if (ios /= 0) goto 52 @@ -172,7 +176,7 @@ C*** Set default reference epoch to Jan. 1, 2010 *** Region 1 is the San Andreas fault in central California *** Region 2 is southern California *** Region 3 is Northern California -*** Region 4 is the Pacific Noerthwest +*** Region 4 is the Pacific Northwest *** Region 5 is western CONUS *** Region 6 is CONUS *** Region 7 is St. Elias, Alaska @@ -185,7 +189,7 @@ C*** Set default reference epoch to Jan. 1, 2010 *** Region 14 is the Juan de Fuca plate *** Region 15 is the Cocos plate *** Region 16 is the Mariana plate -*** REGION 17 is the Philippine Sea plate +*** Region 17 is the Philippine Sea plate IMPLICIT DOUBLE PRECISION (A-H,O-Z) IMPLICIT INTEGER*4 (I-N) @@ -1839,7 +1843,7 @@ C CHARACTER CARD*80 CHARACTER record*120 CHARACTER NAMEF*30,NAME*30,NAMEBB*30, NAMEIF*30 - CHARACTER NAME24 + CHARACTER NAME24*24 CHARACTER BLAB*17 CHARACTER NAMEG*10 CHARACTER TYPE*4 @@ -2275,7 +2279,7 @@ c 6 ' a specified point having a specified velocity.') call MENU1(iopt, frame1) if (iopt .ge. 1 .and. iopt .le. numref) then WRITE(I2,1002) frame1 - 1002 FORMAT('VELOCITIES IN MM/YR RELATIVE TO ', a24 /) + 1002 FORMAT('VELOCITIES IN MM/YR RELATIVE TO ', a24 ) RVN = 0.0D0 RVE = 0.0D0 RVU = 0.0D0 @@ -2320,7 +2324,7 @@ c 6 ' a specified point having a specified velocity.') CALL TOVXYZ(YLAT,GLON,SVN,SVE,SVU,SVX,SVY,SVZ) WRITE(I2,1060) NAME24,LATD,LATM,SLAT,LATDIR,SVN,LOND,LONM, 1 SLON,LONDIR,SVE,EHT,SVU,X,SVX,Y,SVY,Z,SVZ - 1060 FORMAT(/10X,A24,/ + 1060 FORMAT(/A24,/ 1'LATITUDE = ',2I3,F9.5,1X,A1,' NORTH VELOCITY =',F7.2,' mm/yr'/ 2'LONGITUDE = ',2I3,F9.5,1X,A1,' EAST VELOCITY =',F7.2,' mm/yr'/ 3'ELLIPS. HT. = ',F10.3,' m', 6X,'UP VELOCITY =',F7.2,' mm/yr'/ @@ -2415,7 +2419,7 @@ c 6 ' a specified point having a specified velocity.') J = -1 120 J = J + 1 YLAT = XLAT - ON = MINLON + J*JDS + LON = MINLON + J*JDS IF(LON .GT. MAXLON) GO TO 110 YLON = DBLE(LON)/RHOSEC CALL TODMSS(YLON,LOND,LONM,SLON,ISIGN) @@ -2437,16 +2441,16 @@ c 6 ' a specified point having a specified velocity.') TOTVEL = DSQRT(VN*VN + VE*VE) *** code to create vectors to be plotted -c IF (TOTVEL .GE. 5.0D0 .AND. TOTVEL .LE. 50.D0) THEN -c TVE = VE/10.D0 -c TVN = VN/10.D0 -c WRITE(I2,129) ZLON,ZLAT,TVE,TVN,DUMMY,DUMMY,DUMMY, +c IF (TOTVEL .GE. 5.0D0 .AND. TOTVEL .LE. 50.D0) THEN +c TVE = VE/10.D0 +c TVN = VN/10.D0 +c WRITE(I2,129) ZLON,ZLAT,TVE,TVN,DUMMY,DUMMY,DUMMY, c 1 TEMPNA c 129 FORMAT(F10.6, 2x, F9.6, 2X, 5(F7.2, 2x), A8) -c ENDIF +c ENDIF *** code to create contour plots -C WRITE(I2,129) ZLON, ZLAT, TOTVEL +C WRITE(I2,129) ZLON, ZLAT, TOTVEL C 129 FORMAT(F6.2, 1X, F6.2, 1X, F5.1) *** end of temporary code @@ -2813,7 +2817,6 @@ c write (*,*) "FROM XT08 ",RLAT*180.d0/pi,WLON*180.d0/pi,EHT08 parameter (nbbdim = 10000) parameter (rad2deg = 180.d0/3.14159265358979d0) - double precision lat,lon character card*80,namebb*80,nameif*80,name24*80 character record*120 character namef*30 @@ -2821,21 +2824,12 @@ c write (*,*) "FROM XT08 ",RLAT*180.d0/pi,WLON*180.d0/pi,EHT08 character jn*1,jw*1,LATDIR*1,LONDIR*1 character option*1, answer*1, vopt*1 character PID*6,PIDs*6 - character HTDP_version*10 + character HTDP_version*8 LOGICAL FRMXYZ LOGICAL TEST LOGICAL Is_inp_NAD83,Is_out_NAD83 LOGICAL Is_inp_NAD83PAC,Is_out_NAD83PAC LOGICAL Is_inp_NAD83MAR,Is_out_NAD83MAR - logical am_I_in_or_near_AK - logical am_I_in_or_near_AS - logical am_I_in_or_near_CONUS - logical am_I_in_or_near_CQ - logical am_I_in_or_near_Guam - logical am_I_in_or_near_HI - logical am_I_in_or_near_PR - logical am_I_in_or_near_VQ - logical am_I_in_or_near_KW COMMON /CONST/ A, F, E2, EPS, AF, PI, TWOPI, RHOSEC COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6 @@ -2956,48 +2950,6 @@ c write (*,*) "Starting to look at options" 180 call GETPNT(latd, latm, slat, LATDIR, lond, lonm, slon, 1 LONDIR, name24, x, y, z, ylat, ylon, eht) -C Added by JS on 07/22/2015 to limit NAD83 to the US***************************************************** -C This was reversed on 05/11/2017 ***************************************************** - -c if (Is_inp_NAD83 .or. Is_out_NAD83) then -c lat = ylat*rad2deg ; lon = 360.d0 - ylon*rad2deg -c if (lon < 0.d0) lon = lon + 360.d0 -c if(.not. am_I_in_or_near_AK (lat,lon) .and. -c & .not. am_I_in_or_near_CONUS(lat,lon) .and. -c & .not. am_I_in_or_near_PR (lat,lon) .and. -c & .not. am_I_in_or_near_VQ (lat,lon)) then -c write (luout,'(10x,a,8x,a)') -c & "NAD83 (2011) is defined only in US territories",name24 -c stop -c endif -c endif - -c if(.not. am_I_in_or_near_AK (lat,lon) .and. -c & .not. am_I_in_or_near_AS (lat,lon) .and. -c & .not. am_I_in_or_near_CONUS(lat,lon) .and. -c & .not. am_I_in_or_near_CQ (lat,lon) .and. -c & .not. am_I_in_or_near_Guam (lat,lon) .and. -c & .not. am_I_in_or_near_HI (lat,lon) .and. -c & .not. am_I_in_or_near_PR (lat,lon) .and. -c & .not. am_I_in_or_near_VQ (lat,lon) .and. -c & .not. am_I_in_or_near_KW (lat,lon)) then -c write (luout,'(10x,a,8x,a)') -c & "FYI, NAD83 (2011) is defined only in US territories",name24 -c stop -c endif - -c if (Is_inp_NAD83MAR .or. Is_out_NAD83MAR) then -c lat = ylat*rad2deg ; lon = 360.d0 - ylon*rad2deg -c if (lon < 0.d0) lon = lon + 360.d0 -c if(.not. am_I_in_or_near_GUAM (lat,lon) .and. -c & .not. am_I_in_or_near_CQ (lat,lon)) then -c write (luout,'(10x,a,8x,a)') -c & "NAD83 (2011) is defined only in US territories",name24 -c stop -c endif -c endif - -C Until here on 07/22/2015 to limit NAD83 to the US****************************************************** C Added by JS on 09/10/2014 @@ -3125,7 +3077,7 @@ c write (I2,1309) trim(namebb) write (i2,1310) HTDP_version 1310 format (' ***CAUTION: This file was processed using HTDP', - & ' version ',a10, '***') + & ' version ',a8, '***') write (i2,1311) frame2 1311 format (' ***CAUTION: Coordinates in this file are in ', & a24, '***') @@ -3344,28 +3296,8 @@ c 1 name24,0,vxsave,vysave,vzsave,vnsave,vesave,vusave) outlon = ylont*rad2deg call extract_name (name24,iii) -C Added by JS on 07/22/2015 to limit the use in NAD83 to the US -c Reversed on 05/11/2017 - -c if (Is_inp_NAD83 .or. Is_out_NAD83) then -c lat = ylat*rad2deg ; lon = 360.d0 - ylon*rad2deg -c if(.not. am_I_in_or_near_AK (lat,lon) .and. -c & .not. am_I_in_or_near_AS (lat,lon) .and. -c & .not. am_I_in_or_near_CONUS(lat,lon) .and. -c & .not. am_I_in_or_near_CQ (lat,lon) .and. -c & .not. am_I_in_or_near_Guam (lat,lon) .and. -c & .not. am_I_in_or_near_HI (lat,lon) .and. -c & .not. am_I_in_or_near_PR (lat,lon) .and. -c & .not. am_I_in_or_near_VQ (lat,lon) .and. -c & .not. am_I_in_or_near_KW (lat,lon)) then -c write (i2,'(10x,a,8x,a)') -c & "FYI, NAD83(2011) is defined only in US territories",name24 -c endif -c endif write (i2,449) outlat,outlon,ehtnew,name24(1:iii) -C Until here on 07/22/2015 to limit the use in NAD83 to the US - c write (i6,449) outlat,outlon,ehtnew,trim(name24) 449 format (2f16.10,f10.3,4x,a) go to 410 @@ -3462,28 +3394,7 @@ c call PRNTTP(x,y,z,xt,yt,zt,ylat,ylatt,ylon,ylont,eht,ehtnew, c 1 name24,0,vxsave,vysave,vzsave,vnsave,vesave,vusave) call extract_name (name24,iii) -C Added by JS on 07/22/2015 to limit the use in NAD83 to the US - -c if (Is_inp_NAD83 .or. Is_out_NAD83) then -c lat = ylatt*rad2deg ; lon = 360.d0 - ylont*rad2deg -c if(.not. am_I_in_or_near_AK (lat,lon) .and. -c & .not. am_I_in_or_near_AS (lat,lon) .and. -c & .not. am_I_in_or_near_CONUS(lat,lon) .and. -c & .not. am_I_in_or_near_CQ (lat,lon) .and. -c & .not. am_I_in_or_near_Guam (lat,lon) .and. -c & .not. am_I_in_or_near_HI (lat,lon) .and. -c & .not. am_I_in_or_near_PR (lat,lon) .and. -c & .not. am_I_in_or_near_VQ (lat,lon) .and. -c & .not. am_I_in_or_near_KW (lat,lon)) then -c write (i2,'(10x,a,8x,a)') -c & "FYI, NAD83(2011) is defined only in US territories",name24 -c write (i2,1449) xt,yt,zt,name24(1:iii) -c endif -c else - write (i2,1449) xt,yt,zt,name24(1:iii) -c endif - -C Until here on 07/22/2015 to limit the use in NAD83 to the US + write (i2,1449) xt,yt,zt,name24(1:iii) c write (i6,1449) xt,yt,zt,trim(name24) 1449 format (3f20.3,4x,a) @@ -3582,28 +3493,7 @@ c call PRNTTP(x,y,z,xt,yt,zt,ylat,ylatt,ylon,ylont,eht,ehtnew, c 1 name24,0,vxsave,vysave,vzsave,vnsave,vesave,vusave) call extract_name (name24,iii) -C Added by JS on 07/22/2015 to limit the use in NAD83 to the US - -c if (Is_inp_NAD83 .or. Is_out_NAD83) then -c lat = ylatt*rad2deg ; lon = 360.d0 - ylont*rad2deg -c if(.not. am_I_in_or_near_AK (lat,lon) .and. -c & .not. am_I_in_or_near_AS (lat,lon) .and. -c & .not. am_I_in_or_near_CONUS(lat,lon) .and. -c & .not. am_I_in_or_near_CQ (lat,lon) .and. -c & .not. am_I_in_or_near_Guam (lat,lon) .and. -c & .not. am_I_in_or_near_HI (lat,lon) .and. -c & .not. am_I_in_or_near_PR (lat,lon) .and. -c & .not. am_I_in_or_near_VQ (lat,lon) .and. -c & .not. am_I_in_or_near_KW (lat,lon)) then -c write (i2,'(10x,a,8x,a)') -c & "FYI, NAD83(2011) is defined only in US territories",name24 -c write (i2,1449) xt,yt,zt,name24(1:iii) -c endif -c else - write (i2,1449) xt,yt,zt,name24(1:iii) -c endif - -C Until here on 07/22/2015 to limit the use in NAD83 to the US + write (i2,1449) xt,yt,zt,name24(1:iii) c write (i6,1449) xt,yt,zt,trim(name24) go to 511 @@ -3737,11 +3627,11 @@ C*******************************************8*********************************** *** to other reference frames ************************************** C Important note: -C The parameters in common block tranpa are computed using the IGS values of ITRF96==>ITRF97 -C The parameters in common block tranpa1 are computed using the IERS values of ITRF96==>ITRF97 +C The parameters in common block tranpa are computed using the IGS convention ITRF96 <> ITRF97. +C The parameters in common block tranpa1 are computed using the IERS convention ITRF96 = ITRF97. C The latter parameters were added to HTDP in 09/2014. They will be used to transform between C ITRF systems. They will not be used if the transformation involves NAD83 or WGS84 (transit). -C They will be used for the Pacific branches of NAD83. +C However, the NAD 83 Pacific frames use the IERS convention ITRF96 = ITRF97. implicit double precision (a-h, o-z) @@ -3760,12 +3650,12 @@ C They will be used for the Pacific branches of NAD83. & drx1(numref), dry1(numref), drz1(numref), & scale1(numref), dscale1(numref), refepc1(numref) -C Parameters computed with the IGS values of ITRF96==>ITRF97 +C Parameters computed with the IGS convention ITRF96 <> ITRF97 *** From ITRF94 to NAD 83 tx(1) = 0.9910d0 ty(1) = -1.9072d0 - tz(1) = -.5129d0 + tz(1) = -0.5129d0 dtx(1) = 0.d0 dty(1) = 0.d0 dtz(1) = 0.d0 @@ -3773,8 +3663,8 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 ry(1) = 0.46785d-7 rz(1) = 0.56529d-7 drx(1) = 0.00258d-7 - dry(1) = -.03599d-7 - drz(1) = -.00153d-7 + dry(1) = -0.03599d-7 + drz(1) = -0.00153d-7 scale(1) = 0.d0 dscale(1) = 0.0d0 refepc(1) = 1997.0d0 @@ -3782,11 +3672,11 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF88 tx(2) = 0.018d0 ty(2) = 0.000d0 - tz(2) = -.092d0 + tz(2) = -0.092d0 dtx(2) = 0.0d0 dty(2) = 0.0d0 dtz(2) = 0.0d0 - rx(2) = -.0001d0 / rhosec + rx(2) = -0.0001d0 / rhosec ry(2) = 0.0d0 rz(2) = 0.0d0 drx(2) = 0.0d0 @@ -3799,7 +3689,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF89 tx(3) = 0.023d0 ty(3) = 0.036d0 - tz(3) = -.068d0 + tz(3) = -0.068d0 dtx(3) = 0.0d0 dty(3) = 0.0d0 dtz(3) = 0.0d0 @@ -3816,7 +3706,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF90 tx(4) = 0.018d0 ty(4) = 0.012d0 - tz(4) = -.030d0 + tz(4) = -0.030d0 dtx(4) = 0.0d0 dty(4) = 0.0d0 dtz(4) = 0.0d0 @@ -3833,7 +3723,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF91 tx(5) = 0.020d0 ty(5) = 0.016d0 - tz(5) = -.014d0 + tz(5) = -0.014d0 dtx(5) = 0.0d0 dty(5) = 0.0d0 dtz(5) = 0.0d0 @@ -3850,7 +3740,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF92 tx(6) = 0.008d0 ty(6) = 0.002d0 - tz(6) = -.008d0 + tz(6) = -0.008d0 dtx(6) = 0.0d0 dty(6) = 0.0d0 dtz(6) = 0.0d0 @@ -3860,23 +3750,23 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 drx(6) = 0.0d0 dry(6) = 0.0d0 drz(6) = 0.0d0 - scale(6) = -.08d-8 + scale(6) = -0.08d-8 dscale(6) = 0.0d0 refepc(6) = 1988.0d0 *** From ITRF94 to ITRF93 tx(7) = 0.006d0 - ty(7) = -.005d0 - tz(7) = -.015d0 - dtx(7) = -.0029d0 + ty(7) = -0.005d0 + tz(7) = -0.015d0 + dtx(7) = -0.0029d0 dty(7) = 0.0004d0 dtz(7) = 0.0008d0 rx(7) = 0.00039d0 / rhosec - ry(7) = -.00080d0 / rhosec + ry(7) = -0.00080d0 / rhosec rz(7) = 0.00096d0 / rhosec - drx(7) = .00011d0 / rhosec - dry(7) = .00019d0 / rhosec - drz(7) =-.00005d0 / rhosec + drx(7) = 0.00011d0 / rhosec + dry(7) = 0.00019d0 / rhosec + drz(7) = -0.00005d0 / rhosec scale(7) = 0.04d-8 dscale(7) = 0.0d0 refepc(7) = 1988.0d0 @@ -3898,8 +3788,8 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dscale(8) = 0.0d0 refepc(8) = 1996.0d0 -*** From ITRF94 to ITRF97 (based on IGS adopted values) -*** According to IERS: ITRF97 = ITRF96 = ITRF94 +*** From ITRF94 to ITRF97 (also IGS97), based on IGS adopted values +*** According to IGS: ITRF97 <> ITRF96 = ITRF94 tx(9) = 0.00207d0 ty(9) = 0.00021d0 tz(9) = -0.00995d0 @@ -3919,7 +3809,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 *** From ITRF94 to WGS 72 (composition of ITRF94 -> NAD_83 -> WGS_72) tx(10) = 0.9910d0 ty(10) = -1.9072d0 - tz(10) = -.5129d0 - 4.5d0 + tz(10) = -0.5129d0 - 4.5d0 dtx(10) = 0.d0 dty(10) = 0.d0 dtz(10) = 0.d0 @@ -3927,23 +3817,23 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 ry(10) = 0.46785d-7 rz(10) = 0.56529d-7 + 26.85868d-7 drx(10) = 0.00258d-7 - dry(10) = -.03599d-7 - drz(10) = -.00153d-7 + dry(10) = -0.03599d-7 + drz(10) = -0.00153d-7 scale(10) = 0.d0 - 0.2263d-6 dscale(10) = 0.0d0 refepc(10) = 1997.0d0 -*** From ITRF94 to ITRF00 +*** From ITRF94 to ITRF2000 (also IGS00 and IGb00) *** assumes that ITRF94 = ITRF96 and *** uses IGS values for ITRF96 -> ITRF97 -*** and IERS values for ITRF97 -> ITRF00 - tx(11) = -.00463d0 - ty(11) = -.00589d0 - tz(11) = +.00855d0 +*** and IERS values for ITRF97 -> ITRF2000 + tx(11) = -0.00463d0 + ty(11) = -0.00589d0 + tz(11) = 0.00855d0 dtx(11) = -0.00069d0 dty(11) = 0.00070d0 dtz(11) = -0.00046d0 - rx(11) = -.00012467d0 / rhosec + rx(11) = -0.00012467d0 / rhosec ry(11) = 0.00022355d0 / rhosec rz(11) = 0.00006065d0 / rhosec drx(11) = -0.00001347d0 / rhosec @@ -3954,19 +3844,19 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 refepc(11) = 1997.0d0 *** From ITRF94 to PACP00 -*** use PA/ITRF00 rotation rates from Beavan et al., (2002) +*** use PA/ITRF2000 rotation rates from Beavan et al., (2002) tx(12) = 0.9056d0 ty(12) = -2.0200d0 tz(12) = -0.5516d0 - dtx(12) = -.00069d0 + dtx(12) = -0.00069d0 dty(12) = 0.00070d0 dtz(12) = -0.00046d0 rx(12) = 0.027616d0 / rhosec ry(12) = 0.013692d0 / rhosec rz(12) = 0.002773d0 / rhosec - drx(12) = -.000397d0 / rhosec + drx(12) = -0.000397d0 / rhosec dry(12) = 0.001022d0 / rhosec - drz(12) = -.002166d0 / rhosec + drz(12) = -0.002166d0 / rhosec scale(12) = -0.61504d-9 dscale(12) = 0.18201d-9 refepc(12) = 1997.0d0 @@ -3979,8 +3869,8 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dtx(13) = -0.00069d0 dty(13) = 0.00070d0 dtz(13) = -0.00046d0 - rx(13) = .028847d0 / rhosec - ry(13) = .010644d0 / rhosec + rx(13) = 0.028847d0 / rhosec + ry(13) = 0.010644d0 / rhosec rz(13) = 0.008989d0 / rhosec drx(13) = -0.000033d0 / rhosec dry(13) = 0.000120d0 / rhosec @@ -3989,7 +3879,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dscale(13) = 0.18201d-9 refepc(13) = 1997.00d0 -*** From ITRF94 to ITRF2005 +*** From ITRF94 to ITRF2005 (also IGS05) *** assumes that ITRF94 = ITRF96 *** uses IGS values for ITRF96 -> ITRF97 *** uses IERS values for ITRF97 -> ITRF2000 @@ -4000,7 +3890,7 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dtx(14) = -0.00049d0 dty(14) = 0.00060d0 dtz(14) = 0.00134d0 - rx(14) = -.00012467d0 / rhosec + rx(14) = -0.00012467d0 / rhosec ry(14) = 0.00022355d0 / rhosec rz(14) = 0.00006065d0 / rhosec drx(14) = -0.00001347d0 / rhosec @@ -4010,12 +3900,12 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dscale(14) = 0.10201d-9 refepc(14) = 1997.0d0 -*** From ITRF94 to ITRF2008 (also IGS08 and IGB08) +*** From ITRF94 to ITRF2008 (also IGS08 and IGb08) *** assumes that ITRF94 = ITRF96 *** uses IGS values for ITRF96 -> ITRF97 *** uses IERS values for ITRF97 -> ITRF2000 *** uses IERS values for ITRF2000-> ITRF2005 -*** uses IERS values for ITRF2005 -> ITRF2008 (and IGS08 and IGB08) +*** uses IERS values for ITRF2005 -> ITRF2008 tx(15) = -0.00243d0 ty(15) = -0.00389d0 @@ -4023,9 +3913,9 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dtx(15) = -0.00079d0 dty(15) = 0.00060d0 dtz(15) = 0.00134d0 - rx(15) = -0.00012467d0 / rhosec - ry(15) = 0.00022355d0 / rhosec - rz(15) = 0.00006065d0 / rhosec + rx(15) = -0.00012467d0 / rhosec + ry(15) = 0.00022355d0 / rhosec + rz(15) = 0.00006065d0 / rhosec drx(15) = -0.00001347d0 / rhosec dry(15) = 0.00001514d0 / rhosec drz(15) = 0.00001973d0 / rhosec @@ -4033,47 +3923,37 @@ C Parameters computed with the IGS values of ITRF96==>ITRF97 dscale(15) = 0.10201d-9 refepc(15) = 1997.0d0 -*** From ITRF94 to ITRF2014 +*** From ITRF94 to ITRF2014 (also IGS14 and IGb14) *** assumes that ITRF94 = ITRF96 *** uses IGS values for ITRF96 -> ITRF97 *** uses IERS values for ITRF97 -> ITRF2000 *** uses IERS values for ITRF2000-> ITRF2005 -*** uses IERS values for ITRF2005 -> ITRF2008 (and IGS08 and IGB08) +*** uses IERS values for ITRF2005 -> ITRF2008 *** uses IERS values for ITRF2008 -> ITRF2014 -c tx(16) = -0.00403d0 !Differs from ITRF2008 tx(16) = -0.01430d0 -c ty(16) = -0.00579d0 !Differs from ITRF2008 ty(16) = 0.00201d0 -c tz(16) = 0.01125d0 !Differs from ITRF2008 tz(16) = 0.02867d0 - - dtx(16) = -0.00079d0 !Like ITRF2008 - dty(16) = 0.00060d0 !Like ITRF2008 - dtz(16) = 0.00144d0 !Differs from ITRF2008 - -c rx(16) = -0.00012467d0 / rhosec !Like ITRF2008 + dtx(16) = -0.00079d0 + dty(16) = 0.00060d0 + dtz(16) = 0.00144d0 rx(16) = -0.00029978d0 / rhosec -c ry(16) = 0.00022355d0 / rhosec !Like ITRF2008 ry(16) = 0.00042037d0 / rhosec -c rz(16) = 0.00006065d0 / rhosec !Like ITRF2008 rz(16) = 0.00031714d0 / rhosec - drx(16) = -0.00001347d0 / rhosec !Like ITRF2008 - dry(16) = 0.00001514d0 / rhosec !Like ITRF2008 - drz(16) = 0.00001973d0 / rhosec !Like ITRF2008 - -C scale(16) = -1.69504d-9 !Differs from ITRF2008, but only insignificantly + drx(16) = -0.00001347d0 / rhosec + dry(16) = 0.00001514d0 / rhosec + drz(16) = 0.00001973d0 / rhosec scale(16) = -0.36891d-9 - dscale(16) = 0.07201d-9 !Differs from ITRF2008, but only insignificantly - refepc(16) = 2010.0d0 !Differs from ITRF2008 + dscale(16) = 0.07201d-9 + refepc(16) = 2010.0d0 C************************************************************************************************************************* -C Parameters computed with the IERS values of ITRF96==>ITRF97 +C Parameters computed with the IERS convention ITRF96 = ITRF97 *** From ITRF94 to NAD 83 tx1(1) = 0.9910d0 ty1(1) = -1.9072d0 - tz1(1) = -.5129d0 + tz1(1) = -0.5129d0 dtx1(1) = 0.d0 dty1(1) = 0.d0 dtz1(1) = 0.d0 @@ -4081,8 +3961,8 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 ry1(1) = 0.46785d-7 rz1(1) = 0.56529d-7 drx1(1) = 0.00258d-7 - dry1(1) = -.03599d-7 - drz1(1) = -.00153d-7 + dry1(1) = -0.03599d-7 + drz1(1) = -0.00153d-7 scale(1) = 0.d0 dscale1(1) = 0.0d0 refepc1(1) = 1997.0d0 @@ -4090,11 +3970,11 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF88 tx1(2) = 0.018d0 ty1(2) = 0.000d0 - tz1(2) = -.092d0 + tz1(2) = -0.092d0 dtx1(2) = 0.0d0 dty1(2) = 0.0d0 dtz1(2) = 0.0d0 - rx1(2) = -.0001d0 / rhosec + rx1(2) = -0.0001d0 / rhosec ry1(2) = 0.0d0 rz1(2) = 0.0d0 drx1(2) = 0.0d0 @@ -4107,7 +3987,7 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF89 tx1(3) = 0.023d0 ty1(3) = 0.036d0 - tz1(3) = -.068d0 + tz1(3) = -0.068d0 dtx1(3) = 0.0d0 dty1(3) = 0.0d0 dtz1(3) = 0.0d0 @@ -4124,7 +4004,7 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF90 tx1(4) = 0.018d0 ty1(4) = 0.012d0 - tz1(4) = -.030d0 + tz1(4) = -0.030d0 dtx1(4) = 0.0d0 dty1(4) = 0.0d0 dtz1(4) = 0.0d0 @@ -4141,7 +4021,7 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF91 tx1(5) = 0.020d0 ty1(5) = 0.016d0 - tz1(5) = -.014d0 + tz1(5) = -0.014d0 dtx1(5) = 0.0d0 dty1(5) = 0.0d0 dtz1(5) = 0.0d0 @@ -4158,7 +4038,7 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 *** From ITRF94 to ITRF92 tx1(6) = 0.008d0 ty1(6) = 0.002d0 - tz1(6) = -.008d0 + tz1(6) = -0.008d0 dtx1(6) = 0.0d0 dty1(6) = 0.0d0 dtz1(6) = 0.0d0 @@ -4168,28 +4048,29 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 drx1(6) = 0.0d0 dry1(6) = 0.0d0 drz1(6) = 0.0d0 - scale1(6) = -.08d-8 + scale1(6) = -0.08d-8 dscale1(6) = 0.0d0 refepc1(6) = 1988.0d0 *** From ITRF94 to ITRF93 - tx1(7) = 0.006d0 - ty1(7) = -.005d0 - tz1(7) = -.015d0 - dtx1(7) = -.0029d0 - dty1(7) = 0.0004d0 - dtz1(7) = 0.0008d0 - rx1(7) = 0.00039d0 / rhosec - ry1(7) = -.00080d0 / rhosec - rz1(7) = 0.00096d0 / rhosec - drx1(7) = .00011d0 / rhosec - dry1(7) = .00019d0 / rhosec - drz1(7) =-.00005d0 / rhosec - scale1(7) = 0.04d-8 + tx1(7) = 0.006d0 + ty1(7) = -0.005d0 + tz1(7) = -0.015d0 + dtx1(7) = -0.0029d0 + dty1(7) = 0.0004d0 + dtz1(7) = 0.0008d0 + rx1(7) = 0.00039d0 / rhosec + ry1(7) = -0.00080d0 / rhosec + rz1(7) = 0.00096d0 / rhosec + drx1(7) = 0.00011d0 / rhosec + dry1(7) = 0.00019d0 / rhosec + drz1(7) = -0.00005d0 / rhosec + scale1(7) = 0.04d-8 dscale1(7) = 0.0d0 refepc1(7) = 1988.0d0 *** From ITRF94 to ITRF96 +*** According to IERS: ITRF97 = ITRF96 = ITRF94 tx1(8) = 0.d0 ty1(8) = 0.d0 tz1(8) = 0.d0 @@ -4206,20 +4087,20 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 dscale1(8) = 0.0d0 refepc1(8) = 1996.0d0 -*** From ITRF94 to ITRF97 (based on IERS adopted values) +*** From ITRF94 to ITRF97 (also IGS97), based on IERS adopted values *** According to IERS: ITRF97 = ITRF96 = ITRF94 - tx1(9) = 0.00000d0 - ty1(9) = 0.00000d0 - tz1(9) = 0.00000d0 - dtx1(9) = 0.00000d0 - dty1(9) = 0.00000d0 - dtz1(9) = 0.00000d0 - rx1(9) = 0.00000000d0 - ry1(9) = 0.00000000d0 - rz1(9) = 0.00000000d0 - drx1(9) = 0.00000000d0 - dry1(9) = 0.00000000d0 - drz1(9) = 0.00000000d0 + tx1(9) = 0.0d0 + ty1(9) = 0.0d0 + tz1(9) = 0.0d0 + dtx1(9) = 0.0d0 + dty1(9) = 0.0d0 + dtz1(9) = 0.0d0 + rx1(9) = 0.0d0 + ry1(9) = 0.0d0 + rz1(9) = 0.0d0 + drx1(9) = 0.0d0 + dry1(9) = 0.0d0 + drz1(9) = 0.0d0 scale1(9) = 0.d0 dscale1(9) = 0.d0 refepc1(9) = 2000.0d0 @@ -4241,40 +4122,40 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 dscale1(10) = 0.0d0 refepc1(10) = 1997.0d0 -*** From ITRF94 to ITRF00 +*** From ITRF94 to ITRF2000 (also IGS00 and IGb00) *** assumes that ITRF94 = ITRF96 and -*** uses IERS values for ITRF96 -> ITRF97 -*** and IERS values for ITRF97 -> ITRF00 - tx1(11) = -.00670d0 - ty1(11) = -.00430d0 - tz1(11) = +.02270d0 - dtx1(11) = 0.00000d0 - dty1(11) = 0.00060d0 - dtz1(11) = 0.00140d0 - rx1(11) = 0.00000000d0 / rhosec - ry1(11) = 0.00000000d0 / rhosec - rz1(11) = +0.00006000d0 / rhosec - drx1(11) = 0.00000000d0 / rhosec - dry1(11) = 0.00000000d0 / rhosec - drz1(11) = +0.00002000d0 / rhosec +*** uses IERS convention ITRF96 = ITRF97 +*** and IERS values for ITRF97 -> ITRF2000 + tx1(11) = -0.00670d0 + ty1(11) = -0.00430d0 + tz1(11) = 0.02270d0 + dtx1(11) = 0.00000d0 + dty1(11) = 0.00060d0 + dtz1(11) = 0.00140d0 + rx1(11) = 0.0d0 / rhosec + ry1(11) = 0.0d0 / rhosec + rz1(11) = 0.00006000d0 / rhosec + drx1(11) = 0.0d0 / rhosec + dry1(11) = 0.0d0 / rhosec + drz1(11) = 0.00002000d0 / rhosec scale1(11) = -1.58000d-9 dscale1(11) = -0.01000d-9 refepc1(11) = 2000.0d0 *** From ITRF94 to PACP00 -*** use PA/ITRF00 rotation rates from Beavan et al., (2002) - tx1(12) = 0.9035d0 +*** use PA/ITRF2000 rotation rates from Beavan et al., (2002) + tx1(12) = 0.9035d0 ty1(12) = -2.0202d0 tz1(12) = -0.5417d0 - dtx1(12) = 0.00000d0 - dty1(12) = 0.00060d0 - dtz1(12) = 0.0014d0 - rx1(12) = 0.027741d0 / rhosec - ry1(12) = 0.013469d0 / rhosec - rz1(12) = 0.002712d0 / rhosec - drx1(12) = -.000384d0 / rhosec - dry1(12) = 0.001007d0 / rhosec - drz1(12) = -.002166d0 / rhosec + dtx1(12) = 0.0d0 + dty1(12) = 0.00060d0 + dtz1(12) = 0.0014d0 + rx1(12) = 0.027741d0 / rhosec + ry1(12) = 0.013469d0 / rhosec + rz1(12) = 0.002712d0 / rhosec + drx1(12) = -0.000384d0 / rhosec + dry1(12) = 0.001007d0 / rhosec + drz1(12) = -0.002166d0 / rhosec scale1(12) = -1.55000d-9 dscale1(12) = -0.010000d-9 refepc1(12) = 1997.0d0 @@ -4287,8 +4168,8 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 dtx1(13) = -0.00000d0 dty1(13) = 0.00060d0 dtz1(13) = 0.00140d0 - rx1(13) = .028971d0 / rhosec - ry1(13) = .01042d0 / rhosec + rx1(13) = 0.028971d0 / rhosec + ry1(13) = 0.01042d0 / rhosec rz1(13) = 0.008928d0 / rhosec drx1(13) = -0.00002d0 / rhosec dry1(13) = 0.000105d0 / rhosec @@ -4297,9 +4178,9 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 dscale1(13) = -0.01000d-9 refepc1(13) = 1997.00d0 -*** From ITRF94 to ITRF2005 -*** assumes that ITRF94 = ITRF96 -*** uses IERS values for ITRF96 -> ITRF97 +*** From ITRF94 to ITRF2005 (also IGS05) +*** assumes that ITRF94 = ITRF96 +*** uses IERS convention ITRF96 = ITRF97 *** uses IERS values for ITRF97 -> ITRF2000 *** uses IERS values for ITRF2000 -> ITRF2005 tx1(14) = -0.00680d0 @@ -4308,70 +4189,61 @@ C Parameters computed with the IERS values of ITRF96==>ITRF97 dtx1(14) = 0.00020d0 dty1(14) = 0.00050d0 dtz1(14) = 0.00320d0 - rx1(14) = 0.00000000d0 / rhosec - ry1(14) = 0.00000000d0 / rhosec - rz1(14) = +0.00006000d0 / rhosec - drx1(14) = 0.00000000d0 / rhosec - dry1(14) = 0.00000000d0 / rhosec - drz1(14) = +0.00002000d0 / rhosec + rx1(14) = 0.0d0 / rhosec + ry1(14) = 0.0d0 / rhosec + rz1(14) = 0.00006000d0 / rhosec + drx1(14) = 0.0d0 / rhosec + dry1(14) = 0.0d0 / rhosec + drz1(14) = 0.00002000d0 / rhosec scale1(14) = -1.98000d-9 dscale1(14) = -0.09000d-9 refepc1(14) = 2000.0d0 -*** From ITRF94 to ITRF2008 (also IGS08 and IGB08) +*** From ITRF94 to ITRF2008 (also IGS08 and IGb08) *** assumes that ITRF94 = ITRF96 -*** uses IERS values for ITRF96 -> ITRF97 +*** uses IERS convention ITRF96 = ITRF97 *** uses IERS values for ITRF97 -> ITRF2000 -*** uses IERS values for ITRF2000-> ITRF2005 -*** uses IERS values for ITRF2005 -> ITRF2008 (and IGS08 and IGB08) +*** uses IERS values for ITRF2000 -> ITRF2005 +*** uses IERS values for ITRF2005 -> ITRF2008 tx1(15) = -0.00480d0 ty1(15) = -0.00260d0 tz1(15) = 0.03320d0 dtx1(15) = -0.00010d0 dty1(15) = 0.00050d0 dtz1(15) = 0.00320d0 - rx1(15) = 0.00000000d0 / rhosec - ry1(15) = 0.00000000d0 / rhosec - rz1(15) = +0.00006000d0 / rhosec - drx1(15) = -0.00000000d0 / rhosec - dry1(15) = 0.00000000d0 / rhosec - drz1(15) = +0.00002000d0 / rhosec + rx1(15) = 0.0d0 / rhosec + ry1(15) = 0.0d0 / rhosec + rz1(15) = 0.00006000d0 / rhosec + drx1(15) = -0.0d0 / rhosec + dry1(15) = 0.0d0 / rhosec + drz1(15) = 0.00002000d0 / rhosec scale1(15) = -2.92d-9 dscale1(15) = -0.09d-9 refepc1(15) = 2000.0d0 -*** From ITRF94 to ITRF2014 +*** From ITRF94 to ITRF2014 (also IGS14 and IGb14) *** assumes that ITRF94 = ITRF96 -*** uses IERS values for ITRF96 -> ITRF97 +*** uses IERS convention ITRF96 = ITRF97 *** uses IERS values for ITRF97 -> ITRF2000 -*** uses IERS values for ITRF2000-> ITRF2005 -*** uses IERS values for ITRF2005 -> ITRF2008 (and IGS08 and IGB08) +*** uses IERS values for ITRF2000 -> ITRF2005 +*** uses IERS values for ITRF2005 -> ITRF2008 *** uses IERS values for ITRF2008 -> ITRF2014 -c tx1(16) = -0.00640d0 !Differ from ITRF2008 tx1(16) = -0.00740d0 -c ty1(16) = -0.00450d0 !Differ from ITRF2008 ty1(16) = 0.00050d0 -c tz1(16) = 0.03080d0 !Differ from ITRF2008 tz1(16) = 0.06280d0 - - dtx1(16) = -0.00010d0 !Like ITRF2008 - dty1(16) = 0.00050d0 !Like ITRF2008 - dtz1(16) = 0.00330d0 !Differ from ITRF2008 - - rx1(16) = 0.00000000d0 / rhosec !Like ITRF2008 - ry1(16) = 0.00000000d0 / rhosec !Like ITRF2008 -c rz1(16) = +0.00006000d0 / rhosec !Like ITRF2008 - rz1(16) = +0.00026000d0 / rhosec - drx1(16) = -0.00000000d0 / rhosec !Like ITRF2008 - dry1(16) = 0.00000000d0 / rhosec !Like ITRF2008 - drz1(16) = +0.00002000d0 / rhosec !Like ITRF2008 - -c scale1(16) = -2.90d-9 !Differ from ITRF2008, but only insignificantly - scale1(16) = -3.80d-9 - dscale1(16) = -0.12d-9 !Differ from ITRF2008, but only insignificantly - - refepc1(16) = 2010.0d0 !Differ from ITRF2008 + dtx1(16) = -0.00010d0 + dty1(16) = 0.00050d0 + dtz1(16) = 0.00330d0 + rx1(16) = 0.0d0 / rhosec + ry1(16) = 0.0d0 / rhosec + rz1(16) = 0.00026000d0 / rhosec + drx1(16) = 0.0d0 / rhosec + dry1(16) = 0.0d0 / rhosec + drz1(16) = 0.00002000d0 / rhosec + scale1(16) = -3.80d-9 + dscale1(16) = -0.12d-9 + refepc1(16) = 2010.0d0 return end @@ -4656,9 +4528,9 @@ c nframe(10)= 'PNEOS_90 or NEOS_90 ' iframe(22)= 14 nframe(22)= 'ITRF2005 or IGS05 ' iframe(23)= 15 - nframe(23)= 'ITRF2008 or IGS08/IGB08 ' + nframe(23)= 'ITRF2008 or IGS08/IGb08 ' iframe(24)= 16 - nframe(24)= 'ITRF2014 or IGS14 ' + nframe(24)= 'ITRF2014 or IGS14/IGb14 ' write(luout, 100) 100 format( @@ -4674,10 +4546,10 @@ c 1' 4...WGS_72 '/ 1' 9...WGS_84(G1674) (ITRF2008 used) 19...ITRF96 '/ c 1' 10...PNEOS_90 or NEOS_90 (ITRF90 used) 20...ITRF97 or IGS97'/ 1' 10...WGS_84(G1762) (IGb08 used) 20...ITRF97 or IGS97'/ - 1' 11...SIO/MIT_92 (ITRF91 used) 21...ITRF2000 or IGS00/IGb00'/ - 1' 12...ITRF88 22...ITRF2005 or IGS05 '/ - 1' 13...ITRF89 23...ITRF2008 or IGS08/IGb08'/ - 1' 14...ITRF90 or (PNEOS90/NEOS90) 24...ITRF2014 or IGS14 '/ ) + 1' 11...SIO/MIT_92 (ITRF91 used) 21...ITRF2000 or IGS00/IGb00'/ + 1' 12...ITRF88 22...ITRF2005 or IGS05 '/ + 1' 13...ITRF89 23...ITRF2008 or IGS08/IGb08'/ + 1' 14...ITRF90 or (PNEOS90/NEOS90) 24...ITRF2014 or IGS14/IGb14'/) c 1' 14...ITRF90 '/ ) read (luin, *,err=50,iostat=ios) iopt @@ -4712,17 +4584,17 @@ c 1' 14...ITRF90 '/ ) CHARACTER OPT*1,ANSWER*1,BBTYPE*1,VOPT*1 CHARACTER LATDIR*1, LONDIR*1, LATDR*1, LONDR*1 character frame1*24, frame2*24 - character HTDP_version*10 + character HTDP_version*8 LOGICAL TEST LOGICAL Is_iopt_NAD83 COMMON /CONST/ A,F,E2,EPS,AF,PI,TWOPI,RHOSEC COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6 - COMMON /CAUTION/MONTH2,IDAY2,IYEAR2,DATE2,frame1 + COMMON /CAUTION/ MONTH2,IDAY2,IYEAR2,DATE2,frame1 WRITE(LUOUT,20) 20 FORMAT(' ********************************************'/ 1 ' Please enter the time to which the updated'/ - 1 ' positions and/or observations are to correspond.') + 1 ' positions and/or observations are to correspond:') 15 CALL GETMDY(MONTH2,IDAY2,IYEAR2,DATE2,MIN2,TEST) IF(TEST) then write(luout,*) ' Do you wish to re-enter the time? (y/n)' @@ -4735,8 +4607,7 @@ c 1' 14...ITRF90 '/ ) ** Choosing reference frame for positions 35 WRITE(LUOUT,30) 30 FORMAT(' **************************************************'/ - 1 ' Select the reference frame to be used for specifying'/ - 2 ' positions. '/) + 1 ' Specify the reference frame of the input positions:'/) call MENU1( iopt, frame1) IF(IOPT .LT. 1 .OR. IOPT .GT. numref) THEN WRITE(LUOUT,40) @@ -4757,15 +4628,15 @@ c 1' 14...ITRF90 '/ ) 8 ' stations'/ 9 ' and the values for blue book', 9 ' observations.'/ - 9 ' 5...Update positions contained in batch file '/ - 9 ' of delimited records of the form: '/ - 9 ' LAT,LON,EHT,"TEXT" ' / - 9 ' LAT = latitude in degrees (positive north/DBL PREC)'/ - 9 ' LON = longitude in degrees (positive west/DBL PREC)'/ - 9 ' EHT = ellipsoid height in meters (DBL PREC)'/ - 9 ' TEXT = Descriptive text (CHARACTER*24) '/ - 9 ' Example: '/ - 9 ' 40.731671553,112.212671753,34.241,"SALT AIR" '/) + 9 ' 5...Update positions for multiple points contained'/ + 9 ' in a file in Lat-Lon-Ht format:'/ + 9 ' LAT,LON,EHT,TEXT'/ + 9 ' LAT = latitude in degrees (positive north)'/ + 9 ' LON = longitude in degrees (positive west)'/ + 9 ' EHT = ellipsoid height in meters'/ + 9 ' TEXT = Descriptive text (up to 24 characters)'/ + 9 ' Example:'/ + 9 ' 40.731671553,112.212671753,34.241,SALT AIR'/) READ(LUIN,'(A1)',err=501,iostat=ios) OPT if (ios /= 0) goto 501 IF(OPT .eq. '0') THEN @@ -4774,7 +4645,7 @@ c 1' 14...ITRF90 '/ ) WRITE(LUOUT,1020) 1020 FORMAT(' ************************************'/ 1 ' Enter the time', - 2 ' to which the input positions will correspond. ') + 2 ' to which the input positions will correspond: ') 1025 CALL GETMDY(MONTH1,IDAY1,IYEAR1,DATE1,MIN1,TEST) IF(TEST) then write(luout,*) ' Do you wish to re-enter the time? (y/n)' @@ -4902,7 +4773,7 @@ c 1' 14...ITRF90 '/ ) WRITE(LUOUT,122) 122 FORMAT(' *********************************************'/ 1 ' Enter the time', - 1 ' to which the input positions correspond.'/) + 1 ' to which the input positions correspond:'/) 123 CALL GETMDY(MONTH1,IDAY1,IYEAR1,DATE1,MIN1,TEST) IF(TEST) then write(luout,*) ' Do you wish to re-enter the time? (y/n)' @@ -4927,7 +4798,7 @@ c 1' 14...ITRF90 '/ ) write (i2,127) HTDP_version 127 format (' ***CAUTION: This file was processed using HTDP', - & ' version ',a10, '***') + & ' version ', a8, '***') write (i2,128) frame1 128 format (' ***CAUTION: Coordinates in this file are in ', & a24, '***') @@ -4977,8 +4848,8 @@ c ENDIF 634 WRITE(LUOUT, 632) 632 FORMAT(/' ***************************'/ - * ' To what reference frame should the GPS'/ - * ' vectors be transformed?'/ + * ' Specify the reference frame for the updated' + * ' G-File vectors:'// * ' -1...Do not transform GPS vectors.') CALL MENU1(kopt, frame2) IF(KOPT .LT. -1 .OR. KOPT .GT. numref) THEN @@ -4990,7 +4861,7 @@ c ENDIF 1 ' have been transformed to ', a24, ' ***') write( I2, 641) HTDP_version 641 format(' ***CAUTION: Observations were transformed using' - 1 ,' HTDP version ', a10, ' ***') + 1 ,' HTDP version ', a8, ' ***') ENDIF IF(BBTYPE .EQ. '1') THEN @@ -5304,7 +5175,7 @@ C HT(ISN) = HT(ISN) + GH END ****************************************************************** - SUBROUTINE UPBB4(MIN1,MIN2,OPT,IOPT) + SUBROUTINE UPBB4(MIN1,MIN2,OPT) *** Update blue book @@ -5360,7 +5231,7 @@ C*** classical observation records contain only a 2-digit year 1 TYPE .EQ. '*85*' .OR. 1 TYPE .EQ. '*86*' .OR. 1 TYPE .EQ. '*90*') THEN - CONTINUE + CONTINUE ELSEIF (TYPE .EQ. '*12*') THEN IF( CARD(11:14) .EQ. ' ' .OR. @@ -5636,7 +5507,7 @@ c WRITE(LUOUT,500) TYPE END *************************************************************** - SUBROUTINE UPBB5(MIN1,MIN2,OPT,IOPT) + SUBROUTINE UPBB5(MIN1,MIN2,OPT) *** Update 5-digit blue book @@ -6234,10 +6105,10 @@ C DECYR2 = DBLE(IYEAR2) + DBLE(MINO2 - MIN00)/525600.D0 if (ios /= 0) goto 301 110 FORMAT(1X,I4,I2,I2,4X,I4,I2,I2,30X,I2) C CALL TOTIME(IYEAR1,MONTH1,IDAY1,MINO1) -C CALL TOTIME(IYEAR1, 1, 1, MIN00) -C DECYR1 = DBLE(IYEAR1) + DBLE(MINO1 - MIN00)/525600.D0 +C CALL TOTIME(IYEAR1, 1, 1, MIN00) +C DECYR1 = DBLE(IYEAR1) + DBLE(MINO1 - MIN00)/525600.D0 C CALL TOTIME(IYEAR2,MONTH2,IDAY2,MINO2) -C CALL TOTIME(IYEAR2, 1, 1, MIN00) +C CALL TOTIME(IYEAR2, 1, 1, MIN00) CALL IYMDMJ(IYEAR1,MONTH1,IDAY1,MJD1) MINO1 = MJD1 * 24 * 60 CALL IYMDMJ(IYEAR1, 1, 1, MJD0) @@ -6246,7 +6117,7 @@ C CALL TOTIME(IYEAR2, 1, 1, MIN00) MINO2 = MJD2 * 24 * 60 CALL IYMDMJ(IYEAR2, 1, 1, MJD0) DECYR2 = DBLE(IYEAR2) + DBLE(MJD2 - MJD0)/365.D0 -C DECYR2 = DBLE(IYEAR2) + DBLE(MINO2 - MIN00)/525600.D0 +C DECYR2 = DBLE(IYEAR2) + DBLE(MINO2 - MIN00)/525600.D0 MINO = (MINO1 + MINO2) / 2 DECYR = (DECYR1 + DECYR2) / 2.D0 CALL RFCON(IBBREF, JREF) @@ -6385,7 +6256,7 @@ C DECYR2 = DBLE(IYEAR2) + DBLE(MINO2 - MIN00)/525600.D0 IMPLICIT INTEGER*4 (I-N) parameter ( numref = 16 ) - COMMON /REFCON/ IRFCON(36), JRFCON(numref) + COMMON /REFCON/ IRFCON(37), JRFCON(numref) *** From blue book identifier to HTDP indentifier *** WGS 72 Precise @@ -6393,7 +6264,7 @@ c IRFCON(1) = 10 IRFCON(1) = 1 C HTDP no longer supports WGS 72. Hence, if a BlueBook C file contains WGS 72 coordinates, HTDP treats these -C coordinates as if they were NAD 83(2011)coordinates. +C coordinates as if they were NAD 83(2011) coordinates. *** WGS 84 (orig) Precise (set equal to NAD 83) IRFCON(2) = 1 @@ -6456,7 +6327,7 @@ C coordinates as if they were NAD 83(2011)coordinates. *** IGS97 IRFCON(20) = 9 -*** ITRF00 +*** ITRF2000 IRFCON(21) = 11 *** IGS00 @@ -6477,7 +6348,7 @@ C coordinates as if they were NAD 83(2011)coordinates. *** IGS08 IRFCON(27) = 15 -*** IGB08 +*** IGb08 IRFCON(28) = 15 *** ITRF2008 @@ -6492,7 +6363,7 @@ C coordinates as if they were NAD 83(2011)coordinates. *** ITRF2014 IRFCON(32) = 16 -*** IGB14 +*** IGS14 IRFCON(33) = 16 *** NAD83 (2011/2007/CORS96/FBN/HARN) @@ -6504,6 +6375,10 @@ C coordinates as if they were NAD 83(2011)coordinates. *** NAD83 (MA11) IRFCON(36) = 13 +*** IGb14 + IRFCON(37) = 16 + + *** From HTDP identifier to blue book identifier *** NAD 83 (set equal to WGS 84 (transit)) c JRFCON(1) = 2 @@ -6536,7 +6411,7 @@ c JRFCON(1) = 2 *** WGS 72 JRFCON(10) = 1 -*** ITRF00 +*** ITRF2000 JRFCON(11) = 21 *** NAD 83 (PACP00) or NAD 83 (PA11) @@ -6561,12 +6436,15 @@ C coordinates in the Bluebook context *** ITRF2008 or IGS08 JRFCON(15) = 27 -*** IGB08 +*** IGb08 JRFCON(15) = 28 *** ITRF2014 or IGS14 JRFCON(16) = 33 +*** IGb14 + JRFCON(16) = 37 + RETURN END *************************************************** @@ -6578,9 +6456,9 @@ C coordinates in the Bluebook context IMPLICIT INTEGER*4 (I-N) parameter ( numref = 16 ) - COMMON /REFCON/ IRFCON(36), JRFCON(numref) + COMMON /REFCON/ IRFCON(37), JRFCON(numref) - IF (1 .LE. IBBREF .AND. IBBREF .LE. 36) THEN + IF (1 .LE. IBBREF .AND. IBBREF .LE. 37) THEN JREF = IRFCON(IBBREF) ELSE WRITE(6, 10) IBBREF @@ -6601,7 +6479,7 @@ C coordinates in the Bluebook context IMPLICIT INTEGER*4 (I-N) parameter ( numref = 16 ) - COMMON /REFCON/ IRFCON(36), JRFCON(numref) + COMMON /REFCON/ IRFCON(37), JRFCON(numref) IF (JREF .EQ. 0) THEN I = 1 @@ -7250,7 +7128,7 @@ c write (*,*) "From PREDV ",VN, VE, VU 1 16X, ' INPUT VELOCITIES OUTPUT VELOCITIES'/) 130 write( luout, 140) - 140 format( /'**********************************'/ + 140 format(' ************************************************'/ 1 ' Velocities will be transformed at each specified point.'/ 1 ' Please indicate how you wish to input points.'/ 1 ' 0...No more points. Return to main menu.'/ @@ -7491,13 +7369,13 @@ C The parameters in common block tranpa1 are computed using the IERS values of if (iprint .eq. 1) then write( luout, 100) vn1, ve1, vu1, vx1, vy1, vz1 - 100 format( ' ****************************************'/ + 100 format(' ************************************************'/ 1 ' New northward velocity = ', f8.2, ' mm/yr' / 1 ' New eastward velocity = ', f8.2, ' mm/yr'/ 1 ' New upward velocity = ', f8.2, ' mm/yr'/ 1 ' New x velocity = ', f8.2, ' mm/yr'/ 1 ' New y velocity = ', f8.2, ' mm/yr'/ - 1 ' New z velocity = ', f8.2, ' mm/yr'/) + 1 ' New z velocity = ', f8.2, ' mm/yr') endif write( i2, 200) name24, xlat, xlon, @@ -7519,12 +7397,12 @@ C The parameters in common block tranpa1 are computed using the IERS values of IMPLICIT DOUBLE PRECISION (A-H, O-Z) IMPLICIT INTEGER*4 (I-N) - character HTDP_version*10 + character HTDP_version*8 COMMON /FILES/ LUIN, LUOUT, I1, I2, I3, I4, I5, I6 COMMON /VERSION/ HTDP_version WRITE(I2, 10) HTDP_version - 10 FORMAT(' HTDP (VERSION ',a,') OUTPUT' / ) + 10 FORMAT(' HTDP OUTPUT, VERSION ',a / ) RETURN END ********************************************* @@ -8068,7 +7946,7 @@ C------------------------------------------------------------------------------- implicit none - integer*4 i,j,length + integer*4 i,length real*8 x,y,z character name*80,record*120,record1*120,chars*120 character xxxx*80,yyyy*80,zzzz*80 @@ -8183,7 +8061,7 @@ C------------------------------------------------------------------------------- implicit none - integer*4 i,j,length + integer*4 i,length real*8 x,y,z,Vx,Vy,Vz character name*80,record*120,record1*120,chars*120 character xxxx*80,yyyy*80,zzzz*80 @@ -8381,10 +8259,10 @@ C------------------------------------------------------------------------------- implicit none - integer*4 i,j,length - real*8 x,y,z + integer*4 i,length + real*8 x,y character name*24,record*120,record1*120,chars*120 - character xxxx*80,yyyy*80,zzzz*80 + character xxxx*80,yyyy*80 record1 = trim(adjustl(record)) length = len_trim(record1) @@ -8468,7 +8346,7 @@ C------------------------------------------------------------------------------- implicit none - integer*4 i,j,length + integer*4 i,length real*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 @@ -8634,490 +8512,8 @@ C Done return end C----------------------------------------------------------------------------------- - logical function am_I_in_or_near_CONUS (fi,la) - - implicit none - - integer*4 nrows,ncols,xcell,ycell,rec_num - integer*2 sea,code,code_UL,code_UR,code_LL,code_LR - real*8 fi,la,fi_max,la_min,dlamda,dfi,fi_min - real*8 fi1,la1,fi2,la2,la_max - - character grid*80 -c logical am_I_in_or_near_CONUS,CONUS_on_land - logical CONUS_on_land - -C Constants - - parameter (fi_max = 50.d0, - & fi_min = 23.5d0, - & la_min = 235.d0, - & la_max = 295.d0, - & sea = 0) - -C Initialize this function - - am_I_in_or_near_CONUS = .false. - -C If clearly outside Conus - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_CONUS = .true. - return - endif - - return - end -C----------------------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_AK(fi,la) - -C Given your latitude and longitude, this little routine tells you if you -C are in Alaska (AK) or within 10 km of its boundaries. - -C "AK_direct_cell_by_cell_1min" is an unformatted direct access file of nowrs*ncols records -C each record is of length 2 bytes (integer*2), which contains a geographic code defined as follows: - -C 1) The AK code is 278 -C 2) The codes for Canada are 265 to 276 -C 4) The code for sea is 0 - -C Geographic conditions to be in or near AK: - -C 1) The point is on land in AK - -C a) code is 278 - -C 2) In the Canada near the AK: - -C a) code = 265 to 276 -C c) Point is no more than 10km from AK - -C 3) In the sea near AK: - -C a) code = 0 -C c) Point is no more than 10km from AK - - implicit none - - integer*4 nrows,ncols,xcell,ycell,rec_num - integer*2 sea,code,code_UL,code_UR,code_LL,code_LR - real*8 fi,la,fi_max,dlamda,dfi,la_max - real*8 fi1,la1,fi2,la2,la_min,fi_min - - character grid*27,path*46 -c logical am_I_in_or_near_AK,AK_on_land - logical AK_on_land - -C Constants - - parameter (nrows = 1291, - & ncols = 3511, - & dfi = 1.d0/60.d0, - & dlamda = 1.d0/60.d0, - & fi_max = 71.5d0, - & fi_min = 50.d0, - & la_min = 172.d0, - & la_max = 230.5d0, - & sea = 0) - -C Initialize this function - - am_I_in_or_near_AK = .false. - -C When it is clearly not AK - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_AK = .true. - return - endif - -C On 1/22/2016 I disabled the following code following Giovanni's suggestion -C and for simplicity and portability - - return - end -C---------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_HI(fi,la) - -C Given your latitude and longitude, this little routine tells you if you -C are in the Hawaian Islands (Hawaii, Kahoolawe,Kauai,Lanai,Maui,Molokai, Nihau and Oahu) -C or within 5 km of their coastlines. - -C "HI_Hawaii.gmt", "HI_Kahoolawe.gmt", "HI_Kauai.gmt", "HI_Lanai.gmt", "HI_Maui.gmt", "HI_Molokai.gmt", -C "HI_Nihau.gmt" and "HI_Oahu.gmt" are ASCII files that contain coastline polygons for the islands. -C Each polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - implicit none - - integer*4 i,NPC,n_Hawaii,n_Kahoolawe,n_Kauai,n_Lanai,n_Maui - integer*4 n_Molokai,n_Nihau,n_Oahu - - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min -c logical am_I_in_or_near_HI,Hawaii,Kahoolawe,Kauai,Lanai - logical Hawaii,Kahoolawe,Kauai,Lanai - logical Maui,Molokai,Nihau,Oahu - logical UL,UR,LL,LR - -C Constants - - parameter (n_Hawaii = 329, - & n_Kahoolawe = 143, - & n_Kauai = 107, - & n_Lanai = 52, - & n_Maui = 205, - & n_Molokai = 114, - & n_Nihau = 68, - & n_Oahu = 123, - & fi_max = 24.d0, - & fi_min = 18.d0, - & la_min = 199.d0, - & la_max = 207.d0) - - real*8 X1(n_Hawaii) ,Y1(n_Hawaii) !I am not sure the cheap compilers in subversion can handle dynamic memory - real*8 X2(n_Kahoolawe) ,Y2(n_Kahoolawe) - real*8 X3(n_Kauai) ,Y3(n_Kauai) - real*8 X4(n_Lanai) ,Y4(n_Lanai) - real*8 X5(n_Maui) ,Y5(n_Maui) - real*8 X6(n_Molokai) ,Y6(n_Molokai) - real*8 X7(n_Nihau) ,Y7(n_Nihau) - real*8 X8(n_Oahu) ,Y8(n_Oahu) - -C Initialize this function - - am_I_in_or_near_HI = .false. - -C When it is clearly not HI - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_HI = .true. - return - endif - -C On 1/22/2016 I disabled the following code following Giovanni's suggestion -C and to simplify. - - return - end -C------------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_PR(fi,la) - -C Given your latitude and longitude, this little routine tells you if you -C are in Puerto Rico Islands (the main Island, Culebra, Desecheo, Mona or Viequez) -C or within 10 km of their coastlines. - -C "PR_PuertoRico.gmt", "PR_Culebra", "PR_Mona.gmt", "PR_Viequez.gmt", and PR_"Desecheo.gmt" -C are ASCII files that contain coastline polygons for each one of the 5 PR islands. -C Each polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - - implicit none - - integer*4 i,NPC,n_PR,n_Culebra,n_Mona,n_Desecheo,n_Viequez - - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min -c logical am_I_in_or_near_PR,PR - logical PR - logical Culebra,Mona,Desecheo,Viequez - logical UL,UR,LL,LR - -C Constants - - parameter (n_PR = 228, - & n_Culebra = 185, - & n_Mona = 52, - & n_Desecheo = 34, - & n_Viequez = 309, - & fi_max = 19.d0, - & fi_min = 17.d0, - & la_max = 295.d0, - & la_min = 292.d0) - - real*8 X1(n_PR) ,Y1(n_PR) !I am not sure the cheap compilers in subversion can handle dynamic memory - real*8 X2(n_Culebra) ,Y2(n_Culebra) - real*8 X3(n_Mona) ,Y3(n_Mona) - real*8 X4(n_Desecheo),Y4(n_Desecheo) - real*8 X5(n_Viequez) ,Y5(n_Viequez) - -C Initialize this function - - am_I_in_or_near_PR = .false. - -C When it is clearly not PR - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_PR = .true. - return - endif - -C On 1/22/2016 I disabled the following code following Giovanni's suggestion -C and for simplicity. - - return - end -C------------------------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_VQ(fi,la) - -C Given your latitude and longitude, this little routine tells you if you are in -C the Virgin Islands (St. Croix, St John and St Thomas) or within 10 km of their coastlines. - -C "VQ_StCroix.gmt", "VQ_StJohn.gmt" and "VQ_StThomas.gmt" are ASCII files that -C contain coastline polygons for each one of the 3 VQ islands. -C Each polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - - implicit none - - integer*4 i,NPC,n_VQ,n_StCroix,n_StJohn,n_StThomas - - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min -c logical am_I_in_or_near_VQ - logical StCroix,StJohn,StThomas - logical UL,UR,LL,LR - -C Constants - - parameter (n_StCroix = 142, - & n_StJohn = 356, - & n_StThomas = 304, - & fi_max = 18.4d0, - & fi_min = 15.5d0, - & la_max = 295.6d0, - & la_min = 295.d0) - - real*8 X1(n_StCroix) ,Y1(n_StCroix) !I am not sure the cheap compilers in subversion can handle dynamic memory - real*8 X2(n_StJohn ) ,Y2(n_StJohn ) - real*8 X3(n_StThomas),Y3(n_StThomas) - -C Initialize this function - - am_I_in_or_near_VQ = .false. - -C When it is clearly not VQ - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_VQ = .true. - return - endif - -C On 1/22/2016 the following code was disabled following Giovanni's suggestion -C and for simplicity - - return - end -C------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_CQ (fi,la) - -C Given your latitude and longitude, this little routine tells you if you are in -C the North Mariana Islands (Saipan, Tinian and Rota) or within 10 km of their coastlines. - -C "CQ_Rota.gmt", "CQ_Saipan.gmt" and "CQ_Tinian.gmt" are ASCII files that -C contain coastline polygons for each one of the 3 CQ islands. -C Each polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - - implicit none - - integer*4 i,NPC,n_CQ,n_Saipan ,n_Tinian,n_Rota - - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min -c logical am_I_in_or_near_CQ - logical Saipan,Tinian,Rota - logical UL,UR,LL,LR - -C Constants - - parameter (n_Saipan = 122, - & n_Tinian = 70, - & n_Rota = 80, - & fi_max = 15.6d0, - & fi_min = 13.8d0, - & la_max = 146.d0, - & la_min = 144.8d0) - - real*8 X1(n_Saipan ),Y1(n_Saipan ) !I am not sure the cheap compilers in subversion can handle dynamic memory - real*8 X2(n_Tinian ),Y2(n_Tinian ) - real*8 X3(n_Rota ),Y3(n_Rota ) - -C Initialize this function - - am_I_in_or_near_CQ = .false. - -C When it is clearly not CQ - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_CQ = .true. - return - endif - -C On 1/22/2016 the following code was disabled following Giovanni's suggestion -C and for simplicity - - return - end -C------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_AS(fi,la) - -C Given your latitude and longitude, this little routine tells you if you are in -C the American Samoa Islands (Ofu, Rose, Tau and Tutila) or within 10 km of their coastlines. - -C "AS_Ofu.gmt", "AS_Rose.gmt" and " AS_Tau.gmt" and "AS_Tutila.gmt" are ASCII files that -C contain coastline polygons for each one of the 4 AS islands. -C Each polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - - implicit none - - integer*4 i,NPC,n_AS,n_Ofu,n_Rose,n_Tau,n_Tutila - - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min -c logical am_I_in_or_near_AS - logical Ofu,Rose,Tau,Tutila - logical UL,UR,LL,LR - -C Constants - - parameter (n_Ofu = 64, - & n_Rose = 8, - & n_Tau = 76, - & n_Tutila = 239, - & fi_max = -13.7d0, - & fi_min = -14.7d0, - & la_max = 190.8d0, - & la_min = 189.d0) - - real*8 X1(n_Ofu ) ,Y1(n_Ofu ) !not sure the cheap compilers in subversion can handle dynamic memory - real*8 X2(n_Rose) ,Y2(n_Rose) - real*8 X3(n_Tau ) ,Y3(n_Tau ) - real*8 X4(n_Tutila),Y4(n_Tutila) - -C Initialize this function - - am_I_in_or_near_AS = .false. - -C When it is clearly not AS - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_AS = .true. - return - endif - -C On 1/22/2016 the following code was disabled following Giovanni's suggestion -C and for simplicity - - return - end -C------------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_Guam(fi,la) - -C Given your latitude and longitude, this little routine tells you if you are in -C the Guam Island or within 5 km of its coastlines. - -C "GQ_Guam.gmt", is an ASCII file that contains coastline polygons Guam. -C This polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - - implicit none - - integer*4 i,NPC,n_GQ - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min -c logical am_I_in_or_near_Guam - logical GQ - logical UL,UR,LL,LR - -C Constants - - parameter (n_GQ = 300, - & fi_max = 13.7d0, - & fi_min = 13.d0, - & la_max = 145.d0, - & la_min = 144.5d0) - - real*8 X1(n_GQ) ,Y1(n_GQ) !I am not sure the cheap compilers in subversion can handle dynamic memory - -C Initialize this function - - am_I_in_or_near_Guam = .false. - -C When it is clearly not Guam - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_Guam = .true. - return - endif - -C On 1/22/2016 the following code was disabled following Giovanni's suggestion -C and for simplicity - - return - end -C---------------------------------------------------------------------------------------------- - logical function am_I_in_or_near_KW(fi,la) - -C Given your latitude and longitude, this little routine tells you if you are in -C the a little square around Kwajalein of the Marshal islands - -C "Kwajalein.gmt", is an ASCII file that contains a rectangle around Kwajalein. -C This polygon is closed, starts and ends at the same point. -C File format: (East) longitude and latitude in decimal degrees. - - - implicit none - - integer*4 i,NPC,n_KW - - real*8 fi,la,fi1,la1,fi2,la2,fi_max,fi_min,la_max,la_min - -c logical am_I_in_or_near_KW - logical KW - logical UL,UR,LL,LR - -C Constants - - parameter (n_KW = 5, - & fi_max = 10.5d0, - & fi_min = 5.d0 , - & la_max = 174.d0, - & la_min = 165.7d0) - - real*8 X1(n_KW) ,Y1(n_KW) !I am not sure the cheap compilers in subversion can handle dynamic memory - -C Initialize this function - - am_I_in_or_near_KW = .false. - -C When it is clearly not KW - - if (fi>fi_max.or.fi<fi_min.or.la>la_max.or.la<la_min) then - return - else - am_I_in_or_near_KW = .true. - return - endif - -C On 1/22/2016 the following code was disabled following Giovanni's suggestion -C and for simplicity - - return - end +C On 5/10/2020 removed all code and comments entered on 7/22/2015 that limited use +C of HTDP to US and its territories. That functionality was disabled on 5/11/2017. +C The removed code still exists in HTDP v3.2.8 and 3.2.7, commented out to disable. +C But functions remained in those previous versions, causing numerous warnings when +C program compiled. diff --git a/infile2.txt b/infile2.txt deleted file mode 100644 index 74af99f..0000000 --- a/infile2.txt +++ /dev/null @@ -1,2 +0,0 @@ -40.23,120.42,0.0,SALT AIR -35.00,121.00,3.2,test diff --git a/input/.svn/all-wcprops b/input/.svn/all-wcprops new file mode 100644 index 0000000..a14a502 --- /dev/null +++ b/input/.svn/all-wcprops @@ -0,0 +1,35 @@ +K 25 +svn:wc:ra_dav:version-url +V 96 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input +END +test_Mariana.txt +K 25 +svn:wc:ra_dav:version-url +V 113 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input/test_Mariana.txt +END +wa.bfile.txt +K 25 +svn:wc:ra_dav:version-url +V 109 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input/wa.bfile.txt +END +wa.gfile.txt +K 25 +svn:wc:ra_dav:version-url +V 109 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input/wa.gfile.txt +END +test_Pacific.txt +K 25 +svn:wc:ra_dav:version-url +V 113 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input/test_Pacific.txt +END +test_NAmerica.txt +K 25 +svn:wc:ra_dav:version-url +V 114 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input/test_NAmerica.txt +END diff --git a/input/.svn/entries b/input/.svn/entries new file mode 100644 index 0000000..88fe192 --- /dev/null +++ b/input/.svn/entries @@ -0,0 +1,198 @@ +10 + +dir +121793 +https://source.ngs.noaa.gov/svn/repos/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/input +https://source.ngs.noaa.gov/svn/repos + + + +2020-05-18T00:13:26.379714Z +115868 +michael.dennis + + + + + + + + + + + + + + +53ba9092-9722-c17e-94b9-e0e50a2f61ab + +test_Mariana.txt +file + + + + +2021-04-09T01:01:26.242847Z +9a066d287a6d16c8de0bf972f21836cb +2020-05-18T00:13:26.379714Z +115868 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +79 + +wa.bfile.txt +file + + + + +2021-04-09T01:01:26.363843Z +0c483b04973db4dfac02b76754f27048 +2020-05-18T00:13:26.379714Z +115868 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +6916 + +wa.gfile.txt +file + + + + +2021-04-09T01:01:26.384845Z +6175fca73065b218c55eb17466db284d +2020-05-18T00:13:26.379714Z +115868 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +1527 + +test_Pacific.txt +file + + + + +2021-04-09T01:01:26.387838Z +feb966a03604612993b1a251a9596935 +2020-05-18T00:13:26.379714Z +115868 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +91 + +test_NAmerica.txt +file + + + + +2021-04-09T01:01:26.390834Z +d103923b908cfd55d8600ee30559fb82 +2020-05-18T00:13:26.379714Z +115868 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +173 + diff --git a/input/.svn/text-base/test_Mariana.txt.svn-base b/input/.svn/text-base/test_Mariana.txt.svn-base new file mode 100644 index 0000000..9e60389 --- /dev/null +++ b/input/.svn/text-base/test_Mariana.txt.svn-base @@ -0,0 +1,3 @@ +13.4,215.3,0,Guam +15.2,214.25,0,Saipan +0,0,0,Unmodeled location (near Africa)
\ No newline at end of file diff --git a/input/.svn/text-base/test_NAmerica.txt.svn-base b/input/.svn/text-base/test_NAmerica.txt.svn-base new file mode 100644 index 0000000..0673e9a --- /dev/null +++ b/input/.svn/text-base/test_NAmerica.txt.svn-base @@ -0,0 +1,8 @@ +40,100,0,Kansas +37,122,0,California +48,124,0,Washington +45,69,0,Maine +28,81,0,Florida +18.2,66.5,0,Puerto Rico +65,152,0,Alaska +0,0,0,Unmodeled location (near Africa) diff --git a/input/.svn/text-base/test_Pacific.txt.svn-base b/input/.svn/text-base/test_Pacific.txt.svn-base new file mode 100644 index 0000000..4034b59 --- /dev/null +++ b/input/.svn/text-base/test_Pacific.txt.svn-base @@ -0,0 +1,3 @@ +19.5,155.5,0,Hawaii +-14.3,170.7,0,American Samoa +0,0,0,Unmodeled location (near Africa) diff --git a/wa.bfile.txt b/input/.svn/text-base/wa.bfile.txt.svn-base index e178229..e178229 100755..100644 --- a/wa.bfile.txt +++ b/input/.svn/text-base/wa.bfile.txt.svn-base diff --git a/input/.svn/text-base/wa.gfile.txt.svn-base b/input/.svn/text-base/wa.gfile.txt.svn-base new file mode 100644 index 0000000..a2ce076 --- /dev/null +++ b/input/.svn/text-base/wa.gfile.txt.svn-base @@ -0,0 +1,25 @@ +AAA201110 5201110 5SKAGIT AIRPORT PACS SACS +B201110 517 9201110 522 0 1 pages v1109.23IGS 137 1 2 26WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00020006 -152429841 11 38666162 16 -43578164 21 R2781ASEDRR2781ABVSA +D 1 2 9222278 1 3 -8833174 2 3 -9289661 +B201110 517 9201110 522 9 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00060004 143720411 7 -46696601 11 33022280 14 R2781BBVSAR2781B0209 +D 1 2 8734040 1 3 -8117044 2 3 -8983606 +B201110 517 9201110 522 9 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626L1DDFX +Iigs08_1685 IGS 20120422 +C00060005 796935 3 7399405 4 6081745 6 R2781CBVSAR2781CTXY3 +D 1 2 8845538 1 3 -8353980 2 3 -9092949 +B201110 81419201110 82057 1 pages v1109.23IGS 137 1 2 26WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00020006 -152429775 9 38666232 13 -43578178 16 R2811ASEDRR2811ABVSA +D 1 2 9222591 1 3 -8874314 2 3 -9305701 +B201110 81510201110 820 8 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626IFDDPF +Iigs08_1685 IGS 20120422 +C00060004 143720397 10 -46696532 15 33022287 19 R2811BBVSAR2811B0209 +D 1 2 9082499 1 3 -8676777 2 3 -9225080 +B201110 81434201110 82048 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626L1DDFX +Iigs08_1685 IGS 20120422 +C00060005 796872 3 7399339 4 6081703 5 R2811CBVSAR2811CTXY3 +D 1 2 8895522 1 3 -8468069 2 3 -9037687 diff --git a/input/test_Mariana.txt b/input/test_Mariana.txt new file mode 100644 index 0000000..9e60389 --- /dev/null +++ b/input/test_Mariana.txt @@ -0,0 +1,3 @@ +13.4,215.3,0,Guam +15.2,214.25,0,Saipan +0,0,0,Unmodeled location (near Africa)
\ No newline at end of file diff --git a/input/test_NAmerica.txt b/input/test_NAmerica.txt new file mode 100644 index 0000000..0673e9a --- /dev/null +++ b/input/test_NAmerica.txt @@ -0,0 +1,8 @@ +40,100,0,Kansas +37,122,0,California +48,124,0,Washington +45,69,0,Maine +28,81,0,Florida +18.2,66.5,0,Puerto Rico +65,152,0,Alaska +0,0,0,Unmodeled location (near Africa) diff --git a/input/test_Pacific.txt b/input/test_Pacific.txt new file mode 100644 index 0000000..4034b59 --- /dev/null +++ b/input/test_Pacific.txt @@ -0,0 +1,3 @@ +19.5,155.5,0,Hawaii +-14.3,170.7,0,American Samoa +0,0,0,Unmodeled location (near Africa) diff --git a/input/wa.bfile.txt b/input/wa.bfile.txt new file mode 100644 index 0000000..e178229 --- /dev/null +++ b/input/wa.bfile.txt @@ -0,0 +1,85 @@ +000010*AA*HZTLOBS WHPACIW + H PACIFIC, INCORPORATED 20111027 +000020*10*SKAGIT AIRPORT PACS SACS +000030*12*201110201110MBEEDWARDS M B 4WAB0 +000050*25*0006R2781ABVSAMBE00300003 +000060*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000070*26*FILES\ALL\BVSA278A.11O +000080*27*00061110051709 2409 F 500F 2950IN01020 +000090*27*00061110052209 2409 F 570F 2950IN01020 +000100*25*0006R2811ABVSAMBE00100001 +000110*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000120*26*FILES\ALL\BVSA281A.11O +000130*27*00061110081419 2354 F 420F 3010IN00000 +000140*27*00061110082057 2354 F 610F 3010IN00000 +000150*25*0002R2781ASEDRKAR002 002 +000160*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B FILES\ALL\SEDR278A.11O +000170*27*00021110051600 187 +000180*27*00021110052200 187 +000190*25*0002R2811ASEDRKAR002 002 +000200*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B FILES\ALL\SEDR281A.11O +000210*27*00021110081400 187 +000220*27*00021110082200 187 +000230*25*0004R2781B0209KMK00100001 +000240*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000250*26*FILES\ALL\0209278B.11O +000260*27*00041110051658 2354 F 540F 2950IN00020 +000270*27*00041110052210 2354 F 550F 2950IN00020 +000280*25*0004R2811B0209MBE00500005 +000290*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000300*26*FILES\ALL\0209281B.11O +000310*27*00041110081510 2354 F 400F 3020IN00000 +000320*27*00041110082008 2354 F 560F 3010IN00021 +000330*25*0005R2781CTXY3MBE00400004 +000340*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000350*26*FILES\ALL\TXY3278C.11O +000360*27*00051110051642 2409 F 500F 2950IN01020 +000370*27*00051110052212 2409 F 570F 2950IN01020 +000380*25*0005R2811CTXY3MBE00400004 +000390*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000400*26*FILES\ALL\TXY3281C.11O +000410*27*00051110081434 2409 F 470F 3010IN00000 +000420*27*00051110082048 2409 F 610F 3010IN00001 +000430*25*0006R2781BBVSAMBE00300003 +000440*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000450*26*FILES\ALL\BVSA278B.11O +000460*26*DUPLICATE ENTRY SESSION A +000470*27*00061110051709 2409 F 500F 2950IN01020 +000480*27*00061110052209 2409 F 570F 2950IN01020 +000490*25*0006R2811BBVSAMBE00100001 +000500*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000510*26*FILES\ALL\BVSA281B.11O +000520*26*DUPLICATE ENTRY SESSION A +000530*27*00061110081419 2354 F 420F 3010IN00000 +000540*27*00061110082057 2354 F 610F 3010IN00000 +000550*25*0006R2781CBVSAMBE00300003 +000560*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000570*26*FILES\ALL\BVSA278C.11O +000580*26*DUPLICATE ENTRY SESSION A +000590*27*00061110051709 2409 F 500F 2950IN01020 +000600*27*00061110052209 2409 F 570F 2950IN01020 +000610*25*0006R2811CBVSAMBE00100001 +000620*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000630*26*FILES\ALL\BVSA281C.11O +000640*26*DUPLICATE ENTRY SESSION A +000650*27*00061110081419 2354 F 420F 3010IN00000 +000660*27*00061110082057 2354 F 610F 3010IN00000 +000670*70*001055 TRIMBLE R8_GNSS 4624117085 +000680*70*002055 TRIMBLE NETRS 4518249780 +000690*70*003055 TRIMBLE R8_GNSS 4622115269 +000700*70*004055 TRIMBLE R8_GNSS 4622115259 +000710*70*005055 TRIMBLE R8_GNSS 4629119141 +000720*72*001 TRMR8_GNSS NONE 4624117085 +000730*72*002 TRM29659.00 UNAV 0220061148 +000740*72*003 TRMR8_GNSS NONE 4622115269 +000750*72*004 TRMR8_GNSS NONE 4622115259 +000760*72*005 TRMR8_GNSS NONE 4629119141 +000050*80*0002SEDRO WOOLEY DNR CORS ARP 48311759215N122132578964W WABA +000060*86*0002 51565K Y88NGS -21323 30242A41A +000070*80*0004GP 29020 9 48302721580N122134082432W WABA +000080*86*0004 18534A12Y88NGS -21393 -2859A41A +000090*80*0005POT TXY 3 48281493084N122254990300W WABA +000100*86*0005 37845A12Y88NGS -21909 15936A41A +000110*80*0006BVS A 48274570456N122253386127W WABA +000120*86*0006 25035G N88 -21911 3124A41A + *93* 18.240 2.126 +000130*AA* diff --git a/wa.gfile.txt b/input/wa.gfile.txt index 51cb197..a2ce076 100755..100644 --- a/wa.gfile.txt +++ b/input/wa.gfile.txt @@ -1,25 +1,25 @@ AAA201110 5201110 5SKAGIT AIRPORT PACS SACS -B201110 517 9201110 522 0 1 pages v1109.23IGS 128 1 2 26WHPACI2012 626IFDDFX +B201110 517 9201110 522 0 1 pages v1109.23IGS 137 1 2 26WHPACI2012 626IFDDFX Iigs08_1685 IGS 20120422 C00020006 -152429841 11 38666162 16 -43578164 21 R2781ASEDRR2781ABVSA D 1 2 9222278 1 3 -8833174 2 3 -9289661 -B201110 517 9201110 522 9 1 pages v1109.23IGS 128 1 2 25WHPACI2012 626IFDDFX +B201110 517 9201110 522 9 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626IFDDFX Iigs08_1685 IGS 20120422 C00060004 143720411 7 -46696601 11 33022280 14 R2781BBVSAR2781B0209 D 1 2 8734040 1 3 -8117044 2 3 -8983606 -B201110 517 9201110 522 9 1 pages v1109.23IGS 128 1 2 25WHPACI2012 626L1DDFX +B201110 517 9201110 522 9 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626L1DDFX Iigs08_1685 IGS 20120422 C00060005 796935 3 7399405 4 6081745 6 R2781CBVSAR2781CTXY3 D 1 2 8845538 1 3 -8353980 2 3 -9092949 -B201110 81419201110 82057 1 pages v1109.23IGS 128 1 2 26WHPACI2012 626IFDDFX +B201110 81419201110 82057 1 pages v1109.23IGS 137 1 2 26WHPACI2012 626IFDDFX Iigs08_1685 IGS 20120422 C00020006 -152429775 9 38666232 13 -43578178 16 R2811ASEDRR2811ABVSA D 1 2 9222591 1 3 -8874314 2 3 -9305701 -B201110 81510201110 820 8 1 pages v1109.23IGS 128 1 2 25WHPACI2012 626IFDDPF +B201110 81510201110 820 8 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626IFDDPF Iigs08_1685 IGS 20120422 C00060004 143720397 10 -46696532 15 33022287 19 R2811BBVSAR2811B0209 D 1 2 9082499 1 3 -8676777 2 3 -9225080 -B201110 81434201110 82048 1 pages v1109.23IGS 128 1 2 25WHPACI2012 626L1DDFX +B201110 81434201110 82048 1 pages v1109.23IGS 137 1 2 25WHPACI2012 626L1DDFX Iigs08_1685 IGS 20120422 C00060005 796872 3 7399339 4 6081703 5 R2811CBVSAR2811CTXY3 D 1 2 8895522 1 3 -8468069 2 3 -9037687 diff --git a/output/.svn/all-wcprops b/output/.svn/all-wcprops new file mode 100644 index 0000000..29a1b69 --- /dev/null +++ b/output/.svn/all-wcprops @@ -0,0 +1,95 @@ +K 25 +svn:wc:ra_dav:version-url +V 97 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output +END +test_NAmerica_alt.out +K 25 +svn:wc:ra_dav:version-url +V 119 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/test_NAmerica_alt.out +END +test_NAmerica.out +K 25 +svn:wc:ra_dav:version-url +V 115 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/test_NAmerica.out +END +tfile1.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/tfile1.out +END +vfile1.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/vfile1.out +END +wa.bfile.out +K 25 +svn:wc:ra_dav:version-url +V 110 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/wa.bfile.out +END +vfile2.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/vfile2.out +END +vfile3.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/vfile3.out +END +tvfile.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/tvfile.out +END +newfile.out +K 25 +svn:wc:ra_dav:version-url +V 109 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/newfile.out +END +wa.gfile.out +K 25 +svn:wc:ra_dav:version-url +V 110 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/wa.gfile.out +END +test_Pacific.out +K 25 +svn:wc:ra_dav:version-url +V 114 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/test_Pacific.out +END +dfile1.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/dfile1.out +END +vfile.out +K 25 +svn:wc:ra_dav:version-url +V 107 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/vfile.out +END +dfile2.out +K 25 +svn:wc:ra_dav:version-url +V 108 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/dfile2.out +END +test_Mariana.out +K 25 +svn:wc:ra_dav:version-url +V 114 +/svn/repos/!svn/ver/115924/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output/test_Mariana.out +END diff --git a/output/.svn/entries b/output/.svn/entries new file mode 100644 index 0000000..917ad9e --- /dev/null +++ b/output/.svn/entries @@ -0,0 +1,538 @@ +10 + +dir +121793 +https://source.ngs.noaa.gov/svn/repos/NGSIDBAPI/AnalysisApi/htdp/tags/release-3.2.9/web/HTDP-download/output +https://source.ngs.noaa.gov/svn/repos + + + +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + +53ba9092-9722-c17e-94b9-e0e50a2f61ab + +test_NAmerica_alt.out +file + + + + +2021-04-09T01:01:26.682847Z +4b98bb2fb4b0faa30b6939ec953f49eb +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +914 + +test_NAmerica.out +file + + + + +2021-04-09T01:01:26.653839Z +0d6a99831a5c63b5c325793cba490133 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +955 + +tfile1.out +file + + + + +2021-04-09T01:01:26.596844Z +57e57ad071690434d1f17030c41cad55 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +690 + +vfile1.out +file + + + + +2021-04-09T01:01:26.623843Z +ec8d79d64d627e2839b711bba1e8625e +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +5836 + +wa.bfile.out +file + + + + +2021-04-09T01:01:26.632851Z +2cbc12e59666d66d66a96cc2ef4db0ca +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +7201 + +vfile2.out +file + + + + +2021-04-09T01:01:26.629861Z +164095b3e33dd59cfc4fffec8b9beae9 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +529 + +vfile3.out +file + + + + +2021-04-09T01:01:26.635848Z +abdd8f2eddc7f65a2256a5d5acdae676 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +487 + +tvfile.out +file + + + + +2021-04-09T01:01:26.643853Z +0867035b51a672cbc36833a02fbba4c7 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +589 + +newfile.out +file + + + + +2021-04-09T01:01:26.639844Z +f9220c183c06d5eaa46ebb90b6bfe4ab +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +753 + +wa.gfile.out +file + + + + +2021-04-09T01:01:26.647834Z +e8f8cf05315bbe8d779c68f5c3e58d36 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +3312 + +test_Pacific.out +file + + + + +2021-04-09T01:01:26.650834Z +10399f749559c6689988244bef513f5d +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +680 + +dfile1.out +file + + + + +2021-04-09T01:01:26.581844Z +c99e007eba23665a716c5f3527b0baea +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +365 + +vfile.out +file + + + + +2021-04-09T01:01:26.621835Z +3b91b9b4c3325f5e0a34785e2ab4e708 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +904 + +dfile2.out +file + + + + +2021-04-09T01:01:26.618840Z +65d050d4412a11c458fdae8d6ec660d1 +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +365 + +test_Mariana.out +file + + + + +2021-04-09T01:01:26.626842Z +71c75a64527c08bb1990a07a06a40d1f +2020-05-18T00:13:38.007077Z +115869 +michael.dennis + + + + + + + + + + + + + + + + + + + + + +670 + diff --git a/output/.svn/text-base/dfile1.out.svn-base b/output/.svn/text-base/dfile1.out.svn-base new file mode 100644 index 0000000..b455f27 --- /dev/null +++ b/output/.svn/text-base/dfile1.out.svn-base @@ -0,0 +1,8 @@ + HTDP OUTPUT, VERSION 3.2.9 + + DISPLACEMENTS IN METERS RELATIVE TO NAD_83(2011/CORS96/2007) + FROM 01-01-1985 TO 01-01-1995 (month-day-year) + FROM 1985.000 TO 1995.000 (decimal years) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP +beta 36 40 11.28000 N 121 46 19.92000 W 0.445 -0.259 -0.018 diff --git a/output/.svn/text-base/dfile2.out.svn-base b/output/.svn/text-base/dfile2.out.svn-base new file mode 100644 index 0000000..bbb23ff --- /dev/null +++ b/output/.svn/text-base/dfile2.out.svn-base @@ -0,0 +1,8 @@ + HTDP OUTPUT, VERSION 3.2.9 + + DISPLACEMENTS IN METERS RELATIVE TO NAD_83(2011/CORS96/2007) + FROM 10-16-1989 TO 10-18-1989 (month-day-year) + FROM 1989.789 TO 1989.795 (decimal years) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP +beta 36 40 11.28000 N 121 46 19.92000 W 0.074 -0.001 -0.004 diff --git a/output/.svn/text-base/newfile.out.svn-base b/output/.svn/text-base/newfile.out.svn-base new file mode 100644 index 0000000..233e560 --- /dev/null +++ b/output/.svn/text-base/newfile.out.svn-base @@ -0,0 +1,16 @@ + HTDP OUTPUT, VERSION 3.2.9 + + UPDATED POSITIONS IN NAD_83(2011/CORS96/2007) + FROM 5-06-1991 TO 7-04-1995 (month-day-year) + FROM 1991.345 TO 1995.504 (decimal years) + + OLD COORDINATE NEW COORDINATE VELOCITY DISPLACEMENT + + alpha + LATITUDE 38 06 12.96000 N 38 06 12.96502 N 37.19 mm/yr 0.155 m north + LONGITUDE 122 56 7.80000 W 122 56 7.80406 W -23.79 mm/yr -0.099 m east + ELLIP. HT. 0.000 -0.006 -1.37 mm/yr -0.006 m up + X -2732250.837 -2732250.866 -6.90 mm/yr -0.029 m + Y -4217684.424 -4217684.286 33.10 mm/yr 0.138 m + Z 3914499.164 3914499.282 28.42 mm/yr 0.118 m + diff --git a/output/.svn/text-base/test_Mariana.out.svn-base b/output/.svn/text-base/test_Mariana.out.svn-base new file mode 100644 index 0000000..4b605f2 --- /dev/null +++ b/output/.svn/text-base/test_Mariana.out.svn-base @@ -0,0 +1,14 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(MA11/MARP00) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2020=(2020.000) *** + + 13.4000088594 215.3000089897 1.976 Guam + 15.2000082574 214.2500092122 1.960 Saipan + + Unmodeled location (near Africa) is outside of the modeled region. + diff --git a/output/.svn/text-base/test_NAmerica.out.svn-base b/output/.svn/text-base/test_NAmerica.out.svn-base new file mode 100644 index 0000000..c91a06a --- /dev/null +++ b/output/.svn/text-base/test_NAmerica.out.svn-base @@ -0,0 +1,19 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(2011/CORS96/2007) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2020=(2020.000) *** + + 40.0000059056 100.0000131843 -0.965 Kansas + 37.0000054908 122.0000193408 -0.548 California + 48.0000034473 124.0000184377 -0.286 Washington + 45.0000107518 69.0000049046 -1.138 Maine + 28.0000058327 81.0000058995 -1.556 Florida + 18.2000050384 66.4999985203 -1.879 Puerto Rico + 64.9999966235 152.0000304343 0.483 Alaska + + Unmodeled location (near Africa) is outside of the modeled region. + diff --git a/output/.svn/text-base/test_NAmerica_alt.out.svn-base b/output/.svn/text-base/test_NAmerica_alt.out.svn-base new file mode 100644 index 0000000..f55fccf --- /dev/null +++ b/output/.svn/text-base/test_NAmerica_alt.out.svn-base @@ -0,0 +1,17 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(2011/CORS96/2007) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2010 (2010.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2010=(2010.000) *** + + 40.0000062553 100.0000114584 -0.964 Kansas + 37.0000034059 122.0000150307 -0.546 California + 48.0000038522 124.0000179554 -0.285 Washington + 45.0000101759 69.0000029453 -1.137 Maine + 28.0000056127 81.0000046970 -1.554 Florida + 18.2000039077 66.4999994854 -1.879 Puerto Rico + 64.9999985433 152.0000288257 0.484 Alaska + 0.0000050155 359.9999798125 -1.008 Unmodeled location (near Africa) diff --git a/output/.svn/text-base/test_Pacific.out.svn-base b/output/.svn/text-base/test_Pacific.out.svn-base new file mode 100644 index 0000000..520bbb1 --- /dev/null +++ b/output/.svn/text-base/test_Pacific.out.svn-base @@ -0,0 +1,14 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(PA11/PACP00) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2020=(2020.000) *** + + 19.5000122905 155.5000310534 0.174 Hawaii + -14.2999842111 170.7000346593 0.388 American Samoa + + Unmodeled location (near Africa) is outside of the modeled region. + diff --git a/output/.svn/text-base/tfile1.out.svn-base b/output/.svn/text-base/tfile1.out.svn-base new file mode 100644 index 0000000..8a86a3a --- /dev/null +++ b/output/.svn/text-base/tfile1.out.svn-base @@ -0,0 +1,13 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(2011/CORS96/2007) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + test1 + LATITUDE 40 00 0.00000 N 40 00 0.02126 N 0.81 mm/yr north + LONGITUDE 100 00 0.00000 W 100 00 0.04746 W 1.88 mm/yr east + ELLIP. HT. 0.000 -0.965 m -1.14 mm/yr up + X -849609.759 -849610.666 m 2.09 mm/yr + Y -4818376.378 -4818375.039 m 1.05 mm/yr + Z 4077985.572 4077985.454 m -0.11 mm/yr + diff --git a/output/.svn/text-base/tvfile.out.svn-base b/output/.svn/text-base/tvfile.out.svn-base new file mode 100644 index 0000000..9806639 --- /dev/null +++ b/output/.svn/text-base/tvfile.out.svn-base @@ -0,0 +1,15 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING VELOCITIES FROM ITRF2000 or IGS00/IGb00 TO NAD_83(2011/CORS96/2007) + + INPUT VELOCITIES OUTPUT VELOCITIES + + gamma + latitude = 38.000000000 longitude = 123.000000000 + northward velocity -12.00 2.70 mm/yr + eastward velocity -10.00 3.55 mm/yr + upward velocity 2.00 1.34 mm/yr + x velocity -13.27 3.31 mm/yr + y velocity -2.07 -1.42 mm/yr + z velocity -8.22 2.95 mm/yr + diff --git a/output/.svn/text-base/vfile.out.svn-base b/output/.svn/text-base/vfile.out.svn-base new file mode 100644 index 0000000..f0746ea --- /dev/null +++ b/output/.svn/text-base/vfile.out.svn-base @@ -0,0 +1,19 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) + +alpha +LATITUDE = 38 6 12.96000 N NORTH VELOCITY = 37.19 mm/yr +LONGITUDE = 122 56 7.80000 W EAST VELOCITY = -23.79 mm/yr +ELLIPS. HT. = 0.000 m UP VELOCITY = -1.37 mm/yr +X = -2732250.837 m X VELOCITY = -6.90 mm/yr +Y = -4217684.424 m Y VELOCITY = 33.10 mm/yr +Z = 3914499.164 m Z VELOCITY = 28.42 mm/yr + +beta +LATITUDE = 36 40 11.28000 N NORTH VELOCITY = 37.15 mm/yr +LONGITUDE = 121 46 19.92000 W EAST VELOCITY = -25.83 mm/yr +ELLIPS. HT. = 0.000 m UP VELOCITY = -1.33 mm/yr +X = -2696934.816 m X VELOCITY = -9.71 mm/yr +Y = -4354426.684 m Y VELOCITY = 33.37 mm/yr +Z = 3788064.740 m Z VELOCITY = 29.01 mm/yr diff --git a/output/.svn/text-base/vfile1.out.svn-base b/output/.svn/text-base/vfile1.out.svn-base new file mode 100644 index 0000000..f146f46 --- /dev/null +++ b/output/.svn/text-base/vfile1.out.svn-base @@ -0,0 +1,71 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP + +grid1 0 0 34 0 0.00000 N 118 30 0.00000 W 29.66 -24.84 -1.23 +grid1 0 1 34 0 0.00000 N 118 40 0.00000 W 31.41 -25.03 -1.24 +grid1 0 2 34 0 0.00000 N 118 50 0.00000 W 31.99 -25.64 -1.24 +grid1 0 3 34 0 0.00000 N 119 0 0.00000 W 32.92 -26.43 -1.24 +grid1 0 4 34 0 0.00000 N 119 10 0.00000 W 33.91 -27.16 -1.24 +grid1 1 0 34 5 0.00000 N 118 30 0.00000 W 29.10 -24.82 -1.24 +grid1 1 1 34 5 0.00000 N 118 40 0.00000 W 30.16 -25.26 -1.24 +grid1 1 2 34 5 0.00000 N 118 50 0.00000 W 30.88 -25.73 -1.24 +grid1 1 3 34 5 0.00000 N 119 0 0.00000 W 31.79 -26.26 -1.24 +grid1 1 4 34 5 0.00000 N 119 10 0.00000 W 32.80 -26.75 -1.24 +grid1 2 0 34 10 0.00000 N 118 30 0.00000 W 27.87 -24.72 -1.24 +grid1 2 1 34 10 0.00000 N 118 40 0.00000 W 28.90 -25.26 -1.24 +grid1 2 2 34 10 0.00000 N 118 50 0.00000 W 29.76 -25.74 -1.24 +grid1 2 3 34 10 0.00000 N 119 0 0.00000 W 30.69 -26.19 -1.24 +grid1 2 4 34 10 0.00000 N 119 10 0.00000 W 31.67 -26.60 -1.25 +grid1 3 0 34 15 0.00000 N 118 30 0.00000 W 26.61 -24.58 -1.24 +grid1 3 1 34 15 0.00000 N 118 40 0.00000 W 27.67 -25.17 -1.24 +grid1 3 2 34 15 0.00000 N 118 50 0.00000 W 28.63 -25.68 -1.24 +grid1 3 3 34 15 0.00000 N 119 0 0.00000 W 29.60 -26.13 -1.25 +grid1 3 4 34 15 0.00000 N 119 10 0.00000 W 30.55 -26.56 -1.25 +grid1 4 0 34 20 0.00000 N 118 30 0.00000 W 25.31 -24.31 -1.24 +grid1 4 1 34 20 0.00000 N 118 40 0.00000 W 26.40 -24.96 -1.24 +grid1 4 2 34 20 0.00000 N 118 50 0.00000 W 27.48 -25.54 -1.25 +grid1 4 3 34 20 0.00000 N 119 0 0.00000 W 28.50 -26.06 -1.25 +grid1 4 4 34 20 0.00000 N 119 10 0.00000 W 29.45 -26.53 -1.25 +grid1 5 0 34 25 0.00000 N 118 30 0.00000 W 23.87 -24.00 -1.24 +grid1 5 1 34 25 0.00000 N 118 40 0.00000 W 25.07 -24.55 -1.24 +grid1 5 2 34 25 0.00000 N 118 50 0.00000 W 26.26 -25.27 -1.25 +grid1 5 3 34 25 0.00000 N 119 0 0.00000 W 27.38 -25.91 -1.25 +grid1 5 4 34 25 0.00000 N 119 10 0.00000 W 28.35 -26.47 -1.25 +grid1 6 0 34 30 0.00000 N 118 30 0.00000 W 23.00 -23.04 -1.24 +grid1 6 1 34 30 0.00000 N 118 40 0.00000 W 23.58 -23.74 -1.25 +grid1 6 2 34 30 0.00000 N 118 50 0.00000 W 24.94 -24.77 -1.25 +grid1 6 3 34 30 0.00000 N 119 0 0.00000 W 26.20 -25.64 -1.25 +grid1 6 4 34 30 0.00000 N 119 10 0.00000 W 27.23 -26.35 -1.25 +grid1 7 0 34 35 0.00000 N 118 30 0.00000 W 21.94 -20.75 -1.25 +grid1 7 1 34 35 0.00000 N 118 40 0.00000 W 23.09 -23.08 -1.25 +grid1 7 2 34 35 0.00000 N 118 50 0.00000 W 23.32 -23.71 -1.25 +grid1 7 3 34 35 0.00000 N 119 0 0.00000 W 24.88 -25.08 -1.25 +grid1 7 4 34 35 0.00000 N 119 10 0.00000 W 26.02 -26.12 -1.25 +grid1 8 0 34 40 0.00000 N 118 30 0.00000 W 20.41 -17.20 -1.25 +grid1 8 1 34 40 0.00000 N 118 40 0.00000 W 21.75 -20.41 -1.25 +grid1 8 2 34 40 0.00000 N 118 50 0.00000 W 22.28 -22.50 -1.25 +grid1 8 3 34 40 0.00000 N 119 0 0.00000 W 23.27 -23.98 -1.25 +grid1 8 4 34 40 0.00000 N 119 10 0.00000 W 24.47 -25.65 -1.26 +grid1 9 0 34 45 0.00000 N 118 30 0.00000 W 18.68 -13.65 -1.25 +grid1 9 1 34 45 0.00000 N 118 40 0.00000 W 19.93 -16.12 -1.25 +grid1 9 2 34 45 0.00000 N 118 50 0.00000 W 21.55 -20.21 -1.25 +grid1 9 3 34 45 0.00000 N 119 0 0.00000 W 21.12 -21.81 -1.26 +grid1 9 4 34 45 0.00000 N 119 10 0.00000 W 23.88 -22.82 -1.26 +grid1 10 0 34 50 0.00000 N 118 30 0.00000 W 17.38 -11.35 -1.25 +grid1 10 1 34 50 0.00000 N 118 40 0.00000 W 18.19 -12.70 -1.25 +grid1 10 2 34 50 0.00000 N 118 50 0.00000 W 19.49 -14.92 -1.26 +grid1 10 3 34 50 0.00000 N 119 0 0.00000 W 20.04 -17.64 -1.26 +grid1 10 4 34 50 0.00000 N 119 10 0.00000 W 23.06 -19.58 -1.26 +grid1 11 0 34 55 0.00000 N 118 30 0.00000 W 16.45 -9.91 -1.25 +grid1 11 1 34 55 0.00000 N 118 40 0.00000 W 16.58 -10.82 -1.25 +grid1 11 2 34 55 0.00000 N 118 50 0.00000 W 18.23 -12.45 -1.26 +grid1 11 3 34 55 0.00000 N 119 0 0.00000 W 20.78 -14.26 -1.26 +grid1 11 4 34 55 0.00000 N 119 10 0.00000 W 22.01 -16.63 -1.26 +grid1 12 0 35 0 0.00000 N 118 30 0.00000 W 14.59 -9.12 -1.25 +grid1 12 1 35 0 0.00000 N 118 40 0.00000 W 15.97 -10.27 -1.26 +grid1 12 2 35 0 0.00000 N 118 50 0.00000 W 17.56 -11.42 -1.26 +grid1 12 3 35 0 0.00000 N 119 0 0.00000 W 19.22 -12.83 -1.26 +grid1 12 4 35 0 0.00000 N 119 10 0.00000 W 20.42 -14.51 -1.26 diff --git a/output/.svn/text-base/vfile2.out.svn-base b/output/.svn/text-base/vfile2.out.svn-base new file mode 100644 index 0000000..db9907f --- /dev/null +++ b/output/.svn/text-base/vfile2.out.svn-base @@ -0,0 +1,10 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP + +SEDRO WOOLEY DNR CORS AR 48 31 17.59215 N 122 13 25.78964 W 4.90 5.66 -1.55 +GP 29020 9 48 30 27.21580 N 122 13 40.82432 W 4.91 5.67 -1.55 +POT TXY 3 48 28 14.93084 N 122 25 49.90300 W 5.16 5.92 -1.55 +BVS A 48 27 45.70456 N 122 25 33.86127 W 5.16 5.92 -1.55 diff --git a/output/.svn/text-base/vfile3.out.svn-base b/output/.svn/text-base/vfile3.out.svn-base new file mode 100644 index 0000000..89bade9 --- /dev/null +++ b/output/.svn/text-base/vfile3.out.svn-base @@ -0,0 +1,9 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) +NAME OF LINE LATITUDE LONGITUDE NORTH EAST UP + +line1 0 35 17 28.25504 N 120 18 53.31236 W 32.80 -25.08 -1.28 +line1 1 35 17 28.30000 N 120 15 35.43100 W 32.25 -24.70 -1.28 +line1 2 35 17 28.25504 N 120 12 17.54964 W 31.66 -24.28 -1.28 +line1 3 35 17 28.12016 N 120 8 59.66841 W 31.00 -23.82 -1.28 diff --git a/output/.svn/text-base/wa.bfile.out.svn-base b/output/.svn/text-base/wa.bfile.out.svn-base new file mode 100644 index 0000000..5085176 --- /dev/null +++ b/output/.svn/text-base/wa.bfile.out.svn-base @@ -0,0 +1,88 @@ + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in NAD_83(2011/CORS96/2007)*** + ***CAUTION: Coordinates in this file have been updated to 1-01-2010 = (2010.000) *** +000010*AA*HZTLOBS WHPACIW + H PACIFIC, INCORPORATED 20111027 +000020*10*SKAGIT AIRPORT PACS SACS +000030*12*201110201110MBEEDWARDS M B 4WAB0 +000050*25*0006R2781ABVSAMBE00300003 +000060*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000070*26*FILES\ALL\BVSA278A.11O +000080*27*00061110051709 2409 F 500F 2950IN01020 +000090*27*00061110052209 2409 F 570F 2950IN01020 +000100*25*0006R2811ABVSAMBE00100001 +000110*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000120*26*FILES\ALL\BVSA281A.11O +000130*27*00061110081419 2354 F 420F 3010IN00000 +000140*27*00061110082057 2354 F 610F 3010IN00000 +000150*25*0002R2781ASEDRKAR002 002 +000160*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B FILES\ALL\SEDR278A.11O +000170*27*00021110051600 187 +000180*27*00021110052200 187 +000190*25*0002R2811ASEDRKAR002 002 +000200*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B FILES\ALL\SEDR281A.11O +000210*27*00021110081400 187 +000220*27*00021110082200 187 +000230*25*0004R2781B0209KMK00100001 +000240*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000250*26*FILES\ALL\0209278B.11O +000260*27*00041110051658 2354 F 540F 2950IN00020 +000270*27*00041110052210 2354 F 550F 2950IN00020 +000280*25*0004R2811B0209MBE00500005 +000290*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000300*26*FILES\ALL\0209281B.11O +000310*27*00041110081510 2354 F 400F 3020IN00000 +000320*27*00041110082008 2354 F 560F 3010IN00021 +000330*25*0005R2781CTXY3MBE00400004 +000340*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000350*26*FILES\ALL\TXY3278C.11O +000360*27*00051110051642 2409 F 500F 2950IN01020 +000370*27*00051110052212 2409 F 570F 2950IN01020 +000380*25*0005R2811CTXY3MBE00400004 +000390*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000400*26*FILES\ALL\TXY3281C.11O +000410*27*00051110081434 2409 F 470F 3010IN00000 +000420*27*00051110082048 2409 F 610F 3010IN00001 +000430*25*0006R2781BBVSAMBE00300003 +000440*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000450*26*FILES\ALL\BVSA278B.11O +000460*26*DUPLICATE ENTRY SESSION A +000470*27*00061110051709 2409 F 500F 2950IN01020 +000480*27*00061110052209 2409 F 570F 2950IN01020 +000490*25*0006R2811BBVSAMBE00100001 +000500*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000510*26*FILES\ALL\BVSA281B.11O +000520*26*DUPLICATE ENTRY SESSION A +000530*27*00061110081419 2354 F 420F 3010IN00000 +000540*27*00061110082057 2354 F 610F 3010IN00000 +000550*25*0006R2781CBVSAMBE00300003 +000560*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000570*26*FILES\ALL\BVSA278C.11O +000580*26*DUPLICATE ENTRY SESSION A +000590*27*00061110051709 2409 F 500F 2950IN01020 +000600*27*00061110052209 2409 F 570F 2950IN01020 +000610*25*0006R2811CBVSAMBE00100001 +000620*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000630*26*FILES\ALL\BVSA281C.11O +000640*26*DUPLICATE ENTRY SESSION A +000650*27*00061110081419 2354 F 420F 3010IN00000 +000660*27*00061110082057 2354 F 610F 3010IN00000 +000670*70*001055 TRIMBLE R8_GNSS 4624117085 +000680*70*002055 TRIMBLE NETRS 4518249780 +000690*70*003055 TRIMBLE R8_GNSS 4622115269 +000700*70*004055 TRIMBLE R8_GNSS 4622115259 +000710*70*005055 TRIMBLE R8_GNSS 4629119141 +000720*72*001 TRMR8_GNSS NONE 4624117085 +000730*72*002 TRM29659.00 UNAV 0220061148 +000740*72*003 TRMR8_GNSS NONE 4622115269 +000750*72*004 TRMR8_GNSS NONE 4622115259 +000760*72*005 TRMR8_GNSS NONE 4629119141 +000050*80*0002SEDRO WOOLEY DNR CORS ARP 48311759511N122132578449W WABA +000060*86*0002 51565K Y88NGS -21323 30242A41A +000070*80*0004GP 29020 9 48302721877N122134081916W WABA +000080*86*0004 18534A12Y88NGS -21393 -2859A41A +000090*80*0005POT TXY 3 48281493396N122254989763W WABA +000100*86*0005 37845A12Y88NGS -21909 15936A41A +000110*80*0006BVS A 48274570768N122253385590W WABA +000120*86*0006 25035G N88 -21911 3124A41A + *93* 18.240 2.126 +000130*AA* diff --git a/output/.svn/text-base/wa.gfile.out.svn-base b/output/.svn/text-base/wa.gfile.out.svn-base new file mode 100644 index 0000000..e6719c7 --- /dev/null +++ b/output/.svn/text-base/wa.gfile.out.svn-base @@ -0,0 +1,28 @@ + ***CAUTION: Observations in this file have been updated to 1-01-2010 = (2010.000) *** + ***CAUTION: All GPS interstation vectors have been transformed to NAD_83(2011/CORS96/2007) *** + ***CAUTION: Observations were transformed using HTDP version 3.2.9 *** +AAA2010010120100101SKAGIT AIRPORT PACS SACS ZT +B201110 517 9201110 522 0 1 pages v1109.23IGS 134 1 2 26WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00020006 -152429845 11 38666164 16 -43578171 21 R2781ASEDRR2781ABVSA +D 1 2 9222278 1 3 -8833174 2 3 -9289661 +B201110 517 9201110 522 9 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00060004 143720414 7 -46696604 11 33022288 14 R2781BBVSAR2781B0209 +D 1 2 8734040 1 3 -8117044 2 3 -8983606 +B201110 517 9201110 522 9 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626L1DDFX +Iigs08_1685 IGS 20120422 +C00060005 796935 3 7399406 4 6081744 6 R2781CBVSAR2781CTXY3 +D 1 2 8845538 1 3 -8353980 2 3 -9092949 +B201110 81419201110 82057 1 pages v1109.23IGS 134 1 2 26WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00020006 -152429779 9 38666234 13 -43578185 16 R2811ASEDRR2811ABVSA +D 1 2 9222591 1 3 -8874314 2 3 -9305701 +B201110 81510201110 820 8 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626IFDDPF +Iigs08_1685 IGS 20120422 +C00060004 143720400 10 -46696535 15 33022295 19 R2811BBVSAR2811B0209 +D 1 2 9082499 1 3 -8676777 2 3 -9225080 +B201110 81434201110 82048 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626L1DDFX +Iigs08_1685 IGS 20120422 +C00060005 796872 3 7399340 4 6081702 5 R2811CBVSAR2811CTXY3 +D 1 2 8895522 1 3 -8468069 2 3 -9037687 diff --git a/output/dfile1.out b/output/dfile1.out new file mode 100644 index 0000000..b455f27 --- /dev/null +++ b/output/dfile1.out @@ -0,0 +1,8 @@ + HTDP OUTPUT, VERSION 3.2.9 + + DISPLACEMENTS IN METERS RELATIVE TO NAD_83(2011/CORS96/2007) + FROM 01-01-1985 TO 01-01-1995 (month-day-year) + FROM 1985.000 TO 1995.000 (decimal years) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP +beta 36 40 11.28000 N 121 46 19.92000 W 0.445 -0.259 -0.018 diff --git a/output/dfile2.out b/output/dfile2.out new file mode 100644 index 0000000..bbb23ff --- /dev/null +++ b/output/dfile2.out @@ -0,0 +1,8 @@ + HTDP OUTPUT, VERSION 3.2.9 + + DISPLACEMENTS IN METERS RELATIVE TO NAD_83(2011/CORS96/2007) + FROM 10-16-1989 TO 10-18-1989 (month-day-year) + FROM 1989.789 TO 1989.795 (decimal years) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP +beta 36 40 11.28000 N 121 46 19.92000 W 0.074 -0.001 -0.004 diff --git a/output/newfile.out b/output/newfile.out new file mode 100644 index 0000000..233e560 --- /dev/null +++ b/output/newfile.out @@ -0,0 +1,16 @@ + HTDP OUTPUT, VERSION 3.2.9 + + UPDATED POSITIONS IN NAD_83(2011/CORS96/2007) + FROM 5-06-1991 TO 7-04-1995 (month-day-year) + FROM 1991.345 TO 1995.504 (decimal years) + + OLD COORDINATE NEW COORDINATE VELOCITY DISPLACEMENT + + alpha + LATITUDE 38 06 12.96000 N 38 06 12.96502 N 37.19 mm/yr 0.155 m north + LONGITUDE 122 56 7.80000 W 122 56 7.80406 W -23.79 mm/yr -0.099 m east + ELLIP. HT. 0.000 -0.006 -1.37 mm/yr -0.006 m up + X -2732250.837 -2732250.866 -6.90 mm/yr -0.029 m + Y -4217684.424 -4217684.286 33.10 mm/yr 0.138 m + Z 3914499.164 3914499.282 28.42 mm/yr 0.118 m + diff --git a/output/test_Mariana.out b/output/test_Mariana.out new file mode 100644 index 0000000..4b605f2 --- /dev/null +++ b/output/test_Mariana.out @@ -0,0 +1,14 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(MA11/MARP00) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2020=(2020.000) *** + + 13.4000088594 215.3000089897 1.976 Guam + 15.2000082574 214.2500092122 1.960 Saipan + + Unmodeled location (near Africa) is outside of the modeled region. + diff --git a/output/test_NAmerica.out b/output/test_NAmerica.out new file mode 100644 index 0000000..c91a06a --- /dev/null +++ b/output/test_NAmerica.out @@ -0,0 +1,19 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(2011/CORS96/2007) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2020=(2020.000) *** + + 40.0000059056 100.0000131843 -0.965 Kansas + 37.0000054908 122.0000193408 -0.548 California + 48.0000034473 124.0000184377 -0.286 Washington + 45.0000107518 69.0000049046 -1.138 Maine + 28.0000058327 81.0000058995 -1.556 Florida + 18.2000050384 66.4999985203 -1.879 Puerto Rico + 64.9999966235 152.0000304343 0.483 Alaska + + Unmodeled location (near Africa) is outside of the modeled region. + diff --git a/output/test_NAmerica_alt.out b/output/test_NAmerica_alt.out new file mode 100644 index 0000000..f55fccf --- /dev/null +++ b/output/test_NAmerica_alt.out @@ -0,0 +1,17 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(2011/CORS96/2007) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2010 (2010.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2010=(2010.000) *** + + 40.0000062553 100.0000114584 -0.964 Kansas + 37.0000034059 122.0000150307 -0.546 California + 48.0000038522 124.0000179554 -0.285 Washington + 45.0000101759 69.0000029453 -1.137 Maine + 28.0000056127 81.0000046970 -1.554 Florida + 18.2000039077 66.4999994854 -1.879 Puerto Rico + 64.9999985433 152.0000288257 0.484 Alaska + 0.0000050155 359.9999798125 -1.008 Unmodeled location (near Africa) diff --git a/output/test_Pacific.out b/output/test_Pacific.out new file mode 100644 index 0000000..520bbb1 --- /dev/null +++ b/output/test_Pacific.out @@ -0,0 +1,14 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(PA11/PACP00) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in ITRF2014 or IGS14/IGb14 *** + ***CAUTION: Coordinates in this file have been updated to 1-01-2020=(2020.000) *** + + 19.5000122905 155.5000310534 0.174 Hawaii + -14.2999842111 170.7000346593 0.388 American Samoa + + Unmodeled location (near Africa) is outside of the modeled region. + diff --git a/output/tfile1.out b/output/tfile1.out new file mode 100644 index 0000000..8a86a3a --- /dev/null +++ b/output/tfile1.out @@ -0,0 +1,13 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING POSITIONS FROM NAD_83(2011/CORS96/2007) (EPOCH = 01-01-2010 (2010.0000)) + TO ITRF2014 or IGS14/IGb14 (EPOCH = 01-01-2020 (2020.0000)) + + test1 + LATITUDE 40 00 0.00000 N 40 00 0.02126 N 0.81 mm/yr north + LONGITUDE 100 00 0.00000 W 100 00 0.04746 W 1.88 mm/yr east + ELLIP. HT. 0.000 -0.965 m -1.14 mm/yr up + X -849609.759 -849610.666 m 2.09 mm/yr + Y -4818376.378 -4818375.039 m 1.05 mm/yr + Z 4077985.572 4077985.454 m -0.11 mm/yr + diff --git a/output/tvfile.out b/output/tvfile.out new file mode 100644 index 0000000..9806639 --- /dev/null +++ b/output/tvfile.out @@ -0,0 +1,15 @@ + HTDP OUTPUT, VERSION 3.2.9 + + TRANSFORMING VELOCITIES FROM ITRF2000 or IGS00/IGb00 TO NAD_83(2011/CORS96/2007) + + INPUT VELOCITIES OUTPUT VELOCITIES + + gamma + latitude = 38.000000000 longitude = 123.000000000 + northward velocity -12.00 2.70 mm/yr + eastward velocity -10.00 3.55 mm/yr + upward velocity 2.00 1.34 mm/yr + x velocity -13.27 3.31 mm/yr + y velocity -2.07 -1.42 mm/yr + z velocity -8.22 2.95 mm/yr + diff --git a/output/vfile.out b/output/vfile.out new file mode 100644 index 0000000..f0746ea --- /dev/null +++ b/output/vfile.out @@ -0,0 +1,19 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) + +alpha +LATITUDE = 38 6 12.96000 N NORTH VELOCITY = 37.19 mm/yr +LONGITUDE = 122 56 7.80000 W EAST VELOCITY = -23.79 mm/yr +ELLIPS. HT. = 0.000 m UP VELOCITY = -1.37 mm/yr +X = -2732250.837 m X VELOCITY = -6.90 mm/yr +Y = -4217684.424 m Y VELOCITY = 33.10 mm/yr +Z = 3914499.164 m Z VELOCITY = 28.42 mm/yr + +beta +LATITUDE = 36 40 11.28000 N NORTH VELOCITY = 37.15 mm/yr +LONGITUDE = 121 46 19.92000 W EAST VELOCITY = -25.83 mm/yr +ELLIPS. HT. = 0.000 m UP VELOCITY = -1.33 mm/yr +X = -2696934.816 m X VELOCITY = -9.71 mm/yr +Y = -4354426.684 m Y VELOCITY = 33.37 mm/yr +Z = 3788064.740 m Z VELOCITY = 29.01 mm/yr diff --git a/output/vfile1.out b/output/vfile1.out new file mode 100644 index 0000000..f146f46 --- /dev/null +++ b/output/vfile1.out @@ -0,0 +1,71 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP + +grid1 0 0 34 0 0.00000 N 118 30 0.00000 W 29.66 -24.84 -1.23 +grid1 0 1 34 0 0.00000 N 118 40 0.00000 W 31.41 -25.03 -1.24 +grid1 0 2 34 0 0.00000 N 118 50 0.00000 W 31.99 -25.64 -1.24 +grid1 0 3 34 0 0.00000 N 119 0 0.00000 W 32.92 -26.43 -1.24 +grid1 0 4 34 0 0.00000 N 119 10 0.00000 W 33.91 -27.16 -1.24 +grid1 1 0 34 5 0.00000 N 118 30 0.00000 W 29.10 -24.82 -1.24 +grid1 1 1 34 5 0.00000 N 118 40 0.00000 W 30.16 -25.26 -1.24 +grid1 1 2 34 5 0.00000 N 118 50 0.00000 W 30.88 -25.73 -1.24 +grid1 1 3 34 5 0.00000 N 119 0 0.00000 W 31.79 -26.26 -1.24 +grid1 1 4 34 5 0.00000 N 119 10 0.00000 W 32.80 -26.75 -1.24 +grid1 2 0 34 10 0.00000 N 118 30 0.00000 W 27.87 -24.72 -1.24 +grid1 2 1 34 10 0.00000 N 118 40 0.00000 W 28.90 -25.26 -1.24 +grid1 2 2 34 10 0.00000 N 118 50 0.00000 W 29.76 -25.74 -1.24 +grid1 2 3 34 10 0.00000 N 119 0 0.00000 W 30.69 -26.19 -1.24 +grid1 2 4 34 10 0.00000 N 119 10 0.00000 W 31.67 -26.60 -1.25 +grid1 3 0 34 15 0.00000 N 118 30 0.00000 W 26.61 -24.58 -1.24 +grid1 3 1 34 15 0.00000 N 118 40 0.00000 W 27.67 -25.17 -1.24 +grid1 3 2 34 15 0.00000 N 118 50 0.00000 W 28.63 -25.68 -1.24 +grid1 3 3 34 15 0.00000 N 119 0 0.00000 W 29.60 -26.13 -1.25 +grid1 3 4 34 15 0.00000 N 119 10 0.00000 W 30.55 -26.56 -1.25 +grid1 4 0 34 20 0.00000 N 118 30 0.00000 W 25.31 -24.31 -1.24 +grid1 4 1 34 20 0.00000 N 118 40 0.00000 W 26.40 -24.96 -1.24 +grid1 4 2 34 20 0.00000 N 118 50 0.00000 W 27.48 -25.54 -1.25 +grid1 4 3 34 20 0.00000 N 119 0 0.00000 W 28.50 -26.06 -1.25 +grid1 4 4 34 20 0.00000 N 119 10 0.00000 W 29.45 -26.53 -1.25 +grid1 5 0 34 25 0.00000 N 118 30 0.00000 W 23.87 -24.00 -1.24 +grid1 5 1 34 25 0.00000 N 118 40 0.00000 W 25.07 -24.55 -1.24 +grid1 5 2 34 25 0.00000 N 118 50 0.00000 W 26.26 -25.27 -1.25 +grid1 5 3 34 25 0.00000 N 119 0 0.00000 W 27.38 -25.91 -1.25 +grid1 5 4 34 25 0.00000 N 119 10 0.00000 W 28.35 -26.47 -1.25 +grid1 6 0 34 30 0.00000 N 118 30 0.00000 W 23.00 -23.04 -1.24 +grid1 6 1 34 30 0.00000 N 118 40 0.00000 W 23.58 -23.74 -1.25 +grid1 6 2 34 30 0.00000 N 118 50 0.00000 W 24.94 -24.77 -1.25 +grid1 6 3 34 30 0.00000 N 119 0 0.00000 W 26.20 -25.64 -1.25 +grid1 6 4 34 30 0.00000 N 119 10 0.00000 W 27.23 -26.35 -1.25 +grid1 7 0 34 35 0.00000 N 118 30 0.00000 W 21.94 -20.75 -1.25 +grid1 7 1 34 35 0.00000 N 118 40 0.00000 W 23.09 -23.08 -1.25 +grid1 7 2 34 35 0.00000 N 118 50 0.00000 W 23.32 -23.71 -1.25 +grid1 7 3 34 35 0.00000 N 119 0 0.00000 W 24.88 -25.08 -1.25 +grid1 7 4 34 35 0.00000 N 119 10 0.00000 W 26.02 -26.12 -1.25 +grid1 8 0 34 40 0.00000 N 118 30 0.00000 W 20.41 -17.20 -1.25 +grid1 8 1 34 40 0.00000 N 118 40 0.00000 W 21.75 -20.41 -1.25 +grid1 8 2 34 40 0.00000 N 118 50 0.00000 W 22.28 -22.50 -1.25 +grid1 8 3 34 40 0.00000 N 119 0 0.00000 W 23.27 -23.98 -1.25 +grid1 8 4 34 40 0.00000 N 119 10 0.00000 W 24.47 -25.65 -1.26 +grid1 9 0 34 45 0.00000 N 118 30 0.00000 W 18.68 -13.65 -1.25 +grid1 9 1 34 45 0.00000 N 118 40 0.00000 W 19.93 -16.12 -1.25 +grid1 9 2 34 45 0.00000 N 118 50 0.00000 W 21.55 -20.21 -1.25 +grid1 9 3 34 45 0.00000 N 119 0 0.00000 W 21.12 -21.81 -1.26 +grid1 9 4 34 45 0.00000 N 119 10 0.00000 W 23.88 -22.82 -1.26 +grid1 10 0 34 50 0.00000 N 118 30 0.00000 W 17.38 -11.35 -1.25 +grid1 10 1 34 50 0.00000 N 118 40 0.00000 W 18.19 -12.70 -1.25 +grid1 10 2 34 50 0.00000 N 118 50 0.00000 W 19.49 -14.92 -1.26 +grid1 10 3 34 50 0.00000 N 119 0 0.00000 W 20.04 -17.64 -1.26 +grid1 10 4 34 50 0.00000 N 119 10 0.00000 W 23.06 -19.58 -1.26 +grid1 11 0 34 55 0.00000 N 118 30 0.00000 W 16.45 -9.91 -1.25 +grid1 11 1 34 55 0.00000 N 118 40 0.00000 W 16.58 -10.82 -1.25 +grid1 11 2 34 55 0.00000 N 118 50 0.00000 W 18.23 -12.45 -1.26 +grid1 11 3 34 55 0.00000 N 119 0 0.00000 W 20.78 -14.26 -1.26 +grid1 11 4 34 55 0.00000 N 119 10 0.00000 W 22.01 -16.63 -1.26 +grid1 12 0 35 0 0.00000 N 118 30 0.00000 W 14.59 -9.12 -1.25 +grid1 12 1 35 0 0.00000 N 118 40 0.00000 W 15.97 -10.27 -1.26 +grid1 12 2 35 0 0.00000 N 118 50 0.00000 W 17.56 -11.42 -1.26 +grid1 12 3 35 0 0.00000 N 119 0 0.00000 W 19.22 -12.83 -1.26 +grid1 12 4 35 0 0.00000 N 119 10 0.00000 W 20.42 -14.51 -1.26 diff --git a/output/vfile2.out b/output/vfile2.out new file mode 100644 index 0000000..db9907f --- /dev/null +++ b/output/vfile2.out @@ -0,0 +1,10 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) + +NAME OF SITE LATITUDE LONGITUDE NORTH EAST UP + +SEDRO WOOLEY DNR CORS AR 48 31 17.59215 N 122 13 25.78964 W 4.90 5.66 -1.55 +GP 29020 9 48 30 27.21580 N 122 13 40.82432 W 4.91 5.67 -1.55 +POT TXY 3 48 28 14.93084 N 122 25 49.90300 W 5.16 5.92 -1.55 +BVS A 48 27 45.70456 N 122 25 33.86127 W 5.16 5.92 -1.55 diff --git a/output/vfile3.out b/output/vfile3.out new file mode 100644 index 0000000..89bade9 --- /dev/null +++ b/output/vfile3.out @@ -0,0 +1,9 @@ + HTDP OUTPUT, VERSION 3.2.9 + +VELOCITIES IN MM/YR RELATIVE TO NAD_83(2011/CORS96/2007) +NAME OF LINE LATITUDE LONGITUDE NORTH EAST UP + +line1 0 35 17 28.25504 N 120 18 53.31236 W 32.80 -25.08 -1.28 +line1 1 35 17 28.30000 N 120 15 35.43100 W 32.25 -24.70 -1.28 +line1 2 35 17 28.25504 N 120 12 17.54964 W 31.66 -24.28 -1.28 +line1 3 35 17 28.12016 N 120 8 59.66841 W 31.00 -23.82 -1.28 diff --git a/output/wa.bfile.out b/output/wa.bfile.out new file mode 100644 index 0000000..5085176 --- /dev/null +++ b/output/wa.bfile.out @@ -0,0 +1,88 @@ + ***CAUTION: This file was processed using HTDP version 3.2.9 *** + ***CAUTION: Coordinates in this file are in NAD_83(2011/CORS96/2007)*** + ***CAUTION: Coordinates in this file have been updated to 1-01-2010 = (2010.000) *** +000010*AA*HZTLOBS WHPACIW + H PACIFIC, INCORPORATED 20111027 +000020*10*SKAGIT AIRPORT PACS SACS +000030*12*201110201110MBEEDWARDS M B 4WAB0 +000050*25*0006R2781ABVSAMBE00300003 +000060*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000070*26*FILES\ALL\BVSA278A.11O +000080*27*00061110051709 2409 F 500F 2950IN01020 +000090*27*00061110052209 2409 F 570F 2950IN01020 +000100*25*0006R2811ABVSAMBE00100001 +000110*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000120*26*FILES\ALL\BVSA281A.11O +000130*27*00061110081419 2354 F 420F 3010IN00000 +000140*27*00061110082057 2354 F 610F 3010IN00000 +000150*25*0002R2781ASEDRKAR002 002 +000160*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B FILES\ALL\SEDR278A.11O +000170*27*00021110051600 187 +000180*27*00021110052200 187 +000190*25*0002R2811ASEDRKAR002 002 +000200*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B FILES\ALL\SEDR281A.11O +000210*27*00021110081400 187 +000220*27*00021110082200 187 +000230*25*0004R2781B0209KMK00100001 +000240*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000250*26*FILES\ALL\0209278B.11O +000260*27*00041110051658 2354 F 540F 2950IN00020 +000270*27*00041110052210 2354 F 550F 2950IN00020 +000280*25*0004R2811B0209MBE00500005 +000290*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000300*26*FILES\ALL\0209281B.11O +000310*27*00041110081510 2354 F 400F 3020IN00000 +000320*27*00041110082008 2354 F 560F 3010IN00021 +000330*25*0005R2781CTXY3MBE00400004 +000340*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000350*26*FILES\ALL\TXY3278C.11O +000360*27*00051110051642 2409 F 500F 2950IN01020 +000370*27*00051110052212 2409 F 570F 2950IN01020 +000380*25*0005R2811CTXY3MBE00400004 +000390*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000400*26*FILES\ALL\TXY3281C.11O +000410*27*00051110081434 2409 F 470F 3010IN00000 +000420*27*00051110082048 2409 F 610F 3010IN00001 +000430*25*0006R2781BBVSAMBE00300003 +000440*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000450*26*FILES\ALL\BVSA278B.11O +000460*26*DUPLICATE ENTRY SESSION A +000470*27*00061110051709 2409 F 500F 2950IN01020 +000480*27*00061110052209 2409 F 570F 2950IN01020 +000490*25*0006R2811BBVSAMBE00100001 +000500*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\B +000510*26*FILES\ALL\BVSA281B.11O +000520*26*DUPLICATE ENTRY SESSION A +000530*27*00061110081419 2354 F 420F 3010IN00000 +000540*27*00061110082057 2354 F 610F 3010IN00000 +000550*25*0006R2781CBVSAMBE00300003 +000560*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000570*26*FILES\ALL\BVSA278C.11O +000580*26*DUPLICATE ENTRY SESSION A +000590*27*00061110051709 2409 F 500F 2950IN01020 +000600*27*00061110052209 2409 F 570F 2950IN01020 +000610*25*0006R2811CBVSAMBE00100001 +000620*26*FILE: C:\...NTS\SKAGIT AIRPORT 2011 PACS SACS\C +000630*26*FILES\ALL\BVSA281C.11O +000640*26*DUPLICATE ENTRY SESSION A +000650*27*00061110081419 2354 F 420F 3010IN00000 +000660*27*00061110082057 2354 F 610F 3010IN00000 +000670*70*001055 TRIMBLE R8_GNSS 4624117085 +000680*70*002055 TRIMBLE NETRS 4518249780 +000690*70*003055 TRIMBLE R8_GNSS 4622115269 +000700*70*004055 TRIMBLE R8_GNSS 4622115259 +000710*70*005055 TRIMBLE R8_GNSS 4629119141 +000720*72*001 TRMR8_GNSS NONE 4624117085 +000730*72*002 TRM29659.00 UNAV 0220061148 +000740*72*003 TRMR8_GNSS NONE 4622115269 +000750*72*004 TRMR8_GNSS NONE 4622115259 +000760*72*005 TRMR8_GNSS NONE 4629119141 +000050*80*0002SEDRO WOOLEY DNR CORS ARP 48311759511N122132578449W WABA +000060*86*0002 51565K Y88NGS -21323 30242A41A +000070*80*0004GP 29020 9 48302721877N122134081916W WABA +000080*86*0004 18534A12Y88NGS -21393 -2859A41A +000090*80*0005POT TXY 3 48281493396N122254989763W WABA +000100*86*0005 37845A12Y88NGS -21909 15936A41A +000110*80*0006BVS A 48274570768N122253385590W WABA +000120*86*0006 25035G N88 -21911 3124A41A + *93* 18.240 2.126 +000130*AA* diff --git a/output/wa.gfile.out b/output/wa.gfile.out new file mode 100644 index 0000000..e6719c7 --- /dev/null +++ b/output/wa.gfile.out @@ -0,0 +1,28 @@ + ***CAUTION: Observations in this file have been updated to 1-01-2010 = (2010.000) *** + ***CAUTION: All GPS interstation vectors have been transformed to NAD_83(2011/CORS96/2007) *** + ***CAUTION: Observations were transformed using HTDP version 3.2.9 *** +AAA2010010120100101SKAGIT AIRPORT PACS SACS ZT +B201110 517 9201110 522 0 1 pages v1109.23IGS 134 1 2 26WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00020006 -152429845 11 38666164 16 -43578171 21 R2781ASEDRR2781ABVSA +D 1 2 9222278 1 3 -8833174 2 3 -9289661 +B201110 517 9201110 522 9 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00060004 143720414 7 -46696604 11 33022288 14 R2781BBVSAR2781B0209 +D 1 2 8734040 1 3 -8117044 2 3 -8983606 +B201110 517 9201110 522 9 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626L1DDFX +Iigs08_1685 IGS 20120422 +C00060005 796935 3 7399406 4 6081744 6 R2781CBVSAR2781CTXY3 +D 1 2 8845538 1 3 -8353980 2 3 -9092949 +B201110 81419201110 82057 1 pages v1109.23IGS 134 1 2 26WHPACI2012 626IFDDFX +Iigs08_1685 IGS 20120422 +C00020006 -152429779 9 38666234 13 -43578185 16 R2811ASEDRR2811ABVSA +D 1 2 9222591 1 3 -8874314 2 3 -9305701 +B201110 81510201110 820 8 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626IFDDPF +Iigs08_1685 IGS 20120422 +C00060004 143720400 10 -46696535 15 33022295 19 R2811BBVSAR2811B0209 +D 1 2 9082499 1 3 -8676777 2 3 -9225080 +B201110 81434201110 82048 1 pages v1109.23IGS 134 1 2 25WHPACI2012 626L1DDFX +Iigs08_1685 IGS 20120422 +C00060005 796872 3 7399340 4 6081702 5 R2811CBVSAR2811CTXY3 +D 1 2 8895522 1 3 -8468069 2 3 -9037687 diff --git a/pc_build.txt b/pc_build.txt deleted file mode 100644 index be028e1..0000000 --- a/pc_build.txt +++ /dev/null @@ -1,15 +0,0 @@ -Make sure MYS2 is installed. By default help desk install 64-bit. Make sure help desk install 32-bit also. - -So the package you need for 32-bit is: mingw-w64-i686-toolchain - -Launch the following command -C:\MSYS2\msys2_shell.cmd - -Bruce.Tran@NGS-L-004085843 MSYS /c/tmp/htdp -$ export PATH=$PATH:/c/msys2/mingw32/bin - -Bruce.Tran@NGS-L-004085843 MSYS /c/tmp/htdp -$ gfortran -DNGS_PC_ENV -static -m32 -Wall -Wtabs *.f -o htdp - - -After htdp.exe is created copy to htdp directory and commit back into svn. diff --git a/pc_build_readme.txt b/pc_build_readme.txt new file mode 100644 index 0000000..4d2ca66 --- /dev/null +++ b/pc_build_readme.txt @@ -0,0 +1,28 @@ +Install MYS2 (https://www.msys2.org/). + +Install both 64- and 32-bit (https://www.msys2.org/wiki/MSYS2-installation/). + +Install following package for 32-bit: mingw-w64-i686-toolchain + +Details on package management at https://www.msys2.org/wiki/Using-packages/. + + +Assume here that msys2 installed at C:/msys2/ +and mingw32 installed at C:/msys2/mingw2/bin/ + + +Once installed with packages, launch msys2.exe +(at C:/msys2/, or create shortcut). + + +Change path to location of source code (this path just an example): +[user computer] MSYS ~ +$ cd /c/temp/htdp + +Set system path to compiler (must go to where msys2 and mingw32 installed): +[user computer] MSYS /c/temp/htdp +$ export PATH=$PATH:/c/msys2/mingw32/bin + +Run compiler with following command: +[user computer] MSYS /c/temp/htdp +$ gfortran -DNGS_PC_ENV -static -m32 -Wall -Wtabs *.f -o htdp
\ No newline at end of file |