集中载荷 (cload.f)

此子程序用于用户自定义集中载荷，由*CLOAD卡片上的USER参数表征。头部和变量说明如下：
      subroutine cload(xload,kstep,kinc,time,node,idof,coords,vold,
     &  mi,ntrans,trab,inotr,veold)
!
!     用户子程序 cload
!
!
!     输入:
!
!     kstep              步编号
!     kinc               增量编号
!     time(1)            当前步时间
!     time(2)            当前总时间
!     node               节点编号
!     idof               自由度
!     coords(1..3)       节点的全局坐标
!     vold(0..mi(2)
!              ,1..nk)   所有节点的解场（对于模态
!                        动力学：所有请求了输出
!                        或施加了力的节点）
!                        0: 温度
!                        1: 全局x方向的位移
!                        2: 全局y方向的位移
!                        3: 全局z方向的位移
!                        4: 静压力
!     mi(1)              每个单元的最大积分点数（所有
!                        单元的最大值）
!     mi(2)              每个节点的最大自由度数（所有
!                        节点的最大值），用于像v(0:mi(2))...
!     veold(0..3,1..nk)  所有节点解场对时间的导数（对于模态
!                        动力学：所有请求了输出
!                        或施加了力的节点）
!                        0: 温度速率
!                        1: 全局x方向的速度
!                        2: 全局y方向的速度
!                        3: 全局z方向的速度
!     ntrans             转换定义的数目
!     trab(1..6,i)       定义转换i的两点坐标
!     trab(7,i)          -1: 圆柱形转换
!                         1: 矩形转换 
!     inotr(1,j)         应用于节点j的转换编号
!     inotr(2,j)         在节点j中应用了转换的SPC对应于MPC。
!                        inotr(2,j)包含为MPC的非齐次部分生成的
!                        新节点的编号
!
!     输出:
!
!     xload              节点"node"在idof方向上的集中载荷
!                        （全局坐标）
!