Add psfc -> psml routine and lonlat output routine