|
|
ÊÖÅú¤¢¤ê¤¬¤È¤¦¤´¤¶¤¤¤Þ¤¹¡£
Æñ¤·¤¯¤Æ»Í¶ìȬ¶ì¤·¤Æ¤Ð¤«¤ê¤Ç¤¹¡£
¶µ¤¨¤Æ¤¤¤¿¤À¤¤¤¿¥×¥í¥°¥é¥à¤Ç¼Â¹Ô¤¹¤ë¤È¡¤¹³ÎÏfd¤¬¥Þ¥¤¥Ê¥¹¤Ë¤Ê¤Ã¤Æ¤·¤Þ¤¦¤Î¤Ç¤¹¤¬
²¿¤¬¤¤¤±¤Ê¤¤¤Î¤Ç¤·¤ç¤¦¡Ä¹³ÎϤ¬¥Þ¥¤¥Ê¥¹¤Ë¤Ê¤ë¤³¤È¤Ï¤¢¤êÆÀ¤ë¤Î¤Ç¤¹¤«¡©
¤Á¤Ê¤ß¤Ë»ä¤Ï±²ÅÙ¤ò½¬¤Ã¤Æ¤¤¤Ê¤«¤Ã¤¿¤Î¤Ç±²ÅÙ¤ò»È¤ï¤º°Ê²¼¤Î¤è¤¦¤Ê¥×¥í¥°¥é¥à¤Ç²ò¤¤¤Æ¤¤¤Þ¤·¤¿¡£
SUMX=0.
SUMY=0.
DO I=2,MX
DX = 0.5*SQRT((X(I+1,1)-X(I-1,1))**2+(Y(I+1,1)-Y(I-1,1))**2)
UXI=0.5*(U(I+1,1)-U(I-1,1))
VXI=0.5*(V(I+1,1)-V(I-1,1))
UETA=0.5*(U(I,2)-UB1J(I))
VETA=0.5*(V(I,2)-VB1J(I))
c UB1J(I),VB1JʪÂÎÆâ²¾ÁÛή®
YETA=YE(I,1)
YXI=YX(I,1)
XETA=XE(I,1)
XXI=XX(I,1)
AJAC=AJ(I,1)
PXX=(-P(I,J)+2.*(YETA*UXI-YXI*UETA)/AJAC*REI)*HX(I,1)
PYX=((XETA*UXI-XXI*UETA)+(YETA*VXI-YXI*VETA))/AJAC*REI*HY(I,1)
PXY=((XETA*UXI-XXI*UETA)+(YETA*VXI-YXI*VETA))/AJAC*REI*HX(I,1)
PYY=(-P(I,J)+2.*(XXI*VETA+XETA*VXI)/AJAC*REI)*HY(I,1)
c n(HX(I,1),HY(i,2))¢ªË¡Àþ¥Ù¥¯¥È¥ë
SUMX = SUMX + (PXX+PYX)*DX
SUMY = SUMY + (PXY+PYY)*DX
END DO
²¿ÅÙ¤âËÜÅö¤Ë¤¹¤ß¤Þ¤»¤ó¡£
|
|
