FESOM / fdiag

FESOM2 monitoring
1 stars 1 forks source link

XMOC seemingly does not work anymore #17

Open JanStreffing opened 4 weeks ago

JanStreffing commented 4 weeks ago

Might be the same issue as #16 and #15. All seem to have mesh / data dimension mismatch as basis, even in the final error differs.

xmoc:
    Global Ocean:
        levels: [-20, 20, 41]

results in

fdiag reference_vs_remove_PARMS_part2.yml
/work/ab0246/a270092/software/fdiag/fdiag/templates
The names of runs to be processed are:['reference', 'remove_PARMS']
!!! Performing xmoc !!!
Passed unknown parameter: workflow_name
Passed unknown parameter: climatology_path
Passed unknown parameter: climatology_year
Passed unknown parameter: columns
Passed unknown parameter: bbox
Passed unknown parameter: res
Passed unknown parameter: mapproj
Passed unknown parameter: figsize_small
Passed unknown parameter: start_year
Passed unknown parameter: end_year
Passed unknown parameter: start_year_short
Passed unknown parameter: end_year_short
Passed unknown parameter: years_short
Passed unknown parameter: workflow_settings
Passed unknown parameter: ofolder_notebooks
Passed unknown parameter: ofolder_figures
Executing:  47%|█████████████████████████████████████████████████▉                                                         | 7/15 [01:33<01:46, 13.29s/cell]
Traceback (most recent call last):
  File "/work/ab0246/a270092/software/miniconda3/envs/fdiag/bin/fdiag", line 33, in <module>
    sys.exit(load_entry_point('fdiag', 'console_scripts', 'fdiag')())
             ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/work/ab0246/a270092/software/fdiag/fdiag/fdiag.py", line 242, in fdiag
    webpage = analyses[analysis](settings, analysis)
              ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/work/ab0246/a270092/software/fdiag/fdiag/drivers.py", line 318, in drive_xmoc
    pm.execute_notebook(
  File "/work/ab0246/a270092/software/miniconda3/envs/fdiag/lib/python3.12/site-packages/papermill/execute.py", line 131, in execute_notebook
    raise_for_execution_errors(nb, output_path)
  File "/work/ab0246/a270092/software/miniconda3/envs/fdiag/lib/python3.12/site-packages/papermill/execute.py", line 251, in raise_for_execution_errors
    raise error
papermill.exceptions.PapermillExecutionError:
---------------------------------------------------------------------------
Exception encountered at "In [7]":
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Cell In[7], line 5
      3 for key, value in xmoc.items():
      4     stream_func[key] = {}
----> 5     lats, moc = pf.xmoc_data(mesh, value['data'], mask = mask)
      6     stream_func[key]['data'] = moc
      7     stream_func[key]['lats'] = lats

File /work/ab0246/a270092/software/miniconda3/envs/fdiag/lib/python3.12/site-packages/pyfesom2/diagnostics.py:430, in xmoc_data(mesh, data, nlats, mask, return_masked, meshdiag, el_area, nlevels, face_x, face_y)
    427 pos = ((face_y - lats[0]) / dlat).astype("int")
    429 if isinstance(data, xr.DataArray):
--> 430     w = data[:, :].values * mask[:, None]
    431 else:
    432     w = data[:, :] * mask[:, None]

ValueError: operands could not be broadcast together with shapes (48,126858) (126858,1)