/[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.5 by jmc, Mon Aug 25 21:35:12 2014 UTC revision 1.9 by jmc, Sat Sep 10 22:05:35 2016 UTC
# Line 1  Line 1 
1   prefix='dynStD';   prefix='dynStD';
2    %prefix='iceStD';
3   namA={'g02','b02'};   namA={'g02','b02'};
4  %namA={'r17a'};  %namA={'r17a'};
5   Nexp=size(namA,2);   Nexp=size(namA,2);
# Line 7  Line 8 
8  % $Header$  % $Header$
9  % $Name$  % $Name$
10    
   
11  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;  nItMx=1e10*ones(1,Nexp); %nItMx(3)=11;
12  nItMx=2400*ones(1,Nexp);  %nItMx=2400*ones(1,Nexp);
13  namLg=namA ; namLg=strrep(namLg,'_','\_');  namLg=namA ; namLg=strrep(namLg,'_','\_');
14    undef=123456.7;
15  %-----------  %-----------
16  %- test if the variable krd is define :  %- test if the variable krd is define :
17  if size(who('krd'),1) > 0,  if size(who('krd'),1) > 0,
# Line 82  end; Line 83  end;
83  %--  %--
84    
85  list_on=zeros(1,nbV);  list_on=zeros(1,nbV);
86  nbG=8;  nbG=18;
87  nbG=min(nbG,nbV); list_on(1:nbG)=1 ;  nbG=min(nbG,nbV); list_on(1:nbG)=1 ;
88  %if nbG < nbV, list_on=0; list_on(nbG:nbV)=1; end %- to get 2nd set of plots  %if nbG < nbV, list_on=0; list_on(nbG:nbV)=1; end %- to get 2nd set of plots
89  %list_on(1:6)=[1 1 1 1 1 1];  %list_on(1:6)=[1 1 1 1 1 1];
# Line 107  titall='Global Ocean, Cubic-G (32x32) , Line 108  titall='Global Ocean, Cubic-G (32x32) ,
108  %titall='Dyncore test-case 5 (cs-32)' ;  %titall='Dyncore test-case 5 (cs-32)' ;
109    
110  %=========================================================  %=========================================================
111    ng=0; fxb=100; fyb=60; fdx=100; fdy=40; fsc=1.;
112    %fyb=-360; fxb=-2600; fdy=60; fsc=1.5;
113    
 ng=0;  
114  for jv=1:nbV,  for jv=1:nbV,
115  %-------------------  %-------------------
116   flag=list_on(jv); kl=0;   flag=list_on(jv); kl=0;
# Line 125  for jv=1:nbV, Line 127  for jv=1:nbV,
127    ng=ng+1;    ng=ng+1;
128  %- reset "ng" to jv for fix fig number (independent of which one is in list_on):  %- reset "ng" to jv for fix fig number (independent of which one is in list_on):
129    %ng=jv;    %ng=jv;
130    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;
131    var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp);    var=squeeze(vv1(1+kl,:,1,:,:)); dd=zeros(5,Nexp); av=zeros(5,Nexp);
132    for n=1:Nexp,    for n=1:Nexp,
133     dd(:,n)=max(var(isA(n):ieA(n),:,n))-min(var(isA(n):ieA(n),:,n));     tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0;
134     av(:,n)=mean(var(isA(n):ieA(n),:,n));     av(:,n)=sum(tmp);
135       if length(I) < size(tmp,1), av(:,n)=av(:,n)/(size(tmp,1)-length(I)); end
136       tmp(I,:)=NaN; dd(:,n)=max(tmp)-min(tmp);
137    end    end
138      var(find(var==undef))=NaN;
139    for nv=1:4,    for nv=1:4,
140      subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:';      subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:';
141      for n=1:Nexp,      for n=1:Nexp,
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22