hdiff output

r30700/SBM.f 2016-07-06 15:36:00.887567105 +0100 r30699/SBM.f 2016-07-06 15:36:01.255572081 +0100
 63:       character TMPCHAR 63:       character TMPCHAR
 64:       integer TMPINT,nexc,I1,I2 64:       integer TMPINT,nexc,I1,I2
 65:       double precision TMPREAL,concentration 65:       double precision TMPREAL,concentration
 66:       LOGICAL TARR(NATOMS) 66:       LOGICAL TARR(NATOMS)
 67:       INTEGER EXCLUSIONS,NUMOFEXCLUSIONS 67:       INTEGER EXCLUSIONS,NUMOFEXCLUSIONS
 68:       INTEGER MAXEXCLUSIONS,MAXEXCLUSIONSELEC 68:       INTEGER MAXEXCLUSIONS,MAXEXCLUSIONSELEC
 69:       PARAMETER(MAXEXCLUSIONS=20) 69:       PARAMETER(MAXEXCLUSIONS=20)
 70:       PARAMETER(MAXEXCLUSIONSELEC=5) 70:       PARAMETER(MAXEXCLUSIONSELEC=5)
 71:       DIMENSION EXCLUSIONS(NATOMS*MAXEXCLUSIONS),NUMOFEXCLUSIONS(NATOMS) 71:       DIMENSION EXCLUSIONS(NATOMS*MAXEXCLUSIONS),NUMOFEXCLUSIONS(NATOMS)
 72:       INTEGER IOSTATUS,J1,J2,NRIGIDBODY,ATOMID 72:       INTEGER IOSTATUS,J1,J2,NRIGIDBODY,ATOMID
 73:       CHARACTER(LEN=10) CHECK1 73:       CHARACTER CHECK1
 74: !$    INTEGER NTHREADS,OMP_GET_NUM_THREADS,OMP_GET_THREAD_NUM,TID  74: !$    INTEGER NTHREADS,OMP_GET_NUM_THREADS,OMP_GET_THREAD_NUM,TID 
 75: !$OMP PARALLEL 75: !$OMP PARALLEL
 76: !$    NTHREADS = OMP_GET_NUM_THREADS() 76: !$    NTHREADS = OMP_GET_NUM_THREADS()
 77: !$      TID = OMP_GET_THREAD_NUM() 77: !$      TID = OMP_GET_THREAD_NUM()
 78: !$    if(TID .eq. 0)then  78: !$    if(TID .eq. 0)then 
 79: !$      write(*,*) 'OMP enabled. Number of threads:', NTHREADS 79: !$      write(*,*) 'OMP enabled. Number of threads:', NTHREADS
 80: !$    endif 80: !$    endif
 81: !$OMP END PARALLEL 81: !$OMP END PARALLEL
 82:  82: 
 83:       ALLOCATE(XSBM(NATOMS)) 83:       ALLOCATE(XSBM(NATOMS))
148: 148: 
149:         DEALLOCATE(TMPSTT)149:         DEALLOCATE(TMPSTT)
150: 150: 
151: 151: 
152:         ! CONCENTRATION IS THE MONOVALENT ION CONCENTRATION kappa is in units152:         ! CONCENTRATION IS THE MONOVALENT ION CONCENTRATION kappa is in units
153:         ! A^-1153:         ! A^-1
154:         KAPPA=0.32*sqrt(CONCENTRATION)154:         KAPPA=0.32*sqrt(CONCENTRATION)
155:         PREFACTOR=PREFACTOR/DC155:         PREFACTOR=PREFACTOR/DC
156:         read(30,*) ANtemp156:         read(30,*) ANtemp
157:         write(*,*) ANtemp, 'atoms'157:         write(*,*) ANtemp, 'atoms'
158:         if(NATOMS .NE. ANTEMP)THEN158:         if(NATOMS .EQ. ANTEMP)THEN
159:           write(*,*) 'ERROR: Number of atoms in SBM.INP and odata are not consistent'159:           write(*,*) 'Number of atoms in SBM.INP and odata are not consistent'
160:           write(*,*) 'SBM.INP=',ANTEMP, 'odata',NATOMS 
161:           STOP 
162:         ENDIF160:         ENDIF
163:         NUMOFSBMCHARGES=0161:         NUMOFSBMCHARGES=0
164:         ALLOCATE(SBMCHARGE(ANTEMP))162:         ALLOCATE(SBMCHARGE(ANTEMP))
165:         ALLOCATE(SBMCHARGEON(ANTEMP))163:         ALLOCATE(SBMCHARGEON(ANTEMP))
166:         do i=1, ANtemp164:         do i=1, ANtemp
167:           read(30,*) TMPINT,ATYPE,TMPINT,TMPCHAR,TMPCHAR,SBMCHARGE(i),ATMASS(I)165:           read(30,*) TMPINT,ATYPE,TMPINT,TMPCHAR,TMPCHAR,SBMCHARGE(i),ATMASS(I)
168:             if(ATYPE .gt. NATOMTYPES)THEN166:             if(ATYPE .gt. NATOMTYPES)THEN
169:               write(*,*) 'ERROR: Unknown atomtype',ATYPE167:               write(*,*) 'ERROR: Unknown atomtype',ATYPE
170:             endif168:             endif
171:             ATOMTYPES(I)=ATYPE169:             ATOMTYPES(I)=ATYPE
372:        close(30)370:        close(30)
373: 371: 
374: 372: 
375: ! If RIGIDINIT is on, then we will want to keep information about rigid groups373: ! If RIGIDINIT is on, then we will want to keep information about rigid groups
376: ! so that we may ignore their interactions later.374: ! so that we may ignore their interactions later.
377: 375: 
378:       ALLOCATE (SBMRBG(NATOMS))376:       ALLOCATE (SBMRBG(NATOMS))
379:       DO I=1,NATOMS377:       DO I=1,NATOMS
380:         SBMRBG(I)=-I378:         SBMRBG(I)=-I
381:       ENDDO379:       ENDDO
382:         NRIGIDBODY=0380: 
383: ! approach to I/O for rbodyconfig taken from genrigid.f90381: ! approach to I/O for rbodyconfig taken from genrigid.f90
384:       if(RIGIDINIT)THEN382:       if(RIGIDINIT)THEN
385:         OPEN(UNIT=222,FILE='rbodyconfig',status='old')383:         OPEN(UNIT=222,FILE='rbodyconfig',status='old')
386:         DO384:         DO
387:           READ(222,*,IOSTAT=iostatus) CHECK1385:           READ(222,*,IOSTAT=iostatus) CHECK1
388:           IF (iostatus<0) THEN386:           IF (iostatus<0) THEN
389:             CLOSE(222)387:             CLOSE(222)
390:             EXIT388:             EXIT
391:           ELSE IF (TRIM(ADJUSTL(CHECK1)).EQ.'GROUP') then389:           ELSE IF (TRIM(ADJUSTL(CHECK1)).EQ.'GROUP') then
392:             NRIGIDBODY=NRIGIDBODY+1390:             NRIGIDBODY=NRIGIDBODY+1
393:           ENDIF391:           ENDIF
394:         END DO392:         END DO
395:         CLOSE(222)393:         CLOSE(222)
396:   394:    
397:         write(*,*) NRIGIDBODY, ' rigid bodies will be used' 
398:   
399:         OPEN(UNIT=222,FILE='rbodyconfig',status='old')395:         OPEN(UNIT=222,FILE='rbodyconfig',status='old')
400:         DO J1 = 1, NRIGIDBODY396:         DO J1 = 1, NRIGIDBODY
401:           READ(222,*) CHECK1,DUMMY397:           READ(222,*) CHECK1,DUMMY
402:           DO J2 = 1, DUMMY398:           DO J2 = 1, DUMMY
403:             READ(222,*) ATOMID399:             READ(222,*) ATOMID
404:             SBMRBG(ATOMID)=J1400:             SBMRBG(ATOMID)=J1
405:           ENDDO401:           ENDDO
406:         ENDDO402:         ENDDO
407:         CLOSE(222)403:         CLOSE(222)
408:       ENDIF404:       ENDIF


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0