*****Main Program: agedasi How to Compile: -> Makefile "make" or 1) p2c -LTURBO agedasi.pas 2) gcc agedasi.c -lm -lp2c -o agedasi NOTE: When translating by P2C, you may have some warning messages and notes, but they are neglectable. Initialization files: 1) para.dat .002 <- delta 40 <- ndiv .01 <- trend1 .03 <- trend2 .02 <- trend3 100 <- popsize 314 <- mingen 523 <- maxgen 105 <- teachgen .3 <- pcross1 .003 <- pmutation1 .8 <- gap1 .3 <- pcross2 .003 <- pmutation2 .8 <- gap2 100 <- ntime 100 <- nsample detail <- rand_mode, 'write' >> random.dat.*, rout.dat.*,eval.dat 'read' << random.dat.* >> wout.dat.*,fout.dat.*,qout.dat.* 'all' >> rout.dat.*,eval.dat,wout.dat.*,fout.dat.*,qout.dat.* 'detail'<< random.dat.* >> all*.dat.*, rout.dat.*,eval.dat,wout.dat.*,fout.dat.*,qout.dat.*,forecast.dat.* 2) filename.dat result/ <- res_dir 10 <- span_detail (if rand_mode='detail') 1 <- ntime 38 <- numbers of paths 3) loading.dat.* 6 <- nfactor 12 <- nload 1 2 3 5 6 8 9 10 12 13 15 17 <- loading[0,1..nload] -0.5256285 0.5009647 0.06243332 0.03961189 0.2006960 -0.08381246 0.02597303 -0.02328062 0.04314966 0.03699096 0.08974927 -0.0008957278 <- loading[1,1..nload] -0.004053159 0.1347315 0.07956871 0.1302897 -0.08148127 0.003575953 -0.05729144 -0.007092795 -0.03815311 -0.4688561 0.5628648 -0.03801513 -0.03359736 -0.1498163 -0.3578819 0.4885089 0.2497583 0.2265632 0.08001726 0.07951225 0.080272 -0.08378752 -0.06947102 0.01537606 0.08698737 0.1431886 -0.02713363 0.1144817 0.3719848 0.2358101 0.4970994 -0.1008574 -0.08589251 0.03567224 -0.03308006 0.07129067 -0.003221098 0.2047597 -0.01898621 0.004368091 0.3662496 -0.4132485 -0.0884039 -0.05141621 0.3751856 0.02951169 -0.02682375 -0.002225342 0.006387463 0.03059961 0.02125821 0.05153005 0.1263902 0.2070009 0.03141974 0.06762324 0.05679817 0.2283047 0.09568169 0.4938448 0.02398835 0.1076032 0.07406759 -0.02280962 -0.08545549 0.0425681 0.1327297 -0.3411588 0.09180624 0.08521358 0.08750241 -0.08087169 <- loading[nfactor,1..nload] NOTES 1) Result417c uses data.org (old data) instead of data.dat (new data). For result417c, change line 93 in agedasi.pas, assign(dfile,'data.dat'); to assign(dfile,'data.org'); And it uses rate.org (old data) instead of rate.dat (new data). For result417c, change line 42 in data.age, assign(wfile,'rate.dat'); to assign(wfile,'rate.org'); 2) For scenario analysis, change lines below line 66 in "data.age" and change like 1). *****Frequency count program: freq How to Compile: -> Makefile 1) p2c freq.pas 2) gcc freq.c -lm -lp2c -o freq Initialization files: 1) freq.ini 627 Start week of forecast period (teachgen) 686 End week (maxgen) 95 Minimum value of range (minval) 145 Maximum value of range (maxval) 25 Number of frequency intervals (count) 1 Start number of simulation file (start) 50 Number of simulation time (ntime) 50 Maximum size of frequency (zmax) result/ Name of the directory which has rout.dat.* (dir) Output files: freq.dat