渗透系数设置成变量在FLAC3D中到底能不能运行？

飘落

渗透系数设置成随时间变化的变量在FLAC3D中到底能不能运行？这是我用FISH写的函数，怎么调用的时候老是不对？大家看看那里错了？
def constants
    k_wl=1.8e-7
    k_w2=1.2e-7
    a=2.28
    b=0.0406
    d=1.518e-3
    f=1.6e3
    h=1.0e3
    x=60
end
constants
def shentouxishu
       loop t (0,60)
            n=-(b*t)
            q1=(k_w1*f*d)/(h*a)
            q2=(k_w2*f*d)/(h*a)
            s1=q1*exp(n)
            s2=q2*exp(n)
            command
            prop perm s1 poros 0.48 range group lishatop_1 any group lishatop_2 any group lishatop_3 any &                                               group lishabase_1 any group lishabase_2 any group lishabase_3 any
            prop perm s2 poros 0.51 range group ruanshitu_1 any group ruanshitu_2 any group ruanshitu_3 any
            end_command
        end_loop
end
shentouxishu

hxm19920629

mark一下