matlab_class/sample_analysis/plot_one_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='DDtheta'; tt=[53]*336; kk=30; cc=[-2 8];
fld=rdmds2gcmfaces([dirIn nameFld],tt);
fld=fld(:,:,kk); msk=mygrid.hFacC(:,:,kk); 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			%get field:
23			%%%%%%%%%%%
24			nameFld='DDtheta'; tt=[53]*336; kk=30; cc=[-2 8];
25			fld=rdmds2gcmfaces([dirIn nameFld],tt);
26			fld=fld(:,:,kk); msk=mygrid.hFacC(:,:,kk); fld(find(msk==0))=NaN;
27
28			%%%%%%%%%%%%
29			%plot field:
30			%%%%%%%%%%%%
31			figure; set(gcf,'Units','Normalized','Position',[0.1 0.3 0.4 0.6]);
32			[X,Y,FLD]=convert2pcol(mygrid.XC,mygrid.YC,fld);
33			pcolor(X,Y,FLD); axis([-180 180 -90 90]); shading flat; caxis(cc); colorbar;
34
35			figure; set(gcf,'Units','Normalized','Position',[0.5 0.3 0.4 0.6]);
36			FLD=convert2array(fld);
37			imagescnan(FLD','nancolor',[1 1 1]*0.8); axis xy; caxis(cc); colorbar;
38
39			figure; set(gcf,'Units','Normalized','Position',[0.3 0.1 0.4 0.6]);
40			subplot(3,3,7); imagescnan(fld{1}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
41			subplot(3,3,8); imagescnan(fld{2}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
42			subplot(3,3,5); imagescnan(fld{3}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
43			subplot(3,3,6); imagescnan(fld{4}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
44			subplot(3,3,3); imagescnan(fld{5}','nancolor',[1 1 1]*0.8); axis xy; caxis(cc);
45
46