/[MITgcm]/MITgcm_contrib/enderton/Diagnostics/DiagUtility/calc_ZonAv_CS.m
ViewVC logotype

Diff of /MITgcm_contrib/enderton/Diagnostics/DiagUtility/calc_ZonAv_CS.m

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

revision 1.1 by enderton, Mon Jan 31 15:43:28 2005 UTC revision 1.3 by enderton, Tue Oct 18 19:22:15 2005 UTC
# Line 112  end Line 112  end
112  if kpr > 0 & krd > 0,  if kpr > 0 & krd > 0,
113      rac=dir;      rac=dir;
114  %if nBas < 0, hFacC=ones(ncx,nc,dim3); else hFacC=rdmds([rac,'hFacC']); end  %if nBas < 0, hFacC=ones(ncx,nc,dim3); else hFacC=rdmds([rac,'hFacC']); end
115  %nr=size(hFacC,3);      %nr_hfac=size(hFacC,3);
116      if nBas < 0, hFacC=ones(ncx,nc,nr); end%else hFacC=rdmds([rac,'hFacC']); end      if nBas < 0, hFacC=ones(ncx,nc,nr); end%else hFacC=rdmds([rac,'hFacC']); end
117      mskC=min(1,hFacC); mskC=ceil(mskC);      mskC=min(1,hFacC); mskC=ceil(mskC);
118        %hFacC=reshape(hFacC,nPg,nr_hfac); mskC=reshape(mskC,nPg,nr_hfac);
119        if nr == 1, hFacC = hFacC(:,:,1); mskC = mskC(:,:,1); end % Cludge by Daniel.
120      hFacC=reshape(hFacC,nPg,nr); mskC=reshape(mskC,nPg,nr);      hFacC=reshape(hFacC,nPg,nr); mskC=reshape(mskC,nPg,nr);
121      a6c=reshape(arc,nPg,1);      a6c=reshape(arc,nPg,1);
122      if nBas > 0,      if nBas > 0,
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22