anchor.diff:     DEALLOCATE(Pe_elem)
      DEALLOCATE(Qe_elem)
      DEALLOCATE(Se_elem)
  +   
  +   DEALLOCATE(Fi_elem)
  +   DEALLOCATE(Ki_elem)
  +   DEALLOCATE(Mi_elem)
  +   DEALLOCATE(Gi_elem)
      
      9999 IF(allo_stat/=0) THEN
               IF(ALLOCATED(Fc_elem)) DEALLOCATE(Fc_elem)
  
positive.code:    ALLOCATE(Mi_elem(dof_node,dof_node,node_elem),STAT = allo_stat)
     IF(allo_stat/=0) GOTO 9999
     Mi_elem = 0.0D0
     
     DO i=1,node_elem
         E10 = 0.0D0
         E1 = 0.0D0
  
negative.code:             IF(ALLOCATED(Pe_elem)) DEALLOCATE(Pe_elem)
              IF(ALLOCATED(Qe_elem)) DEALLOCATE(Qe_elem)
              IF(ALLOCATED(Se_elem)) DEALLOCATE(Se_elem)
              
              IF(ALLOCATED(Fi_elem)) DEALLOCATE(Fi_elem)
              IF(ALLOCATED(Ki_elem)) DEALLOCATE(Ki_elem)
              IF(ALLOCATED(Mi_elem)) DEALLOCATE(Mi_elem)
              IF(ALLOCATED(Gi_elem)) DEALLOCATE(Gi_elem)
          ENDIF
  
     END SUBROUTINE ElementMatrix
  