/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/Breakdown_FRM.py:362: RuntimeWarning: divide by zero encountered in divide
return (LAMP_mag * np.mean(sample_PANEL, axis=0)) / (np.piS110)
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 1785, in stationsEnsemblesReflectance
if not ProcessL2.ensemblesReflectance(node, sasGroup, referenceGroup, ancGroup,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/ProcessL2.py", line 1413, in ensemblesReflectance
x_unc = sensor.FRML2(pds, rho_scalar, rho_vec, rho_unc, wavelengths, x_slice, x_breakdown_unc)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/BaseInstrument.py", line 584, in FRML2
sample_Lw = UNC_Obj_FRM.MCP.run_samples(Propagate.Lw_FRM, [sample_LT, sample_rho, sample_LI])
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/home/gossn/.conda/envs/hypercp/lib/python3.11/site-packages/punpy/mc/mc_propagation.py", line 1484, in run_samples
MC_y = list(map(func, *[x[sli] for x in MC_x]))
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/Uncertainty_Analysis.py", line 540, in Lw_FRM
return lt - (rho * li)
~~~~^~~~
ValueError: operands could not be broadcast together with shapes (189,) (199,)
/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/Breakdown_FRM.py:362: RuntimeWarning: divide by zero encountered in divide
return (LAMP_mag * np.mean(sample_PANEL, axis=0)) / (np.piS110)
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 1785, in stationsEnsemblesReflectance
if not ProcessL2.ensemblesReflectance(node, sasGroup, referenceGroup, ancGroup,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/ProcessL2.py", line 1413, in ensemblesReflectance
x_unc = sensor.FRML2(pds, rho_scalar, rho_vec, rho_unc, wavelengths, x_slice, x_breakdown_unc)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/BaseInstrument.py", line 584, in FRML2
sample_Lw = UNC_Obj_FRM.MCP.run_samples(Propagate.Lw_FRM, [sample_LT, sample_rho, sample_LI])
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/home/gossn/.conda/envs/hypercp/lib/python3.11/site-packages/punpy/mc/mc_propagation.py", line 1484, in run_samples
MC_y = list(map(func, *[x[sli] for x in MC_x]))
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/tcenas/s3optcal/workdata/ocean_colour/HyperCP/HyperCP/Source/PIU/Uncertainty_Analysis.py", line 540, in Lw_FRM
return lt - (rho * li)
~~~~^~~~
ValueError: operands could not be broadcast together with shapes (189,) (199,)