sampi0.for

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

*   SAMPI0.FOR
*
      SUBROUTINE SAMPI0
*
*--initialises all the sampling variables
*
         INCLUDE 'dsmc0.h'
*
      NCOL=0
      NSMP=0
      MOVT=0.
      SELT=0.
      SEPT=0.
      DO 100 L=1,MNSP
        DO 50 N=1,MNC
          CS(1,N,L)=1.E-6
          DO 20 M=2,7
            CS(M,N,L)=0.
20        CONTINUE
50      CONTINUE
100   CONTINUE
      DO 200 M=1,MNSP
        DO 150 N=1,MNSP
          COL(M,N)=0.
150     CONTINUE
200   CONTINUE
 
	DO K=1,2
	 DO I=1,3
	   SMDDIS(K,I)=0.0D0
	 ENDDO
	   SMDTIME(K)=0.0D0
	 DO I=1,9
	   SMCSS(I,K,1)=0.0D0
	 ENDDO
	ENDDO

      RETURN
      END