21 |
|
|
22 |
%load paths and grid: |
%load paths and grid: |
23 |
%-------------------- |
%-------------------- |
24 |
global mygrid; if isempty(mygrid); grid_load('GRID/',5,'compact'); end; |
gcmfaces_global; if isempty(mygrid); grid_load; end; |
25 |
|
|
26 |
%set more optional paramaters to default values |
%set more optional paramaters to default values |
27 |
global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
global myparms; year0=myparms.yearInAve(1); year1=myparms.yearInAve(2); |
95 |
|
|
96 |
if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end; |
if isempty(vecCaxis); cc=3*sqrt(nanmean(misfit_map.^2)); else; cc=vecCaxis(vv,kkk); end; |
97 |
subplot(nk,length(listVar),vv+(kkk-1)*length(listVar)); |
subplot(nk,length(listVar),vv+(kkk-1)*length(listVar)); |
98 |
m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc}); |
m_map_gcmfaces(misfit_map,1,{'myCaxis',[-1 1]*cc},{'myShading','flat'}); |
99 |
title(sprintf('%s misfits at %dm (%d to %d)',varCur,dep,year0,year1)); |
title(sprintf('%s misfits at %dm (%d to %d)',varCur,dep,year0,year1)); |
100 |
end; |
end; |
101 |
end; |
end; |