В ответ на:

---

Несколько раз вызывать эту функцию не вариант?

---

code:

    subplot(1,2,1);axes('View',[0 90])
    subplot(1,2,1);hold on
    subplot(1,2,1);mesh(psi(:,1:n_green),lmbd(:,1:n_green),S(:,1:n_green));colormap(CS_green)
    subplot(1,2,1);mesh(psi(:,n_green:n_red),lmbd(:,n_green:n_red),S(:,n_green:n_red));colormap(CS_red)
  

В ответ на:

---

Если нет, то можно просто задать разные карты высот для разных диапазонов и соответственно сместить значения в необходимой области.

---

code:
function h = multimap(X1, Y1, Z1, M1, X2, Y2, Z2, M2)
%MULTIMAP: make two meshes with different colormaps 
%
% X1, Y1, Z1, M1: the first mesh data and colormap
% X2, Y2, Z2, M2: the second mesh data and colormap
%
% Usage 
% [X1,Y1] = meshgrid(-3:.125:3);
% Z1 = peaks(X1,Y1) + 10;
% M1 = hot(64);
% [X2,Y2] = meshgrid(-3:.125:3);
% Z2 = - peaks(X2,Y2) - 10;
% M2 = cool(64);
% h = multimap(X1, Y1, M1, X2, Y2, M2)
%% plot separately for comparison
% figure
% mesh(X1, Y1, Z1); colormap(M1);
% figure
% mesh(X2, Y2, Z2); colormap(M2);


cindex1 = getci(Z1, M1);
cindex2 = getci(Z2, M2) + size(M1, 1);

h1 = mesh(X1,Y1,Z1, cindex1);
set(h1, 'CDataMapping','direct');
hold on
h2 = mesh(X2,Y2,Z2, cindex2);
set(h2, 'CDataMapping','direct');
colormap( [M1; M2]);
h = [h1 h2];

function cindex = getci(Z, M)
%GETCI: helper function to get color index
% 
cmin=min(Z(:));
cmax=max(Z(:));
m = size(M, 1);
cindex = fix((Z-cmin)/(cmax-cmin)*m);