fluid_supply.f90

visual fortran 源码

subroutine supply()
use varible
use non_common

ENTRY UGRID()
XU(2)=0
DX=XL/(L1-2)
do I=3,L1
XU(I)=XU(I-1)+DX
end do
YV(2)=0
DY=YL/(M1-2)
do J=3,M1
YV(J)=YV(J-1)+DY
end do
return

ENTRY OTHERS()
do I=2,L1
if(I==2)then
F(2,2,3)=0
else
F(I,2,3)=F(I-1,2,3)-RHO(I-1,1)*V(I-1,2)*R(1)*XCV(I-1)
end if
  do J=3,M1
  RHOM=FX(I)*RHO(I,J-1)+FXM(I)*RHO(I-1,J-1)
  F(I,J,3)=F(I,J-1,3)+RHOM*U(I,J-1)*ARX(J-1)
  end do
continue
end do

do J=2,M2
P(1,J)=(P(2,J)*XCVS(3)-P(3,J)*XDIF(2))/XDIF(3)
P(L1,J)=(P(L2,J)*XCVS(L2)-P(L3,J)*XDIF(L1))/XDIF(L2)
end do
do I=2,L2
P(I,1)=(P(I,2)*YCVS(3)-P(I,3)*YDIF(2))/YDIF(3)
P(I,M1)=(P(I,M2)*YCVS(M2)-P(I,M3)*YDIF(M1))/YDIF(M2)
end do
P(1,1)=P(2,1)+P(1,2)-P(2,2)
P(L1,1)=P(L2,1)+P(L1,2)-P(L2,2)
P(1,M1)=P(2,M1)+P(1,M2)-P(2,M2)
P(L1,M1)=P(L2,M1)+P(L1,M2)-P(L2,M2)
PREF=P(IPREF,JPREF)
do J=1,M1
   do I=1,L1
   P(I,J)=P(I,J)-PREF
   end do
end do

end subroutine supply