Hey @Ara
The issue you are facing arises because you are using two datasets (ROTI and S4) but applying only one colormap at a time and overwriting it during plotting. This causes confusion in the colorbar because the last “colormap()” and “colorbar()” command wins, and both datasets might appear as though they share the same scale and colors, even though they represent different phenomena.
You need to create two separate colorbars, each corresponding only to its dataset. Right now, both datasets are plotted on the same axes (ax), so they share a single colormap context. To separate them visually kindly refer to the following steps:
- Keep your main plotting on a single axes (ax).
- Use invisible dummy axes to separately handle each colormap + colorbar:
ax_roti for ROTI using jet
ax_s4 for S4 using parula
- Plot ROTI and S4 as usual, but apply colormaps only to their respective dummy axes.
Modified Code Snippet:
Replace your ROTI + S4 plotting and colorbar code with this:
% ----------- ROTI Plot -----------
scatterm(lat_ROTI, lon_ROTI, 30, ROTI, 'filled');
% Create dummy axes for ROTI colormap and colorbar
ax_roti = axes('Position', ax.Position, 'Color', 'none', 'Visible', 'off');
colormap(ax_roti, jet);
cb_roti = colorbar(ax_roti, 'southoutside');
cb_roti.Label.String = 'ROTI';
cb_roti.Position = [0.35 0.08 0.3 0.02]; % Adjust as needed
caxis(ax_roti, [0 0.5]);
% ----------- S4 Plot -----------
for station_num = 1:29
if isfield(station_data(station_num), 'data')
…
…
…
end
end
end
% Create dummy axes for S4 colormap and colorbar
ax_s4 = axes('Position', ax.Position, 'Color', 'none', 'Visible', 'off');
colormap(ax_s4, parula);
cb_s4 = colorbar(ax_s4, 'eastoutside');
cb_s4.Label.String = 'S4';
cb_s4.Position = [0.88 0.25 0.02 0.5]; % Adjust as needed to avoid overlap
caxis(ax_s4, [0 0.3]);
For more information on “colormap” and “colorbar”, kindly refer the following MATLAB documentations:
