% 2 - Plot grid
% Make some plots showing the 16 fields in absolute and in relative
% magnitude
% csg=load_cs_grid(510,510,510,'t510');
% csg=make_cs_segments_and_patches(csg);
% bathy=load_cs_bathy(csg,'t510','bathy.bin','float32');
csg=load_cs_grid(32,32,32,'t32');
csg=make_cs_segments_and_patches(csg);
bathy=load_cs_bathy(csg,'t32','bathy.bin','float32');

gridarr=csg.gridarr;
 xcpos=1;   ycpos=2;  dxfpos=3;  dyfpos=4;
  rapos=5;   xgpos=6;   ygpos=7;  dxvpos=8;  
 dyupos=9; razpos=10; dxcpos=11; dycpos=12; 
rawpos=13; raspos=14; dxgpos=15; dygpos=16;
clf

subplot(4,4,1);
imagesc(flipud(gridarr(:,:,xcpos)'));colorbar;title('XC (lon)');axis equal

subplot(4,4,2);
imagesc(flipud(gridarr(:,:,ycpos)'));colorbar;title('YC (lat)');axis equal

subplot(4,4,3);
phi=gridarr(:,:,dxfpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DXF (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,4);
phi=gridarr(:,:,dyfpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DYF (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,5);
phi=gridarr(:,:,rapos)./1e6;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('RA (km^2, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,6);
imagesc(flipud(gridarr(:,:,xgpos)'));colorbar;title('XG (lon)');axis equal

subplot(4,4,7);
imagesc(flipud(gridarr(:,:,ygpos)'));colorbar;title('YG (lat)');axis equal

subplot(4,4,8);
phi=gridarr(:,:,dxvpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DXV (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,9);
phi=gridarr(:,:,dyupos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DYU (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,10);
phi=gridarr(:,:,razpos)./1e6;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('RAZ (km^2, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,11);
phi=gridarr(:,:,dxcpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DXC (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,12);
phi=gridarr(:,:,dycpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DYC (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,13);
phi=gridarr(:,:,rawpos)./1e6;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('RAW (km^2, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,14);
phi=gridarr(:,:,raspos)./1e6;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('RAS (km^2, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,15);
phi=gridarr(:,:,dxgpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DXG (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;

subplot(4,4,16);
phi=gridarr(:,:,dygpos)./1e3;
phimin=min(phi(find(phi~=NaN & phi~=0)));
phimax=max(phi(find(phi~=NaN & phi~=0)));
ttl=sprintf('DYG (km, min=%0.5g, max=%0.5g)',phimin,phimax);
imagesc(flipud(phi'));title(ttl,'FontSize',8);axis equal;caxis([phimin phimax]);colorbar;