Fixed a bug when the longitude lies in a small sector between the end and start

src/AVISOLOCAL.cpp

@@ -183,6 +183,10 @@ template<class DataType> AVISOLOCALData::Data AVISOLOCALData::ReadVarRaw(const N
  real lat1 = p->latb;
  real lat2 = p->late;
 
+ // Special case when the longitude lies in a small sector between the end and the start
+ if(lon1 < lonb) lon1 = lone;
+ if(lon2 > lone) lon2 = lonb;
+
  auto yb = static_cast<size_t>(Floor((lat1 - latb) / latstep));
  auto ye = static_cast<size_t>(Ceil((lat2 - latb) / latstep));
  if(ye > dn.ny - 1) ye = dn.ny - 1;

src/layereddata.cpp

@@ -142,6 +142,10 @@ std::pair<const BaseParameters*, MString> LayeredData::Parameters(michlib_intern
   real lat1 = args.at("latb").ToReal();
   real lat2 = args.at("late").ToReal();
 
+  // Special case when the longitude lies in a small sector between the end and the start
+  if(lon1 < lonb) lon1 = lone;
+  if(lon2 > lone) lon2 = lonb;
+
   ppar->yb = static_cast<size_t>(Floor((lat1 - latb) / latstep));
   ppar->ye = static_cast<size_t>(Ceil((lat2 - latb) / latstep));
   if(ppar->ye > dname.ny - 1) ppar->ye = dname.ny - 1;