/[MITgcm]/MITgcm_contrib/high_res_cube/eddy_flux/run_calcUVBave.m
ViewVC logotype

Annotation of /MITgcm_contrib/high_res_cube/eddy_flux/run_calcUVBave.m

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


Revision 1.4 - (hide annotations) (download)
Sat Aug 14 19:31:09 2004 UTC (20 years, 11 months ago) by edhill
Branch: MAIN
CVS Tags: HEAD
Changes since 1.3: +28 -4 lines
 o now the averaging scales nicely across 6 nodes (2yrs of data per node)

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.4
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 edhill 1.1
64 edhill 1.4 /scratch/edhill/matlab_r13/bin/matlab -nojvm -nodisplay
65    
66     clear all ; close all
67 edhill 1.3 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 edhill 1.4 !/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 edhill 1.3
82    
83    

  ViewVC Help
Powered by ViewVC 1.1.22