; --- model geometry (hollow cylinder - 9 degree wedge) ---
gen zone cshell p0 0. 0 1. p1 1.e2 0 1. &
p2 0. 0 0. p4 1.e2 0 0. &
p3 0.98769e2 0.15643e2 1. &
p5 0.98769e2 0.15643e2 0. &
dim 1 1 1 1 size n1 n2 1 ratio 1.1 1 1
range name in cyl end1 0 0 -100 end2 0 0 100 rad rin
range name out cyl end1 0 0 -100 end2 0 0 100 rad 99. not
range name xl y -.001 .001 z -.001 .001
; --- fish函数 ---
def sol_pp
tab1 = tab1 + 2
tab2 = tab2 + 2
pnt = gp_head
loop while pnt # null
rad = sqrt(gp_ypos(pnt)^2 + gp_zpos(pnt)^2)
if rad < 1.e-4 then
x = gp_xpos(pnt)
if x > 0.99 then
table(tab1,x) = gp_pp(pnt) / p0
ft = fltime
e_val = ss * x * x / ft
val = exp_int
table(tab2,x) = - val * c_e + 1.
end_if
end_if
pnt = gp_next(pnt)
end_loop
end
只是个fish 函数 我不理解的是“rad = sqrt(gp_ypos(pnt)^2 + gp_zpos(pnt)^2)”表示什么意思 而且rad 的值怎么会小于1.e-4 ?
还有gp_head 是第一个网格的位置吗??