求助lingo高手，为什么一直报错63模型定义不正确实在找不到问题

MODEL:
sets:
quyu/1..13/:Qs;
huishou/1..13/:Ri;
zhuanyun/1..13/:Rj;
chaijie/1..13/:Rk;
ershou/1..1/;
linjian/1..1/;
yuancailiao/1..1/;
feiliao/1..1/;
lk1(chaijie,ershou):Qkh;
lk2(chaijie,linjian):Qkl;
lk3(chaijie,yuancailiao):Qkm;
lk4(huishou,zhuanyun):Qij;
lk5(quyu,zhuanyun):Qsj;
lk6(zhuanyun,chaijie):Qjk;
lk7(huishou,zhuanyun):Dij;
lk8(zhuanyun,chaijie):Djk;
lk9(chaijie,ershou):Dkh;
lk10(chaijie,linjian):Dkl;
lk11(chaijie,yuancailiao):Dkm;
lk12(chaijie,feiliao):Dkn;
lk13(chaijie,feiliao):Qkn;
endsets
data:
Ri=1 1 1 1 1 1 1 1 1 1 1 1 1;
GS=2416439;SQ=150;phk=700;pkl=450;pkm=236.5;
cr=120;cd=20;cb1=10;cb2=30;
Fw=50000;Fi=30000;Fj=150000;Fk=900000;
Ow=2;Oi=5;Oj=10;Ok=70;
Cij=1;Cjk=1;Ckh=1;W=200;Ckl=0.0001;Ckm=0.0001;Ckn=0.0001;Uj=350000;Uk=500000;
Dij=
2    3.7    2.6    6.8    12    4.2    61    30    17.5    43    47    41    86    
3.7    2    2.5    4.7    14    5.5    70    34    14    46    43    52    82    
2.6    2.5    2    5.1    14    4    63    32    15    45    46    45    85    
6.8    4.7    5.1    2    16    8.4    68    36    11    48    40    46    84    
12    14    14    16    2    17    54    23    28    36    58    31    108    
4.2    5.5    4    8.4    17    2    59    33    18    39    33    47    96    
61    70    63    68    54    59    2    41    78    25    109    80    139    
30    34    32    36    23    33    41    2    47    29    75    47    120    
17.5    14    15    11    28    18    78    47    2    59    34    53    86    
43    46    45    48    36    39    25    29    59    2    90    63    116    
47    43    46    40    58    53    109    75    34    90    2    82    43    
41    52    45    46    31    47    80    47    53    63    82    2    132    
86    82    85    84    108    96    139    120    86    116    43    132    2    
;
Djk=
2    3.7    2.6    6.8    12    4.2    61    30    17.5    43    47    41    86    
3.7    2    2.5    4.7    14    5.5    70    34    14    46    43    52    82    
2.6    2.5    2    5.1    14    4    63    32    15    45    46    45    85    
6.8    4.7    5.1    2    16    8.4    68    36    11    48    40    46    84    
12    14    14    16    2    17    54    23    28    36    58    31    108    
4.2    5.5    4    8.4    17    2    59    33    18    39    33    47    96    
61    70    63    68    54    59    2    41    78    25    109    80    139    
30    34    32    36    23    33    41    2    47    29    75    47    120    
17.5    14    15    11    28    18    78    47    2    59    34    53    86    
43    46    45    48    36    39    25    29    59    2    90    63    116    
47    43    46    40    58    53    109    75    34    90    2    82    43    
41    52    45    46    31    47    80    47    53    63    82    2    132    
86    82    85    84    108    96    139    120    86    116    43    132    2    
;
Dkh=
6.8    
4.7    
5.1    
2    
16    
8.4    
68    
36    
11    
48    
40    
46    
84    
;
Dkl=
4.2    
5.5    
4    
8.4    
17    
2    
59    
33    
18    
39    
53    
47    
96    
;
Dkm=
4.2    
5.5    
4    
8.4    
17    
2    
59    
33    
18    
39    
53    
47    
96    
;
Dkn=
12    
14    
14    
16    
2    
17    
54    
23    
28    
36    
58    
31    
108    
;
Qs=
137819    174340    210168    352033    180187    279197    67135    178788    314712    88314    141463    149140    143143    
;
enddata
max=U-C1-C2-C3-C41-C42-C43-C44-C45-C46;
U=GS*SQ+@sum(chaijie(k):@sum(ershou(h):Qkh(k,h)*pkh))+@sum(chaijie(k):@sum(linjian(l):Qkl(k,l)*pkl))+@sum(chaijie(k):@sum(yuancailiao(m):Qkm(k,m)*pkm));
C1=(cr+cd)*rw*SQ+(cr-cb1)*ep*SQ+(cr-cb2)*u*SQ-Fw+@sum(huishou(i):Ri(i)*Fi)+@sum(zhuanyun(j):Rj(j)*Fj)+@sum(chaijie(k):Rk(k)*Fk);
C2=Fw+@sum(huishou(i):Ri(i)*Fi)+@sum(zhuanyun(j):Rj(j)*Fj)+@sum(chaijie(k):Rk(k)*Fk);
C3=Ow*SQ+(ep+rw)*Oi*SQ+Oj*(@sum(quyu(s):@sum(zhuanyun(j):Qsj(s,j)*Rj(j)))+@sum(huishou(i):@sum(zhuanyun(j):Qij(i,j)*Rj(j))))
   +Ok*(@sum(zhuanyun(j):@sum(chaijie(k):Qjk(j,k)*Rj(j)*Rk(k))));
C41=Cij*(@sum(quyu(i):@sum(zhuanyun(j):Qij(i,j))*Dij(i,j)*Rj(j)));
C42=Cjk*(@sum(zhuanyun(v):@sum(chaijie(k):Qjk(v,k)*Djk(v,k)*Rj(v)*Rk(k))));
C43=Ckh*(@sum(ershou(h):@sum(chaijie(k):Qkh(k,h)*Dkh(k,h)*Rk(k))));
C44=W*Ckl*(@sum(linjian(l):@sum(chaijie(k):Qkl(k,l)*Dkl(k,l)*Rk(k))));
C45=W*Ckm*(@sum(yuancailiao(m):@sum(chaijie(k):Qkm(k,m)*Dkm(k,m)*Rk(k))));
C46=W*Ckn*(@sum(feiliao(n):@sum(chaijie(k):Qkn(k,n)*Dkn(k,n)*Rk(k))));
@for(huishou(i):@sum(zhuanyun(z):Qij(i,z))=(ep+rw)*Qs(i));
@for(quyu(s):@sum(zhuanyun(z):Qsj(s,z))=u*Qs(s));
ep+rw+u=1;
@for(zhuanyun(z):@sum(huishou(i):Qij(i,z))+@sum(quyu(s):Qsj(s,z))=@sum(chaijie(k):Qjk(z,k)));
@for(chaijie(k):@sum(ershou(h):Qkh(k,h))=0.45*@sum(zhuanyun(z):Qjk(z,k)));
@for(chaijie(k):@sum(linjian(l):Qkl(k,l))=0.25*@sum(zhuanyun(z):Qjk(z,k)));
@for(chaijie(k):@sum(yuancailiao(m):Qkm(k,m))=0.2*@sum(zhuanyun(z):Qjk(z,k)));
@for(chaijie(k):@sum(feiliao(n):Qkn(k,n))=0.1*@sum(zhuanyun(z):Qjk(z,k)));
@for(zhuanyun(z):@sum(huishou(i):Qij(i,z))+@sum(quyu(s):Qsj(s,z))<=Uj*Rj(z));
@for(chaijie(k):@sum(zhuanyun(z):Qjk(z,k))<=Uk*Rk(k));
@sum(zhuanyun(z):Rj(z))>=1;@sum(zhuanyun(z):Rj(z))<=13;
@sum(chaijie(k):Rk(k))>=1;@sum(chaijie(k):Rk(k))<=13;
@for(zhuanyun(z):@bin(Rj(z)));
@for(chaijie(k):@bin(Rk(k)));