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)

  
#=======  split into 4 nodes (3 each)  =======
cat ilist_199[2-4].txt | awk '{print $1}' > ilist_92--94
cat ilist_199[5-7].txt | awk '{print $1}' > ilist_95--97
cat ilist_199[89].txt ilist_2000.txt | awk '{print $1}' > ilist_98--00
cat ilist_200[1-3].txt | awk '{print $1}' > ilist_01--03
wc -l ilist_??--??
scp ilist_??--?? ta54:

#=======  split into 6 nodes (2 each)  =======
cat ilist_1992.txt ilist_1993.txt  | awk '{print $1}' > ilist_40
cat ilist_1994.txt ilist_1995.txt  | awk '{print $1}' > ilist_41
cat ilist_1996.txt ilist_1997.txt  | awk '{print $1}' > ilist_42
cat ilist_1998.txt ilist_1999.txt  | awk '{print $1}' > ilist_43
cat ilist_2000.txt ilist_2001.txt  | awk '{print $1}' > ilist_44
cat ilist_2002.txt ilist_2003.txt  | awk '{print $1}' > ilist_45
wc -l ilist_4[0-5]
scp ilist_4[0-5] ta54:


/scratch/edhill/matlab_r13/bin/matlab -nojvm -nodisplay

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

  !/sbin/ifconfig eth0 | grep 'inet addr:' | sed -e 's|\.| |g' | awk '{print $5}' > myip
  load myip
  ips = [ myip  201 200 200 200 ]

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


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
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	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	load ilist
40
41	calcUVBave(Spat, Tpat, Upat, Vpat, ilist, nlay, ne)
42
43
44
45	#======= split into 4 nodes (3 each) =======
46	cat ilist_199[2-4].txt \| awk '{print $1}' > ilist_92--94
47	cat ilist_199[5-7].txt \| awk '{print $1}' > ilist_95--97
48	cat ilist_199[89].txt ilist_2000.txt \| awk '{print $1}' > ilist_98--00
49	cat ilist_200[1-3].txt \| awk '{print $1}' > ilist_01--03
50	wc -l ilist_??--??
51	scp ilist_??--?? ta54:
52
53	#======= split into 6 nodes (2 each) =======
54	cat ilist_1992.txt ilist_1993.txt \| awk '{print $1}' > ilist_40
55	cat ilist_1994.txt ilist_1995.txt \| awk '{print $1}' > ilist_41
56	cat ilist_1996.txt ilist_1997.txt \| awk '{print $1}' > ilist_42
57	cat ilist_1998.txt ilist_1999.txt \| awk '{print $1}' > ilist_43
58	cat ilist_2000.txt ilist_2001.txt \| awk '{print $1}' > ilist_44
59	cat ilist_2002.txt ilist_2003.txt \| awk '{print $1}' > ilist_45
60	wc -l ilist_4[0-5]
61	scp ilist_4[0-5] ta54:
62
63
64	/scratch/edhill/matlab_r13/bin/matlab -nojvm -nodisplay
65
66	clear all ; close all
67	nlay = 50;
68	ne = 510;
69
70	Spat = 'c5l/Stave/data/Stave.%10.10d.data';
71	Tpat = 'c5l/Ttave/data/Ttave.%10.10d.data';
72	Upat = 'c5l/_Veltave/data/uVeltave.%10.10d.data';
73	Vpat = 'c5l/_Veltave/data/vVeltave.%10.10d.data';
74	load ilist
75
76	!/sbin/ifconfig eth0 \| grep 'inet addr:' \| sed -e 's\|\.\| \|g' \| awk '{print $5}' > myip
77	load myip
78	ips = [ myip 201 200 200 200 ]
79
80	calcUVBave_ttcp(ips, Spat, Tpat, Upat, Vpat, ilist, nlay, ne)
81
82
83