1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
function ppp_ustar2h (Ustar,name)
## usage: ppp_Ustar2h (Ustar[,name])
##
##
if nargin<2
name = "Ustar";
endif
[N,N_U] = size(Ustar);
## Open the file
filename = sprintf("%s.h", name);
fid = fopen(filename,"w");
## Header
header = sprintf("/*\n File %s generated by ppp_ustar2h on %s */\n", \
filename, ctime(time));
def = sprintf("#define N_U %i\n#define N_T %i\n", N_U, N);
fprintf(fid, "%s%sdouble %s[N_T][N_U] = {\n",header,def,name);
for i=1:N
fprintf(fid, "{");
for j=1:N_U
if j<N_U
comma = ",";
|
|
>
>
>
>
|
>
>
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
function ppp_ustar2h (Ustar,DT,name)
## usage: ppp_Ustar2h (Ustar[,name])
##
##
if nargin<2
DT = 1;
endif
if nargin<3
name = "Ustar";
endif
[N,N_U] = size(Ustar);
## Open the file
filename = sprintf("%s.h", name);
fid = fopen(filename,"w");
## Header
header = sprintf("/*\n File %s generated by ppp_ustar2h on %s */\n", \
filename, ctime(time));
def = sprintf("#define N_U %i\n#define N_T %i\n#define DT %g\n", \
N_U, N, DT);
def = sprintf("%sdouble U[N_U];\n",def);
fprintf(fid, "%s%sdouble %s[N_T][N_U] = {\n",header,def,name);
for i=1:N
fprintf(fid, "{");
for j=1:N_U
if j<N_U
comma = ",";
|