anchor.diff:        if ( C .LT. 0.) c = c+cm
         RNDM = RNDM - c
         if ( RNDM .LT. 0. ) RNDM = RNDM+1.
  -C
  +!
         return
         END
  -C.**************************************************************************
  -C
  -C.**************************************************************************
  +!.**************************************************************************
  +!
  +!.**************************************************************************
         subroutine RNDMST(NA1,NA2,NA3,NB1)
  -      REAL*4 u,c,cd,cm
  +      REAL u,c,cd,cm
         COMMON/RANDOM/ u(97),c,cd,cm,i,j
  -C
  +!
         ma1 = na1
         ma2 = na2
         ma3 = na3
         mb1 = nb1
         i   = 97
         j   = 33
  -C
  +!
         do ii2 = 1,97
           s = 0.0
           t = 0.5
  
positive.code:       A1= SQRT(-2.0*ALOG(Y))
        A = A1*SIN(X)
        B = A1*COS(X)
  C
        RETURN
        END
  C
        real*4 function RNDM(rdummy)
        REAL*4          u,c,cd,cm, rrdummy
        COMMON/RANDOM/  u(97),c,cd,cm,i,j
  C
        rrdummy = rdummy
  C
        RNDM = u(i)-u(j)
        if (RNDM .LT. 0.) RNDM = RNDM+1.
        U(i) = RNDM
  
negative.code:       ichoix=i
  
  
  !-------- Do the interaction ----------
  
  !      ichoix=3
  
        if (ichoix.eq.1) then
  
  