anchor.diff:      llims(3,2) = minval(lngm, mask)
       ulims(3,2) = maxval(lngm, mask)
   
  -    !$omp section
       mask(:,:) = angdist(latm, lngm, lat2, lng1) .le. (dmax / radearth)
       call growmask(mask)
       llims(4,1) = minval(latm, mask)
       ulims(4,1) = maxval(latm, mask)
       llims(4,2) = minval(lngm, mask)
       ulims(4,2) = maxval(lngm, mask)
  -
  -    !$omp end parallel sections
   
       llat = minval(llims(:,1))
       llng = minval(llims(:,2))
  
positive.code:     llims(1,2) = minval(lngm, mask)
      ulims(1,2) = maxval(lngm, mask)
  
      !$omp section
      mask(:,:) = angdist(latm, lngm, lat1, lng2) .le. (dmax / radearth)
      call growmask(mask)
      llims(2,1) = minval(latm, mask)
  
negative.code:       lngm(i,j) = (lng1 + lng2) / 2 + real(j - 1, dp) / p - 180
      end forall
  
      mask(:,:) = angdist(latm, lngm, lat1, lng1) .le. (dmax / radearth)
      call growmask(mask)
      llims(1,1) = minval(latm, mask)
  