anchor.diff:  ! be interrogated by a routine using an instance of the type. If you add a field
   ! here, be sure to change the table of parameters used to size and index fieldmask above.
   
  -      REAL(ReKi), ALLOCATABLE :: Force(:,:)              ! Field: Force vectors (3,NNodes)
  -      REAL(ReKi), ALLOCATABLE :: Moment(:,:)             ! Field: Moment vectors (3,NNodes)
  -      REAL(R8Ki), ALLOCATABLE :: Orientation(:,:,:)      ! Field: Direction Cosine Matrix (DCM) (3,3,NNodes)
  -      REAL(R8Ki), ALLOCATABLE :: TranslationDisp(:,:)    ! Field: Translational displacements (3,NNodes)
  -      REAL(ReKi), ALLOCATABLE :: RotationVel(:,:)        ! Field: Rotational velocities (3,NNodes)
  -      REAL(ReKi), ALLOCATABLE :: TranslationVel(:,:)     ! Field: Translational velocities (3,NNodes)
  -      REAL(ReKi), ALLOCATABLE :: RotationAcc(:,:)        ! Field: Rotational accelerations (3,NNodes)
  -      REAL(ReKi), ALLOCATABLE :: TranslationAcc(:,:)     ! Field: Translational accelerations (3,NNodes)
  -      REAL(ReKi), ALLOCATABLE :: Scalars(:,:)            ! Scalars (nScalars,NNodes)
  -      INTEGER                 :: nScalars                ! Stores value of nScalars when created
  +      REAL(ReKi), ALLOCATABLE :: Force(:,:)              !< Field: Force vectors (3,NNodes)
  +      REAL(ReKi), ALLOCATABLE :: Moment(:,:)             !< Field: Moment vectors (3,NNodes)
  +      REAL(R8Ki), ALLOCATABLE :: Orientation(:,:,:)      !< Field: Direction Cosine Matrix (DCM) (3,3,NNodes)
  +      REAL(R8Ki), ALLOCATABLE :: TranslationDisp(:,:)    !< Field: Translational displacements (3,NNodes)
  +      REAL(ReKi), ALLOCATABLE :: RotationVel(:,:)        !< Field: Rotational velocities (3,NNodes)
  +      REAL(ReKi), ALLOCATABLE :: TranslationVel(:,:)     !< Field: Translational velocities (3,NNodes)
  +      REAL(ReKi), ALLOCATABLE :: RotationAcc(:,:)        !< Field: Rotational accelerations (3,NNodes)
  +      REAL(ReKi), ALLOCATABLE :: TranslationAcc(:,:)     !< Field: Translational accelerations (3,NNodes)
  +      REAL(ReKi), ALLOCATABLE :: Scalars(:,:)            !< Scalars (nScalars,NNodes)
  +      INTEGER                 :: nScalars                !< Stores value of nScalars when created
   !bjj: to be added later:       REAL(ReKi), ALLOCATABLE :: ElementScalars(nElementScalars,nelemlist)       ! Scalars associated with elements
   !bjj: to be added later:       INTEGER             :: nElementScalars               ! Stores value of nElementScalars when created
   
        ! Keeping track of siblings:
  -      TYPE(MeshType),POINTER :: SiblingMesh      => NULL() ! Pointer to mesh's (only) sibling
  +      TYPE(MeshType),POINTER :: SiblingMesh      => NULL() !< Pointer to mesh's (only) sibling
   
      END TYPE MeshType
   
   CONTAINS
  +!> This function returns the number of nodes in a given type of element.
      INTEGER FUNCTION NumNodes( Xelement )
  -     INTEGER, INTENT(IN) :: Xelement
  +     INTEGER, INTENT(IN) :: Xelement !< type of element
   
        SELECT CASE ( Xelement )
          CASE ( ELEMENT_POINT )
  
positive.code:      END SELECT
     END FUNCTION NumNodes
  
     LOGICAL FUNCTION HasMotionFields(Mesh)
  
        TYPE(MeshType), INTENT(IN) :: Mesh
  
        IF (     Mesh%FieldMask(MASKID_TRANSLATIONDISP)    &
            .or. Mesh%FieldMask(MASKID_ORIENTATION)        &
  
negative.code: 
     END FUNCTION HasLoadFields
  
  !> This subroutine copies the element record data from one ElemRecType data structure to another. It calls the Fortran 2003 
  !! intrinsic MOVE_ALLOC routine to move the address of the Src\%ElemNodes array to the Dest\%ElemNodes array without physically
  !! copying any of the array. On exist Src\%ElemNodes will be deallocated. 
     SUBROUTINE Mesh_MoveAlloc_ElemRecType( Src, Dest )
        TYPE(ElemRecType), INTENT(INOUT) :: Src   !< mesh containing ElemNodes to be moved
        TYPE(ElemRecType), INTENT(INOUT) :: Dest  !< mesh that will receive the ElemNodes array from Src 
     
        Dest%Xelement   = Src%Xelement
        Dest%Nneighbors = Src%Nneighbors
  