[Enh] enable selcting axis during variogram fitting.

Fitting a directional variogram currently only works when provided with all directions with the same bins.

We could add an `axis` argument to the fitting routine to only fit the variogram along the specified axis. This would also help with spatio-temporal latlon data when fitting along the time axis.