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

Diff of /MITgcm_contrib/jmc_script/grph_StD.m

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

revision 1.1 by jmc, Thu May 29 23:47:13 2008 UTC revision 1.2 by jmc, Wed Oct 8 18:18:04 2008 UTC
# Line 1  Line 1 
1  Nexp=3;   prefix='stDiag';
2    Nc=4; nc=Nc*ones(1,Nexp); prefix='stDiag';   namA={'mG_2','vmUpW'};
3   namA(1,1:Nc)='r06a';  %namA=char{'r525','r526','r527','r528'};
4   namA(2,1:Nc)='r07a';   Nexp=size(namA,2);
  namA(3,1:Nc)='r07b';  
5  %-  %-
6    %namA(2)={'r528'};
7    
8  % $Header$  % $Header$
9  % $Name$  % $Name$
10    
11    
12  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;
13  %nItMx=360*ones(1,Nexp);  %nItMx=60*ones(1,Nexp);
14  namLg=namA ; for n=1:Nexp, namLg(n,find(namLg(n,:)=='_'))='-'; end  namLg=namA ; namLg=strrep(namLg,'_','\_');
15  %-----------  %-----------
16  %- ngEn = Nb of Energy plot: = 2 or = 4  ;  %- ngEn = Nb of Energy plot: = 2 or = 4  ;
17  ngEn=2;  ngEn=2;
# Line 25  if krd > 0, Line 25  if krd > 0,
25  %- define list of fields to read in:  %- define list of fields to read in:
26  %listV={'Eta','U','V','W','T','S','DETADT2','RELHUM','Phi'};  %listV={'Eta','U','V','W','T','S','DETADT2','RELHUM','Phi'};
27  %listV={'Eta','U','V','W','T','S','CONVADJ','DETADT2'};  %listV={'Eta','U','V','W','T','S','CONVADJ','DETADT2'};
28    %listV={'Eta','UE_VEL_C','VN_VEL_C','W','T','DETADT2','Phi'};
29  %- or take all them:  %- or take all them:
30  clear listV ; listV='all_flds';   clear listV ; listV='all_flds';
31  %-----------  %-----------
32    
33  %- start to read the longest record:  %- start to read the longest record:
34    n=1; rf=-1; if strcmp(char(listV),'all_flds'), rf=0; end    n=1; rf=-1; if strcmp(char(listV),'all_flds'), rf=0; end
35    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
36      read_StD(prefix,namA(n,1:nc(n)),listV);      read_StD(prefix,namA(n),listV);
37    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
38  %- set global dims: & load vvA --> vvB  %- set global dims: & load vvA --> vvB
39    nbV=size(listV,2);    nbV=size(listV,2);
# Line 42  clear listV ; listV='all_flds'; Line 43  clear listV ; listV='all_flds';
43  %----  %----
44   for n=2:Nexp,   for n=2:Nexp,
45    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
46      read_StD(prefix,namA(n,1:nc(n)),listV);      read_StD(prefix,namA(n),listV);
47    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
48    if (nrec < nIt),    if (nrec < nIt),
49      fprintf('\n');      fprintf('\n');
# Line 66  elseif krd < 0, Line 67  elseif krd < 0,
67  end  end
68  if krd ~= 0,  if krd ~= 0,
69   ttA=squeeze(tiA(:,2,:));   ttA=squeeze(tiA(:,2,:));
70   ttA=ttA/86400; titT='days'; ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';   ttA=ttA/86400; titT='days'; %ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';
71  end  end
72  %=========================================================  %=========================================================
73    
74  ttax1=0 ; ttax2=0 ; ttay=zeros(nbV,2);  ttax1=0 ; ttax2=0 ; ttay=zeros(nbV,2);
75  %-- fixed time axis bound :  %-- fixed time axis bound :
76  % ttax1=3; ttax2=3.3;  % ttax1=15.; ttax2=20.;
77  %-- fixed Y axis bound :  %-- fixed Y axis bound :
78  % ttay(4,:)=[0 0.6];  % ttay(4,:)=[0 0.6];
79  %-----------  %-----------
# Line 104  linA(6,:)='c-'; Line 105  linA(6,:)='c-';
105  ieA=min(ieA,nItMx);  ieA=min(ieA,nItMx);
106  %titall='AIM , Cubic-G (32x32) , cpl-FM Forcing' ;  %titall='AIM , Cubic-G (32x32) , cpl-FM Forcing' ;
107  titall='Global Ocean, Cubic-G (32x32) , NCEP Forc (2)' ;  titall='Global Ocean, Cubic-G (32x32) , NCEP Forc (2)' ;
108    %titall='Dyncore test-case 5 (cs-32)' ;
109    
110  %=========================================================  %=========================================================
111    
# Line 134  for ng=1:nbV, Line 136  for ng=1:nbV,
136        ttmn=sprintf([ttmn,' %2.1e ;'],dd(nv,n));        ttmn=sprintf([ttmn,' %2.1e ;'],dd(nv,n));
137        ttav=sprintf([ttav,' %3.2e ;'],av(nv,n));        ttav=sprintf([ttav,' %3.2e ;'],av(nv,n));
138      end ; hold off ;      end ; hold off ;
     if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;  
139      AA=axis ; dAA=AA(4)-AA(3);      AA=axis ; dAA=AA(4)-AA(3);
140      if AA(3)*AA(4) <= 0, AA(3)=min(AA(3),-dAA/10); AA(4)=max(AA(4),dAA/10); end      if AA(3)*AA(4) <= 0, AA(3)=min(AA(3),-dAA/10); AA(4)=max(AA(4),dAA/10); end
141        if ttax1 < ttax2, AA(1)=ttax1; AA(2)=ttax2; end;
142      axis(AA); grid ;      axis(AA); grid ;
    %AA=axis; text(AA(1)*.9+AA(2)*.1,AA(3)*0.1+AA(4)*0.9,ttmn);  
143      if nv == 1, title(['Avr ',titv,'  ',ttmn]); end      if nv == 1, title(['Avr ',titv,'  ',ttmn]); end
144      if nv == 2, title(['Std-Dev ',titv,'  ',ttav]); end      if nv == 2, title(['Std-Dev ',titv,'  ',ttav]); end
145      if nv == 3, title(['min ',titv,'  ',ttav]); legend(namLg(1:Nexp,:),0); end      if nv == 3, title(['min ',titv,'  ',ttav]); legend(namLg,0); end
146      if nv == 4, title(['Max ',titv,'  ',ttav]); end      if nv == 4, title(['Max ',titv,'  ',ttav]); end
    %if nv == 2, title(['Del-2 ',titv,'  ',ttav]); end  
147    end ; xlabel(titT);    end ; xlabel(titT);
148  %--  %--
149    axes('position',[.01,.01,.99,.99],'Visible','off');    axes('position',[.01,.01,.99,.99],'Visible','off');
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22