loop i="i" (1,izones)什么意思？if state(i,j) > 0中的state(i,j)可以直接赋值吗？ 能帮我解释一下吗？ [复制链接]

def cfi
  loop i (1,izones)
    loop j (1,jzones)
       if state(i,j) > 0 then
           h_mm=hb_mmr
           h_ss=hb_ssr
       else
           h_mm=hb_mmi
           h_ss=hb_ssi
       end_if
       effsxx = sxx(i,j)  + pp(i,j)
       effsyy = syy(i,j)  + pp(i,j)
       effszz = szz(i,j)  + pp(i,j)
       tension(i,j)=0.5*hb_sc*(sqrt(h_mm^2+4*h_ss)-h_mm)
       temp1=-0.5*(effsxx+effsyy)
       temp2=sqrt(sxy(i,j)^2+0.25*(effsxx-effsyy)^2)
       s3=min(temp1-temp2,-effszz)
       if s3<0.0 then
         s3=0.0
       end_if
       h_lam=1+0.5*h_mm*hb_sc/sqrt(h_mm*hb_sc*s3+hb_sc^2*h_ss)
       fi=2*atan(sqrt(h_lam))/degrad-90
       friction(i,j)=fi
       sci=s3*(1-h_lam)+sqrt(h_mm*s3*hb_sc+hb_sc^2*h_ss)
       coe=0.5*sci/sqrt(h_lam)
       cohesion(i,j)=coe
       psi = hoek_psi
       dilation(i,j)=psi
     end_loop
  end_loop
end

有点不太明白。loop i="i" (1,izones)什么意思？if state(i,j) > 0中的state(i,j)可以直接赋值吗？ 班主能帮我解释一下吗？

  loop i (1,izones)，是i从1到izones遍历，调用的时候会定义izones的值，一般应该都是所有单元的个数。state(i,j)可以赋值，可以在if前面state(i,j)=1。

呵呵，加油!