p['extractors'] = [
        {
            'name': 'tooth1',
            'nodes': [-10.20362, -17.46838, -229.9061],
            'variables': [ 'sigvm' ]
        },
        {
            'name': 'tooth2',
            'nodes': [-11.92466, 26.3042, -229.5354],
            'variables': [ 'sigxx', 'sigyy', 'sigzz', 'sigvm' ]
        },
        {
            'name': 'Lmuscle',
            'nodes': f'{path}/Lmuscle.stl',
            'variables': [ 'x0', 'y0', 'z0', 'sigvm' ]
        }
    ]

extracting results...
        internal energy = 95.219 N.mm
        contact_pts_Fx = 17.915 N
        [...]
        tooth1_sigvm = 4.507 MPa
        tooth2_sigxx = -1.000 MPa
        tooth2_sigyy = -1.025 MPa
        tooth2_sigzz = -0.687 MPa
        tooth2_sigvm = 3.803 MPa
        Lmuscle_x0 = (95 values in "Lmuscle_x0.tsv")
        Lmuscle_y0 = (95 values in "Lmuscle_y0.tsv")
        Lmuscle_z0 = (95 values in "Lmuscle_z0.tsv")
        Lmuscle_sigvm = (95 values in "Lmuscle_sigvm.tsv")
Total Execution time: 2.2s