   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
