fields.f90

出场插值程序

!-----------------------------------------------------------
!             二维流场插值，吴志根编写，版权所有!
!----------------------------------------------------------
	PARAMETER(L1=102)						  ! 网格
	PARAMETER(M1=102)
	PARAMETER(NN=10)
	PARAMETER(NI=NN-1,NJ=1)
!-----------------------------------------------------------		
	REAL U(L1,M1,NN)
	REAL V(L1,M1,NN)
	REAL P(L1,M1,NN)
!---------------------------------------------------------	
	REAL X(NI)
	REAL Y(NI)
	REAL B(NJ)
	REAL A(NI,NJ)
	REAL S(NJ)
	REAL TC
!----------------------------------------//插值系数求解
	DATA X/1,1.25,1.5,1.75,2,2.25,2.5,2.75,3/		 !  Log10(Re)
	DATA B/80/										 !  Re
		  
	DO J=1,NJ
	S(J)=1.0
	DO I=1,NI
	S(J)=S(J)*(B(J)-X(I))
	ENDDO
	ENDDO

	TC=0.0
	OPEN(55,FILE='COF.DAT') 
	DO K=1,NJ
	DO J=1,NI
	Y(J)=1.0
	DO I=1,NI
	IF(I.NE.J) Y(J)=Y(J)*(X(J)-X(I))
	ENDDO
	A(J,K)=S(k)/(B(K)-X(J))/Y(J)
	TC=TC+ABS(A(J,K))
	WRITE(*,*) A(J,K)
	WRITE(55,*) A(J,K)
	ENDDO
	WRITE(*,*)
	ENDDO
	WRITE(*,*)TC
	CLOSE(55) 
!--------------------------------------------------------------------
		OPEN(11,FILE='1.DAT')		!已计算流场
		OPEN(12,FILE='1.25.DAT') 
		OPEN(13,FILE='1.5.DAT')
		OPEN(14,FILE='1.75.DAT')
		OPEN(15,FILE='2.DAT')
		OPEN(16,FILE='2.25.DAT')
		OPEN(17,FILE='2.5.DAT')
		OPEN(18,FILE='2.75.DAT') 						
		OPEN(19,FILE='3.DAT')

		OPEN(20,FILE='OUT.DAT')			
!--------------------------------------------------------------------
		DO I=1,L1
		DO J=1,M1
		READ(11,*)U(I,J,1),V(I,J,1),P(I,J,1)
		READ(12,*)U(I,J,2),V(I,J,2),P(I,J,2)
		READ(13,*)U(I,J,3),V(I,J,3),P(I,J,3)
		READ(14,*)U(I,J,4),V(I,J,4),P(I,J,4)
		READ(15,*)U(I,J,5),V(I,J,5),P(I,J,5)
		READ(16,*)U(I,J,6),V(I,J,6),P(I,J,6)
		READ(17,*)U(I,J,7),V(I,J,7),P(I,J,7)
		READ(18,*)U(I,J,8),V(I,J,8),P(I,J,8)
		READ(19,*)U(I,J,9),V(I,J,9),P(I,J,9) 
!--------------------------------------
		DO JJ=1,NJ
		U(I,J,NI+JJ)=0.0
		V(I,J,NI+JJ)=0.0
		P(I,J,NI+JJ)=0.0
		DO II=1,NI
		U(I,J,NI+JJ)=U(I,J,NI+JJ)+A(II,JJ)*U(I,J,II)
		V(I,J,NI+JJ)=V(I,J,NI+JJ)+A(II,JJ)*V(I,J,II)
		P(I,J,NI+JJ)=P(I,J,NI+JJ)+A(II,JJ)*P(I,J,II)
		ENDDO
		ENDDO
!---------------------------------------
		WRITE(20,*)U(I,J,10),V(I,J,10),P(I,J,10)
!---------------------------------------
		ENDDO
		ENDDO

		CLOSE(11)
		CLOSE(12)
		CLOSE(13)
		CLOSE(14)
		CLOSE(15)
		CLOSE(16)
		CLOSE(17)
		CLOSE(18)
		CLOSE(19) 
		CLOSE(20) 
!-------------------------------------------------------------------------
	END