clcl;

sFile = 'DentonEmic17_ds02';

if( exist([sFile '.mat'],'file')==2 )
  load([sFile '.mat']);
else
  aa = dlmread([sFile '.csv'],',',16,0);
  save([sFile '.mat'],'aa','-v7.3');
end

%%

t = 90;
ifld = 3;

aac = aa(aa(:,4)==t,:);

q = unique(aac(:,1));
nq = numel(q);
r = unique(aac(:,3));
nr = numel(r);
f = reshape(aac(:,6+ifld),[nq,nr]);

maxAbsVal = maxAll( abs(f) );  %q0
zcolorM = [ -maxAbsVal ; 0 ; maxAbsVal ];
colorsM = [ 0 0 .35 ; 1 1 1 ;  .35 0 0 ];
colorM = colormap_colors(zcolorM,colorsM);

imagesc(r,q,f,[-maxAbsVal maxAbsVal]);
colormap(colorM);
axis xy
xlabel('r')
ylabel('q')
title('Should be like Figure 1Cb')
