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

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

?? 2.txt

?? 基于CIP的海嘯波浪數值模擬
?? TXT
?? 第 1 頁 / 共 5 頁
字號: 
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){
/***** 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 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.){
printf("(sound_velocity) INVALID iphi value (=%16.8e) at (%d,%d)\n",iphi(phi[i][j]),i,j);
exit(-1);
}
c2=(1.-iphi(phi[i][j]))*cair2 + iphi(phi[i][j])*cwater2;
if(c2<=0.){
printf("Negative Sound Speed =%16.8e at (%2d,%2d), iphi=%16.8e, cair2=%16.8e, cwater2=%16.8e\n",
c2, i,j,iphi(phi[i][j]),cair2, cwater2);
exit(-1);
}
return(c2);
}
int bc_p(double pn[IMAX][JMAX]){
int i,j;
/********** B.C. for Pressure & Sound speed *****************************/
/* LEFT side */
/*... outlet ...*/
for(j=2; j<=OUTLET+1; j++){
pn[0 ][j]=pn[2 ][j];
pn[1 ][j]=pn[2 ][j];
}
/*... wall ...*/
for(j=OUTLET+2; j<=JMAX-3; j++){
pn[0 ][j]=pn[3 ][j];
pn[1 ][j]=pn[2 ][j];

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
午夜国产不卡在线观看视频| 成人动漫精品一区二区| 国产一区二区不卡在线| 成人免费毛片a| 欧美一区二区三区小说| 国产精品不卡视频| 蜜桃一区二区三区四区| 欧洲精品视频在线观看| 久久午夜老司机| 五月综合激情婷婷六月色窝| 成人一区二区在线观看| 日韩欧美成人激情| 亚洲一区二区美女| 99re在线视频这里只有精品| 日韩欧美国产一区二区在线播放| 亚洲欧美电影院| 国产成人精品免费在线| 欧美成人vps| 毛片不卡一区二区| 欧美网站大全在线观看| 亚洲男人的天堂在线观看| 国产精品456| 久久综合五月天婷婷伊人| 三级欧美在线一区| 欧美三级视频在线观看| 亚洲欧美激情一区二区| 色综合亚洲欧洲| 日韩欧美国产高清| 一区二区激情小说| 色综合天天综合网国产成人综合天| 国产蜜臀av在线一区二区三区 | 欧美放荡的少妇| 亚洲一区免费在线观看| 色老汉av一区二区三区| 亚洲免费观看高清在线观看| av不卡免费电影| 中文字幕在线免费不卡| 99精品桃花视频在线观看| 国产日产欧美精品一区二区三区| 国产乱码字幕精品高清av| 久久这里只有精品视频网| 久久er99精品| 久久久精品国产99久久精品芒果| 国产又黄又大久久| 国产视频911| 北岛玲一区二区三区四区| 亚洲欧洲性图库| 色香蕉久久蜜桃| 亚洲午夜精品网| 欧美一三区三区四区免费在线看 | 一区二区三区中文在线观看| 一本久道久久综合中文字幕| 一区二区三区高清不卡| 91精品国产美女浴室洗澡无遮挡| 久久99久久精品| 久久久精品国产免大香伊| 高清beeg欧美| 一片黄亚洲嫩模| 欧美一区二区三区在| 国产aⅴ综合色| 亚洲码国产岛国毛片在线| 欧美日韩大陆在线| 国产伦理精品不卡| 综合分类小说区另类春色亚洲小说欧美| 波多野结衣在线一区| 亚洲第一在线综合网站| 欧美xxxx在线观看| 91在线视频免费91| 日本在线播放一区二区三区| 国产色产综合产在线视频| 一道本成人在线| 久久国内精品视频| 自拍偷拍国产亚洲| 日韩欧美国产综合在线一区二区三区| 国产宾馆实践打屁股91| 亚洲香蕉伊在人在线观| www一区二区| 色噜噜久久综合| 国产主播一区二区| 亚洲成精国产精品女| 精品国内片67194| 色综合天天综合| 韩国毛片一区二区三区| 亚洲综合激情小说| 日本一区二区三区在线观看| 欧美理论片在线| 成人免费黄色大片| 日本亚洲电影天堂| 亚洲精品国产成人久久av盗摄| 日韩一级在线观看| 色综合天天性综合| 国产传媒一区在线| 首页国产丝袜综合| 亚洲色图欧美偷拍| 欧美极品另类videosde| 日韩女同互慰一区二区| 色屁屁一区二区| 北条麻妃一区二区三区| 韩国成人福利片在线播放| 日韩福利电影在线观看| 亚洲免费电影在线| 欧美国产精品一区二区三区| 日韩免费观看高清完整版在线观看| 91蜜桃网址入口| 成人国产精品免费观看| 国产麻豆精品在线| 国产尤物一区二区| 九九国产精品视频| 美女高潮久久久| 麻豆精品精品国产自在97香蕉| 亚洲已满18点击进入久久| 亚洲欧美影音先锋| 国产精品久久久久久久第一福利 | 欧美丝袜丝nylons| 北条麻妃一区二区三区| 成人app软件下载大全免费| 懂色av一区二区三区蜜臀| 国产精品一线二线三线精华| 韩国精品一区二区| 国产suv精品一区二区6| 成人免费看黄yyy456| 99久久精品免费看国产免费软件| 成人精品免费看| 97久久精品人人做人人爽| 91视频com| 欧美揉bbbbb揉bbbbb| 欧美一区二区视频在线观看| 日韩欧美亚洲国产另类| 久久精品日韩一区二区三区| 欧美国产在线观看| 成人免费小视频| 亚洲一区精品在线| 青青青伊人色综合久久| 国产精品一区在线观看你懂的| 国产福利一区在线观看| 99精品在线观看视频| 在线观看免费亚洲| 欧美日韩国产高清一区二区三区| 欧美三级电影一区| 9191成人精品久久| 2020国产精品自拍| 亚洲国产精品国自产拍av| 日本一二三不卡| 亚洲精品国产品国语在线app| 性做久久久久久久久| 日日骚欧美日韩| 捆绑紧缚一区二区三区视频| 麻豆一区二区99久久久久| a亚洲天堂av| 欧美日韩精品一区二区三区四区 | 国产成人午夜视频| 色综合天天综合网天天看片| 欧美色视频一区| 久久精品免视看| 中文字幕日本不卡| 亚洲综合一区二区精品导航| 日韩精品一级中文字幕精品视频免费观看| 日韩制服丝袜av| 国产真实乱对白精彩久久| 色哟哟欧美精品| 欧美日韩亚洲综合在线| 精品卡一卡二卡三卡四在线| 久久蜜桃av一区精品变态类天堂| 伊人色综合久久天天| 日本欧美在线观看| 国产不卡免费视频| 欧美三级三级三级爽爽爽| 国产精品女人毛片| 五月婷婷综合激情| 国产suv精品一区二区6| 欧美一级在线观看| 中文字幕在线观看不卡视频| 三级一区在线视频先锋| 岛国精品一区二区| 欧美videossexotv100| 久久综合成人精品亚洲另类欧美 | 国产精品视频一二| 奇米色一区二区三区四区| 色综合欧美在线视频区| 日韩精品影音先锋| 亚洲免费在线观看视频| 丁香桃色午夜亚洲一区二区三区| 欧美三级乱人伦电影| 国产视频一区在线观看| 三级不卡在线观看| 91福利视频网站| 中文字幕精品在线不卡| 捆绑调教一区二区三区| 岛国精品在线播放| 国产欧美va欧美不卡在线| 蜜臀99久久精品久久久久久软件| 色综合av在线| 亚洲另类色综合网站| 国产成人午夜精品5599| 欧美一级艳片视频免费观看| 亚洲精品中文字幕在线观看| 成人av动漫网站| 久久久精品欧美丰满| 蜜臀av性久久久久蜜臀aⅴ四虎| 国产91对白在线观看九色|