anchor.diff:     REAL(BDKi)                  :: rot0_temp(3)
      REAL(BDKi)                  :: rotu_temp(3)
      REAL(BDKi)                  :: rot_temp(3)
  +   REAL(BDKi)                  :: R0_temp(3,3)
      INTEGER(IntKi)              :: inode
      INTEGER(IntKi)              :: temp_id
      INTEGER(IntKi)              :: temp_id2
  
positive.code:    DO i=1,3
         uu0(i+3) = rot_temp(i)
     ENDDO
  
  END SUBROUTINE BD_GaussPointDataAt0
  !-----------------------------------------------------------------------------------------------------------------------------------
  
negative.code:        DO i=1,3
             uu0(i) = uu0(i) + hhi*Nuu0(temp_id+i)
             uu0(i+3) = uu0(i+3) + hhi*Nrr0(temp_id2+i)
  !           E10(i) = E10(i) + hpi*Nuu0(temp_id+i)
         ENDDO
     ENDDO
  
  