Zhang_Rho_Uncertainty Elapsed Time: 75.4 s
/tcenas/home/gossn/.conda/envs/hypercp/lib/python3.11/site-packages/comet_maths/linear_algebra/matrix_calculation.py:274: UserWarning: One of the correlation matrices is not positive definite. It has been slightly changed (maximum difference of 6.661338147750939e-16) to accomodate our method.
warnings.warn(
Traceback (most recent call last):
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Main.py", line 641, in multi2Clicked
self.processMulti()
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Main.py", line 634, in processMulti
Controller.processFilesMultiLevel(
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/Controller.py", line 852, in processFilesMultiLevel
Controller.processSingleLevel(pathOut, fp, calibrationMap, 'L2')
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/Controller.py", line 772, in processSingleLevel
root = Controller.processL2(root,outFilePath)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/Controller.py", line 470, in processL2
node = ProcessL2.processL2(root,station)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/ProcessL2.py", line 1877, in processL2
if not ProcessL2.stationsEnsemblesReflectance(node, root,station):
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/ProcessL2.py", line 1769, in stationsEnsemblesReflectance
if not ProcessL2.ensemblesReflectance(node, sasGroup, referenceGroup, ancGroup,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/ProcessL2.py", line 1401, in ensemblesReflectance
x_unc, l2_bd = sensor.ClassBasedL2(node, uncGroup, rho_scalar, rho_vec, rho_unc, wavelengths.tolist(),
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/BaseInstrument.py", line 272, in ClassBasedL2
rho = utils.interp_common_wvls(np.array(list(rhoVec.values()), dtype=float),
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/utils.py", line 62, in interp_common_wvls
new_y = np.interp(newWaveBands, x, y) # InterpolatedUnivariateSpline(x, y, k=3)(newWavebands)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/home/gossn/.conda/envs/hypercp/lib/python3.11/site-packages/numpy/lib/function_base.py", line 1599, in interp
return interp_func(x, xp, fp, left, right)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
ValueError: fp and xp are not of the same length.