Hi Alessandro,
The problem is the order of the operations.
You do a rotation of surfaces that define volume 1.
That’s why they are kept.
The rotation is properly handled if you write it as follows:
tmp[] = Physical Surface {102,103};
Rotate {{1, 0, 0}, {-0.5, 0.65, -0.25}, Pi/6} {
Surface {tmp[]};
}
Surface Loop(5) = {28, 23, 25, 27, 26, 24};
Surface Loop(6) = {9, 4, 6, 8, 7, 10, 1, 22, 18, 21, 20, 17, 19};
Surface Loop(7) = {15, 11, 13, 12, 16, 14};
Volume(1) = {5, 6, 7};
Physical Volume(999) = {1};
Best regards,
Ruth
—
Prof. Ruth V. Sabariego
KU Leuven
Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
On 20 Oct 2017, at 10:46, Alessandro Vicini
<[email protected]<mailto:[email protected]>> wrote:
Hullo everybody. In the last two lines of the attached geo file I retrieve the
list of elementary surfaces ID from their physical IDs, and then apply a
rotation to the list. However, only the edges are rotated, and the "old"
surfaces and edges are not deleted after the operation. Why...? Thank you.
Alessandro
<Test.geo>_______________________________________________
gmsh mailing list
[email protected]<mailto:[email protected]>
http://onelab.info/mailman/listinfo/gmsh
_______________________________________________
gmsh mailing list
[email protected]
http://onelab.info/mailman/listinfo/gmsh