亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? 555.cpp

?? 基于CIP的海嘯波浪數值模擬
?? CPP
?? 第 1 頁 / 共 5 頁
字號: 
exit(1);
}
fprintf(fpdat,"t=%16.8e\n",time);
fprintf(fpdat,"%16s %16s %16s %16s %16s %16s %16s\n","x","y","pressure","density","u","w","phi");
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fpdat,"%16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e\n",
(i-2)*DX,(j-2)*DY,p[i][j],rho[i][j],u[i][j],v[i][j],iphi(phi[i][j]));
}
}
fclose(fpdat);
return(0);
}
int output1(double f[IMAX][JMAX], int count){
FILE *fp;
int i, j;
int stp;
char svname[25]= "./Data/ot1????.dat";
stp=count/OUTPUT;
svname[10] = '0' + stp/1000;
svname[11] = '0' + (stp%1000)/100;
svname[12] = '0' + ((stp%1000)%100)/10;
svname[13] = '0' + ((stp%1000)%100)%10;
if ((fp=fopen (svname,"w+")) == NULL){
printf("%s was not opened!\n",svname);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
printf("in %d steps, data outputted to %s...........",count, svname);
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fp,"%2d %2d %30.22e\n",i,j,f[i][j]);
}
}
printf("Done\n");
fclose(fp);
return(0);
}
int output2(double f[IMAX][JMAX], int count){
FILE *fp;
int i, j;
int stp;
char svname[25]= "./Data/ot2????.dat";
stp=count/OUTPUT;
svname[10] = '0' + stp/1000;
svname[11] = '0' + (stp%1000)/100;
svname[12] = '0' + ((stp%1000)%100)/10;
svname[13] = '0' + ((stp%1000)%100)%10;
if ((fp=fopen (svname,"w+")) == NULL){
printf("%s was not opened!\n",svname);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
printf("in %d steps, data outputted to %s...........",count, svname);
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fp,"%2d %2d %30.22e\n",i,j,f[i][j]);
}
}
printf("Done\n");
fclose(fp);
return(0);
}
int output3(double f[IMAX][JMAX], int count){
FILE *fp;
int i, j;
int stp;
char svname[25]= "./Data/ot3????.dat";
stp=count/OUTPUT;
svname[10] = '0' + stp/1000;
svname[11] = '0' + (stp%1000)/100;
svname[12] = '0' + ((stp%1000)%100)/10;
svname[13] = '0' + ((stp%1000)%100)%10;
if ((fp=fopen (svname,"w+")) == NULL){
printf("%s was not opened!\n",svname);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
printf("in %d steps, data outputted to %s...........",count, svname);
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fp,"%2d %2d %30.22e\n",i,j,f[i][j]);
}
}
printf("Done\n");
fclose(fp);
return(0);
}
int output4(double f[IMAX][JMAX], int count){
FILE *fp;
int i, j;
int stp;
char svname[25]= "./Data/ot4????.dat";
stp=count/OUTPUT;
svname[10] = '0' + stp/1000;
svname[11] = '0' + (stp%1000)/100;
svname[12] = '0' + ((stp%1000)%100)/10;
svname[13] = '0' + ((stp%1000)%100)%10;
if ((fp=fopen (svname,"w+")) == NULL){
printf("%s was not opened!\n",svname);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
printf("in %d steps, data outputted to %s...........",count, svname);
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fp,"%2d %2d %30.22e\n",i,j,f[i][j]);
}
}
printf("Done\n");
fclose(fp);
return(0);
}
int output5(double f[IMAX][JMAX], int count){
FILE *fp;
int i, j;
int stp;
char svname[25]= "./Data/ot5????.dat";
stp=count/OUTPUT;
svname[10] = '0' + stp/1000;
svname[11] = '0' + (stp%1000)/100;
svname[12] = '0' + ((stp%1000)%100)/10;
svname[13] = '0' + ((stp%1000)%100)%10;
if ((fp=fopen (svname,"w+")) == NULL){
printf("%s was not opened!\n",svname);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
printf("in %d steps, data outputted to %s...........",count, svname);
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fp,"%2d %2d %30.22e\n",i,j,f[i][j]);
}
}
printf("Done\n");
fclose(fp);
return(0);
}
int output1d(double f[IMAX][JMAX], int count){
FILE *fp;
int i, j;
int stp;
char svname[25]= "./Data/1dp????.dat";
stp=count/OUTPUT;
svname[10] = '0' + stp/1000;
svname[11] = '0' + (stp%1000)/100;
svname[12] = '0' + ((stp%1000)%100)/10;
svname[13] = '0' + ((stp%1000)%100)%10;
if ((fp=fopen (svname,"w+")) == NULL){
printf("%s was not opened!\n",svname);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
printf("in %d steps, data outputted to %s...........",count, svname);
for(j=0; j<JMAX; j++){
fprintf(fp,"%16.8e\n",f[IMAX-4][j]);
}
fclose(fp);
printf("Done\n");
return(0);
}
int output1d_cav(double u[IMAX][JMAX], double v[IMAX][JMAX], int count){
FILE *fpu, *fpv;
int i, j;
int stp;
char dcau[25]="./Data/cau????.dat";
char dcav[25]="./Data/cav????.dat";
stp=count/OUTPUT_CAV;
dcau[10] = '0' + stp/1000;
dcau[11] = '0' + (stp%1000)/100;
dcau[12] = '0' + ((stp%1000)%100)/10;
dcau[13] = '0' + ((stp%1000)%100)%10;
dcav[10] = '0' + stp/1000;
dcav[11] = '0' + (stp%1000)/100;
dcav[12] = '0' + ((stp%1000)%100)/10;
dcav[13] = '0' + ((stp%1000)%100)%10;
if ((fpu=fopen (dcau,"w+")) == NULL){
printf("%s was not opened!\n",dcau);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
if ((fpv=fopen (dcav,"w+")) == NULL){
printf("%s was not opened!\n",dcav);
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
fclose(fpu);
fclose(fpv);
return(0);
}
int set_initial_condition(double rho[IMAX][JMAX], double gxrho[IMAX][JMAX], double gyrho[IMAX][JMAX],
							double rhon[IMAX][JMAX], double gxrhon[IMAX][JMAX], double gyrhon[IMAX][JMAX],
							double u[IMAX][JMAX], double gxu[IMAX][JMAX], double gyu[IMAX][JMAX],
							double un[IMAX][JMAX], double gxun[IMAX][JMAX], double gyun[IMAX][JMAX],
							double v[IMAX][JMAX], double gxv[IMAX][JMAX], double gyv[IMAX][JMAX],
							double vn[IMAX][JMAX], double gxvn[IMAX][JMAX], double gyvn[IMAX][JMAX],
							double phi[IMAX][JMAX], double gxphi[IMAX][JMAX], double gyphi[IMAX][JMAX],
							double phin[IMAX][JMAX], double gxphin[IMAX][JMAX], double gyphin[IMAX][JMAX],
							double p[IMAX][JMAX], double gxp[IMAX][JMAX], double gyp[IMAX][JMAX],
							double pn[IMAX][JMAX], double gxpn[IMAX][JMAX], double gypn[IMAX][JMAX]){
		int i,j;
		double zz, dlt;
		dlt=DD0R*D0;
		for(i=0; i<IMAX; i++){
			for(j=0; j<JMAX; j++){
			/* water area */
			// if(j<=INLET-1)   //slx.add 由于D0=INLET*DY,表明有INLET－1個縱向格子為水,
			// 但從下面的壓力計算式來看，好像是認為INLET+1,但最下面兩個格子似乎是虛擬的水
				if(j<=INLET+1){
					/***** for phi *****/
					phi[i][j]=dphi(IN_PHI); /* water */
					phin[i][j]=dphi(IN_PHI); /* water */
					/***** for rho *****/
					rho[i][j]=IN_RHO;
					rhon[i][j]=IN_RHO;
					/***** for pressure *****/
					p[i][j]=RHOwater*G*(INLET-(j-1.5))*DY+Patm;
					pn[i][j]=p[i][j];
					u[i][j]=IN_U;
					v[i][j]=IN_V;
					un[i][j]=u[i][j];
					vn[i][j]=v[i][j];
				}
				else{
					/* air area */
					/***** for phi *****/
					phi[i][j]=dphi(0.); /* air */
					phin[i][j]=dphi(0.); /* air */
					/***** for rho *****/
					rho[i][j]=RHOair;
					rhon[i][j]=RHOair;
					/***** for pressure *****/
					p[i][j]=Patm;
					pn[i][j]=p[i][j];
					/***** for velocity *****/
					u[i][j]=IN_U;
					v[i][j]=0.;
					un[i][j]=u[i][j];
					vn[i][j]=v[i][j];
				}
			}//for.j  slx.add
		}
		/***** for gradient *****/
		for(i=0; i<IMAX; i++){
		for(j=0; j<JMAX; j++){
		gxrho[i][j]=0.;
		gyrho[i][j]=0.;
		gxp[i][j]=0.;
		gyp[i][j]=0.;
		gxphi[i][j]=0.;
		gyphi[i][j]=0.;
		gxu[i][j]=0.;
		gyu[i][j]=0.;
		gxv[i][j]=0.;
		gyv[i][j]=0.;
		gxrhon[i][j]=0.;
		gyrhon[i][j]=0.;
		gxpn[i][j]=0.;
		gypn[i][j]=0.;
		gxphin[i][j]=0.;
		gyphin[i][j]=0.;
		gxun[i][j]=0.;
		gyun[i][j]=0.;
		gxvn[i][j]=0.;
		gyvn[i][j]=0.;
		}
		}
		/***** for gradient *****/
		for(i=2; i<IMAX-1; i++){
			for(j=2; j<JMAX-1; j++){
				ITVLM+=iphi(phi[i][j]);
			}
		}
		return(0);
}
int re_calc_output(double t, int count, double rho[IMAX][JMAX], double gxrho[IMAX][JMAX], double gyrho[IMAX][JMAX],
double rhon[IMAX][JMAX], double gxrhon[IMAX][JMAX], double gyrhon[IMAX][JMAX],
double u[IMAX][JMAX], double gxu[IMAX][JMAX], double gyu[IMAX][JMAX],
double un[IMAX][JMAX], double gxun[IMAX][JMAX], double gyun[IMAX][JMAX],
double v[IMAX][JMAX], double gxv[IMAX][JMAX], double gyv[IMAX][JMAX],
double vn[IMAX][JMAX], double gxvn[IMAX][JMAX], double gyvn[IMAX][JMAX],
double phi[IMAX][JMAX], double gxphi[IMAX][JMAX], double gyphi[IMAX][JMAX],
double phin[IMAX][JMAX], double gxphin[IMAX][JMAX], double gyphin[IMAX][JMAX],
double p[IMAX][JMAX], double gxp[IMAX][JMAX], double gyp[IMAX][JMAX],
double pn[IMAX][JMAX], double gxpn[IMAX][JMAX], double gypn[IMAX][JMAX]){
FILE *fp;
int i,j;
if ((fp=fopen ("./Data/ContCalcData.dat","w+")) == NULL){
printf("’./Data/ContCalcData.dat’ was not opened!\n");
printf("’Maybe you should make ’Data’ directory....\n");
exit(1);
}
fprintf(fp,"%10d\n",count);
fprintf(fp,"%16.8e\n",t);
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
fprintf(fp,"%16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e\n",
rho[i][j], gxrho[i][j], gyrho[i][j], rhon[i][j], gxrhon[i][j], gyrhon[i][j],
u[i][j], gxu[i][j], gyu[i][j], un[i][j]);
fprintf(fp,"%16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e\n",
gxun[i][j], gyun[i][j], v[i][j], gxv[i][j], gyv[i][j],
vn[i][j], gxvn[i][j], gyvn[i][j], phi[i][j], gxphi[i][j]);
fprintf(fp,"%16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e %16.8e\n",
gyphi[i][j], phin[i][j], gxphin[i][j], gyphin[i][j], p[i][j],
gxp[i][j], gyp[i][j], pn[i][j], gxpn[i][j], gypn[i][j]);
}
}
fclose(fp);
return(0);
}
int update(double rho[IMAX][JMAX], double rhon[IMAX][JMAX],
double u[IMAX][JMAX], double un[IMAX][JMAX],
double v[IMAX][JMAX], double vn[IMAX][JMAX],
double p[IMAX][JMAX], double pn[IMAX][JMAX],
double phi[IMAX][JMAX], double phin[IMAX][JMAX]){
int i,j;
double oldp[IMAX][JMAX];
/* -- for printf old pressure -- */
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
oldp[i][j]=p[i][j];
}
}
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
rho[i][j]=rhon[i][j];
if(rho[i][j]<=RHOair){
rho[i][j]=RHOair;
phi[i][j]=dphi(0.);
}
if(rho[i][j]>=RHOwater){
rho[i][j]=RHOwater;
phi[i][j]=dphi(1.);
}
phi[i][j]=phin[i][j];
if(iphi(phi[i][j])<0.){
rho[i][j]=RHOair;
phi[i][j]=dphi(0.);
}
if(iphi(phi[i][j])>1.){
rho[i][j]=RHOwater;
phi[i][j]=dphi(1.);
}
u[i][j]=un[i][j];
v[i][j]=vn[i][j];
p[i][j]=pn[i][j];
if(p[i][j]<=0.){
printf("(update) Negative Pressure %16.3e oldp=%16.3e at (%d,%d)\n",p[i][j],oldp[i][j],i,j);
printf("%16.3e %16.3e %16.3e\n",p[i-1][j+1],p[i][j+1],p[i+1][j+1]);
printf("%16.3e %16.3e %16.3e\n",p[i-1][j ],p[i][j ],p[i+1][j ]);
printf("%16.3e %16.3e %16.3e\n",p[i-1][j-1],p[i][j-1],p[i+1][j-1]);
printf("--- oldp ---\n");
printf("%16.3e %16.3e %16.3e\n",oldp[i-1][j+1],oldp[i][j+1],oldp[i+1][j+1]);
printf("%16.3e %16.3e %16.3e\n",oldp[i-1][j ],oldp[i][j ],oldp[i+1][j ]);
printf("%16.3e %16.3e %16.3e\n",oldp[i-1][j-1],oldp[i][j-1],oldp[i+1][j-1]);
return(-1);
}
}
}
return(0);
}
int check(double rho[IMAX][JMAX], double p[IMAX][JMAX], double tmp1[IMAX][JMAX], double tmp2[IMAX][JMAX]){
int i,j;
for(i=0; i<IMAX; i++){
for(j=0; j<JMAX; j++){
if(rho[i][j]<=0.){
printf("(check) Negative density %16.8e at (%d,%d)\n",rho[i][j],i,j);
printf(" %16d %16d %16d %16d %16d\n", i-2,i-1,i,i+1,i+2);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j+2,rho[i-2][j+2],rho[i-1][j+2],rho[i][j+2],rho[i+1][j+2],rho[i+2][j+2]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j+1,rho[i-2][j+1],rho[i-1][j+1],rho[i][j+1],rho[i+1][j+1],rho[i+2][j+1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j ,rho[i-2][j ],rho[i-1][j ],rho[i][j ],rho[i+1][j ],rho[i+2][j ]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j-1,rho[i-2][j-1],rho[i-1][j-1],rho[i][j-1],rho[i+1][j-1],rho[i+2][j-1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j-2,rho[i-2][j-2],rho[i-1][j-2],rho[i][j-2],rho[i+1][j-2],rho[i+2][j-2]);
printf("----- pressure -----\n");
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j+2,p[i-2][j+2],p[i-1][j+2],p[i][j+2],p[i+1][j+2],p[i+2][j+2]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j+1,p[i-2][j+1],p[i-1][j+1],p[i][j+1],p[i+1][j+1],p[i+2][j+1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j ,p[i-2][j ],p[i-1][j ],p[i][j ],p[i+1][j ],p[i+2][j ]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j-1,p[i-2][j-1],p[i-1][j-1],p[i][j-1],p[i+1][j-1],p[i+2][j-1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j-2,p[i-2][j-2],p[i-1][j-2],p[i][j-2],p[i+1][j-2],p[i+2][j-2]);
printf("----- tmp1 -----\n");
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j+2,tmp1[i-2][j+2],tmp1[i-1][j+2],tmp1[i][j+2],tmp1[i+1][j+2],tmp1[i+2][j+2]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j+1,tmp1[i-2][j+1],tmp1[i-1][j+1],tmp1[i][j+1],tmp1[i+1][j+1],tmp1[i+2][j+1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j ,tmp1[i-2][j ],tmp1[i-1][j ],tmp1[i][j ],tmp1[i+1][j ],tmp1[i+2][j ]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j-1,tmp1[i-2][j-1],tmp1[i-1][j-1],tmp1[i][j-1],tmp1[i+1][j-1],tmp1[i+2][j-1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j-2,tmp1[i-2][j-2],tmp1[i-1][j-2],tmp1[i][j-2],tmp1[i+1][j-2],tmp1[i+2][j-2]);
printf("----- tmp2 -----e\n");
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j+2,tmp2[i-2][j+2],tmp2[i-1][j+2],tmp2[i][j+2],tmp2[i+1][j+2],tmp2[i+2][j+2]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j+1,tmp2[i-2][j+1],tmp2[i-1][j+1],tmp2[i][j+1],tmp2[i+1][j+1],tmp2[i+2][j+1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j ,tmp2[i-2][j ],tmp2[i-1][j ],tmp2[i][j ],tmp2[i+1][j ],tmp2[i+2][j ]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j-1,tmp2[i-2][j-1],tmp2[i-1][j-1],tmp2[i][j-1],tmp2[i+1][j-1],tmp2[i+2][j-1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",
j-2,tmp2[i-2][j-2],tmp2[i-1][j-2],tmp2[i][j-2],tmp2[i+1][j-2],tmp2[i+2][j-2]);
exit(-1);
}
if(p[i][j]<=0.){
printf("(check) Negative Pressure %16.3e at (%d,%d)\n",p[i][j],i,j);
printf(" %16d %16d %16d %16d %16d\n", i-2,i-1,i,i+1,i+2);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j+2,p[i-2][j+2],p[i-1][j+2],p[i][j+2],p[i+1][j+2],p[i+2][j+2]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j+1,p[i-2][j+1],p[i-1][j+1],p[i][j+1],p[i+1][j+1],p[i+2][j+1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j ,p[i-2][j ],p[i-1][j ],p[i][j ],p[i+1][j ],p[i+2][j ]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j-1,p[i-2][j-1],p[i-1][j-1],p[i][j-1],p[i+1][j-1],p[i+2][j-1]);
printf("%3d : %16.8e %16.8e %16.8e %16.8e %16.8e\n",j-2,p[i-2][j-2],p[i-1][j-2],p[i][j-2],p[i+1][j-2],p[i+2][j-2]);
exit(-1);
}
}
}
return(0);
}
double sound_velocity(int i, int j, double rho[IMAX][JMAX], double p[IMAX][JMAX], double phi[IMAX][JMAX]){
	double cair2, cwater2, c2, csrho;
	//----- original -----
	// cair2=GAMMA*p[i][j]/rho[i][j];
	// cwater2=KAPPA*(p[i][j]+B)/rho[i][j];
	//----- improved -----
	cair2=GAMMA*p[i][j]/RHOair;
	cwater2=KAPPA*(p[i][j]+B)/RHOwater;
	//--------------------
	if(iphi(phi[i][j])<0. || iphi(phi[i][j])>1.){

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
国产精品资源网| 久久婷婷综合激情| 欧美精品一区二区三| 中文字幕亚洲区| 免费精品99久久国产综合精品| 国产成人亚洲综合a∨婷婷| 欧美日韩aaa| ...av二区三区久久精品| 国内久久精品视频| 日韩欧美一级二级三级久久久| 亚洲图片另类小说| 国产一区二区免费视频| 日韩视频免费观看高清完整版 | 成人黄色在线看| 日韩一区二区三| 亚洲va国产天堂va久久en| 不卡的av网站| 欧美国产在线观看| 国产在线播放一区二区三区| 欧美一区二区视频观看视频| 亚洲夂夂婷婷色拍ww47| aaa国产一区| 亚洲欧洲精品天堂一级| 丰满白嫩尤物一区二区| 国产婷婷色一区二区三区在线| 久久99精品久久久久婷婷| 欧美日韩精品欧美日韩精品一 | 亚洲激情欧美激情| 97精品视频在线观看自产线路二| 久久九九99视频| 国产精品一区二区久久不卡| 精品国产乱码久久| 国内外成人在线| 国产人伦精品一区二区| 国产suv精品一区二区6| 亚洲国产精品ⅴa在线观看| 国产成人自拍高清视频在线免费播放| 精品国产乱码久久久久久影片| 蜜乳av一区二区三区| 日韩精品一区二区三区在线播放 | 日韩精品影音先锋| 黄网站免费久久| 国产三级欧美三级日产三级99| 国产毛片一区二区| 中文字幕一区二区三区色视频| caoporm超碰国产精品| 一级特黄大欧美久久久| 欧美三级中文字幕在线观看| 午夜精品成人在线视频| 日韩欧美一区在线观看| 国产黄色精品视频| 亚洲欧洲美洲综合色网| 色噜噜狠狠成人中文综合| 亚洲无线码一区二区三区| 51精品国自产在线| 国产一区二区三区免费| ...xxx性欧美| 日韩一区二区三区电影| 国产成人丝袜美腿| 亚洲图片欧美色图| 精品国产伦理网| 色婷婷综合久久久中文字幕| 亚洲一区av在线| 精品成人a区在线观看| 波多野结衣中文字幕一区二区三区| 亚洲免费在线观看| 日韩免费看的电影| 色香色香欲天天天影视综合网| 日韩和欧美的一区| 国产精品美女久久久久久久久久久| 欧美亚洲另类激情小说| 国产中文字幕一区| 亚洲最色的网站| 久久日韩粉嫩一区二区三区| 91麻豆国产自产在线观看| 久久成人麻豆午夜电影| 亚洲综合小说图片| 精品盗摄一区二区三区| 欧美日韩高清一区| 不卡av电影在线播放| 久久99久久久久久久久久久| 亚洲你懂的在线视频| 精品蜜桃在线看| 在线免费观看日韩欧美| 成人一区二区三区在线观看| 日本在线不卡一区| 亚洲激情五月婷婷| 国产欧美日韩综合精品一区二区| 91麻豆精品国产91久久久久久| av中文字幕不卡| 国产激情精品久久久第一区二区| 日韩va亚洲va欧美va久久| 一区二区三区四区av| 中文字幕亚洲区| 欧美激情一区二区三区| 日韩欧美激情在线| 91精品视频网| 欧美日韩另类国产亚洲欧美一级| 972aa.com艺术欧美| 国产成人精品aa毛片| 久久国产尿小便嘘嘘尿| 三级在线观看一区二区| 亚洲高清中文字幕| 亚洲一区二区av电影| 亚洲图片自拍偷拍| 亚洲午夜一区二区三区| 亚洲精品高清在线观看| 中文字幕中文乱码欧美一区二区| 国产亚洲va综合人人澡精品| 精品成人一区二区三区| 久久先锋影音av| 国产日韩v精品一区二区| 国产欧美精品一区二区三区四区 | 91精品国产91综合久久蜜臀| 欧美色综合久久| 欧美在线播放高清精品| 91九色最新地址| 欧美日本乱大交xxxxx| 欧美另类z0zxhd电影| 91精品国产麻豆| 精品乱码亚洲一区二区不卡| 久久久久久久久久久久久久久99| 精品国产一区二区三区久久久蜜月| 欧美成人bangbros| 日韩欧美国产一区二区三区| 精品国产乱码久久久久久牛牛| 2014亚洲片线观看视频免费| 国产视频在线观看一区二区三区| 国产喂奶挤奶一区二区三区| 中文字幕不卡的av| 亚洲三级免费电影| 亚洲小说春色综合另类电影| 日韩成人av影视| 国产黄色精品网站| 日本精品视频一区二区| 欧美日韩国产欧美日美国产精品| 日韩一区二区精品| 欧美成人伊人久久综合网| 日韩视频不卡中文| 国产精品网曝门| 亚洲电影在线播放| 久久99精品视频| 91丨九色丨蝌蚪丨老版| 欧美日韩精品免费观看视频| 日韩视频免费观看高清完整版| 国产精品视频你懂的| 亚洲一二三四在线| 国产一区二区三区免费播放| 91啪亚洲精品| 26uuu亚洲综合色| 一区二区三区免费在线观看| 久久爱www久久做| 色老汉av一区二区三区| 精品国产91亚洲一区二区三区婷婷| 国产精品污www在线观看| 午夜电影一区二区| www.综合网.com| 日韩欧美一级在线播放| 亚洲人吸女人奶水| 狠狠色丁香婷婷综合| 欧美日精品一区视频| 国产欧美一区二区在线| 爽爽淫人综合网网站| 成人av在线影院| 精品久久人人做人人爱| 亚洲女同一区二区| 国产成人精品三级麻豆| 欧美精品视频www在线观看| 中文字幕中文字幕在线一区| 久久99精品网久久| 欧美精品第一页| 亚洲自拍与偷拍| 99天天综合性| 国产欧美一区视频| 久久精品久久99精品久久| 欧美亚洲国产一区在线观看网站| 久久精品人人做人人爽人人| 青草国产精品久久久久久| 在线中文字幕一区| 综合网在线视频| 成人免费毛片嘿嘿连载视频| 欧美电影免费观看高清完整版| 五月天一区二区三区| 91美女片黄在线| 中文字幕日韩av资源站| 国产美女娇喘av呻吟久久| 日韩一本二本av| 日韩高清国产一区在线| 欧美精品自拍偷拍动漫精品| 亚洲一区二区三区视频在线播放| 99国产精品久久久久久久久久 | 一区二区三区在线观看国产| 风间由美一区二区av101| 精品久久久久香蕉网| 日本女优在线视频一区二区| 91精选在线观看| 青娱乐精品视频在线| 精品女同一区二区| 欧美色国产精品| 亚洲欧美日韩久久|