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

Diff of /MITgcm_contrib/jmc_script/grph_MON.m

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

revision 1.1 by jmc, Mon Nov 9 17:36:04 2020 UTC revision 1.2 by jmc, Tue Feb 8 17:06:00 2022 UTC
# Line 227  for ng=1:size(list_on,2) Line 227  for ng=1:size(list_on,2)
227    for nv=1:ngEx-ngEk,    for nv=1:ngEx-ngEk,
228      vvM=mean(vvA(max(2,isA(n)):ieA(n),nv,:),1);      vvM=mean(vvA(max(2,isA(n)):ieA(n),nv,:),1);
229      subplot(100*ngEx+10+nv+2); ttmn=' Mx-mn:'; ttav=' Av:';      subplot(100*ngEx+10+nv+2); ttmn=' Mx-mn:'; ttav=' Av:';
230        if nv == 1, titv2=' \theta'; else titv2=' S '; end
231        if list_log(ng) > 0,
232          var=zeros(nrec,Nexp); for n=1:Nexp, var(isA(n):ieA(n),n)=vvA(isA(n):ieA(n),nv,n); end
233          if min(var(:))*max(var(:)) >= 0, msk=var; ttyax='log';
234            var(find(msk==0))=1; var=log10(abs(var)); var(find(msk==0))=NaN;
235          end
236        else var=squeeze(vvA(:,nv,:)); end
237      for n=1:Nexp,      for n=1:Nexp,
238        plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));        plot(ttA(isA(n):ieA(n),n),var(isA(n):ieA(n),n),char(linA(n)));
239        if n == 1, hold on ; end ;        if n == 1, hold on ; end ;
# Line 235  for ng=1:size(list_on,2) Line 242  for ng=1:size(list_on,2)
242      end ; hold off ;      end ; hold off ;
243      if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;      if ttax1 < ttax2, AA=axis; axis([ttax1 ttax2 AA(3:4)]); end;
244      grid ;      grid ;
245      if nv == 1, title(['mean ',titv,' \theta']); end     %title(['mean ',titv1,titv2]);
246      if nv == 2, title(['mean ',titv,' S ']); end      title(['mean ',titv1,titv2,' ; ',ttav]);
247  %--  %--
248      titX='   '; for n=1:Nexp, titX=[titX,sprintf(' %4.3e',vvM(n))]; end      titX='   '; for n=1:Nexp, titX=[titX,sprintf(' %4.3e,',vvM(n))]; end
249      AA=axis; TM=text(AA(2),AA(3)+0.1*(AA(4)-AA(3)),titX);      AA=axis; TM=text(AA(2),AA(3)+0.1*(AA(4)-AA(3)),titX(1:end-1));
250      set(TM,'HorizontalAlignment','right');      set(TM,'HorizontalAlignment','right');
251  %--  %--
252    end ;    end ;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22