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

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

?? 555.cpp.bak

?? 基于CIP的海嘯波浪數(shù)值模擬
?? BAK
?? 第 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一区二区三区免费野_久草精品视频
国产欧美精品一区aⅴ影院| 亚洲精品国产一区二区精华液| 国产日韩欧美高清| 亚洲制服丝袜一区| 国产精品18久久久久久久久久久久| 色婷婷亚洲综合| 中文字幕乱码亚洲精品一区| 青青草原综合久久大伊人精品 | 欧美成人aa大片| 亚洲欧美国产77777| 国产一区二区免费看| 欧美嫩在线观看| 亚洲精品视频在线看| 成人动漫中文字幕| 欧美zozozo| 青青草原综合久久大伊人精品 | 日韩精品影音先锋| 午夜欧美大尺度福利影院在线看| 成人精品小蝌蚪| www国产成人| 美女视频黄久久| 91精品国产乱码| 五月天丁香久久| 在线播放欧美女士性生活| 亚洲欧美日韩国产综合| 成人性视频免费网站| 国产午夜精品一区二区| 精品一区二区三区在线播放| 日韩一区二区视频| 奇米精品一区二区三区在线观看 | 色综合咪咪久久| 日韩一区欧美一区| 91欧美激情一区二区三区成人| 国产精品久久久一本精品 | 国产一本一道久久香蕉| 精品国产乱码久久久久久1区2区| 日韩精品视频网| 欧美电影免费观看高清完整版在线 | 国产99精品视频| 国产亚洲精品久| 成人手机电影网| 欧美国产综合一区二区| 成人精品一区二区三区四区| 中文字幕不卡的av| 99久久伊人网影院| 亚洲欧美国产毛片在线| 欧美日韩小视频| 蜜桃免费网站一区二区三区| 久久嫩草精品久久久精品| 国产精品 欧美精品| 中文字幕日韩av资源站| 欧美日韩综合在线| 美日韩一级片在线观看| 国产校园另类小说区| 91亚洲精华国产精华精华液| 亚洲一区二区三区小说| 欧美成人精品1314www| 成人一区二区视频| 亚洲影院在线观看| 亚洲精品在线电影| 91麻豆精品视频| 视频一区视频二区中文| 久久久久一区二区三区四区| 99精品欧美一区| 蜜桃91丨九色丨蝌蚪91桃色| 国产精品美女久久久久aⅴ| 色天天综合久久久久综合片| 欧美aaaaaa午夜精品| 国产精品久久久久久久岛一牛影视 | 一本大道久久a久久精品综合| 水野朝阳av一区二区三区| 国产性天天综合网| 欧美日韩国产片| 成人av在线网站| 奇米影视7777精品一区二区| 日韩美女精品在线| 26uuu国产在线精品一区二区| 色综合中文字幕| 国产精品99久久久| 青椒成人免费视频| 亚洲一卡二卡三卡四卡| 久久精品亚洲精品国产欧美| 欧美欧美欧美欧美首页| 波多野结衣中文字幕一区| 麻豆精品视频在线| 亚洲国产中文字幕在线视频综合| 久久精品视频在线看| 欧美一区二区三区四区视频| 一本久久a久久精品亚洲| 国产一区中文字幕| 美女精品一区二区| 亚洲va天堂va国产va久| 亚洲欧洲精品天堂一级| 久久香蕉国产线看观看99| 678五月天丁香亚洲综合网| aaa欧美日韩| 夫妻av一区二区| 精东粉嫩av免费一区二区三区| 性做久久久久久久久| 亚洲精品视频观看| 亚洲色图欧美在线| 国产精品毛片大码女人| 久久免费看少妇高潮| 欧美电影免费观看高清完整版在线| 欧美在线高清视频| 一本大道久久a久久综合| av影院午夜一区| 国产在线国偷精品产拍免费yy| 日本亚洲三级在线| 日本大胆欧美人术艺术动态 | 国产在线不卡一区| 久久国产精品99久久人人澡| 日韩一区精品视频| 日精品一区二区| 日本va欧美va瓶| 美国精品在线观看| 国产在线看一区| 国产精品69久久久久水密桃| 国产乱码精品一区二区三区av| 国产精品伊人色| 成人看片黄a免费看在线| 国产成人免费视频网站高清观看视频| 国产一二三精品| 成人一区二区三区视频在线观看 | 欧美一区二区三区视频在线观看 | 欧美日韩另类国产亚洲欧美一级| 欧美日韩在线综合| 宅男噜噜噜66一区二区66| 91精品国产综合久久国产大片| 日韩一区二区三区精品视频| 日韩区在线观看| 中文字幕第一区| 一区二区久久久| 日韩黄色一级片| 国产精品一区在线观看乱码| av电影在线观看一区| 欧美日韩一卡二卡| 精品国产亚洲在线| 国产精品不卡一区二区三区| 亚洲一区二三区| 韩国视频一区二区| 色综合久久中文综合久久牛| 这里只有精品免费| 欧美激情中文不卡| 天天影视涩香欲综合网| 国产精品一区二区你懂的| 91一区一区三区| 欧美成人猛片aaaaaaa| 1000部国产精品成人观看| 五月婷婷久久综合| 懂色中文一区二区在线播放| 欧洲激情一区二区| 久久久午夜电影| 亚洲va欧美va国产va天堂影院| 狠狠色综合日日| 在线观看免费视频综合| 久久众筹精品私拍模特| 一区二区三区不卡视频| 精品一区二区在线看| 欧洲人成人精品| 中文字幕在线视频一区| 美女视频免费一区| 欧美艳星brazzers| 国产日韩av一区二区| 日本三级亚洲精品| 在线免费一区三区| 亚洲国产精品成人综合 | 亚洲777理论| 国产91丝袜在线播放九色| 777午夜精品免费视频| 亚洲图片你懂的| 国产成人精品在线看| 日韩视频一区二区三区在线播放 | 一色屋精品亚洲香蕉网站| 激情成人综合网| 日韩欧美国产一区二区三区| 亚洲综合免费观看高清完整版 | 日本高清成人免费播放| 国产午夜一区二区三区| 久久福利资源站| 欧美一区二区三区色| 亚洲午夜在线观看视频在线| 91免费看`日韩一区二区| 国产日韩三级在线| 韩国欧美一区二区| 日韩精品一区国产麻豆| 青青国产91久久久久久| 欧美一区二区三区视频在线观看| 亚洲午夜久久久久| 91福利精品视频| 专区另类欧美日韩| 91啪在线观看| 亚洲精品菠萝久久久久久久| 97久久超碰国产精品| **欧美大码日韩| 91视频.com| 亚洲风情在线资源站| 91黄色激情网站| 亚洲电影一区二区三区| 欧美喷潮久久久xxxxx|