Error $141 - Can't figure out what is wrong

GAMS
1

Hello,

I need your help with the following error:
Error 141 - Symbol declared, but no values have been assigned.

Would you be able to help me fix or point out what I can do to fix my code below?
I greatly appreciate you assistance with this issue.


Sets bus /124/,slack(bus) /13/,Gen /g1g12/,t /t1t24/, genD(gen) /g1g2,g5,g7g11/,
genN(gen) /g3,g4,g6,g12/;
scalars Sbase /100/, VOLL /10000/, VOLW /50/;
alias(bus,node);
Sets gn NODES /Anderlues, Antwerpen, Arlon, Berneau, Blaregnies, Brugge, Dudzele,
Gent, Liege, Loenhout, Mons, Namur, Petange, Peronnes, Sinsin,
Voeren, Wanze, Warnand, Zeebrugge, Zomergem/
a PIPES /L1L24/;
Alias (gn,gm);
set Pnm(a,gn,gm) ‘arc description’
/
L1 .Zeebrugge .Dudzele
L2 .Zeebrugge .Dudzele
L3 .Dudzele .Brugge
L4 .Dudzele .Brugge
L5 .Brugge .Zomergem
L6 .Loenhout .Antwerpen
L7 .Antwerpen .Gent
L8 .Gent .Zomergem
L9 .Zomergem .Peronnes
L10 .Voeren .Berneau
L11 .Voeren .Berneau
L12 .Berneau .Liege
L13 .Berneau .Liege
L14 .Liege .Warnand
L15 .Liege .Warnand
L16 .Warnand .Namur
L17 .Namur .Anderlues
L18 .Anderlues .Peronnes
L19 .Peronnes .Mons
L20 .Mons .Blaregnies
L21 .Warnand .Wanze
L22 .Wanze .Sinsin
L23 .Sinsin .Arlon
L24 .Arlon .Petange
/;

table Ndata(gn,*) ‘Node Data’
slo sup Sd plo pup c
Anderlues 0 1.2 0 0 66.2 0
Antwerpen 0 0 4.03 1.25 80 0
Arlon 0 0 0.22 0 66.2 0
Berneau 0 0 0 0 66.2 0
Blaregnies 0 0 15.62 2.08 66.2 0
Brugge 0 0 3.92 1.25 80 0
Dudzele 0 8.4 0 0 77 2.28
Gent 0 0 5.26 1.25 80 0
Liege 0 0 6.39 1.25 66.2 0
Loenhout 0 4.8 0 0 77 2.28
Mons 0 0 6.85 0 66.2 0
Namur 0 0 2.12 0 66.2 0
Petange 0 0 1.92 1.04 66.2 0
Peronnes 0 0.96 1.92 1.04 66.2 1.68
Sinsin 0 0 0 0 63 0
Voeren 0 22.01 0 2.08 66.2 1.68
Wanze 0 0 0 0 66.2 0
Warnand 0 0 0 0 66.2 0
Zeebrugge 0 11.59 0 0 77 2.28
Zomergem 0 0 0 0 80 0;

set GElink(gn,gen)
/Loenhout . g12
Voeren . g6
Sinsin . g3
Petange . g4/;

table AData(a,gn,gm,*) ‘Arc Data’
act C2mn
L1. Zeebrugge . Dudzele 0 9.07027
L2. Zeebrugge . Dudzele 0 9.07027
L3. Dudzele . Brugge 0 6.04685
L4. Dudzele . Brugge 0 6.04685
L5. Brugge . Zomergem 0 1.39543
L6. Loenhout . Antwerpen 0 0.10025
L7. Antwerpen . Gent 0 0.14865
L8. Gent . Zomergem 0 0.22689
L9. Zomergem . Peronnes 0 0.65965
L10. Voeren . Berneau 1 7.25622
L11. Voeren . Berneau 1 0.10803
L12. Berneau . Liege 0 1.81405
L13. Berneau . Liege 0 0.02701
L14. Liege . Warnand 0 1.45124
L15. Liege . Warnand 0 0.02161
L16. Warnand . Namur 0 0.86384
L17. Namur . Anderlues 0 0.90703
L18. Anderlues . Peronnes 0 7.25622
L19. Peronnes . Mons 0 3.62811
L20. Mons . Blaregnies 0 1.45124
L21. Warnand . Wanze 0 0.05144
L22. Wanze . Sinsin 1 0.00642
L23. Sinsin . Arlon 0 0.00170
L24. Arlon . Petange 0 0.02782;

table GD(Gen,*) ‘Generating units characteristics’
Pmax Pmin b CostsD costst RU RD SU SD UT DT uini U0 So
g1 400 100 5.47 0 0 47 47 105 108 1 1 1 5 0
g2 400 100 5.47 0 0 47 47 106 112 1 1 1 6 0
g3 152 30.4 13.32 1430.4 1430.4 14 14 43 45 8 4 1 2 0
g4 152 30.4 13.32 1430.4 1430.4 14 14 44 57 8 4 1 2 0
g5 155 54.25 16 0 0 21 21 65 77 8 8 0 0 2
g6 155 54.25 10.52 312 312 21 21 66 73 8 8 1 10 0
g7 310 108.5 10.52 624 624 21 21 112 125 8 8 1 10 0
g8 350 140 10.89 2298 2298 28 28 154 162 8 8 1 5 0
g9 350 75 20.7 1725 1725 49 49 77 80 8 8 0 0 2
g10 591 206.85 20.93 3056.7 3056.7 21 21 213 228 12 10 0 0 8
g11 60 12 26.11 437 437 7 7 19 31 4 2 0 0 1
g12 300 0 0 0 0 35 35 315 326 0 0 1 2 0;


set GB(bus, Gen) ‘connectivity index of each generating unit to each bus’
/
18.g1
21.g2

  1. g3
  2. g4
    15.g5
    16.g6
    23.g7
    23.g8
  3. g9
    13.g10
    15.g11
    22.g12
    /;
    Table BusData(bus,*) ‘Demands of each bus in MW’
    Pd Qd
    1 108 22
    2 97 20
    3 180 37
    4 74 15
    5 71 14
    6 136 28
    7 125 25
    8 171 35
    9 175 36
    10 195 40
    13 265 54
    14 194 39
    15 317 64
    16 100 20
    18 333 68
    19 181 37
    20 128 26;

Table branch(bus, node,*) ‘Network technical’
r x b limit
1 .2 0.0026 0.0139 0.4611 175
1 .3 0.0546 0.2112 0.0572 175
1 .5 0.0218 0.0845 0.0229 175
2 .4 0.0328 0.1267 0.0343 175
2 .6 0.0497 0.192 0.052 175
3 .9 0.0308 0.119 0.0322 175
3 .24 0.0023 0.0839 0 400
4 .9 0.0268 0.1037 0.0281 175
5 .10 0.0228 0.0883 0.0239 175
6 .10 0.0139 0.0605 2.459 175
7 .8 0.0159 0.0614 0.0166 175
8 .9 0.0427 0.1651 0.0447 175
8 .10 0.0427 0.1651 0.0447 175
9 .11 0.0023 0.0839 0 400
9 .12 0.0023 0.0839 0 400
10.11 0.0023 0.0839 0 400
10.12 0.0023 0.0839 0 400
11.13 0.0061 0.0476 0.0999 500
11.14 0.0054 0.0418 0.0879 500
12.13 0.0061 0.0476 0.0999 500
12.23 0.0124 0.0966 0.203 500
13.23 0.0111 0.0865 0.1818 500
14.16 0.005 0.0389 0.0818 500
15.16 0.0022 0.0173 0.0364 500
15.21 0.00315 0.0245 0.206 1000
15.24 0.0067 0.0519 0.1091 500
16.17 0.0033 0.0259 0.0545 500
16.19 0.003 0.0231 0.0485 500
17.18 0.0018 0.0144 0.0303 500
17.22 0.0135 0.1053 0.2212 500
18.21 0.00165 0.01295 0.109 1000
19.20 0.00255 0.0198 0.1666 1000
20.23 0.0014 0.0108 0.091 1000
21.22 0.0087 0.0678 0.1424 500 ;

Table DataWDL(t,*)
w d g
t1 0.05 0.68 0.86
t2 0.08 0.65 0.88
t3 0.10 0.61 0.87
t4 0.28 0.60 0.82
t5 0.37 0.58 0.79
t6 0.58 0.60 0.77
t7 0.67 0.63 0.74
t8 0.58 0.65 0.71
t9 0.49 0.70 0.70
t10 0.57 0.78 0.69
t11 0.77 0.74 0.67
t12 0.71 0.86 0.67
t13 0.84 0.88 0.68
t14 0.93 0.84 0.69
t15 0.95 0.82 0.70
t16 0.97 0.82 0.70
t17 0.85 0.87 0.73
t18 0.61 1.00 0.76
t19 0.60 0.98 0.80
t20 0.57 0.93 0.84
t21 0.57 0.89 0.90
t22 0.57 0.81 0.93
t23 0.71 0.74 0.97
t24 0.78 0.73 1.00
;
Parameters Wcap(bus),conex(bus,node),SD(gn);
branch(bus, node, ‘bij’)branch(bus,node, 'Limit')=1/branch(bus,node, 'x'); conex(bus, node)(branch(bus, node, ‘limit’)and branch(node, bus, ‘limit’))=1;
conex(bus, node)$(conex(node,bus))=1;
Variables f(a,gn,gm,t),sg(gn,t),pressure(gn,t),
EC,Pij(bus,node,t),Pg(Gen,t),delta(bus,t),lsh(bus,t),Pw(bus,t),pc(bus,t),Gc,OF;
Pnm(a,gn,gm)$adata(a,gn,gm,‘c2mn’)=yes;


Equations const1,const2,const3,const4,const5,const6,CG1,CG2,CG3,CG4,Objective;

const1(bus,node,t)$conex(bus,node)…Pij(bus,node,t)=e=branch(bus,node,‘bij’)*
(delta(bus,t)-delta(node,t));

const2(bus,t)…lsh(bus,t)$BusData(bus,‘pd’)+Pw(bus,t)$Wcap(bus)+sum(Gen$GB(bus,Gen),
Pg(Gen,t))-DataWDL(t,‘d’)*BusData(bus,‘pd’)/Sbase=e=+sum(node$conex(node,bus),
Pij(bus,node,t));

const3…EC=e=sum((bus,GenD,t)$GB(bus,GenD),Pg(GenD,t)*GD(GenD,‘b’)Sbase)
+sum((bus,t),VOLLlsh(bus,t)Sbase$BusData(bus,‘pd’)+VOLWPc(bus,t)*sbase$Wcap(bus));

const4(gen,t)…pg(gen,t+1)-pg(gen,t)=l=GD(gen,‘RU’)/Sbase;

const5(gen,t)…pg(gen,t-1)-pg(gen,t)=l=GD(gen,‘RD’)/Sbase;

const6(bus,t)$Wcap(bus)…pc(bus,t)=e=DataWDL(t,‘w’)*Wcap(bus)/Sbase-pw(bus,t);

Pg.lo(Gen,t)=GD(Gen,‘Pmin’)/Sbase;
Pg.up(Gen,t)=GD(Gen,‘Pmax’)/Sbase;
delta.up(bus,t)=pi/2;
delta.lo(bus,t)=-pi/2;
delta.fx(slack,t)=0;
Pij.up(bus,node,t)((conex(bus,node)))=1*branch(bus,node,'Limit')/Sbase; Pij.lo(bus,node,t)((conex(bus,node)))=-1*branch(bus,node,‘Limit’)/Sbase;
lsh.up(bus,t)=DataWDL(t,‘d’)*BusData(bus,‘pd’)/Sbase;
lsh.lo(bus,t)=0;

Pw.up(bus,t)=DataWDL(t,‘w’)*Wcap(bus)/Sbase;
Pw.lo(bus,t)=0;
Pc.up(bus,t)=DataWDL(t,‘w’)*Wcap(bus)/Sbase;
Pc.lo(bus,t)=0;
SD(gn)=Ndata(gn,‘SD’);
CG1(gn,t)…sum(Pnm(a,gn,gm),f(Pnm,t))=e=sum(Pnm(a,gm,gn),f(Pnm,t))
+sg(gn,t)$(Ndata(gn,‘Sup’)>0)-DataWDL(t,‘G’)*SD(gn)
-sum((GenN)$Gelink(gn,GenN),Pg(GenN,t)*GD(GenN,‘b’)*Sbase/35315);

CG2(Pnm(a,gn,gm),t)$(AData(a,gn,gm,‘C2mn’) AND AData(a,gn,gm,‘ACT’)=0)
…signpower(f(Pnm,t),2)=e= AData(Pnm,‘C2mn’)*(pressure(gn,t)-pressure(gm,t));

CG3(Pnm(a,gn,gm),t)$(AData(a,gn,gm,‘C2mn’) AND AData(a,gn,gm,‘ACT’)=1)
…-sqr(f(Pnm,t))=g= AData(Pnm,‘C2mn’)*(pressure(gn,t)-pressure(gm,t));

CG4…Gc=e=sum((gn,t),35315*ndata(gn,‘c’)*sg(gn,t)$Ndata(gn,‘Sup’));

Objective…OF=e=EC+Gc;
sg.lo(gn,t)=0;
sg.up(gn,t)=ndata(gn,‘sup’);

pressure.lo(gn,t)=sqr(ndata(gn,‘plo’));
pressure.up(gn,t)=sqr(ndata(gn,‘pup’));
f.lo(Pnm(a,gn,gm),t)(AData(a,gn,gm,'C2mn'))= -sqrt(AData(a,gn,gm,'C2mn')*(pressure.up(gn,t)-pressure.lo(gn,t))); f.up(Pnm(a,gn,gm),t)(AData(a,gn,gm,‘C2mn’))=
sqrt(AData(a,gn,gm,‘C2mn’)*(pressure.up(gn,t)-pressure.lo(gn,t)));
f.lo (Pnm(a,gn,gm),t)$(AData(a,gn,gm,‘C2mn’)AND AData(a,gn,gm,‘ACT’)=1)=0;
Model overall /all/;
Solve overall using nlp min OF;




Thank you.

2

Hi, please follow thew forum rules. Don’t just paste program in your text. At least use a code section or better attach the file. All the table formatting is broken if you just paste in the text!

The error is pretty clear

 254  Pw.up(bus,t)=DataWDL(t,'w')*Wcap(bus)/Sbase;
****                                 $141
**** 141  Symbol declared but no values have been assigned. Check for missing
****         data definition, assignment, data loading or implicit assignment
****         via a solve statement.
****         A wild shot: You may have spurious commas in the explanatory
****         text of a declaration. Check symbol reference list.

You have declared Wcap but never filled it with data. So GAMS complains because it doesn’t know what data to use.

-Michael