/[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.5 by jmc, Sun Sep 6 18:15:09 2015 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;
128    %fxb=-2600; fyb=160; %fxb=100;
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+100*ng fyb+40*ng 500 700]);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));
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22