def gas_fvf(z, temp, pressure):
  """
  Gas FVF calculated in oilfield unit, result in res ft3/scf
  inputs temp in Fahrenheit, pressure in psia or psig
  """
  temp = temp + 459.67
  Bg = 0.0282793 * z * temp / pressure 
  return(Bg)