/[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.6 by jmc, Thu Mar 5 20:57:26 2015 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 32  if krd > 0, Line 33  if krd > 0,
33    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
34      read_StD(prefix,namA(n),listV);      read_StD(prefix,namA(n),listV);
35    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
   vv1(find(vv1==undef))=NaN;  
36  %- set global dims: & load vvA --> vvB  %- set global dims: & load vvA --> vvB
37    nbV=size(listV,2);    nbV=size(listV,2);
38    nrec=nIt; n3d=nk; nReg=nRg;    nrec=nIt; n3d=nk; nReg=nRg;
# Line 43  if krd > 0, Line 43  if krd > 0,
43    [ntA(n),rList,tim,vv1,listV] = ...    [ntA(n),rList,tim,vv1,listV] = ...
44      read_StD(prefix,namA(n),listV);      read_StD(prefix,namA(n),listV);
45    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);    nIt=ntA(n); nk=size(vv1,1); nRg=size(vv1,3);
   vv1(find(vv1==undef))=NaN;  
46    if (nrec < nIt),    if (nrec < nIt),
47      fprintf('\n');      fprintf('\n');
48      error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);      error([' Nb of records=',int2str(nIt),' exceeds nrec=',int2str(nrec)]);
# Line 130  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.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22