anchor.diff:  *
         RETURN
         END
  -      SUBROUTINE DIA0TRM(ITASK,LUIN,LUOUT,VECIN,VECOUT,FACTOR)
  -*
  -*
  -* Direct calculation of diagonal terms for perturbation operator
  -*
  -* Itask = 1 : VEC2 = (Diag + Factor )    VEC1
  -* Itask = 2 : VEC2 = (Diag + factor ) -1 VEC1
  -* 
  -* if ICISTR.NE.1 VECOUT is LUOUT, VECIN is LUIN
  +      SUBROUTINE DIABLK(IASM,IATP,IBSM,IBTP,IFULL,DIAG)
  +*
  +* Obtain diagonal block of determinant block 
  +* IASM,IATP,IBSM,IBTP
  +*
  +*
  +*. Type of operator is taken from OPER
  +*. Type of CI space is taken from CANDS
   *
   c      IMPLICIT REAL*8(A-H,O-Z)
   * =====
  
positive.code: CNW   DIMENSION FIFA(*)
        integer FIFA
  *
        INCLUDE 'orbinp.inc'
        INCLUDE 'glbbas.inc'
        INCLUDE 'lucinp.inc'
  
negative.code:       ELSE
          KLSVST = 1
        END IF
        MAXA = IMNMX(int_mb(KNSTSO(IAGP)),NSMST*NOCTPA,2)
        CALL MEMMAN(KLRJKA,MAXA,'ADDL  ',2,'KLRJKA')
  *
  ** Info on block structure of internal state
  