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
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
#include "libggcm.h"
#include <string.h>
#include <stdlib.h>
/** Our fields are stored in an MHDdata record, in no particular order
* (actually it depends on the input order, which we don't rely on). This
* scans through and copies the data for the desired field out of the
* structure. */
float *ggcm_field(MHDdata * data, int nx, int ny, int nz,
char * const field_name)
{
float *ret = (float*)malloc(sizeof(float) * (nx*ny*nz));
MHDdata *dat = data;
#if 0
int cnt=0;
float *cur = ret;
#endif
while(dat != NULL) {
if(strcmp(dat->field_name, field_name) == 0) {
#if 0
for(cnt = 0; cnt < (nx*ny*nz); ++cnt) {
*cur++ = dat->data[cnt];
}
#else
memcpy(ret, dat->data, sizeof(float) * (nx*ny*nz));
#endif
break;
}
dat = dat->next;
}
return ret;
}
float ***ggcm_field_matrix(const MHDdata * data, int nx, int ny, int nz,
char field_name[128])
{
float ***fmat = m_alloc(nx, ny, nz);
int i,j,k, count=0;
MHDdata *dat = (MHDdata*) data;
while (dat != NULL) {
if (strcmp(dat->field_name, field_name) == 0) {
for (k = 0; k < nz; ++k) {
for (j = 0; j < ny; ++j) {
for (i = 0; i < nx; ++i) {
fmat[i][j][k] = dat->data[count++];
}
}
}
}
dat = dat->next;
}
return fmat;
}
/** alias for ggcm_field_matrix */
__deprecated float ***ggcmFieldMatrix(const MHDdata *data,
int nx, int ny, int nz,
char field_name[128]) {
return ggcm_field_matrix(data, nx, ny, nz, field_name);
}
|
