I searched the help-glpk archive. I found out that I have this problem before,
http://lists.gnu.org/archive/html/help-glpk/2004-09/msg00006.html And the answer is just putting a default 0 to the parameter according to Andrew. http://lists.gnu.org/archive/html/help-glpk/2004-09/msg00008.html Answer to my present question is, putting default 0; param Area_Stand{i in I,j in J,k in K,l in L} default 0; instead of just this, param Area_Stand{i in I,j in J,k in K,l in L}; Thanks, Noli On 12/17/09, Noli Sicad <[email protected]> wrote: > Hi, > > I am trying to read a dbf file (*.dbf) from shapefile (i.e. ESRI > ArcGIS) or GIS file (*.shp) and writing a dbf table. > > I am wondering how can I solve this problem - "no value" problem as > you can see in the console output below, it is not part of the data > set from dbf file. > > ~~~~~~~~~~~~~~~ > Reading tab8... > Display statement at line 30 > Area_Stand[A,R,1,3] = 131.507 > Area_Stand[A,R,2,3] = 122.955 > Area_Stand[A,'',1,1] = 57.899 > Area_Stand[A,'',2,1] = 15.573 > Area_Stand[A,'',3,1] = 78.544 > Area_Stand[A,'',4,1] = 67.906 > Area_Stand[A,'',5,1] = 145.412 > Area_Stand[A,R,1,2] = 155.223 > Area_Stand[A,'',2,2] = 64.358 > Area_Stand[A,'',3,2] = 35.798 > Area_Stand[A,R,4,2] = 124.597 > Area_Stand[B,'',1,3] = 46.048 > Area_Stand[B,'',2,3] = 73.329 > Area_Stand[B,'',3,3] = 40.189 > Area_Stand[B,'',1,1] = 4.184 > Area_Stand[B,'',2,1] = 15.175 > Area_Stand[B,'',3,1] = 52.294 > Area_Stand[B,'',1,2] = 78.506 > Area_Stand[B,'',2,2] = 32.534 > Area_Stand[B,'',3,2] = 22.573 > Area_Stand[B,'',4,2] = 6.566 > Area_Stand[B,'',5,2] = 58.724 > Area_Stand[B,'',4,3] = 69.777 > Area_Stand[C,'',1,3] = 96.944 > Area_Stand[C,'',2,3] = 23.426 > Area_Stand[C,'',3,3] = 43.471 > Area_Stand[C,'',1,1] = 208.631 > Area_Stand[C,'',2,1] = 245.654 > Area_Stand[C,'',3,1] = 43.55 > Area_Stand[C,'',4,2] = 6.934 > Area_Stand[C,'',5,2] = 163.692 > Area_Stand[C,'',6,2] = 72.195 > Area_Stand[C,'',7,2] = 23.371 > Area_Stand[C,'',8,2] = 72.455 > Generating Area... > ReadingDbfGIS_7.mod:41: no value for Area_Stand[A,R,1,1] > MathProg model processing error > ~~~~~~~~~~~~~~~~~~~~~ > > My problem are in here: > > # Problem here > param A{i in I}:= sum {j in J, k in K, l in L} Area_Stand[i,j,k,l]; > > > # Problem here as well. > subject to Area {i in I}: > sum {j in J,k in K,l in L} XForestLand[i,j,k,l] = sum {j in J,k in > K,l in L} Area_Stand[i,j,k,l]; > > How can I put conditions on these two equation e.g if statements in > these to equation. > > Please see the mathprog script and console output below. > > Thanks in advance, Noli > ~~~~~~~~~~~~~~~~~~~~~ > > Mathprog script > ~~~~~~ > #Reading dbf (*,dbf) from shapefile (*..shp), GIS file i.e. ArcGIS > file and writing dbf table. > > /* > Noli Sicad <[email protected]> > > stands8.shp > stands8.dbf > stands8.shx > stabds8.prj > > */ > > set I; > table tab IN "xBASE" "vegtype.dbf": I <- [VEG_TYPE]; > display I; > > set J; > table tab2 IN "xBASE" "vegtype2.dbf": J <- [VEG_TYPE2]; > display J; > > set K; > table tab3 IN "xBASE" "climatec.dbf": K <- [CLIMATE_C]; > display K; > > set L; > table tab4 IN "xBASE" "MGT.dbf": L <- [MANAGEMENT]; > display L; > > param Area_Stand{i in I,j in J,k in K,l in L}; > > set S, dimen 4; > table tab8 IN "xBASE" "stands8.dbf": S <- [VEG_TYPE, > VEG_TYPE2,CLIMATE_C,MANAGEMENT], Area_Stand ~ACRES; > display Area_Stand; > > var XForestLand{i in I,j in J,k in K,l in L}; > > # Problem here > param A{i in I}:= sum {j in J, k in K, l in L} Area_Stand[i,j,k,l]; > > # Problem here as well. > subject to Area {i in I}: > sum {j in J,k in K,l in L} XForestLand[i,j,k,l] = sum {j in J,k in > K,l in L} Area_Stand[i,j,k,l]; > > ~~~~~~~~~~~~~~~ > Console output > > ~~~~~~~~~~~~~~~ > Parameter(s) specified in the command line: > --cover --clique --gomory --mir -m ReadingDbfGIS_7.mod > Reading model section from ReadingDbfGIS_7.mod... > ReadingDbfGIS_7.mod:78: warning: unexpected end of file; missing end > statement inserted > 78 lines were read > Reading tab... > Display statement at line 12 > I: > A > B > C > D > Reading tab2... > Display statement at line 16 > J: > R > '' > Reading tab3... > Display statement at line 20 > K: > 1 > 2 > 3 > 4 > 5 > 6 > 7 > 8 > Reading tab4... > Display statement at line 24 > L: > 1 > 2 > 3 > Reading tab8... > Display statement at line 30 > Area_Stand[A,R,1,3] = 131.507 > Area_Stand[A,R,2,3] = 122.955 > Area_Stand[A,'',1,1] = 57.899 > Area_Stand[A,'',2,1] = 15.573 > Area_Stand[A,'',3,1] = 78.544 > Area_Stand[A,'',4,1] = 67.906 > Area_Stand[A,'',5,1] = 145.412 > Area_Stand[A,R,1,2] = 155.223 > Area_Stand[A,'',2,2] = 64.358 > Area_Stand[A,'',3,2] = 35.798 > Area_Stand[A,R,4,2] = 124.597 > Area_Stand[B,'',1,3] = 46.048 > Area_Stand[B,'',2,3] = 73.329 > Area_Stand[B,'',3,3] = 40.189 > Area_Stand[B,'',1,1] = 4.184 > Area_Stand[B,'',2,1] = 15.175 > Area_Stand[B,'',3,1] = 52.294 > Area_Stand[B,'',1,2] = 78.506 > Area_Stand[B,'',2,2] = 32.534 > Area_Stand[B,'',3,2] = 22.573 > Area_Stand[B,'',4,2] = 6.566 > Area_Stand[B,'',5,2] = 58.724 > Area_Stand[B,'',4,3] = 69.777 > Area_Stand[C,'',1,3] = 96.944 > Area_Stand[C,'',2,3] = 23.426 > Area_Stand[C,'',3,3] = 43.471 > Area_Stand[C,'',1,1] = 208.631 > Area_Stand[C,'',2,1] = 245.654 > Area_Stand[C,'',3,1] = 43.55 > Area_Stand[C,'',4,2] = 6.934 > Area_Stand[C,'',5,2] = 163.692 > Area_Stand[C,'',6,2] = 72.195 > Area_Stand[C,'',7,2] = 23.371 > Area_Stand[C,'',8,2] = 72.455 > Generating Area... > ReadingDbfGIS_7.mod:41: no value for Area_Stand[A,R,1,1] > MathProg model processing error >>Exit code: 1 Time: 0.209 > _______________________________________________ Help-glpk mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-glpk
