/[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.7 by jmc, Sat Mar 7 14:45:49 2015 UTC revision 1.8 by jmc, Sun Sep 6 18:15:09 2015 UTC
# Line 9  Line 9 
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;  undef=123456.7;
15  %-----------  %-----------
# Line 83  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 108  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;
112    %fyb=160; fxb=-2600;
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 126  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+100*ng fyb+40*ng 500 700]);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     tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0;     tmp=var(isA(n):ieA(n),:,n); [I]=find(tmp(:,1)==undef); tmp(I,:)=0;
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22