test

jswhit · Jul 2, 2024 · b09798a · b09798a
1 parent 155a0f2
commit b09798a
Showing 1 changed file with 8 additions and 5 deletions.
diff --git a/src/pygrib/_pygrib.pyx b/src/pygrib/_pygrib.pyx
@@ -1321,8 +1321,9 @@ cdef class gribmessage(object):
             else:
                 missval = 1.e30
             if self.expand_reduced:
-                lonsperlat = self['pl'].astype(np.int64)
-                nx = lonsperlat.max()
+                lonsperlat = self['pl']
+                #nx = lonsperlat.max()
+                nx = 2*ny
                 datarr = redtoreg(datarr, lonsperlat, missval=missval)
             else:
                 nx = None
@@ -1559,8 +1560,9 @@ cdef class gribmessage(object):
                 lats = self['distinctLatitudes']
                 if lat2 < lat1 and lats[-1] > lats[0]: lats = lats[::-1]
                 ny = self['Nj']
-                lonsperlat = self['pl'].astype(np.int64)
-                nx = lonsperlat.max()
+                lonsperlat = self['pl']
+                #nx = lonsperlat.max()
+                nx = 2*ny
                 lon1 = self['longitudeOfFirstGridPointInDegrees']
                 lon2 = self['longitudeOfLastGridPointInDegrees']
                 lons = np.linspace(lon1,lon2,nx)
@@ -1572,7 +1574,8 @@ cdef class gribmessage(object):
             if self.expand_reduced:
                 ny = self['Nj']
                 lonsperlat = self['pl'].astype(np.int64)
-                nx = lonsperlat.max()
+                #nx = lonsperlat.max()
+                nx = 2*ny
                 lat1 = self['latitudeOfFirstGridPointInDegrees']
                 lat2 = self['latitudeOfLastGridPointInDegrees']
                 lon1 = self['longitudeOfFirstGridPointInDegrees']