main()
{
printf("ok!\n");
double I=0.08;
double L=0.02;
double Uavg=14;
double x[10000][3],v[3],x0[3];
double dt;
double dk;
int kmax;
int nt;
int GridNumber;
FILE *fd0;
ne=500; //! spectral sample size
dt=0.001;
kmax=500;//(int)(1/(dt));
dk=1;//0.5;
genspec_(kmax,dk);
nt=2000; //! number of time steps
fd0=fopen("grid.txt","r");
GridNumber=0;
do
{
double t0,t1,t2;
fscanf(fd0,"%lf,%lf,%lf",&t0,&t1,&t2);
x[GridNumber][0]=t0;
x[GridNumber][1]=t1;
x[GridNumber][2]=t2;
GridNumber++;
} while(!feof(fd0));
fclose(fd0);
for(int i=0;i
FILE *fd;
fd=fopen("nnn.txt","w");
fprintf(fd,"%d,%d\n",GridNumber,nt);
for(int it=0;it
double t=it*dt;
printf("%f\n",t);
//fprintf(fd,"%lf",t);
for(int ix=0;ix
x0[0]=x[ix][0]-0.004700;
x0[1]=x[ix][1];
x0[2]=x[ix][2];
I=Turb_Itensity(x0[2]);
L=Turb_LengthScale(x0[2]);
Uavg=VelocityProfile(x0[2]);
genvec_2(t,x0,v,kmax,dk,I,L,Uavg);
fprintf(fd,"%lf,%lf,%lf\n",v[0],v[1],v[2]);
}
//fprintf(fd,"\n");
}
fclose(fd);
delspec_();
return 0;
}