If the levelsarray in the input yaml file does not include all level values found in the control file then the data load fails when the time data is added to the dataset. The issue is that the several variables associated with levels (level, level_yaxis_z, level_assim, level_nassim) need to all be reduced if the actual levels requested is not all available. This is already done for Channel so it should be a simple addition to the subset_coordinates method in mon_data_space.py.
If the
levels
array in the input yaml file does not include all level values found in the control file then the data load fails when the time data is added to the dataset. The issue is that the several variables associated with levels (level
,level_yaxis_z
,level_assim
,level_nassim
) need to all be reduced if the actual levels requested is not all available. This is already done forChannel
so it should be a simple addition to thesubset_coordinates
method inmon_data_space.py
.