29_41_1.f

cfd经典PISO算法的FORTRAN程序

      parameter( nx =42, ny= 42, nxy=nx*ny)
      parameter( nxm=nx-1, nym=ny-1, nxm2=nx-2, nym2=ny-2)
      parameter( itimmax=1, itermax=450)
      parameter (ipr=2, jpr=2)

      common/fieldn/u(nxy),v(nxy),t(nxy),p(nxy),pp(nxy)
      common/fieldo/uo(nxy),vo(nxy),to(nxy)
      common/coeffs/ap(nxy),apu(nxy),apv(nxy),su(nxy),sv(nxy) 
      common/coeff1/ae1(nxy),aw1(nxy),as1(nxy),an1(nxy) 
      common/coeff2/ae2(nxy),aw2(nxy),as2(nxy),an2(nxy) 
      common/coeff3/ae3(nxy),aw3(nxy),as3(nxy),an3(nxy) 
      common/piso2/ uprime(nxy),vprime(nxy),utilde(nxy),vtilde(nxy)
      common/pgrad/dpx(nxy),dpy(nxy)
      common/grids/dx,dy,vol,vdydx,vdxdy,adydx,adxdy,apt,li(nx)
      common/mflux/flxe(nxy),flxw(nxy),flxs(nxy),flxn(nxy) 
      common/mflux2/flxe2(nxy),flxw2(nxy),flxs2(nxy),flxn2(nxy)

      common/init1/rho,vis,prm,grav,beta,th,tc,tref,ulid,xl,yl,dt
      common/init2/re,rei,gr,ra,grre,prre
      common/convg/resu(4),rest(4),urf(4),srmax,srmin,fnxy
      common/monitor/iter, itim, time
      common/flags/lcalcv,lcalcp,lcalct
      logical lcalcv,lcalcp,lcalct,lmoni,lunst