real function fahrenheit(c)
     real :: c    ! temperature, in Celsius degrees
     fahrenheit = (9.0/5.0)*c + 32.0  ! Convert Celsius to fahrenheit
end function fahrenheit