Line | |
---|
1 | function [ output_args ] = make_bck_twin( totaltime) |
---|
2 | |
---|
3 | outfile='../obs_float/vel_bck_all.dat'; |
---|
4 | format long |
---|
5 | ufile='../obs_float/uzero.dat'; |
---|
6 | vfile='../obs_float/vzero.dat'; |
---|
7 | fid=fopen(outfile,'w'); |
---|
8 | |
---|
9 | for i=1:totaltime |
---|
10 | |
---|
11 | |
---|
12 | time_step=i; |
---|
13 | ubck=load(ufile); |
---|
14 | vbck=load(vfile); |
---|
15 | |
---|
16 | |
---|
17 | vel_bck=zeros(size(ubck,1),5); |
---|
18 | |
---|
19 | vel_bck(:,1)=time_step; |
---|
20 | vel_bck(:,2)=ubck(:,2); |
---|
21 | vel_bck(:,3)=ubck(:,3); |
---|
22 | vel_bck(:,4)=ubck(:,4); |
---|
23 | vel_bck(:,5)=vbck(:,4); |
---|
24 | |
---|
25 | |
---|
26 | |
---|
27 | fprintf(fid,'%d %d %d %f %f\n',vel_bck'); |
---|
28 | |
---|
29 | |
---|
30 | |
---|
31 | end |
---|
32 | |
---|
33 | fclose(fid); |
---|
34 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.