脚本的第一部分必须从文本文件中读取风数据..文件的第一个数字是行的编号,每行包含6个数据:纬度,经度,高度,模块,高程和方位角。 当我在调试模式下在VC ++上运行我的脚本时,我注意到 dat 不是wind_dat的向量,而只是* wind_dat *的类型元素。但是,如果我在 fscanf 的 循环中插入断点,我确定所有变量都被读取,但未保存在 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);
}
答案 0 :(得分:0)
在“观察”窗口中,键入dat,100
以查看数组dat
指向的前100个元素。
调试器无法从变量声明中知道dat
实际指向多个对象。所以它需要一个提示。