Fix a transposition error in field calculation.

CHANGELOG

 Changelog
 =========
+2021-11-11 v1.1.3
+-----------------
+  * Fix a transpositon error in field calculation 
+
 2021-10-26 v1.1.2
 -----------------
   * Adapt grad functions to scalar input 

pymagglobal/__init__.py

@@ -46,4 +46,4 @@ from pymagglobal.core import local_curve, dipole_series, file2splines, \
     Model
 from pymagglobal import utils
 
-__version__ = '1.1.2'
+__version__ = '1.1.3'

pymagglobal/core.py

@@ -577,6 +577,7 @@ def field(z_at, splines, cov_splines=None, field_type='nez', check=True,
                               field_cov,
                               rot_mats)
 
+    field_cov = field_cov.T
     if field_type == 'dif':
         grad_D = grad_d(*field)
         cov_D = np.einsum('ij, ijk, ik -> i',
@@ -599,4 +600,4 @@ def field(z_at, splines, cov_splines=None, field_type='nez', check=True,
         return np.array(nez2dif(*field)), \
             np.sqrt(np.array([cov_D, cov_I, cov_F]))
 
-    return field, np.sqrt(field_cov).T
+    return field, np.sqrt(field_cov)