anchor.diff:  * 
         IF(I_DO_COMHAM.EQ.1) THEN
            WRITE(6,*) ' Complete Hamiltonian matrix will be constructed '
  +CBERT VECTOR DOUBLES!!! NOT USED
            CALL MEMMAN(KLHMAT,NVAR*NVAR,'ADDL  ',2,'HMAT  ')
   C COMHAM(H,NVAR,NBLOCK,LBLOCK,VEC1,VEC2)
            ECOREL = 0.0D0
  
positive.code:       END IF
  C
        RETURN
        END
  
negative.code: C?       WRITE(6,*) ' IH0_CSF = ', IH0_CSF
           CALL CSDIAG(WORK(KPVEC2),WORK(KPDIA_SD),
       &        NCONF_PER_OPEN(1,ISM),MAXOP,ISM,
       &        int_mb(KSDREO_I(ISM)),NPCMCNF,NPCSCNF,IPRCSF,
       &        ICNFBAT,NCOCCLS_ACT,dbl_mb(KCIOCCLS_ACT),
       &        int_mb(KCLBT),LUDIA,LUSC52)
  C        CSDIAG(CSFDIA,DETDIA,NCNFTP,MAXOP,ISM,
  C    &         ICTSDT,NPCMCNF,NPCSCNF,IPRCSF,
  C    &         ICNFBAT,NOCCLS_ACT,IOCCLS_ACT,
  