area.f90

使用WENO格式求解流体力学中的激波问题

FUNCTION AREA(X1,Y1,X2,Y2,X3,Y3)
IMPLICIT NONE
DOUBLE PRECISION AREA,X1,X2,X3,Y1,Y2,Y3
AREA=ABS(X1*Y2-Y1*X2+X2*Y3-X3*Y2+X3*Y1-X1*Y3)
area=area/2.0
END FUNCTION AREA