anchor.diff:  *. open orbitals
         INTEGER KZ_PTDT(*), KREO_PTDT(*)
   *. The work array containing used for WORK(KZ_PTDET()),WORK(KREO_PTDT())
  -      DIMENSION WORK(*)
  +c     DIMENSION WORK(*)
   *. Specific input
         DIMENSION IBLOCK(8,NBLOCK)
   *. Scratch space 
  
positive.code: *. Address of this spinprojection pattern   
  C  IZNUM_PTDT(IAB,NOPEN,NALPHA,Z,NEWORD,IREORD)
               IPTDT = IZNUM_PTDT(IDET_VC,NOPEN,NOPEN_AL,
       &                WORK(KZ_PTDT(NOPEN+1)),WORK(KREO_PTDT(NOPEN+1)),
       &              1)
               ISIGNP = 1
               IF(IPTDT.EQ.0) THEN
  C?            IF(IRESTR2 .EQ. 1) THEN
  
negative.code:       RETURN
        END
        SUBROUTINE REO_GASDET_S(IREO,NSSOA,NSSOB,NOCTPA,NOCTPB,
       &                 IOCTPA,IOCTPB,NBLOCK,IBLOCK,
       &                 NAEL,NBEL,
       &                 IASTR,IBSTR,IBLTP,NSMST,
       &                 NELFSPGP,NOCCLS,NGAS,
  