dvm_combination.f90

圆柱绕流问题的涡方法数值模拟（fortran源码）

subroutine combination(k0,A,B)
use varible
real dis
real ::A(nlimit,knp),B(nlimit,knp)
integer :: k0,k,L,k1,L1
 do L=1,knp
    do k=1,k0
	if(ndp(k,l)==0)then
	 continue
	else
	  do L1=1,knp
	     do k1=1,k0
		    if (ndp(k1,L1)==0)then
		      continue 
		    else
		      if (L==L1.and.k==k1)then
		        continue 
		      else 
		        dis=sqrt((A(k,L)-A(k1,L1))**2+(B(k,L)-B(k1,L1))**2)
		       !taosum=tao(k,L)+tao(k1,L1)
		        if(dis<sigma(k1))then 
		          !A(k1,L1)=(tao(k1,L1)*A(k1,L1)+tao(k,L)*A(k,L))/taosum
		          !B(k1,L1)=(tao(k1,L1)*B(k1,L1)+tao(k,L)*B(k,L))/taosum
		          tao(k1,L1)=tao(k1,L1)+tao(k,L)
		          A(k,L)=0.0
		          B(k,L)=0.0
		          tao(k,L)=0.0
		          ndp(k,L)=0.0
			      continue
			    else 
			      continue
			    end if
		      end if
		    end if
	     end do
	  end do
	end if
    end do
 end do
 return
end subroutine