fluid_module.f90

visual fortran 源码

!-----------------------常量模块-----------------------------------------
module const
integer,parameter::NI=7,NJ=12,NIJ=12,L1=7,L2=6,L3=5,M1=12,M2=11,M3=10,NFMAX=10,NFX3=13
end module const

!-----------------------变量声明模块-----------------------------------------
module varible
use const 
!T,EQUIVLENCE(F,T)
!F(NI,NJ,NFMAX-->F(NI,NJ,NFX3
!
!
!
real,SAVE:: AIP(NIJ,NIJ),AIM(NIJ,NIJ),AJP(NIJ,NIJ),AJM(NIJ,NIJ),AP(NIJ,NIJ)
real,SAVE :: X(NIJ),Y(NIJ),XU(NIJ),YV(NIJ),XDIF(NIJ),YDIF(NIJ),XCV(NIJ),YCV(NIJ),XCVS(NIJ),YCVS(NIJ)
real,SAVE :: F(NIJ,NIJ,NFX3),P(NIJ,NIJ),RHO(NIJ,NIJ),GAM(NIJ,NIJ),CON(NIJ,NIJ)
real,SAVE :: YCVR(NIJ),YCVRS(NIJ),ARX(NIJ),ARXJ(NIJ),ARXJP(NIJ),R(NIJ),RMN(NIJ),SX(NIJ),SXMN(NIJ)
real,SAVE :: XCVI(NIJ),XCVIP(NIJ),DU(NIJ,NIJ),DV(NIJ,NIJ),FV(NIJ),FVP(NIJ),FX(NIJ),FXM(NIJ)
real,SAVE :: FY(NIJ),FYM(NIJ),PT(NIJ),QT(NIJ)
real,SAVE :: NF,NP,NRHO,NGAM,RELAX(NFX3),temp
integer,SAVE :: IST,JST,ITER,LAST,IPREF,JPREF,MODE
logical,SAVE :: LSOLVE(NFX3),LBLK(NFX3),LSTOP
real,SAVE ::TIME,DT,XL,YL,NTIMES(NFX3),RHOCON
real,SAVE:: SMAX,SSUM,FLOW,DIFF,ACOF
real,SAVE ::U(NIJ,NIJ),V(NIJ,NIJ),PC(NIJ,NIJ),COF(NIJ,NIJ,6),COFU(NIJ,NIJ,6),COFV(NIJ,NIJ,6),COFP(NIJ,NIJ,6)
!==========================================
real,save::T(NIJ,NIJ),w(NIJ,NIJ)
equivalence (F(1,1,4),w(1,1)),(F(1,1,5),T(1,1))
!==========================================
equivalence (F(1,1,1),U(1,1)),(F(1,1,2),V(1,1)),(F(1,1,3),PC(1,1))
equivalence (COF(1,1,1),CON(1,1))
data LSTOP,LSOLVE,LBLK/.FALSE.,NFX3*.FALSE.,NFX3*.FALSE./
data RELAX,NTIMES/NFX3*1.0,NFX3*100/
data DT,IPREF,JPREF,RHOCON/1.0E30,1,1,1.0/
end module varible

module non_common
use const
real :: APT,VOL,GM,GMM,SXT,SXB,DX,RHOM
integer :: I,J,K
end module non_common