FIXME : Add documentation
pure complex function favg_x (fl, fr) complex, intent (in) :: fl, fr favg_x = 0.5*(fl+fr) end function favg_x