/[MITgcm]/MITgcm/verification/global_ocean.90x40x15/diags_matlab/mit_meridflux.m
ViewVC logotype

Contents of /MITgcm/verification/global_ocean.90x40x15/diags_matlab/mit_meridflux.m

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

Revision 1.1.2.1 - (show annotations) (download)
Wed Oct 23 18:26:35 2002 UTC (21 years, 6 months ago) by mlosch
Branch: release1
CVS Tags: release1_p12, release1_p10, release1_p16, release1_p15, release1_p11, release1_p14, release1_p13, release1_p17, release1_p8, release1_p9, release1_p6, release1_p7, release1_p13_pre, release1_p12_pre
Branch point for: release1_50yr
Changes since 1.1: +0 -0 lines 
o fixed the verification/global_ocean.90x40x15 experiment:
 - new bathymetry (the world according to A., JMC, and M.)
 - new initial fields and forcing fields (*.bin files)
 - new POLY3.COEFFS (for the next release one should switch to a full
   equation of state)
 - fixed several errors and redundancies in the data file
 - experiment uses looped cells
 - added matlab directory with diagnostic scripts for plotting of output
1 function T = mit_meridflux(surface_flux,dx,dy)
2 % function T = mit_meridflux(surface_flux,dx,dy)
3
4 if nargin == 2
5 area = dx;
6 elseif nargin == 3
7 area = dx.*dy;
8 else
9 error('need two or three input parameters')
10 end
11 [nx ny nt] = size(surface_flux);
12 for kt=1:nt
13 fdxdy(:,:,kt) = change(surface_flux(:,:,kt),'==',NaN,0).*area;
14 end
15
16 T = squeeze(sum(cumsum(fdxdy,2),1));
17
18 return
  ViewVC Help
Powered by ViewVC 1.1.22