program main implicit none integer,parameter :: MX=41 integer,parameter :: MY=41 integer,parameter :: MZ=1 real,parameter :: RE=10000.0 real,parameter :: DT=0.001 real,parameter :: XX=1.0 real,parameter :: YY=1.0 real,parameter :: EPS=1.0d-5 integer,parameter :: NMAX=50000 integer,parameter :: NPMAX=100 real,dimension(MX,MY) :: U,V,UT,VT,P,RHS real,dimension(MX) :: X,XX1,XX2,XX3 real,dimension(MY) :: Y,YY1,YY2,YY3 real DX,DY,PP,err,A,B,UN,VN,Pmin integer i,j,k,n,np !Irregular intervals Grid B=0.95 !0