update last commit

Author: Mark Wieczorek

pyshtools/shclasses/shcoeffs.py

@@ -2351,7 +2351,7 @@ def expand(self, grid='DH2', lat=None, colat=None, lon=None, degrees=True,
             if lmax is None:
                 lmax = self.lmax
             if lmax_calc is None:
-                lmax_calc = self.lmax
+                lmax_calc = min(self.lmax, lmax)
             if backend is None:
                 backend = preferred_backend()
             if name is None:
@@ -2432,7 +2432,7 @@ def gradient(self, grid='DH2', lmax=None, lmax_calc=None, units=None,
         if lmax is None:
             lmax = self.lmax
         if lmax_calc is None:
-            lmax_calc = self.lmax
+            lmax_calc = min(self.lmax, lmax)
         if backend is None:
             backend = preferred_backend()
         if name is None:

pyshtools/shclasses/shgravcoeffs.py

@@ -2534,7 +2534,7 @@ def expand(self, a=None, f=None, r=None, colat=None, lat=None, lon=None,
             if lmax is None:
                 lmax = self.lmax
             if lmax_calc is None:
-                lmax_calc = self.lmax
+                lmax_calc = min(self.lmax, lmax)
             if name is None:
                 name = self.name
 
@@ -2647,7 +2647,7 @@ class instance. The components are
         if lmax is None:
             lmax = self.lmax
         if lmax_calc is None:
-            lmax_calc = self.lmax
+            lmax_calc = min(self.lmax, lmax)
         if name is None:
             name = self.name
 
@@ -2741,7 +2741,7 @@ class instance.
         if lmax is None:
             lmax = self.lmax
         if lmax_calc is None:
-            lmax_calc = self.lmax
+            lmax_calc = min(self.lmax, lmax)
         if name is None:
             name = self.name
 

pyshtools/shclasses/shmagcoeffs.py

@@ -2041,7 +2041,7 @@ def expand(self, a=None, f=None, r=None, lat=None, colat=None, lon=None,
             if lmax is None:
                 lmax = self.lmax
             if lmax_calc is None:
-                lmax_calc = self.lmax
+                lmax_calc = min(self.lmax, lmax)
             if name is None:
                 name = self.name
 
@@ -2150,7 +2150,7 @@ class instance. The components are
         if lmax is None:
             lmax = self.lmax
         if lmax_calc is None:
-            lmax_calc = self.lmax
+            lmax_calc = min(self.lmax, lmax)
         if name is None:
             name = self.name