So this code:
Code: Select all
j = 0
forall {
for i=1, n3d()-1 { //COSTAS: Need to start from ZERO; Not sure... i_membrane doesn't start at 0
j = j + 1 // global segment index
//Use indexes which chop up the segment into appropriate fractions; e.g. 1/8, 3/8, 5/8, 7/8
i_segment.x[j] = i_membrane(((i-1)*2+1)/(2*(n3d()-1)))*area(((i-1)*2+1)/(2*(n3d()-1)))
fprint("%f\n", i_segment.x[j]) //current segment saved
}
}
wopen()
Code: Select all
...
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.018351
0.001033
0.001033
0.001033
0.001033
0.000085
0.000085
0.000085
0.000085
0.000085
0.000005
0.000005
0.000005
0.000005
...
Hoc syntax is driving me bananas!
Thanks in advance for any help.
Cat