anchor.diff:         CALL BD_GravityForce(mmm,mEta,gravity,Fg)
          Fg(:) = Fg(:) + DistrLoad_GL(:,igp)
   
  -!WRITE(*,*) 'Fi'
  -!WRITE(*,*) Fi
  -!WRITE(*,*) 'Fb'
  -!WRITE(*,*) Fb
  -!WRITE(*,*) 'Fg'
  -!WRITE(*,*) Fg
          DO i=1,node_elem
              DO j=1,dof_node
                  temp_id1 = (i-1) * dof_node+j
  
positive.code:        rho(1:3,1:3) = EMass0_GL(4:6,4:6,igp)
         CALL BD_GaussPointDataMass(hhx,hpx,Nvvv,Naaa,RR0,node_elem,dof_node,vvv,aaa,vvp,mmm,mEta,rho)
         CALL BD_ElasticForce(E1,RR0,kapa,Stif,cet,Fc,Fd,Oe,Pe,Qe)
  !WRITE(*,*) 'Fc'
  !WRITE(*,*) Fc
  !WRITE(*,*) 'Fd'
  !WRITE(*,*) Fd
  !WRITE(*,*) 'aaa'
  !WRITE(*,*) aaa
         IF(damp_flag .EQ. 1) THEN
             CALL BD_DissipativeForce(beta,Stif,vvv,vvp,E1,Fc,Fd,Sd,Od,Pd,Qd,betaC,Gd,Xd,Yd)
         ENDIF
  
negative.code: 
     CALL BD_ReadInput(InitInp%InputFile,InputFileData,InitInp%RootName,ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        IF( ErrStat >= AbortErrLev ) RETURN
     CALL BD_ValidateInputData( InputFileData, ErrStat2, ErrMsg2 )
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
        IF( ErrStat >= AbortErrLev ) RETURN
  
     p%analysis_type  = InputFileData%analysis_type
     p%damp_flag  = InputFileData%InpBl%damp_flag
  