&fortin

  xsep = 0.1d0
  rho0 = 1.0d0
  rho1 = 1.000001d0
  pl   = 1000.d0
  pr   = 0.01d0

/