function []=demo_grids();
% DEMO_GRIDS llustrates the use of the various grids in this approach.

listSol={'ECCOv4','ECCOv3','ECCOv2','ECCO2','OCCA','SOSE'};
listNfc=[5 1 1 6 1 1];
listFrm={'compact','straight','straight','cube','straight','straight'};

for sol=1:length(listSol);
gcmfaces_global; mygrid=[]; 
grid_load([listSol{sol} filesep],listNfc(sol),listFrm{sol});
fld=read_bin([listSol{sol} filesep 'someTHETA.data'],1,0,64);
figure; gcmfaces_sphere(fld,[-2 20],10,'S',1);
gcmfaces_caption(['Some temperture field at 300m depth for: ' sol]);
end;


