anchor.diff:  
      INTEGER(IntKi):: i,j
   
  -   eee = 0.0D0 
  +   eee = 0.0D0
      DO i=1,3
          eee(i) = E1(i) - RR0(i,1)
          eee(i+3) = kapa(i)
  
positive.code:                                   EStif0_GL,EMass0_GL,gravity,DistrLoad_GL,&
                                    ngp,node_elem,dof_node,damp_flag,beta,&
                                    elf,elm)
                             
     !-------------------------------------------------------------------------------
     ! This subroutine total element forces and mass matrices
     !-------------------------------------------------------------------------------
  
negative.code:            u%PointLoad%Force(1:3,i) = temp6(1:3)
             u%PointLoad%Moment(1:3,i) = temp6(4:6)
         ENDDO
  
         DO i=1,p%ngp * p%elem_total + 2
             temp6(:) = 0.0D0
             temp6(1:3) = u%DistrLoad%Force(1:3,i)
  