ub1.f90

这是一个用fortran编写的“第一类边界条件的有限元算法”

!***************************************************************************!
!                SUBROUTINE UB1(nd1,nb1,u1,nd,iw,sk,u)                      !
!    功能    在定带宽储存的总体系数矩阵和右侧列向量上加上第一类边界条件。   !
!    参数说明                                                               !
!      nd1---整型变量,输入参数,第一类边界条件的节点数。                     !
!      nb1---nd1个元的一维整数组，输入参数，存放第一类边界条件的节点号。    !
!      u1----nd1个元的一维实数组，输入参数，存放第一类边界节点的场值。      !
!      nd----整型变量,输入参数,节点总数。                                   !
!      iw----整型变量,输入参数,半带宽。                                     !
!      sk----nd*iw的二维整数组,输入、输出参数,定带宽存放总体系数矩阵，输出  !
!                时,第一类边界条件已带入。                                  !
!      u-----nd个元的一维实数组，输出参数，存放加入第一类边界条件后的右端列 !
!                向量。                                                     !
!***************************************************************************!
 subroutine UB1(nd1,nb1,u1,nd,iw,sk,u)

   dimension nb1(nd1),u1(nd1),sk(nd,iw),u(nd)

     do 10 i=1,nd
10	   u(i)=0
!     end do
	 do 20 i=1,nd1
	   j=nb1(i)
	   sk(j,iw)=sk(j,iw)*1.E10
20     u(j)=sk(j,iw)*u1(i)
  
   return
 end subroutine UB1