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

Diff of /MITgcm_contrib/jmc_script/plot_StD.m

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

revision 1.1 by jmc, Wed Oct 8 18:18:04 2008 UTC revision 1.2 by jmc, Fri Sep 27 23:11:26 2013 UTC
# Line 1  Line 1 
1   prefix='stDiag';   prefix='dynStD';
2   namA='r527';   namA='r527';
3   Nexp=1; Nc=size(namA,2);   Nexp=1; Nc=size(namA,2);
4  %--  %--
5    
# Line 64  elseif krd < 0, Line 64  elseif krd < 0,
64  end  end
65  if krd ~= 0,  if krd ~= 0,
66   ttA=squeeze(tiA(:,2,:));   ttA=squeeze(tiA(:,2,:));
67   ttA=ttA/86400; titT='days'; %ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';   ttA=ttA/3600; titT='hrs';  ttA=ttA/24; titT='days';
68     ttA=ttA/30 ; titT='month'; ttA=ttA/12 ; titT='year';
69  end  end
70  %=========================================================  %=========================================================
71    
# Line 165  for ng=1:nbV, Line 166  for ng=1:nbV,
166      subplot(410+nv);      subplot(410+nv);
167      var=squeeze(vv1(2:nk,:,1,nv,n))';      var=squeeze(vv1(2:nk,:,1,nv,n))';
168      mnV=min(var(:)); MxV=max(var(:)); ccv=c_levs(mnV,MxV,-12);      mnV=min(var(:)); MxV=max(var(:)); ccv=c_levs(mnV,MxV,-12);
169      [cs,h]=contour(ttA(isA(n):ieA(n),n),[1:nk-1],var(isA(n):ieA(n),:)',ccv);     %- pressure coordinate, from bottom (k=1) to top (k=nr=nk-1)
170       %[cs,h]=contour(ttA(isA(n):ieA(n),n),[1:nk-1],var(isA(n):ieA(n),:)',ccv);
171       %- z coordinate, from bottom (k=nr=nk-1) to top (k=1)
172        [cs,h]=contour(ttA(isA(n):ieA(n),n),-[1:nk-1],var(isA(n):ieA(n),:)',ccv);
173     %clabel(cs);isoline0(h);     %clabel(cs);isoline0(h);
174      colorbar;      colorbar;
175      if nv == 1, title(['Avr ',titv]); end      if nv == 1, title(['Avr ',titv]); end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22