/[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.3 by jmc, Tue Sep 30 22:10:26 2014 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 11  Line 14 
14  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;
15  %nItMx=60*ones(1,Nexp);  %nItMx=60*ones(1,Nexp);
16  namLg=namA ; namLg=strrep(namLg,'_','\_');  namLg=namA ; namLg=strrep(namLg,'_','\_');
17    undef=123456.7;
18  %-----------  %-----------
19  %- test if the variable krd is define :  %- test if the variable krd is define :
20  if size(who('krd'),1) > 0,  if size(who('krd'),1) > 0,
# Line 32  if krd > 0, Line 36  if krd > 0,
36    [ntA(n),rList,tim,vv1,listV,kList] = ...    [ntA(n),rList,tim,vv1,listV,kList] = ...
37      read_StD(prefix,namA(n,:),listV);      read_StD(prefix,namA(n,:),listV);
38    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
39      vv1(find(vv1==undef))=NaN;
40  %- set global dims: & load vvA --> vvB  %- set global dims: & load vvA --> vvB
41    nbV=size(listV,2);    nbV=size(listV,2);
42    nrec=nIt; n3d=nk; nReg=nRg;    nrec=nIt; n3d=nk; nReg=nRg;
# Line 42  if krd > 0, Line 47  if krd > 0,
47    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
48      read_StD(prefix,namA(n,:),listV);      read_StD(prefix,namA(n,:),listV);
49    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
50      vv1(find(vv1==undef))=NaN;
51    if (nrec < nIt),    if (nrec < nIt),
52      fprintf('\n');      fprintf('\n');
53      error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);      error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);
# Line 66  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 79  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 110  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 135  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.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22