sample0.for

给予DSMC的二维CFD模拟FORTRAN代码

*   SAMPLE0.FOR
*
      SUBROUTINE SAMPLE0
*
*--samples the molecules in the flow
*
         INCLUDE 'dsmc0.h'
*
      NSMP=NSMP+1
      DO 100 NN=1,MNSG
        DO 50 N=1,MNC
          L=IC(2,N,NN)
          IF (L.GT.0) THEN
            DO 10 J=1,L
              K=IC(1,N,NN)+J
              M=IR(K)
              I=IPS(M)
              CS(1,N,I)=CS(1,N,I)+1
              DO 5 LL=1,3
                CS(LL+1,N,I)=CS(LL+1,N,I)+PV(LL,M)
                CS(LL+4,N,I)=CS(LL+4,N,I)+PV(LL,M)**2
5             CONTINUE
10          CONTINUE
          END IF
50      CONTINUE
100   CONTINUE
      RETURN
      END