DO NN=1,ML%N_RESIDUE(NR)
         IF (ABS(ML%NU_RESIDUE(NN,NR))>TWO_EPSILON_EB) THEN
            IF (ML%ADJUST_H) &
               MATL_COEF_MATRIX(REAC_COUNTER,MATL_MATRIX_POINTER(ML%RESIDUE_MATL_INDEX(NN,NR))) = -ML%NU_RESIDUE(NN,NR)   <--
            CALL INTERPOLATE1D_UNIFORM(0,MATERIAL(ML%RESIDUE_MATL_INDEX(NN,NR))%H,ML%TMP_REF(NR),ANS)
            MATL_COEF_VECTOR(REAC_COUNTER) = MATL_COEF_VECTOR(REAC_COUNTER )+ ML%NU_RESIDUE(NN,NR)*ANS
         ENDIF
      ENDDO