&fortin

  rho0 = 1.0

  ! p0 = rho0 u_phi**2/(gamma * M**2) - 1/2  -- specify M to get p0

  ! M ~ 0.001
  p0 = 599999.5

  ! M ~ 0.01
  ! p0 = 5999.5  

  ! M ~ 0.1
  ! p0 = 59.5

  ! t_r = 1.0


/

&tagging

  denerr = 3
  dengrad = 0.01
  max_denerr_lev = 3
  max_dengrad_lev = 3

  presserr = 3
  pressgrad = 0.01
  max_presserr_lev = 3
  max_pressgrad_lev = 3

/

&extern

  eos_assume_neutral = T

/
