anchor.diff:          CALL APRBLM2(WORK(KFI),NTOOBS,NTOOBS,NSMOB,1)
         END IF
         CALL FA_FROM_INIINT
  -     &(WORK(KFA),WORK(KMOMO),WORK(KMOMO),WORK(KRHO1),1)
  +     &(WORK(KFA),WORK(KMOMO),WORK(KMOMO),dbl_mb(KRHO1),1)
         KINT2_A(IE2ARR_F)  = KINT2_FSAVE
   *
         ONE = 1.0D0
  
positive.code:              CALL COPVCD(LUSC1,LUHC,dbl_mb(KVEC1),1,LBLK)
            END IF
            IF(IDENSI.GE.1)
       &    CALL DENSI2(IDENSI,WORK(KRHO1),WORK(KRHO2),
       &          dbl_mb(KVEC1),dbl_mb(KVEC2),LUSC1,LUHC,EXPS2,
       &          ISPNDEN,WORK(KSRHO1),WORK(KRHO2AA),WORK(KRHO2AB),
       &          WORK(KRHO2BB),1 )
  
negative.code:       CALL APRBLM2(WORK(KFI),NTOOBS,NTOOBS,NMSOB,1)
        END IF
        CALL FA_FROM_INIINT
       &(WORK(KFA),WORK(KMOMO),WORK(KMOMO),dbl_mb(KRHO1),1)
        KINT2_A(IE2ARR_F)  = KINT2_FSAVE
  *
        ONE = 1.0D0
  