0

スクリプトの最初の部分は、テキスト ファイルから風データを読み取る必要があります。ファイルの最初の数値は行の番号で、各行には 6 つのデータが含まれています: 緯度、経度、高さ、モジュール、標高、方位角。スクリプトを VC++ でデバッグ モードで実行すると、datが wi​​nd_dat のベクトルではなく、*wind_dat* の型要素であることに気付きました。ただし、 fscanfのforサイクルにブレークポイントを挿入すると、すべての変数が読み取られると確信していますが、datには保存されず、最初の行のみが保存されます..

ここに、スクリプトの関連コードを投稿します。

typedef struct wind {
double xlat;
double ylong;
double zheight;
double mod;
double elev;
double az;} wind_dat;

static wind_dat *  dataprocess( const char *, int*);
static void   pol2cart( double *, double, double, double);

wind_dat       *dat;
int         num;


int main(int argc, char **argv)
{   

dat=dataprocess(filename,&num);
...
}

static wind_dat* dataprocess (const char *filename, int *num)
{
int it;
double coord[3], obs_cart[3], a[3][3];
FILE *ptrfile;

if ((fopen_s(&ptrfile,filename, "r"))!=0)
{printf("\nError while opening file.\n");
exit(EXIT_FAILURE);}
fscanf_s(ptrfile,"%d",num);

if ((dat=(wind_dat*)malloc((*num)*sizeof(wind_dat)))==NULL)
    exit(EXIT_FAILURE);

pol2cart(obs_cart,obs[0],obs[1],obs[2]+6371000);
a[0][0]=cos(az)*cos(obs[1]+PI)+sin(az)*sin(PI+obs[1])*cos(PI-obs[0]);
a[0][1]=cos(az)*sin(obs[1]+PI)-sin(az)*cos(obs[1]+PI)*cos(PI-obs[0]);
a[0][2]=-sin(az)*sin(PI-obs[0]);
a[1][0]=sin(az)*cos(obs[1]+PI)-cos(az)*sin(obs[1]+PI)*cos(PI-obs[0]);
a[1][1]=sin(az)*sin(obs[1]+PI)+cos(az)*cos(obs[1]+PI)*cos(PI-obs[0]);
a[1][2]=cos(az)*sin(PI-obs[0]);
a[2][0]=sin(obs[1]+PI)*sin(PI-obs[0]);
a[2][1]=-sin(PI-obs[0])*cos(obs[1]+PI);
a[2][2]=cos(PI-obs[0]);
for (it=0;it<(*num);it++)
    {fscanf_s(ptrfile,"%lf %lf %lf %lf %lf %lf",&dat[it].xlat, &dat[it].ylong, &dat[it].zheight, &dat[it].mod, &dat[it].elev, &dat[it].az);

    pol2cart(coord, dat[it].xlat, dat[it].ylong, dat[it].zheight+6371000);
    coord[0]=coord[0]-obs_cart[0];
    coord[1]=coord[1]-obs_cart[1];
    coord[2]=coord[2]-obs_cart[2];
    dat[it].xlat =  coord[0]*a[0][0]+coord[1]*a[0][1]+coord[2]*a[0][2];
    dat[it].ylong = coord[0]*a[1][0]+coord[1]*a[1][1]+coord[2]*a[1][2];
    dat[it].zheight=coord[0]*a[2][0]+coord[1]*a[2][1]+coord[2]*a[2][2];
    dat[it].az=dat[it].az-az/PI*180;
    dat[it].mod=dat[it].mod*arrow_factor;
    }
fclose(ptrfile);
return(dat);
}

static void pol2cart(double *dat, double x, double y, double z)
{

dat[0]=z*cos(x*180/PI)*cos(y*180/PI);
dat[1]=z*cos(x*180/PI)*sin(y*180/PI);
dat[2]=z*sin(x*180/PI);
}
4

1 に答える 1