hdiff output

r30462/commons.f90 2016-05-17 11:30:10.483311003 +0100 r30461/commons.f90 2016-05-17 11:30:11.383322816 +0100
 81:       INTEGER, ALLOCATABLE :: NPERMSIZE(:), PERMGROUP(:) 81:       INTEGER, ALLOCATABLE :: NPERMSIZE(:), PERMGROUP(:)
 82:       INTEGER, ALLOCATABLE :: NSETS(:) 82:       INTEGER, ALLOCATABLE :: NSETS(:)
 83:       INTEGER, ALLOCATABLE :: SETS(:,:) 83:       INTEGER, ALLOCATABLE :: SETS(:,:)
 84:       INTEGER, ALLOCATABLE :: TAGNUM(:) 84:       INTEGER, ALLOCATABLE :: TAGNUM(:)
 85:       LOGICAL, ALLOCATABLE, DIMENSION(:) :: FROZEN, INTFROZEN 85:       LOGICAL, ALLOCATABLE, DIMENSION(:) :: FROZEN, INTFROZEN
 86:  86: 
 87:       INTEGER NFSTART, NFFINISH, NINTS, NCONNMIN, CONNECTMIN1, CONNECTMIN2, NFREEZE, NPATHS, PARALLEL, NRANDOMMETRIC, & 87:       INTEGER NFSTART, NFFINISH, NINTS, NCONNMIN, CONNECTMIN1, CONNECTMIN2, NFREEZE, NPATHS, PARALLEL, NRANDOMMETRIC, &
 88:      &        METRICUPAIR, METMATMAX, NPRUNE   88:      &        METRICUPAIR, METMATMAX, NPRUNE  
 89:  89: 
 90:       LOGICAL YESNO, TEST1, TEST2, DEBUG, PRINTT, ADDPT, ADDPT2, ADDPT3, TWOD, BULKT, ANGLEAXIS, TAGT, & 90:       LOGICAL YESNO, TEST1, TEST2, DEBUG, PRINTT, ADDPT, ADDPT2, ADDPT3, TWOD, BULKT, ANGLEAXIS, TAGT, &
 91:      &        CHARMMT, AMBERT, STARTFROMPATH, RELATIVEET, EXTRACTTSFILET, & 91:      &        CHARMMT, AMBERT, STARTFROMPATH, RELATIVEET, &
 92:      &        KMCT, UNRST, KMCCOMMITT, REGROUPT, REGROUPRATET, REGROUPPET, NOPOINTS, ADDPATH, NGTT, GTT, GT2T, & 92:      &        KMCT, UNRST, KMCCOMMITT, REGROUPT, REGROUPRATET, REGROUPPET, NOPOINTS, ADDPATH, NGTT, GTT, GT2T, &
 93:      &        DIJKSTRAT, DIJPAIRT, DIJINITT, EXTRACTMINT, EXTRACTMINFILET, EXTRACTTST, DIJKSTRAWAITT, UNTRAPMETRICT, & 93:      &        DIJKSTRAT, DIJPAIRT, DIJINITT, EXTRACTMINT, EXTRACTTST, DIJKSTRAWAITT, UNTRAPMETRICT, &
 94:      &        EXPCOSTFUNCTION, COPYOPTIMT, CALCORDERT, CONNECTREGIONT, SHORTCUTT, MERGEDBT, UNTRAPT, AMHT,  AMHALLATOMMINT, & 94:      &        EXPCOSTFUNCTION, COPYOPTIMT, CALCORDERT, CONNECTREGIONT, SHORTCUTT, MERGEDBT, UNTRAPT, AMHT,  AMHALLATOMMINT, &
 95:      &        CHECKCONNECTIONST, AMHALLATOMTST, AMHQT,AMHQENGMINT, AMHQCONTT ,AMHRMSDT, AMHRELQT, AMH_RELCOT, DIAGT, ARNOLDIT, & 95:      &        CHECKCONNECTIONST, AMHALLATOMTST, AMHQT,AMHQENGMINT, AMHQCONTT ,AMHRMSDT, AMHRELQT, AMH_RELCOT, DIAGT, ARNOLDIT, &
 96:      &        GT2Sparse, GT2Switch, GT2AltPbb, GT2Rescale, GT2Normalise, GT2DisconnectSources, BARRIERSORT, & 96:      &        GT2Sparse, GT2Switch, GT2AltPbb, GT2Rescale, GT2Normalise, GT2DisconnectSources, BARRIERSORT, &
 97:      &        PERMDIST, PERMISOMER, RIGIDBODY, DIJINITSTARTT, DIJINITCONTT, RETAINSP, REMOVESP, NOFRQS, & 97:      &        PERMDIST, PERMISOMER, RIGIDBODY, DIJINITSTARTT, DIJINITCONTT, RETAINSP, REMOVESP, NOFRQS, &
 98:      &        BARRIERSHORT, FREEZE, RATESHORT, DUMMYRUNT, REWEIGHTT, REGROUPFREET, RFMULTIT, REGROUPFREEABT, READMINT, & 98:      &        BARRIERSHORT, FREEZE, RATESHORT, DUMMYRUNT, REWEIGHTT, REGROUPFREET, RFMULTIT, REGROUPFREEABT, READMINT, &
 99:      &        DUMPGROUPST, FREEPAIRT, KSHORTESTPATHST, KSHORT_FULL_PRINTT, DIJINITFLYT, BHINTERPT, ICINTERPT, & 99:      &        DUMPGROUPST, FREEPAIRT, KSHORTESTPATHST, KSHORT_FULL_PRINTT, DIJINITFLYT, BHINTERPT, ICINTERPT, &
100:      &        DUMMYTST, DOCKT, DSTAGE(6), USEPAIRST, LOWESTFRQT, BISECTT, NGTDISCONNECTALL, ANGLEAXIS2, TFOLDT, &100:      &        DUMMYTST, DOCKT, DSTAGE(6), USEPAIRST, LOWESTFRQT, BISECTT, NGTDISCONNECTALL, ANGLEAXIS2, TFOLDT, &
101:      &        SLURMT, INDEXCOSTFUNCTION, CVT, CVMINIMAT, DOST, IMFRQT, CLOSEFILEST, PULLT, FRICTIONT, ATOMMATCHFULL, &101:      &        SLURMT, INDEXCOSTFUNCTION, CVT, CVMINIMAT, DOST, IMFRQT, CLOSEFILEST, PULLT, FRICTIONT, ATOMMATCHFULL, &
102:      &        INTCONSTRAINTT, CHECKCONINT, INTLJT, INTERPCOSTFUNCTION, REMOVEUNCONNECTEDT, ATOMMATCHDIST, &102:      &        INTCONSTRAINTT, CHECKCONINT, INTLJT, INTERPCOSTFUNCTION, REMOVEUNCONNECTEDT, ATOMMATCHDIST, &
103:      &        DBPT, DBPTDT, DMBLPYT, EFIELDT, MSSTOCKT, NTIPT, PAHAT, PAPT, PATCHYDT, STOCKAAT, RBAAT, RBSYMT, TRAPT, SILANET, &103:      &        DBPT, DBPTDT, DMBLPYT, EFIELDT, MSSTOCKT, NTIPT, PAHAT, PAPT, PATCHYDT, STOCKAAT, RBAAT, RBSYMT, TRAPT, SILANET, &


r30462/keywords.f 2016-05-17 11:30:10.715314052 +0100 r30461/keywords.f 2016-05-17 11:30:11.603325700 +0100
 47:       MSSTOCKT = .FALSE. 47:       MSSTOCKT = .FALSE.
 48:       NTIPT    = .FALSE. 48:       NTIPT    = .FALSE.
 49:       PAHAT    = .FALSE. 49:       PAHAT    = .FALSE.
 50:       PAPT     = .FALSE. 50:       PAPT     = .FALSE.
 51:       PATCHYDT = .FALSE. 51:       PATCHYDT = .FALSE.
 52:       SILANET  = .FALSE. 52:       SILANET  = .FALSE.
 53:       STOCKAAT = .FALSE. 53:       STOCKAAT = .FALSE.
 54:       EFIELDT  = .FALSE. 54:       EFIELDT  = .FALSE.
 55:       DEBUG=.FALSE. 55:       DEBUG=.FALSE.
 56:       EXTRACTMINT=.FALSE. 56:       EXTRACTMINT=.FALSE.
 57:       EXTRACTMINFILET=.FALSE. 
 58:       EXTRACTTST=.FALSE. 57:       EXTRACTTST=.FALSE.
 59:       EXTRACTTSFILET=.FALSE. 
 60:       ADDPT=.FALSE. 58:       ADDPT=.FALSE.
 61:       ADDPT2=.FALSE. 59:       ADDPT2=.FALSE.
 62:       ADDPT3=.FALSE. 60:       ADDPT3=.FALSE.
 63:       DUMPGROUPST=.FALSE. 61:       DUMPGROUPST=.FALSE.
 64:       DSCALE=3.0D0 62:       DSCALE=3.0D0
 65:       PSCALE=0.5D0 63:       PSCALE=0.5D0
 66:       CONNECTIONS=0 64:       CONNECTIONS=0
 67:       MAXTSATTEMPTS=10 65:       MAXTSATTEMPTS=10
 68:       ISEED=1 66:       ISEED=1
 69:       NATTEMPT=0 67:       NATTEMPT=0
1050: C  GMIN executable.1048: C  GMIN executable.
1051: C1049: C
1052:       ELSE IF (WORD.EQ.'EXECGMIN') THEN1050:       ELSE IF (WORD.EQ.'EXECGMIN') THEN
1053:          CALL READA(EXECGMIN)1051:          CALL READA(EXECGMIN)
1054: C1052: C
1055: C  Write the coordinates of minimum WHICHMIN to file extractedmin and stop.1053: C  Write the coordinates of minimum WHICHMIN to file extractedmin and stop.
1056: C1054: C
1057:       ELSE IF (WORD.EQ.'EXTRACTMIN') THEN1055:       ELSE IF (WORD.EQ.'EXTRACTMIN') THEN
1058:          CALL READI(WHICHMIN)1056:          CALL READI(WHICHMIN)
1059:          EXTRACTMINT=.TRUE.1057:          EXTRACTMINT=.TRUE.
1060:       ELSE IF (WORD.EQ.'EXTRACTMINFILE') THEN 
1061:          EXTRACTMINT=.TRUE. 
1062:          EXTRACTMINFILET=.TRUE. 
1063: C1058: C
1064: C  Write the coordinates of minimum WHICHTS to file extractedts and stop.1059: C  Write the coordinates of minimum WHICHTS to file extractedts and stop.
1065: C1060: C
1066:       ELSE IF (WORD.EQ.'EXTRACTTS') THEN1061:       ELSE IF (WORD.EQ.'EXTRACTTS') THEN
1067:          EXTRACTTST=.TRUE.1062:          EXTRACTTST=.TRUE.
1068:          CALL READI(WHICHTS)1063:          CALL READI(WHICHTS)
1069:       ELSE IF (WORD.EQ.'EXTRACTTSFILE') THEN 
1070:          EXTRACTTST=.TRUE. 
1071:          EXTRACTTSFILET=.TRUE. 
1072: C1064: C
1073: C Choose connection pairs based on free energy barriers between regrouped minima1065: C Choose connection pairs based on free energy barriers between regrouped minima
1074: C and product minima.1066: C and product minima.
1075: C1067: C
1076:       ELSE IF (WORD.EQ.'FREEPAIRS') THEN1068:       ELSE IF (WORD.EQ.'FREEPAIRS') THEN
1077:          FREEPAIRT=.TRUE.1069:          FREEPAIRT=.TRUE.
1078:          CALL READF(REGROUPFREETHRESH)1070:          CALL READF(REGROUPFREETHRESH)
1079:          CALL READF(FREETHRESH)1071:          CALL READF(FREETHRESH)
1080:          CALL READF(EINC)1072:          CALL READF(EINC)
1081:          IF (EINC.LE.0.0D0) THEN1073:          IF (EINC.LE.0.0D0) THEN


r30462/mergedb.f90 2016-05-17 11:30:10.939316990 +0100 r30461/mergedb.f90 2016-05-17 11:30:12.427336510 +0100
 21: !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 21: !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
 22: ! 22: !
 23: !  This subroutine merges the database information in directory PATHNAME. 23: !  This subroutine merges the database information in directory PATHNAME.
 24: ! 24: !
 25: SUBROUTINE MERGEDB 25: SUBROUTINE MERGEDB
 26: USE COMMONS,ONLY : NATOMS, IYTS, IZTS, UTSDATA, UTS, HORDERMIN, TOPPOINTER, HORDERTS, PLUS, MINUS, GPFOLD, & 26: USE COMMONS,ONLY : NATOMS, IYTS, IZTS, UTSDATA, UTS, HORDERMIN, TOPPOINTER, HORDERTS, PLUS, MINUS, GPFOLD, &
 27:    &              MAXMIN, MAXTS, FVIBTS, EMIN, FVIBMIN, IXMIN, IYMIN, IZMIN, NEGEIG, PATHNAME, ETS, DEBUG, & 27:    &              MAXMIN, MAXTS, FVIBTS, EMIN, FVIBMIN, IXMIN, IYMIN, IZMIN, NEGEIG, PATHNAME, ETS, DEBUG, &
 28:    &              NMIN, UNRST, CHARMMT, IDIFFTOL, EDIFFTOL, UMINDATA, UMIN, NTS, IXTS, NMINA, NMINB, & 28:    &              NMIN, UNRST, CHARMMT, IDIFFTOL, EDIFFTOL, UMINDATA, UMIN, NTS, IXTS, NMINA, NMINB, &
 29:    &              LOCATIONA, LOCATIONB, ANGLEAXIS, PERMDIST, BOXLX, BOXLY, BOXLZ, GEOMDIFFTOL, TWOD, & 29:    &              LOCATIONA, LOCATIONB, ANGLEAXIS, PERMDIST, BOXLX, BOXLY, BOXLZ, GEOMDIFFTOL, TWOD, &
 30:    &              RIGIDBODY, BULKT, ZSYM, PERMISOMER, IMFRQT, CLOSEFILEST, AMHT, PAIRDISTMAX, DIJINITT, & 30:    &              RIGIDBODY, BULKT, ZSYM, PERMISOMER, IMFRQT, CLOSEFILEST, AMHT, PAIRDISTMAX, DIJINITT, &
 31:    &              DIJINITCONTT, ALLTST, ADDPT2, ADDPT3, NOPT, LPERMDIST 31:    &              DIJINITCONTT, ALLTST, ADDPT2, ADDPT3, NOPT
 32: USE UTILS,ONLY : GETUNIT 32: USE UTILS,ONLY : GETUNIT
 33:  33: 
 34: USE PORFUNCS 34: USE PORFUNCS
 35: IMPLICIT NONE 35: IMPLICIT NONE
 36:  36: 
 37: INTEGER J1, J2, ISTAT, NMINOLD, NTSOLD, NMINDB, NDUMMY, J3, LUNIT, J4, LPLUNIT, LPDUNIT, PLUNIT, PDUNIT 37: INTEGER J1, J2, ISTAT, NMINOLD, NTSOLD, NMINDB, NDUMMY, J3, LUNIT, J4, LPLUNIT, LPDUNIT, PLUNIT, PDUNIT
 38: DOUBLE PRECISION LOCALPOINTS(NOPT), NEWEMIN, NEWETS 38: DOUBLE PRECISION LOCALPOINTS(NOPT), NEWEMIN, NEWETS
 39: DOUBLE PRECISION NEWFVIBMIN, NEWFVIBTS, NEWPOINTSMIN(NOPT), NEWNEGEIG, LPLUS(NOPT), LMINUS(NOPT), & 39: DOUBLE PRECISION NEWFVIBMIN, NEWFVIBTS, NEWPOINTSMIN(NOPT), NEWNEGEIG, LPLUS(NOPT), LMINUS(NOPT), &
 40:   &  NEWPOINTSTS(NOPT), NEWIXMIN,  NEWIYMIN, NEWIZMIN, LPAIRDIST(PAIRDISTMAX), & 40:   &  NEWPOINTSTS(NOPT), NEWIXMIN,  NEWIYMIN, NEWIZMIN, LPAIRDIST(PAIRDISTMAX), &
 41:   &  NEWIXTS,  NEWIYTS, NEWIZTS, DISTANCE, DIST2, LOCALPOINTS2(NOPT), RMAT(3,3) 41:   &  NEWIXTS,  NEWIYTS, NEWIZTS, DISTANCE, DIST2, LOCALPOINTS2(NOPT), RMAT(3,3)
151: !151: !
152: !  DO J2=1,NMIN152: !  DO J2=1,NMIN
153:    DO J2=1,NMIN+NEWMIN153:    DO J2=1,NMIN+NEWMIN
154:       IF (ABS(NEWEMIN-EMIN(J2)).LT.EDIFFTOL) THEN154:       IF (ABS(NEWEMIN-EMIN(J2)).LT.EDIFFTOL) THEN
155:          DISTANCE=1.0D100155:          DISTANCE=1.0D100
156:          READ(UMIN,REC=J2) (LOCALPOINTS2(J3),J3=1,NOPT)156:          READ(UMIN,REC=J2) (LOCALPOINTS2(J3),J3=1,NOPT)
157:          CALL MINPERMDIST(LOCALPOINTS,LOCALPOINTS2,NATOMS,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT,TWOD,DISTANCE,DIST2,RIGIDBODY, &157:          CALL MINPERMDIST(LOCALPOINTS,LOCALPOINTS2,NATOMS,DEBUG,BOXLX,BOXLY,BOXLZ,BULKT,TWOD,DISTANCE,DIST2,RIGIDBODY, &
158:   &                       RMAT,.FALSE.)158:   &                       RMAT,.FALSE.)
159:          IF (.NOT.(ADDPT2.OR.ADDPT3)) THEN159:          IF (.NOT.(ADDPT2.OR.ADDPT3)) THEN
160:             IF ((ABS(NEWEMIN-EMIN(J2)).LT.EDIFFTOL).AND.(DISTANCE.GT.GEOMDIFFTOL)) THEN160:             IF ((ABS(NEWEMIN-EMIN(J2)).LT.EDIFFTOL).AND.(DISTANCE.GT.GEOMDIFFTOL)) THEN
161:                IF (PERMDIST.OR.LPERMDIST) THEN161:                PRINT '(A)',' likely error?'
162:                   PRINT '(A)',' likely error?'162:                PRINT '(A,3G20.10)',' NEWEMIN,EMIN(J2),DISTANCE=',NEWEMIN,EMIN(J2),DISTANCE
163:                   PRINT '(A,3G20.10)',' NEWEMIN,EMIN(J2),DISTANCE=',NEWEMIN,EMIN(J2),DISTANCE163:                PRINT '(A,2I10)',' J2,MINDB=',J2,MINDB
164:                   PRINT '(A,2I10)',' J2,MINDB=',J2,MINDB164:                IF(DEBUG) THEN
165:                   IF(DEBUG) THEN165:                    PRINT '(A)','LOCALPOINTS:'
166:                       PRINT '(A)','LOCALPOINTS:'166:                    PRINT '(3F20.10)',LOCALPOINTS(1:NOPT)
167:                       PRINT '(3F20.10)',LOCALPOINTS(1:NOPT)167:                    PRINT '(A)','LOCALPOINTS2:'
168:                       PRINT '(A)','LOCALPOINTS2:'168:                    PRINT '(3F20.10)',LOCALPOINTS2(1:NOPT)
169:                       PRINT '(3F20.10)',LOCALPOINTS2(1:NOPT) 
170:                   ENDIF 
171: !                 STOP !!! DJW 
172:                ENDIF169:                ENDIF
173: !170: !
174: ! csw34> Added CHARMM structure dumping for debugging purposes171: ! csw34> Added CHARMM structure dumping for debugging purposes
175: !172: !
176:                IF (CHARMMT) THEN173:                IF (CHARMMT) THEN
177:                   PRINT '(A)','Dumping both sets of coordinates to .crd files'174:                   PRINT '(A)','Dumping both sets of coordinates to .crd files'
178:                   CALL CHARMMDUMP(LOCALPOINTS,'localpoints.crd')175:                   CALL CHARMMDUMP(LOCALPOINTS,'localpoints.crd')
179:                   CALL CHARMMDUMP(LOCALPOINTS2,'localpoints2.crd')176:                   CALL CHARMMDUMP(LOCALPOINTS2,'localpoints2.crd')
180:                ENDIF 177:                ENDIF 
 178: !           STOP !!! DJW
181:             ENDIF179:             ENDIF
182:          ENDIF180:          ENDIF
183:          181:          
184:          IF (DISTANCE.LT.GEOMDIFFTOL) THEN182:          IF (DISTANCE.LT.GEOMDIFFTOL) THEN
185:             IF (DEBUG) PRINT '(2(A,I10))','mergedb> minimum ',MINDB,' is database minimum ',J2183:             IF (DEBUG) PRINT '(2(A,I10))','mergedb> minimum ',MINDB,' is database minimum ',J2
186:             IF (ABS(NEWFVIBMIN-FVIBMIN(J2))/FVIBMIN(J2).GT.1.0D-3) THEN184:             IF (ABS(NEWFVIBMIN-FVIBMIN(J2))/FVIBMIN(J2).GT.1.0D-3) THEN
187:                WRITE(*,'(A,F15.5,A,F15.5)') 'mergedb> WARNING, NEWFVIBMIN=',NEWFVIBMIN,' should be ',FVIBMIN(J2)185:                WRITE(*,'(A,F15.5,A,F15.5)') 'mergedb> WARNING, NEWFVIBMIN=',NEWFVIBMIN,' should be ',FVIBMIN(J2)
188:             ENDIF186:             ENDIF
189:             IF (NEWHORDERMIN.NE.HORDERMIN(J2)) THEN187:             IF (NEWHORDERMIN.NE.HORDERMIN(J2)) THEN
190:                WRITE(*,'(A,I6,A,I6)') 'mergedb> ERROR, NEWHORDERMIN=',NEWHORDERMIN,' should be ',HORDERMIN(J2)188:                WRITE(*,'(A,I6,A,I6)') 'mergedb> ERROR, NEWHORDERMIN=',NEWHORDERMIN,' should be ',HORDERMIN(J2)


r30462/setup.f 2016-05-17 11:30:11.163319937 +0100 r30461/setup.f 2016-05-17 11:30:12.639339291 +0100
 88: 776      NDUMMY=NDUMMY-1 88: 776      NDUMMY=NDUMMY-1
 89:          WRITE(*,'(A,I10)')'setup> number of minima written=',NDUMMY 89:          WRITE(*,'(A,I10)')'setup> number of minima written=',NDUMMY
 90:          CLOSE(LUNIT) 90:          CLOSE(LUNIT)
 91:          CLOSE(LUNIT2) 91:          CLOSE(LUNIT2)
 92:          CLOSE(1) 92:          CLOSE(1)
 93:          STOP 93:          STOP
 94:       ENDIF 94:       ENDIF
 95:  95: 
 96:       IF (EXTRACTMINT) THEN 96:       IF (EXTRACTMINT) THEN
 97:          OPEN(UNIT=1,FILE='extractedmin', STATUS='UNKNOWN') 97:          OPEN(UNIT=1,FILE='extractedmin', STATUS='UNKNOWN')
 98:          IF (EXTRACTMINFILET) THEN 98:          IF (WHICHMIN.EQ.-123) THEN
 99:             LUNIT=GETUNIT() 
100:             OPEN(LUNIT,FILE='extractmin',STATUS='OLD') 
101:             DO 
102:               READ(LUNIT,*,END=177) WHICHMIN 
103:               READ(UMIN,REC=WHICHMIN) (LOCALPOINTS(J2),J2=1,NOPT) 
104:               IF (MLP3T.OR.MLPB3T) THEN 
105:                   WRITE(1,'(F25.15)') (LOCALPOINTS(J2),J2=1,NOPT) 
106:                ELSE 
107:                   WRITE(1,'(3F25.15)') (LOCALPOINTS(J2),J2=1,NOPT) 
108:                ENDIF 
109:             ENDDO 
110: 177         CLOSE(LUNIT) 
111:          ELSEIF (WHICHMIN.EQ.-123) THEN 
112:             NDUMMY=0 99:             NDUMMY=0
113:             PRINT '(A)', 'setup> rewriting binary points.min file from extractedmin file'100:             PRINT '(A)', 'setup> rewriting binary points.min file from extractedmin file'
114:             LUNIT=GETUNIT()101:             LUNIT=GETUNIT()
115:             OPEN(UNIT=LUNIT,FILE='points.min.new',ACCESS='DIRECT',FORM='UNFORMATTED',STATUS='UNKNOWN',RECL=8*NOPT) 102:             OPEN(UNIT=LUNIT,FILE='points.min.new',ACCESS='DIRECT',FORM='UNFORMATTED',STATUS='UNKNOWN',RECL=8*NOPT) 
116:             DO 103:             DO 
117:                NDUMMY=NDUMMY+1104:                NDUMMY=NDUMMY+1
118:                READ(1,*,END=777) (LOCALPOINTS(J2),J2=1,NOPT)105:                READ(1,*,END=777) (LOCALPOINTS(J2),J2=1,NOPT)
119:                WRITE(LUNIT,REC=NDUMMY) (LOCALPOINTS(J2),J2=1,NOPT)106:                WRITE(LUNIT,REC=NDUMMY) (LOCALPOINTS(J2),J2=1,NOPT)
120:             ENDDO107:             ENDDO
121: 777         NDUMMY=NDUMMY-1108: 777         NDUMMY=NDUMMY-1
200: 187: 
201:            ENDIF188:            ENDIF
202: 189: 
203:          ENDIF190:          ENDIF
204:          CLOSE(1)191:          CLOSE(1)
205:          STOP192:          STOP
206:       ENDIF193:       ENDIF
207: 194: 
208:       IF (EXTRACTTST) THEN195:       IF (EXTRACTTST) THEN
209:          OPEN(UNIT=1,FILE='extractedts', STATUS='UNKNOWN')196:          OPEN(UNIT=1,FILE='extractedts', STATUS='UNKNOWN')
210:          IF (EXTRACTTSFILET) THEN197:          IF (WHICHTS.EQ.-123) THEN
211:             LUNIT=GETUNIT() 
212:             OPEN(LUNIT,FILE='extractts',STATUS='OLD') 
213:             DO 
214:               READ(LUNIT,*,END=166) WHICHTS 
215:               READ(UTS,REC=WHICHTS) (LOCALPOINTS(J2),J2=1,NOPT) 
216:               IF (MLP3T.OR.MLPB3T) THEN 
217:                   WRITE(1,'(F25.15)') (LOCALPOINTS(J2),J2=1,NOPT) 
218:                ELSE 
219:                   WRITE(1,'(3F25.15)') (LOCALPOINTS(J2),J2=1,NOPT) 
220:                ENDIF 
221:             ENDDO 
222: 166         CLOSE(LUNIT) 
223:  
224:          ELSEIF (WHICHTS.EQ.-123) THEN 
225:             NDUMMY=0198:             NDUMMY=0
226:             PRINT '(A)', 'setup> rewriting binary points.ts file from extractedts file'199:             PRINT '(A)', 'setup> rewriting binary points.ts file from extractedts file'
227:             LUNIT=GETUNIT()200:             LUNIT=GETUNIT()
228:             OPEN(UNIT=LUNIT,FILE='points.ts.new',ACCESS='DIRECT',FORM='UNFORMATTED',STATUS='UNKNOWN',RECL=8*NOPT) 201:             OPEN(UNIT=LUNIT,FILE='points.ts.new',ACCESS='DIRECT',FORM='UNFORMATTED',STATUS='UNKNOWN',RECL=8*NOPT) 
229:             DO202:             DO
230:                NDUMMY=NDUMMY+1203:                NDUMMY=NDUMMY+1
231:                READ(1,*,END=778) (LOCALPOINTS(J2),J2=1,NOPT)204:                READ(1,*,END=778) (LOCALPOINTS(J2),J2=1,NOPT)
232:                WRITE(LUNIT,REC=NDUMMY) (LOCALPOINTS(J2),J2=1,NOPT)205:                WRITE(LUNIT,REC=NDUMMY) (LOCALPOINTS(J2),J2=1,NOPT)
233:             ENDDO206:             ENDDO
234: 778         NDUMMY=NDUMMY-1207: 778         NDUMMY=NDUMMY-1


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0