hires module works well in Ocean_Day conditions. moving development back to core

ancillary.pyx

@@ -10,6 +10,7 @@ cdef extern void snow_mask(char *, unsigned char)
 cdef extern float cithr(int, float, float)
 cdef extern void check_reg_uniformity(int, int, int, int, int, unsigned char, unsigned char *,
                                       float *, float *, unsigned char *, int *, int *, int *, int *)
+#cdef extern int get_b1_thresholds(float *, float *, float *, float *)
 
 import cython
 from cython.view cimport array as cvarray
@@ -176,3 +177,29 @@ def py_check_reg_uniformity(np.ndarray[unsigned char, ndim=2] eco_type, np.ndarr
     return scene_uniformity
 
 
+#@cython.boundscheck(False)
+#@cython.wraparound(False)
+#@cython.initializedcheck(False)
+#def py_get_b1_thresholds(np.ndarray[float, ndim=2] locut, np.ndarray[float, ndim=2] midpt,
+#                         np.ndarray[float, ndim=2] hicut, np.ndarray[float, ndim=2] power):
+#
+#    cdef int i, j
+#
+#    cdef np.ndarray_locut_r = np.zeros((), order='C', dtype=np.float32)
+#
+#    cdef float[:, ::1] locut_mv = locut
+#    cdef float[:, ::1] midpt_mv = midpt
+#    cdef float[:, ::1] hicut_mv = hicut
+#    cdef float[:, ::1] power_mv = power
+#    lines = 0
+#    eles = 0
+#
+#    for i in range(lines):
+#        for j in range(eles):
+#
+#            get_b1_thresholds(&locut[0, 0], &midpt[0, 0], &hicut[0, 0], &power[0, 0])
+#
+#            locut_mv[i][j] = locut
+#            midpt_mv[i][j] = midpt
+#            hicut_mv[i][j] = hicut
+#            power_mv[i][j] = power

hires_module.py

@@ -183,9 +183,9 @@ def main(*,
                     longitude=(['number_of_lines', 'number_of_pixels'], viirs_data.longitude.values),
                     )
                 )
-
-    # out_xr.to_netcdf(f'hires_MVCM.{date}.{time}.nc')
-    out_xr.to_netcdf(out_fname)
+    out_path = '/ships19/hercules/pveglio/mvcm_viirs_hires/outputs'
+    out_xr.to_netcdf(f'{out_path}/hires_MVCM.{date}.{time}.nc')
+    # out_xr.to_netcdf(out_fname)
 
 #    # np.savez(f'{data_path}/pyMVCM_{date}.{time}',
 #    np.savez(f'{data_path}/outputs/pyMVCM_{date}.{time}',

setup.py

@@ -31,7 +31,7 @@ include_dirs = ['include',
                 numpy.get_include(), ]
 
 library_dirs = ['/opt/hdfeos2/2.20-gcc-8.3/lib',
-                '/opt/netcdf4/4.7.0fcc-8.3/lib']
+                '/opt/netcdf4/4.7.0-gcc-8.3/lib']
 
 extensions = [Extension('ancillary_data', sourcefiles,
                         include_dirs=include_dirs,