anchor.diff:          type(hecmwST_local_mesh),target :: hecMESH
           type(fstr_solid)                :: fstrSOLID
   
  -        integer :: i, j, ng, isect, ndof, id, nn
  +        integer :: i, j, ng, isect, ndof, id, nn            
   
           if( hecMESH%n_elem <=0 ) then
                stop "no element defined!"
  
positive.code:            enddo
             deallocate(fstrSOLID%output_ctrl)
          endif
  
          if( associated(fstrSOLID%GL) ) then
              deallocate(fstrSOLID%GL               ,STAT=ierror)
              if( ierror /= 0 ) then
  
negative.code:               call hecmw_abort( hecmw_comm_get_comm())
              end if
          endif
  
  end subroutine
  
  !> Initial setting of heat analysis
  