matlab_class/sample_analysis/plot_std_field.m


%%%%%%%%%%%%%%%%%
%load parameters:
%%%%%%%%%%%%%%%%%
choiceV3orV4='v4'

if strcmp(choiceV3orV4,'v3')&isempty(whos('mygrid'));
  gcmfaces_path;
  global mygrid; mygrid=[]; clear all;
  dirGrid='/net/altix3700/raid4/gforget/mysetups/ecco_v3/RUNS/fwd3y_fg_as_v4_ncep/';
  nF=1; grid_load(dirGrid,nF); global mygrid;
  dirIn='/net/altix3700/raid4/gforget/mysetups/ecco_v3/RUNS/fwd3y_fg_as_v4_ncep/DIAGfiles/';
elseif strcmp(choiceV3orV4,'v4')&isempty(whos('mygrid'));
  gcmfaces_path;
  global mygrid; mygrid=[]; clear all;
  dirGrid='/net/altix3700/raid5/gforget/mysetups/ecco_v4/RUNS/GRIDev4/';
  nF=5; grid_load(dirGrid,nF); global mygrid;
  dirIn='/net/altix3700/raid4/gforget/mysetups/ecco_v4/RUNS/runs3yVISC/fwd3y_ncep_evolution4/DIAGfiles/';
end;


%%%%%%%%%%%
%get field:
%%%%%%%%%%%
nameFld='DDetan'; tt=[53:78]*336; cc=[0 0.10];
fld=rdmds2gcmfaces([dirIn nameFld],tt);
fld=std(fld,[],3); msk=mygrid.hFacC(:,:,1); fld(find(msk==0))=NaN;
 
%%%%%%%%%%%%
%plot field:
%%%%%%%%%%%%
figure; set(gcf,'Units','Normalized','Position',[0.1 0.3 0.4 0.6]);
[X,Y,FLD]=convert2pcol(mygrid.XC,mygrid.YC,fld); 
pcolor(X,Y,FLD); axis([-180 180 -90 90]); shading flat; caxis(cc); colorbar; 

figure; set(gcf,'Units','Normalized','Position',[0.5 0.3 0.4 0.6]);
FLD=convert2array(fld);
imagescnan(FLD','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); colorbar;

figure; set(gcf,'Units','Normalized','Position',[0.3 0.1 0.4 0.6]);
subplot(3,3,7); imagescnan(fld{1}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
subplot(3,3,8); imagescnan(fld{2}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); 
subplot(3,3,5); imagescnan(fld{3}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); 
subplot(3,3,6); imagescnan(fld{4}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); 
subplot(3,3,3); imagescnan(fld{5}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); 


1	gforget	1.1
2			%%%%%%%%%%%%%%%%%
3			%load parameters:
4			%%%%%%%%%%%%%%%%%
5			choiceV3orV4='v4'
6
7			if strcmp(choiceV3orV4,'v3')&isempty(whos('mygrid'));
8			gcmfaces_path;
9			global mygrid; mygrid=[]; clear all;
10			dirGrid='/net/altix3700/raid4/gforget/mysetups/ecco_v3/RUNS/fwd3y_fg_as_v4_ncep/';
11			nF=1; grid_load(dirGrid,nF); global mygrid;
12			dirIn='/net/altix3700/raid4/gforget/mysetups/ecco_v3/RUNS/fwd3y_fg_as_v4_ncep/DIAGfiles/';
13			elseif strcmp(choiceV3orV4,'v4')&isempty(whos('mygrid'));
14			gcmfaces_path;
15			global mygrid; mygrid=[]; clear all;
16			dirGrid='/net/altix3700/raid5/gforget/mysetups/ecco_v4/RUNS/GRIDev4/';
17			nF=5; grid_load(dirGrid,nF); global mygrid;
18			dirIn='/net/altix3700/raid4/gforget/mysetups/ecco_v4/RUNS/runs3yVISC/fwd3y_ncep_evolution4/DIAGfiles/';
19			end;
20
21
22			%%%%%%%%%%%
23			%get field:
24			%%%%%%%%%%%
25			nameFld='DDetan'; tt=[53:78]*336; cc=[0 0.10];
26			fld=rdmds2gcmfaces([dirIn nameFld],tt);
27			fld=std(fld,[],3); msk=mygrid.hFacC(:,:,1); fld(find(msk==0))=NaN;
28
29			%%%%%%%%%%%%
30			%plot field:
31			%%%%%%%%%%%%
32			figure; set(gcf,'Units','Normalized','Position',[0.1 0.3 0.4 0.6]);
33			[X,Y,FLD]=convert2pcol(mygrid.XC,mygrid.YC,fld);
34			pcolor(X,Y,FLD); axis([-180 180 -90 90]); shading flat; caxis(cc); colorbar;
35
36			figure; set(gcf,'Units','Normalized','Position',[0.5 0.3 0.4 0.6]);
37			FLD=convert2array(fld);
38			imagescnan(FLD','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); colorbar;
39
40			figure; set(gcf,'Units','Normalized','Position',[0.3 0.1 0.4 0.6]);
41			subplot(3,3,7); imagescnan(fld{1}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
42			subplot(3,3,8); imagescnan(fld{2}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
43			subplot(3,3,5); imagescnan(fld{3}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
44			subplot(3,3,6); imagescnan(fld{4}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
45			subplot(3,3,3); imagescnan(fld{5}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
46
47