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