diff --git a/Source/LK1/L1D/GET_GRID_6X6_MASS.f90 b/Source/LK1/L1D/GET_GRID_6X6_MASS.f90
index e3b99824..8575549d 100644
--- a/Source/LK1/L1D/GET_GRID_6X6_MASS.f90
+++ b/Source/LK1/L1D/GET_GRID_6X6_MASS.f90
@@ -37,6 +37,7 @@ SUBROUTINE GET_GRID_6X6_MASS ( AGRID, IGRID, FOUND, GRID_MGG )
       USE CONSTANTS_1, ONLY           :  ZERO
       USE DOF_TABLES, ONLY            :  TDOF
       USE SPARSE_MATRICES, ONLY       :  I2_MGG, J_MGG, MGG
+      USE MODEL_STUF, ONLY            :  GRID_SEQ
  
       USE GET_GRID_6X6_MASS_USE_IFs
 
@@ -85,7 +86,7 @@ SUBROUTINE GET_GRID_6X6_MASS ( AGRID, IGRID, FOUND, GRID_MGG )
          ENDDO
       ENDDO
 
-      IGRID_DOF_NUM = 6*(IGRID - 1) + 1
+      IGRID_DOF_NUM = 6*(GRID_SEQ(IGRID) - 1) + 1
 k_do: DO K=1,NTERM_MGG
 
          IF ((I2_MGG(K) >= IGRID_DOF_NUM) .AND. (I2_MGG(K) <= IGRID_DOF_NUM + 5)) THEN
diff --git a/Source/LK1/L1E/MGGC_MASS_MATRIX.f90 b/Source/LK1/L1E/MGGC_MASS_MATRIX.f90
index 2c7f514a..a0eebcdc 100644
--- a/Source/LK1/L1E/MGGC_MASS_MATRIX.f90
+++ b/Source/LK1/L1E/MGGC_MASS_MATRIX.f90
@@ -84,8 +84,7 @@ SUBROUTINE MGGC_MASS_MATRIX
       IROW_START = 1
 i_do1:DO I=1,NGRID
 
-!xx      GRID_NUM = GRID_ID(INV_GRID_SEQ(I))               ! GRID_NUM's are in TDOFI order (internal DOF order)
-         GRID_NUM = GRID_ID(I)
+         GRID_NUM = GRID_ID(INV_GRID_SEQ(I))               ! GRID_NUM's are in TDOFI order (internal DOF order)
          CALL MGG_CONM2_PROC ( I, GRID_NUM, MGG_CONM2, MGG_CONM2_NONZERO )
          CALL GET_GRID_NUM_COMPS ( GRID_NUM, NUM_COMPS, SUBR_NAME )