hdiff output

r32889/genrigid.f90 2017-06-30 16:30:20.066749540 +0100 r32888/genrigid.f90 2017-06-30 16:30:20.362753487 +0100
379:   CLOSE(222)379:   CLOSE(222)
380: 380: 
381: !  vr274> Calling function for allocation to make more general (setup rigid bodies from code)381: !  vr274> Calling function for allocation to make more general (setup rigid bodies from code)
382:   CALL GENRIGID_ALLOCATE(NRIGIDBODY,MAXSITE)382:   CALL GENRIGID_ALLOCATE(NRIGIDBODY,MAXSITE)
383: 383: 
384: ! hk286 > initialise SITESRIGIDBODY, RIGIDGROUPS, RIGIDSINGLES 384: ! hk286 > initialise SITESRIGIDBODY, RIGIDGROUPS, RIGIDSINGLES 
385:   OPEN(UNIT=222,FILE='rbodyconfig',status='unknown')385:   OPEN(UNIT=222,FILE='rbodyconfig',status='unknown')
386:   DO J1 = 1, NRIGIDBODY386:   DO J1 = 1, NRIGIDBODY
387:      READ(222,*) CHECK1, NSITEPERBODY(J1)387:      READ(222,*) CHECK1, NSITEPERBODY(J1)
388:      DO J2 = 1, NSITEPERBODY(J1)388:      DO J2 = 1, NSITEPERBODY(J1)
 389:         write(*,*) J1, J2, NSITEPERBODY(J1), CHECK1
389:         READ(222,*) RIGIDGROUPS(J2,J1)390:         READ(222,*) RIGIDGROUPS(J2,J1)
390: ! csw34> check to make sure the current atom is not already in a rigid body391: ! csw34> check to make sure the current atom is not already in a rigid body
391:         IF (RIGIDISRIGID(RIGIDGROUPS(J2,J1))) THEN392:         IF (RIGIDISRIGID(RIGIDGROUPS(J2,J1))) THEN
392:             PRINT *," genrigid> ERROR: atom ",RIGIDGROUPS(J2,J1)," is in multiple rigid bodies! Stopping."393:             PRINT *," genrigid> ERROR: atom ",RIGIDGROUPS(J2,J1)," is in multiple rigid bodies! Stopping."
393:             STOP394:             STOP
394:         ELSE       395:         ELSE       
395: ! csw34> if not, flag the current atom396: ! csw34> if not, flag the current atom
396:             RIGIDISRIGID(RIGIDGROUPS(J2,J1))=.TRUE.397:             RIGIDISRIGID(RIGIDGROUPS(J2,J1))=.TRUE.
397:             RB_BY_ATOM(RIGIDGROUPS(J2,J1)) = J1398:             RB_BY_ATOM(RIGIDGROUPS(J2,J1)) = J1
398:         ENDIF399:         ENDIF


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0