FEA way, Mustaine3, thank you very much for your quick answers.
Indeed, the problem was the surface definition, a different method had to be used. If you use the *.rpy command starightaway, you will get the following command:
side1Faces1 = s1.getSequenceFromMask(mask=('[#4 ]', ), )
Which is...