high_res_cube/eddy_flux/run_calcUVBave.m

%
%  Ed Hill
%  Mon Aug  2 23:43:41 EDT 2004
%

matlab -nojvm -nodisplay

  clear all ; close all

  %  Get the 1-yr list of days
  load ilist
  deltatT = 1200;                % model time step size (s)
  tavefreq = 259200;             % averaging period (s)
  startDate = datenum(1992,1,1); % model integration starting date
  for iy = 1992:2004
    fid = fopen(sprintf('ilist_%d.txt',iy),'w');
    for it = 1:length(ilist)
      iter = ilist(it);
      dv = datevec(startDate + (iter*deltatT-tavefreq/2)/60/60/24);
      newyear = dv(1);
      if newyear == iy
        a = sprintf('%10.10d %d',iter,dv(1));
        fprintf(fid,'%s\n',a);
      end
    end
    fclose(fid);
  end
  

  clear all ; close all
  nlay = 50;
  ne = 510;
  
  Spat = 'Stave/data/Stave.%10.10d.data';
  Tpat = 'Ttave/data/Ttave.%10.10d.data';
  Upat = '_Veltave/data/uVeltave.%10.10d.data';
  Vpat = '_Veltave/data/vVeltave.%10.10d.data';
  
  load ilist

  calcUVBave(Spat, Tpat, Upat, Vpat, ilist, nlay, ne)


  clear all ; close all
  nlay = 50;
  ne = 510;
  
  Spat = 'c5l/Stave/data/Stave.%10.10d.data';
  Tpat = 'c5l/Ttave/data/Ttave.%10.10d.data';
  Upat = 'c5l/_Veltave/data/uVeltave.%10.10d.data';
  Vpat = 'c5l/_Veltave/data/vVeltave.%10.10d.data';
  
  load ilist

  calcUVBave_ttcp(Spat, Tpat, Upat, Vpat, ilist, nlay, ne)


1	edhill	1.1	%
2			% Ed Hill
3			% Mon Aug 2 23:43:41 EDT 2004
4			%
5
6			matlab -nojvm -nodisplay
7
8			clear all ; close all
9	edhill	1.2
10			% Get the 1-yr list of days
11			load ilist
12			deltatT = 1200; % model time step size (s)
13			tavefreq = 259200; % averaging period (s)
14			startDate = datenum(1992,1,1); % model integration starting date
15			for iy = 1992:2004
16			fid = fopen(sprintf('ilist_%d.txt',iy),'w');
17			for it = 1:length(ilist)
18			iter = ilist(it);
19			dv = datevec(startDate + (iter*deltatT-tavefreq/2)/60/60/24);
20			newyear = dv(1);
21			if newyear == iy
22			a = sprintf('%10.10d %d',iter,dv(1));
23			fprintf(fid,'%s\n',a);
24			end
25			end
26			fclose(fid);
27			end
28
29
30			clear all ; close all
31	edhill	1.1	nlay = 50;
32			ne = 510;
33
34			Spat = 'Stave/data/Stave.%10.10d.data';
35			Tpat = 'Ttave/data/Ttave.%10.10d.data';
36			Upat = '_Veltave/data/uVeltave.%10.10d.data';
37			Vpat = '_Veltave/data/vVeltave.%10.10d.data';
38
39	edhill	1.2	load ilist
40
41			calcUVBave(Spat, Tpat, Upat, Vpat, ilist, nlay, ne)
42
43	edhill	1.1
44
45	edhill	1.3	clear all ; close all
46			nlay = 50;
47			ne = 510;
48
49			Spat = 'c5l/Stave/data/Stave.%10.10d.data';
50			Tpat = 'c5l/Ttave/data/Ttave.%10.10d.data';
51			Upat = 'c5l/_Veltave/data/uVeltave.%10.10d.data';
52			Vpat = 'c5l/_Veltave/data/vVeltave.%10.10d.data';
53
54			load ilist
55
56			calcUVBave_ttcp(Spat, Tpat, Upat, Vpat, ilist, nlay, ne)
57
58
59