anchor.diff:  
     open(99, FILE='ccp.out')
   
  -  write(99, FMT=020) prob, n, np, mi, f, feas, f - PEN * (plim - p), &
  +  write(99, FMT=020) prob, n, np, mi, f, feas, &
  +       f - PEN * max(0.0D0, plim - p) ** 2.0D0, &
          p, fcnt, flag
   
     close(99)
  
positive.code: 
    ! np = INT(n * npfrac)
  
    ! me = 0
  
    ! mi = 10
  
  
negative.code: 
      end if
  
      f = PEN * max(0.0D0, plim - f) ** 2.0D0
  
      ! Quadratic term
  
  