&fortin
  T_l = 1.d9
  u_l = 0.0d0
  rho_l = 1.0d6

  T_r = 1.d6
  u_r = 0.0d0
  rho_r = 1.d6

  use_Tinit = 1

  idir = 1
  frac = 0.5d0
/