/[MITgcm]/MITgcm_contrib/high_res_cube/matlab-grid-generator/bin/meanprofile.m
ViewVC logotype

Annotation of /MITgcm_contrib/high_res_cube/matlab-grid-generator/bin/meanprofile.m

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

Revision 1.1.1.1 - (hide annotations) (download) (vendor branch)
Tue Nov 11 18:08:08 2003 UTC (23 years ago) by cnh
Branch: MAIN, initial
CVS Tags: baseline, HEAD
Changes since 1.1: +0 -0 lines 
Checking in work done with Dimitri on high-resolution cube gridding and parallel 
communications. 
   o code is in a contrib experiment for now so we can continue collaborating
     on it. However most code is general and will be moved into main branch once 
     it is fully hardened.
   o There are README files in the contrib root and in the subdirectories that
     explain the contents
1 cnh 1.1 function [Tm,Tsd] = meanprofile(GRID,T1,varargin);
2    
3     % Remove NaN's
4     T1( find(isnan(T1)) )=0;
5    
6     if nargin==3
7     dTm=(mean(T1,4)-mean(varargin{1},4)).*GRID.mskc;
8     mapmsk=GRID.mskc(:,:,1);
9     elseif nargin==4
10     dTm=(mean(T1,4)-mean(varargin{1},4)).*GRID.mskc;
11     mapmsk=varargin{2};
12     else
13     dTm=mean(T1,4).*GRID.mskc;
14     mapmsk=GRID.mskc(:,:,1);
15     end
16    
17     nxy=prod(size(GRID.rac));
18     nr=prod(size(GRID.mskc,3));
19     ra=reshape((GRID.rac(:).*mapmsk(:))*ones(1,nr),size(GRID.mskc)).*GRID.mskc;
20     dTm=reshape(dTm,[nxy nr]);
21     ra=reshape(ra,[nxy nr]);
22     Tm=squeeze( sum(dTm.*ra) ./ sum(ra) );
23     Tsd=(dTm'-Tm'*ones(1,nxy))';
24     Tsd=sqrt( squeeze( sum((Tsd.^2).*ra) ./ sum(ra) ) );
25    
26     function [F] = impliedflux(GRID,Q);
27     F=mean(Q,3).*GRID.mskc(:,:,1).*GRID.rac;
28     F=-cumsum([0 squeeze( sum(F,1) )]);
  ViewVC Help
Powered by ViewVC 1.1.22