   0          ityp ( 0 - uniform trian,  >0 - component of solution)
   1           ndim
   0           ifv ( 0 -cell vertex, 1- cell centered)
 0.02          pos   (=positivity)
  500          numel (= $c$)
  20.         epsilon1
   1.          p
