Context Navigation

takecorrectionasbck.m @ 160

Last change on this file since 160 was 129, checked in by jbrlod, 10 years ago
last version of Varanth
File size: 1.3 KB

Line
1	outfile='obs_float/ucorr.dat';
2	load obs_float/uv_total.dat;
3	load obs_float/uzero.dat
4
5	umean(:,1)=uzero(:,1);
6	umean(:,2)=uzero(:,2);
7	umean(:,3)=uzero(:,3);
8	U1=zeros(5046,1);
9	V1=zeros(5046,1);
10	e=1;
11	for r=1:1:117
12
13	Unew=zeros(1,5046);
14	Vnew=zeros(1,5046);
15	for count=1:5046
16	Unew(count)=uv_total(e,end-1);
17	Vnew(count)=uv_total(e,end);
18	e=e+1;
19	end
20	U1=U1+Unew';
21	V1=V1+Vnew';
22
23
24	end
25
26	U1=U1./117;
27	V1=V1./117;
28	%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
29	%%% the uv_total.
30	U1=U1;
31	V1=V1;
32	umean(:,4)=U1;
33
34
35	fid=fopen(outfile,'w');
36	fprintf(fid,'%d %d %d %f\n',umean');
37	fclose(fid);
38
39	outfile='obs_float/vcorr.dat';
40	load obs_float/uv_total.dat;
41	load obs_float/uzero.dat
42
43	vmean(:,1)=uzero(:,1);
44	vmean(:,2)=uzero(:,2);
45	vmean(:,3)=uzero(:,3);
46	U1=zeros(5046,1);
47	V1=zeros(5046,1);
48	e=1;
49	for r=1:1:117
50
51	Unew=zeros(1,5046);
52	Vnew=zeros(1,5046);
53	for count=1:5046
54	Unew(count)=uv_total(e,end-1);
55	Vnew(count)=uv_total(e,end);
56	e=e+1;
57	end
58	U1=U1+Unew';
59	V1=V1+Vnew';
60
61
62	end
63
64	U1=U1./117;
65	V1=V1./117;
66	%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
67	%%% the uv_total.
68	U1=U1;
69	V1=V1;
70	vmean(:,4)=V1;
71
72
73	fid=fopen(outfile,'w');
74	fprintf(fid,'%d %d %d %f\n',vmean');
75	fclose(fid);

Note: See TracBrowser for help on using the repository browser.