compare.models function

Compare gravity models