/[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.7 by jmc, Sat Mar 7 14:45:49 2015 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 128  for jv=1:nbV, Line 129  for jv=1:nbV,
129    figure(ng); set(ng,'position',[100+100*ng 60+40*ng 500 700]);clf;    figure(ng); set(ng,'position',[100+100*ng 60+40*ng 500 700]);clf;
130    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);
131    for n=1:Nexp,    for n=1:Nexp,
132     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;
133     av(:,n)=mean(var(isA(n):ieA(n),:,n));     av(:,n)=sum(tmp);
134       if length(I) < size(tmp,1), av(:,n)=av(:,n)/(size(tmp,1)-length(I)); end
135       tmp(I,:)=NaN; dd(:,n)=max(tmp)-min(tmp);
136    end    end
137      var(find(var==undef))=NaN;
138    for nv=1:4,    for nv=1:4,
139      subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:';      subplot(410+nv); ttmn=' Mx-mn:'; ttav=' Av:';
140      for n=1:Nexp,      for n=1:Nexp,
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22