边坡在降雨工况下的安全系数求解问题，代码写得零零散散，求大神帮忙改正呀！！！
model restore 'initial.sav'（前面初始平衡）
model configure fluid
model fluid active on
model mechanical active on
zone fluid cmodel assign isotropic（设置各项同性模型）
zone fluid property permeability 8.45e-8 range group '1'（设置渗透率）
zone fluid property porosity 0.5 range group '1'（设置孔隙率）
zone fluid property permeability 9.26e-8 range group '2'
zone fluid property porosity 0.5 range group '2'
zone fluid biot off（比奥关闭）
zone water density 1e3（设置水的密度 ）
zone gridpoint initialize pore-pressure 0 gradient 0 0 1e4（初始化孔隙压力）
zone gridpoint initialize saturation 0.27 range group '1'
zone gridpoint initialize saturation 0.27 range group '2'（设置饱和度）
zone gridpoint initialize fluid-modulus 2e9
zone gridpoint initialize fluid-tension -1e10（设置水的弹性模量和抗拉强度）
zone initialize density 2740 range group '1'
zone initialize density 2200 range group '2'（设置岩土干密度）
zone water plane normal （设置地下水面）这里是否需要设置？
zone gridpoint fix pore-pressure range group 'Bottom'（固定底面孔隙水压力）
zone fluid cmodel assign null range group 'Bottom'（底面设为不透水边界，四周不管）还是四周也要设置什么边界条件？ 这两句命令是不是重复了？
zone interface '1' create by-face range position-x 0.1 1161.8 position-y 0.1 699.9 position-z 120 390（顶面建立接触面）
apply pwell 0.01 range 'interface'（施加流量到顶面接触面，但命令错误）假设我100mm/d的降雨量下两天，是否需要设置时间步？
model solve（求解）
model save 'jiangyu.sav'
model fluid active off
model mechanical active off（关闭渗流）
zone face apply velocity-normal 0 range group 'South' or 'North'
zone face apply velocity 0 0 0 range group 'Bottom'
zone face apply velocity-normal 0 range group 'East' or 'West'（静力学边界）
model factor-of-safety ratio-local 1e-3 bracket 1.20 2.20 filename 'lu40slop.sav'（安全系数折减）
地下是我看的一个网上将降雨转为流量的代码 但是我不理解什么意思。
fish define liuliang
igp_add=interface.node.head(interface.head)
loop while igp_add # null
xx=interface.node.pos(igp_add,1)
yy=interface.node.pos(igp_add,2)
zz=interface.node.pos(igp_add,3)
p_gp=gp.near(xx,yy,zz)
nm_gp=gp.id(p_gp)
nodeare=interface.node.area(igp_add)（降雨转流量）
