subsection Boundary velocity model
  subsection Function
    set Variable names      = x,y,z
    set Function constants  = cm=0.01, year=1
    set Function expression =  if (x<56e3 && y<=48e3 | x<72e3 && y>48e3,-1*cm/year,1*cm/year);0;0
  end
end
