/[MITgcm]/MITgcm_contrib/jmc_script/plot_StD.m
ViewVC logotype

Diff of /MITgcm_contrib/jmc_script/plot_StD.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by jmc, Thu Mar 5 20:57:26 2015 UTC revision 1.6 by jmc, Wed Aug 23 22:17:40 2017 UTC
# Line 1  Line 1 
1   prefix='dynStD';   prefix='dynStD';
2   prefix='oceStD';  %prefix='oceStD';
3   pCoords=0;   pCoords=0;
4   namA='c06';   namA='c06';
5   Nexp=1; Nc=size(namA,2);   Nexp=1; Nc=size(namA,2);
6     nAvr=1;
7    %- to plot annual mean instead of 10.d aver:
8    %nAvr=36;
9  %--  %--
10    
11  % $Header$  % $Header$
# Line 69  if krd ~= 0, Line 72  if krd ~= 0,
72   ttA=squeeze(tiA(:,2,:));   ttA=squeeze(tiA(:,2,:));
73   ttA=ttA/3600; titT='hrs';  ttA=ttA/24 ; titT='days';   ttA=ttA/3600; titT='hrs';  ttA=ttA/24 ; titT='days';
74   ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';   ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';
75    %- change to plot annual mean:
76     if nAvr > 1 & rem(nrec,nAvr)==0,
77      nrec=nrec/nAvr; ntA=ntA/nAvr;
78      vvA=reshape(   vvA,[n3d nAvr nrec nReg 5 nbV Nexp]);
79      vvA=reshape(mean(vvA,2),[n3d nrec nReg 5 nbV Nexp]);
80      ttA=reshape(   ttA,[nAvr nrec Nexp]);
81      ttA=reshape(mean(ttA,1),[nrec Nexp]);
82     end
83  end  end
84  %=========================================================  %=========================================================
85    
# Line 82  fprintf('Total length: ntA=');fprintf(' Line 93  fprintf('Total length: ntA=');fprintf('
93  for n=1:Nexp,  for n=1:Nexp,
94  fprintf(' exp %i : time(d):%10.2f ->%10.2f \n', n,ttA(1,n),ttA(ntA(n),n) );  fprintf(' exp %i : time(d):%10.2f ->%10.2f \n', n,ttA(1,n),ttA(ntA(n),n) );
95  end;  end;
96  %--  %-
97    
98  list_on=zeros(1,nbV);  list_on=zeros(1,nbV);
99  nbG=9;  nbG=10;
100  nbG=min(nbG,nbV); list_on(1:nbG)=1 ;  nbG=min(nbG,nbV); list_on(1:nbG)=1 ;
101  %list_on(1:6)=[1 1 1 1 1 1];  %list_on(1:6)=[1 1 1 1 1 1];
102  %list_on(5:7)=0;  %list_on(5:7)=0;
# Line 113  dxRed=0; dyRed=0.03; dxB=0.02; dyB=0.9; Line 124  dxRed=0; dyRed=0.03; dxB=0.02; dyB=0.9;
124  xyP(:,2)=xyP(:,2)+0.010;  xyP(:,2)=xyP(:,2)+0.010;
125  xyB(:,2)=xyB(:,2)+0.010;  xyB(:,2)=xyB(:,2)+0.010;
126    
127    fxb=100; fyb=60; fdx=100; fdy=40; fsc=1.;
128    %fyb=-360; fxb=-2600; fdy=60; fsc=1.5;
129    
130  for ng=1:nbV,  for ng=1:nbV,
131  %-------------------  %-------------------
132   yax=[1:nk-1]; if pCoords == 0, yax=-[1:nk-1]; end   yax=[1:nk-1]; if pCoords == 0, yax=-[1:nk-1]; end
# Line 138  for ng=1:nbV, Line 152  for ng=1:nbV,
152    
153   if flag == 1   if flag == 1
154  %--  %--
155    figure(ng); set(ng,'position',[100+100*ng 60+40*ng 500 700]);clf;    figure(ng); set(ng,'position',[fxb+fdx*ng fyb+fdy*ng [500 700]*fsc]);clf;
156      colormap jet
157    if kList(ng) == 1,    if kList(ng) == 1,
158     var=squeeze(vv1(1,:,1,:,:));     var=squeeze(vv1(1,:,1,:,:));
159     dd=squeeze(max(var)-min(var)); av=squeeze(mean(var));     dd=squeeze(max(var)-min(var)); av=squeeze(mean(var));
# Line 157  for ng=1:nbV, Line 172  for ng=1:nbV,
172      axis(AA); grid ;      axis(AA); grid ;
173      if nv == 1, title(['Avr ',titv,'  ',ttmn]); xlabel(titT); end      if nv == 1, title(['Avr ',titv,'  ',ttmn]); xlabel(titT); end
174      if nv == 2, title(['Std-Dev ',titv,'  ',ttav]); end      if nv == 2, title(['Std-Dev ',titv,'  ',ttav]); end
175      if nv == 3, title(['min ',titv,'  ',ttav]); legend(namLg(1:Nexp,:),0); end      if nv == 3, title(['min ',titv,'  ',ttav]);
176                    legend(namLg(1:Nexp,:),'Location','best'); end
177      if nv == 4, title(['Max ',titv,'  ',ttav]); end      if nv == 4, title(['Max ',titv,'  ',ttav]); end
178     %if nv == 2, title(['Del-2 ',titv,'  ',ttav]); end     %if nv == 2, title(['Del-2 ',titv,'  ',ttav]); end
179     end ; %xlabel(titT);     end ; %xlabel(titT);
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22