Fixes to extract_timeseries

* 'format' arg now called 'fmt'
* Now returns passed lat/lon instead of found lat/lon in order to
  ensure that the same number of locations are returned even if
  some of them overlap to the same grid point in the data (needs a
  better solution eventually)