def METHOD_NAME(self, cases, ID, nodes, elements, regions, loads):
    #print "regions**** = ",regions
    #nNodes = self.nnodes
    #nElements = self.nelements
    #result_names = ['Cp', 'Mach', 'U', 'V', 'W', 'E', 'rho',
                                  #'rhoU', 'rhoV', 'rhoW', 'rhoE']
    #nelements, three = elements.shape
    #print regions
    icase = 0
    itime = 0
    region_res = GuiResult(ID, header='Region', title='Region',
                           location='centroid', scalar=regions)
    cases[icase] = region_res(nid_res, (itime, 'Region'))
    icase += 1
    # centroidal
    Xc = zeros(len(elements), 'float64')
    Yc = zeros(len(elements), 'float64')
    Zc = zeros(len(elements), 'float64')
    area = zeros(len(elements), 'float64')
    Xn = zeros(len(nodes), 'float64')
    Yn = zeros(len(nodes), 'float64')
    Zn = zeros(len(nodes), 'float64')
    for i, element in enumerate(elements):
        p1, p2, p3, p4 = element
        P1 = array(nodes[p1])
        P2 = array(nodes[p2])
        P3 = array(nodes[p3])
        P4 = array(nodes[p4])
        a = P3 - P1
        b = P4 - P2
        A = 0.5 * norm(cross(a, b))
        x, y, z = (P1 + P2 + P3 + P4) / 4.0
        Xc[i] = x
        Yc[i] = y
        Zc[i] = z
        area[i] = A
    for i, node in enumerate(nodes):
        Xn[i] = node[0]
        Yn[i] = node[1]
        Zn[i] = node[2]
    cases[(ID, icase, 'node_x', 1, 'node', '%.2f', '')] = Xn
    cases[(ID, icase + 1, 'node_y', 1, 'node', '%.2f', '')] = Yn
    cases[(ID, icase + 2, 'node_z', 1, 'node', '%.2f', '')] = Zn
    cases[(ID, icase + 3, 'centroid_x', 1, 'centroid', '%.2f', '')] = Xc
    cases[(ID, icase + 4, 'centroid_y', 1, 'centroid', '%.2f', '')] = Yc
    cases[(ID, icase + 5, 'centroid_z', 1, 'centroid', '%.2f', '')] = Zc
    cases[(ID, icase + 6, 'Area', 1, 'centroid', '%.2f', '')] = area
    icase += 7
    #print("load.keys() = ", loads.keys())
    #for key in result_names:
        #if key in loads:
            #nodal_data = loads[key]
            #cases[(ID, key, 1, 'node', '%.3f')] = nodal_data
    return cases