Igeometry=2 disabled in current version about spec HOT 4 CLOSED

I see what you mean. And I agree that adding an if check would fix the problem, provided that these variables (Dxyz, Nxyz) are not used in fixed boundary (I think indeed they are not because they describe the computational boundary). Can you first check that these are indeed never used in fixed-boundary? If so, then I would suggest that you add this IF statement in preset.h in a new temporary branch, so we can check that SPEC runs fine in fixed/free boundary test cases.

from spec.

Done as suggested. A simple grep shows Dxyz, Nxyz and Jxyz only appear in casing.h and bnorm.h, and therefore should be irrelevant to fixed-boundary.

from spec.

Great. I am going to test if the branch compiles and runs well on my end, and then you can merge the branch into master. Also, while doing that, can you (first) add a test input file for fixed-boundary, cylindrical geometry? That can be done in the folder InputFiles/TestCases/ and with the format G2VxxLyFi.001.sp where xx is the number of volumes and y is the Lconstraint number. Also one should add a file G2VxxLyFi.001.info with some description of the test case.

from spec.

Works for me. I wrote some review on the pull request. I think you can close this issue after merging.

from spec.