请教：初始化地下水位上下的密度

diwurou

我做了一个简单是算例，模型有2种材料（2个group），想达到这个目的：如果是group 1则密度为2700，如果是group 2则在水面上密度为2500，水面下为2600。
命令流如下：
gen zone brick size 4 4 4 p0 0 0 0 p1 20 0 0 p2 0 20 0 p3 0 0 20
gen zone brick size 4 4 4 p0 20 0 0 p1 40 0 0 p2 20 20 0 p3 20 0 20
group 1
gen zone brick size 4 4 4 p0 20 0 20 p1 40 0 20 p2 20 20 20 p3 30 0 40 &
                          p4 40 20 20 p5 30 20 40 p6 40 0 40 p7 40 20 40
group 2 range group 1 not
m m
pro bulk 2e8 shear 1e8 fric 45 coh 1e5 ten 1e5
set grav 0 0 -9.81
water dens 1000
;形成水位面
water table face 20,0,20 20,20,20 25,0,28
water table face 25,0,28 20,20,20 25,20,28
water table face 25,0,28 25,20,28 35,0,38
water table face 35,0,38 25,20,28 35,20,38
water table face 35,0,38 35,20,38 40,0,40
water table face 40,0,40 35,20,38 40,20,40
water table face 0,0,20  0,20,20  20,0,20
water table face 20,0,20 0,20,20  20,20,20

def ini_dens
   pnt=zone_head
   loop while pnt#null
if z_group(pnt)=1 then
  z_density(pnt)=2700
else
   if z_pp(pnt)#0 then
      z_density(pnt)=2600
   else
      z_density(pnt)=2500
   endif  
endif
  pnt=z_next(pnt)
endloop
end
ini_dens
模型请见附件，本人对FISH不是很懂，请问上面的函数哪里出了问题。

kunpeng0920

我也正是这个问题，求解答

肖邦不懂我

好像if z_group(pnt)=1 then得改成if z_group(pnt)=‘1’ then。  因为组名是字符串，我自己试了下，加上就没错误了，希望对你有所帮助

肖邦不懂我

另外，bad type conversion 错误主要指的就是变量所赋值与其类型不符