def METHOD_NAME(load_stack_superimpose):
    """
    Comparison of RMSF values computed with Biotite with results 
    obtained from GROMACS 2021.5.     
    """
    stack, supimp = load_stack_superimpose
    ca_mask = ((stack[0].atom_name == "CA") & (stack[0].element == "C"))
    rmsf = struc.rmsf(struc.average(supimp[:, ca_mask]), supimp[:, ca_mask])/10
    # Gromacs RMSF:
    # echo "C-alpha" | gmx rmsf -s 1l2y.gro -f 1l2y.xtc -o rmsf.xvg -res                    
    rmsf_gmx = np.array([
        0.1379, 0.036, 0.0261, 0.0255, 0.029, 0.0204, 0.0199, 
        0.0317, 0.0365, 0.0249, 0.0269, 0.032, 0.0356, 0.0446, 
        0.059, 0.037, 0.0331, 0.0392, 0.0403, 0.0954
    ])
    assert np.allclose(rmsf, rmsf_gmx, atol=1e-02