各位高手,下面是我的一段程序,计算地震液化。里面有一个额外变量zextra_1。另外三个变量yl_1, yl_2, c_ppr 计算结束时,我用pri fish命令显示fish程序内变量值,得到yl_2=1.5e6 yl_1=1.3e6 c_ppr=0.13.但是我在flac命令行中输入pl contour zextra_1 时,显示的云图数据竟然是100多~这让我很困惑,我程序想法是让z_extra 1的值等于c_ppr的值想(小于1),但是结果却很大。请高手帮我看看,有什么不对的地方,不胜感激!!!
config zextra 1
def get_old_stress
p_z = zone_head
loop while p_z # null
yl_1 = (z_sig1(p_z) + z_sig2(p_z) + z_sig3(p_z)) / 3.0 + z_pp(p_z)
p_z = z_next(p_z)
endloop
end
get_old_stress
def get_ppr
whilestepping
p_z = zone_head
loop while p_z # null
yl_2 = (z_sig1(p_z) + z_sig2(p_z) + z_sig3(p_z)) / 3.0 + z_pp(p_z)
c_ppr =1 - (yl_1)/(yl_2)
z_extra(p_z,1) = c_ppr
p_z = z_next(p_z)
endloop
end
get_ppr