2016-12-08 10:25 PM
goto "MasterEnd"in the master script:
!area calculation for listing ! area calculation polyArea = 0 for q = 1 to nc - 1 idx = 2 * q - 1 x1 = field_poly[1][idx] y1 = field_poly[1][idx + 1] x2 = field_poly[1][idx + 2] y2 = field_poly[1][idx + 3] polyArea = polyArea + (x2 + x1) * (y2 - y1) next q polyArea = abs (polyArea) / 2 holeArea = 0 for wo = 1 to nHoles for q = 1 to n_hpoly[wo] - 1 idx = 2 * q - 1 x1 = hole_poly[wo][idx] y1 = hole_poly[wo][idx + 1] x2 = hole_poly[wo][idx + 2] y2 = hole_poly[wo][idx + 3] holeArea = holeArea + (x2 + x1) * (y2 - y1) next q next wo holeArea = abs (holeArea) / 2 polyArea = polyArea - holeArea Ceiling_area=polyArea Parameters Ceiling_area=Ceiling_areaThe calculation is copied from the 2d script - and does not interfere with the calculation done there.
line2 sectpand add a line below:[1], sectp [2], sectp[w + 1][1], sectp[w + 1][2]
hotline2 sectpIt is basically the same code as above but adding a hotline in the same place the gridline is. I did not add the IDs for the hotlines so anchoring dimensions to them may be risky (no time to analyze the whole code).[1], sectp [2], sectp[w + 1][1], sectp[w + 1][2]