求帮忙啊，一个程序的小问题啊！

相思枕畔

MODEL:
    title economic pipe;
    sets:
         tubulation/1..9/,Q,A,EK;
         caliber/1,2,3,4,5/;
         LINK(tubulation,caliber):C,X,D;
    ENDSETS
    DATA:
    !各段管长;
    L=700,700,700,500,500,500,500,500,500;
    !各段流量;
    Q=0.135,0.060,0.015,0.015,0.045,0.015,0.015,0.015,0.015;
    !节点服务水压标高;
    EK=47,46,44,46,44,43,45,46,43;
    !水源到节点的管段数;
    A=1,2,3,2,2,3,3,3,3;
    !标准管径单价;
    C=42.0 55.2 64.0 83.2 103.2   
      38.7 42.0 55.2 64.0 0.00
      25.0 30.8 0.00 0.00 0.00
      25.0 30.8 0.00 0.00 0.00
      30.8 38.7 42.0 55.2 0.00
      25.0 30.8 0.00 0.00 0.00
      25.0 30.8 0.00 0.00 0.00
      25.0 30.8 0.00 0.00 0.00
      25.0 30.8 0.00 0.00 0.00;
    !标准管径;
    D=0.25 0.30 0.35 0.40 0.50  
      0.20 0.25 0.30 0.35 0.00
      0.10 0.15 0.00 0.00 0.00
      0.10 0.15 0.00 0.00 0.00
      0.15 0.20 0.25 0.30 0.00
      0.10 0.15 0.00 0.00 0.00
      0.10 0.15 0.00 0.00 0.00
      0.10 0.15 0.00 0.00 0.00
      0.10 0.15 0.00 0.00 0.00;
    ENDDATA
    !OBJECTIVE FUNCTION(目标函数);
    [OBJ]MIN=@SUM(LINK(I,J):C(I,J)*X(I,J));
    !pipe range CONSTRAINTE(管长约束);
    @FOR(tubulation(I):[RANGE_CON]@SUM(caliber(J):X(I,J))=L(I););
    !Stress constraint(压力约束);
    @FOR(caliber(J)@SUM(tubulation(I)0.135^2*X(1,J)/D(1,J)^5.342))*(1.942*10^-3))<(60-47));
    @FOR(CALIBER(J)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.060^2*X(2,J)/D(2,J)^5.342))*(1.942*10^-3))<(60-46));
    @for(caliber(j)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.060^2*X(2,J)/D(2,J)^5.342+0.015^2*X(3,J)/D(3,J)^5.342))*(1.942*10^-3))<(60-44));
    @for(caliber(j)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.015^2*X(4,J)/D(4,J)^5.342))*(1.942*10^-3))<(60-46));
    @FOR(CALIBER(J)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.045^2*X(5,J)/D(5,J)^5.342))*(1.942*10^-3))<(60-44));
    @for(caliber(j)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.045^2*X(5,J)/D(5,J)^5.342+0.015^2*X(6,J)/D(6,J)^5.342))*(1.942*10^-3))<(60-43));
    @for(caliber(j)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.045^2*X(5,J)/D(5,J)^5.342+0.015^2*X(7,J)/D(7,J)^5.342))*(1.942*10^-3))<(60-45));
    @for(caliber(j)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.060^2*X(2,J)/D(2,J)^5.342+0.015^2*X(8,J)/D(8,J)^5.342))*(1.942*10^-3))<(60-46));
    @for(caliber(j)@SUM(TUBULATION(I)0.135^2*X(1,J)/D(1,J)^5.342+0.060^2*X(2,J)/D(2,J)^5.342+0.015^2*X(9,J)/D(9,J)^5.342))*(1.942*10^-3))<(60-43));
    !速度要求;
    !非负假定;
    @FOR(LINK(I,J):X(I,J)>0;);
    END

这个程序能运行，但是有些问题，能不能把数组d中的0剔除，下面d是做分母的，或者其他的方法，也可以，求大神指导！！