anchor.diff:     
      
         ! Define system output initializations (set up and initialize input meshes) here:
  -   call Init_u(temp_L2, InitInp, p, u, ErrStat2, ErrMsg2)
  +   call Init_u(InitInp, p, u, ErrStat2, ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
               
         if (ErrStat >= AbortErrLev) then
  
positive.code:       call cleanup()
        return
     end if
     p%Shp(:,:) = 0.0D0
     p%Der(:,:) = 0.0D0
     p%Jacobian(:,:) = 0.0D0
  
     CALL BD_InitShpDerJaco(p%quadrature,p%GL,GLL,p%uuN0,&
             p%node_elem,p%elem_total,p%dof_node,p%ngp,    &
  
negative.code: ! This subroutine total element forces and mass matrices
  !-------------------------------------------------------------------------------
  
     REAL(BDKi),INTENT(IN   )    :: NGPpos(:,:)
     REAL(BDKi),INTENT(IN   )    :: EMass0_GL(:,:,:) ! Nodal material properties for each element
     REAL(BDKi),INTENT(  OUT)    :: elem_mass  ! Total element force (Fd, Fc, Fb)
     REAL(BDKi),INTENT(  OUT)    :: elem_CG(:)
     REAL(BDKi),INTENT(  OUT)    :: elem_IN(:,:)
     INTEGER(IntKi),INTENT(IN   ):: ngp ! Number of Gauss points
     REAL(BDKi),    INTENT(IN   ):: gw(:)
     REAL(BDKi),    INTENT(IN   ):: Jaco(:)
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
     REAL(BDKi)                  :: mmm
     INTEGER(IntKi)              :: igp
  !   INTEGER(IntKi)              :: ErrStat2                     ! Temporary Error status
  !   CHARACTER(ErrMsgLen)        :: ErrMsg2                      ! Temporary Error message
  