We observed inconsistent behaviour in bottom-related plotting using dfm_tools, where variables plotted near the seabed appear patchy for certain model configurations (e.g. latest DCSM configuration). After debugging and comparing model runs, we found that this behaviour is directly related to the MDU parameter: KeepzLayeringAtBed. The issue occurs when KeepzLayeringAtBed = 1. With KeepzLayeringAtBed = 2, the plotting behaves as expected (smooth fields without patchiness).
Test cases:
- KeepzLayeringAtBed = 2 and bottom plotting works: p:\11212771-004-dcsm-fm-sed-waq\models\dflowfm3d-noordzee_0_5nm-j22_6-v1a\computations\2012-2016_coarse\A02_IM_WQ_solvedInit_nothermo\
- KeepzLayeringAtBed = 1 and bottom plotting becomes patchy: p:\11212771-004-dcsm-fm-sed-waq\models\dflowfm3d-noordzee_0_5nm-j22_6-v1a\computations\2012-2016_coarse\A02_IM_WQ_init\
Minimal code example:
import dfm_tools as dfmt
ds = dfmt.open_partitioned_dataset()
ds_sel = ds.sel(time='2012-12-06', method='nearest')
ds_sel = dfmt.get_Dataset_atdepths(ds_sel, depths=1, reference='bedlevel')
ds_sel = ds_sel['mesh2d_tem1']
ds_sel .ugrid.plot()
Trying a range of depths above the bedlevel all gave us the same issue in the new mdu set-up.
Example figures:
With the first filepath, this creates the figure below (KeepzLayeringAtBed = 2):
With the second filepath, it created the following figure (KeepzLayeringAtBed = 1):
We observed inconsistent behaviour in bottom-related plotting using dfm_tools, where variables plotted near the seabed appear patchy for certain model configurations (e.g. latest DCSM configuration). After debugging and comparing model runs, we found that this behaviour is directly related to the MDU parameter: KeepzLayeringAtBed. The issue occurs when KeepzLayeringAtBed = 1. With KeepzLayeringAtBed = 2, the plotting behaves as expected (smooth fields without patchiness).
Test cases:
Minimal code example:
import dfm_tools as dfmt
ds = dfmt.open_partitioned_dataset()
ds_sel = ds.sel(time='2012-12-06', method='nearest')
ds_sel = dfmt.get_Dataset_atdepths(ds_sel, depths=1, reference='bedlevel')
ds_sel = ds_sel['mesh2d_tem1']
ds_sel .ugrid.plot()
Trying a range of depths above the bedlevel all gave us the same issue in the new mdu set-up.
Example figures:
With the first filepath, this creates the figure below (KeepzLayeringAtBed = 2):
With the second filepath, it created the following figure (KeepzLayeringAtBed = 1):
