1PAGE 1 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND * * * * * S O U R C E L I S T I N G * * * * * 0001 TITLE WORLD POPULATION SECTOR IN NDTRAN VERSION TWO 00000010 0002 * EULER 00000050 0003 * SYMBOL 00000060 0004 * NOSTATS 00000070 0005 * CHECK 00000080 0006 NOTE LIMITS TO GROWTH BY DONELLA MEADOWS ET.AL. 00000020 0007 NOTE THIS IS THE VERSION USED IN THE TECHNICAL REPORT 00000030 0008 NOTE LAST REVISED DECEMBER 14, L972 00000040 0009 NOTE 00000090 0010 NOTE FOUR POPULATION LEVELS 00000100 0011 NOTE 00000110 0012 A POP.K=P1.K+P2.K+P3.K+P4.K 00000120 0013 L P1.K=INTGRL(B.JK-D1.JK-MAT1.JK) 00000130 0014 N P1=P1I 00000140 0015 C P1I=65E7 00000150 0016 R D1.KL=AD1.K 00000160 0017 A AD1.K=P1.K*M1.K 00000170 0018 A M1.K=TABHL(M1T,LE.K,20,70,10) 00000180 0019 T M1T=.0567,.0366,.0243,.0155,.0082,.0023 00000190 0020 R MAT1.KL=(P1.K)(1-M1.K)/15 00000200 0021 L P2.K=INTGRL(MAT1.JK-D2.JK-MAT2.JK) 00000210 0022 N P2=P2I 00000220 0023 C P2I=70E7 00000230 0024 R D2.KL=AD2.K 00000240 0025 A AD2.K=P2.K*M2.K 00000250 0026 A M2.K=TABHL(M2T,LE.K,20,70,10) 00000260 0027 T M2T=.0266,.0171,.0110,.0065,.0040,.0016 00000270 0028 R MAT2.KL=(P2.K)(1-M2.K)/30 00000280 0029 L P3.K=INTGRL(MAT2.JK-D3.JK-MAT3.JK) 00000290 0030 N P3=P3I 00000300 0031 C P3I=19E7 00000310 0032 R D3.KL=AD3.K 00000320 0033 A AD3.K=P3.K*M3.K 00000330 0034 A M3.K=TABHL(M3T,LE.K,20,70,10) 00000340 0035 T M3T=.0562,.0373,.0252,.0171,.0118,.0083 00000350 0036 R MAT3.KL=(P3.K)(1-M3.K)/20 00000360 0037 L P4.K=INTGRL(MAT3.JK-D4.JK) 00000370 0038 N P4=P4I 00000380 0039 C P4I=6E7 00000390 0040 R D4.KL=AD4.K 00000400 0041 A AD4.K=P4.K*M4.K 00000410 0042 A M4.K=TABHL(M4T,LE.K,20,70,10) 00000420 0043 T M4T=.13,.11,.09,.07,.07,.07 00000430 0044 NOTE 00000440 0045 NOTE DEATH RATE EQUATIONS 00000450 0046 NOTE 00000460 0047 A AD.K=AD1.K+AD2.K+AD3.K+AD4.K 00000470 0048 S CDR.K=1000*AD.K/POP.K 00000480 $ 1) ***** W A R N I N G ***** ND0588 0049 A LE.K=LEN*LMF.K*LMHS.K*LMP.K*LMC.K 00000490 0050 C LEN=28 00000500 0051 A LMF.K=TABHL(LMFT,FPC.K/SFPC,0,5,1) 00000510 0052 T LMFT=0,1,1.2,1.3,1.35,1.4 00000520 1PAGE 2 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND 0053 C SFPC=230 00000530 0054 A LMHS.K=CLIP(LMHS2.K,LMHS1.K,TIME.K,IPHST) 00000540 0055 C IPHST=40 00000550 0056 A LMHS1.K=TABHL(LMHS1T,EHSPC.K,0,100,20) 00000560 0057 T LMHS1T=1,1.1,1.4,1.6,1.7,1.8 00000570 0058 A LMHS2.K=TABHL(LMHS2T,EHSPC.K,0,100,20) 00000580 0059 T LMHS2T=1,1.4,1.6,1.8,1.95,2.0 00000590 0060 EXPND SMOOTH(EHSPC,HSAPC,HSID) 00000600 0061+ L $L11.K=INTGRL(HSAPC.JK-$R11.JK) 0062+ N $L11=HSAPC*HSID 0063+ A EHSPC.K=$L11.K/HSID 0064+ R $R11.KL=EHSPC.K 0065+ MEND 0066 C HSID=20 00000610 0067 R HSAPC.KL=HSPC.K 00000620 0068 A HSPC.K=TABHL(HSAPCT,SOPC.K,0,2000,250) 00000630 0069 T HSAPCT=0,20,50,95,140,175,200,220,230 00000640 0070 A LMP.K=TABHL(LMPT,PPOLX.K,0,100,10) 00000650 0071 T LMPT=1.0,.99,.97,.95,.90,.85,.75,.65,.55,.40,.20 00000660 0072 A FPU.K=TABHL(FPUT,POP.K,0,16E9,2E9) 00000670 0073 T FPUT=0,.2,.4,.5,.58,.65,.72,.78,.80 00000680 0074 A LMC.K=1-(CMI.K*FPU.K) 00000690 0075 A CMI.K=TABHL(CMIT,IOPC1.K,0,1600,200) 00000700 0076 T CMIT=.5,.05,-.1,-.08,-.02,.05,.1,.15,.2 00000710 0077 NOTE 00000720 0078 NOTE BIRTH RATE EQUATIONS 00000730 0079 NOTE 00000740 0080 R B.KL=CLIP(AD.K,(TF.K*P2.K*0.5/RLT),TIME.K,PET) 00000750 0081 C RLT=30 00000760 0082 C PET=4000 00000770 0083 S CBR.K=1000*B.JK/POP.K 00000780 0084 A TF.K=MIN(MTF.K,(MTF.K*(1-FCE.K)+DTF.K*FCE.K)) 00000790 0085 A MTF.K=MTFN*FM.K 00000800 0086 C MTFN=12 00000810 0087 A FM.K=TABHL(FMT,LE.K,0,80,10) 00000820 0088 T FMT=0,.2,.4,.6,.8,.9,1,1.05,1.1 00000830 0089 A DTF.K=DCFS.K*CMPLE.K 00000840 0090 A DCFS.K=CLIP(2,DCFSN*FRSN.K,TIME.K,ZPGT) 00000850 0091 C ZPGT=4000 00000860 0092 C DCFSN=4 00000870 0093 A FRSN.K=TABHL(FRSNT,FIE.K,-.2,.2,.1) 00000880 0094 T FRSNT=.5,.6,.7,.85,1 00000890 0095 A FIE.K=(IOPC1.K/AIOPC.K) 00000900 0096 EXPND SMOOTH(AIOPC,IOPC,IEAT) 00000910 0097+ L $L12.K=INTGRL(IOPC.JK-$R12.JK) 0098+ N $L12=IOPC*IEAT 0099+ A AIOPC.K=$L12.K/IEAT 0100+ R $R12.KL=AIOPC.K 0101+ MEND 0102 C IEAT=3 00000920 0103 A SFSN.K=TABHL(SFSNT,DIOPC1.K,0,800,200) 00000930 $ 1) ***** W A R N I N G ***** ND0588 0104 T SFSNT=1.25,1,.9,.8,.75 00000940 0105 EXPND DLINF3(DIOPC1,IOPC1,SAD) 00000950 0106+ R $R13.KL=(IOPC1.K-$L13.K)/(SAD/3) 1PAGE 3 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND 0107+ L $L13.K=INTGRL($R13.JK) 0108+ N $L13=IOPC1 0109+ R $R23.KL=($L13.K-$L23.K)/(SAD/3) 0110+ L $L23.K=INTGRL($R23.JK) 0111+ N $L23=IOPC1 0112+ R $R33.KL=($L23.K-$L33.K)/(SAD/3) 0113+ L $L33.K=INTGRL($R33.JK) 0114+ N $L33=IOPC1 0115+ A DIOPC1.K=$L33.K 0116+ MEND 0117 C SAD=20 00000960 0118 A CMPLE.K=TABHL(CMPLET,PLE.K,0,80,10) 00000970 0119 T CMPLET=3,2.1,1.6,1.4,1.3,1.2,1.1,1.05,1 00000980 0120 EXPND DLINF3(PLE,LE,LPD) 00000990 0121+ R $R14.KL=(LE.K-$L14.K)/(LPD/3) 0122+ L $L14.K=INTGRL($R14.JK) 0123+ N $L14=LE 0124+ R $R24.KL=($L14.K-$L24.K)/(LPD/3) 0125+ L $L24.K=INTGRL($R24.JK) 0126+ N $L24=LE 0127+ R $R34.KL=($L24.K-$L34.K)/(LPD/3) 0128+ L $L34.K=INTGRL($R34.JK) 0129+ N $L34=LE 0130+ A PLE.K=$L34.K 0131+ MEND 0132 C LPD=20 00001000 0133 A FCE.K=CLIP(1.0,(TABHL(FCET,FCFPC.K,0,3,.5)),TIME.K,FCEST) 00001010 0134 C FCEST=4000 00001020 0135 T FCET=.75,.85,.9,.95,.98,.99,1 00001030 0136 EXPND DLINF3(FCFPC,FCAPC,HSID) 00001040 0137+ R $R15.KL=(FCAPC.K-$L15.K)/(HSID/3) 0138+ L $L15.K=INTGRL($R15.JK) 0139+ N $L15=FCAPC 0140+ R $R25.KL=($L15.K-$L25.K)/(HSID/3) 0141+ L $L25.K=INTGRL($R25.JK) 0142+ N $L25=FCAPC 0143+ R $R35.KL=($L25.K-$L35.K)/(HSID/3) 0144+ L $L35.K=INTGRL($R35.JK) 0145+ N $L35=FCAPC 0146+ A FCFPC.K=$L35.K 0147+ MEND 0148 A FCAPC.K=FSAFC.K*SOPC.K 00001050 0149 A FSAFC.K=TABHL(FSAFCT,NFC.K,0,10,2) 00001060 0150 T FSAFCT=0,.005,.015,.025,.03,.035 00001070 0151 A NFC.K=(MTF.K/DTF.K)-1 00001080 0152 NOTE 00001090 0153 NOTE INTERFACE AND PROGRAM VARIABLES 00001100 0154 NOTE 00001110 0155 A IO.K=CLIP(IO2,IO1.K,TIME.K,LT) 00001120 0156 C LT=500 00001130 0157 A IO1.K=CLIP(IO12.K,IO11.K,TIME.K,LT2) 00001140 0158 C LT2=500 00001150 0159 A IO11.K=.7E11*EXP(TIME.K*.037) 00001160 0160 A IO12.K=POP.K*CIO 00001170 0161 C CIO=100 00001180 0162 C IO2=30E11 00001190 0163 R IOPC.KL=IOPC1.K 00001200 1PAGE 4 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND 0164 A IOPC1.K=IO.K/POP.K 00001210 0165 A PPOLX.K=1+RAMP(PS,PT) 00001220 0166 C PS=0 00001230 0167 C PT=10 00001240 0168 A SO.K=CLIP(SO2,SO1.K,TIME.K,LT) 00001250 0169 A SO1.K=CLIP(SO12.K,SO11.K,TIME.K,LT2) 00001260 0170 A SO11.K=1.5E11*EXP(TIME.K*.030) 00001270 0171 A SO12.K=POP.K*CSO 00001280 0172 C CSO=150 00001290 0173 C SO2=32E11 00001300 0174 A SOPC.K=SO.K/POP.K 00001310 0175 A F.K=CLIP(F2,F1.K,TIME.K,LT) 00001320 0176 A F1.K=CLIP(F12.K,F11.K,TIME.K,LT2) 00001330 0177 A F11.K=4.0E11*EXP(TIME.K*.020) 00001340 0178 A F12.K=POP.K*CFOOD 00001350 0179 C CFOOD=250 00001360 0180 C F2=24E11 00001370 0181 A FPC.K=F.K/POP.K 00001380 0182 A POP1.K=P1.K/POP.K 00001390 $ 1) ***** W A R N I N G ***** ND0588 0183 A POP2.K=P2.K/POP.K 00001400 $ 1) ***** W A R N I N G ***** ND0588 0184 A POP3.K=P3.K/POP.K 00001410 $ 1) ***** W A R N I N G ***** ND0588 0185 A POP4.K=P4.K/POP.K 00001420 $ 1) ***** W A R N I N G ***** ND0588 0186 PARM DT=1 00001430 0187 PARM START=1900 00001440 0188 PARM STOP=2100 00001450 0189 PARM PRTPER=5 00001460 0190 PRINT SOPC,IOPC,FPC,POP,CBR 00001470 1PAGE 5 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND * * * * * * S Y M B O L T A B L E * * * * * * VARIABLE NAME VARIABLE TYPE OUTPUT REQUEST VARIABLE NUMBER $L11 LEVEL 96 $L12 LEVEL 135 $L13 LEVEL 150 $L14 LEVEL 180 $L15 LEVEL 212 $L23 LEVEL 160 $L24 LEVEL 190 $L25 LEVEL 221 $L33 LEVEL 169 $L34 LEVEL 199 $L35 LEVEL 230 $R11 RATE 98 $R12 RATE 137 $R13 RATE 149 $R14 RATE 179 $R15 RATE 210 $R23 RATE 159 $R24 RATE 189 $R25 RATE 220 $R33 RATE 168 $R34 RATE 198 $R35 RATE 229 AD AUXILIARY 80 AD1 AUXILIARY 38 AD2 AUXILIARY 52 AD3 AUXILIARY 65 AD4 AUXILIARY 77 AIOPC AUXILIARY 134 B RATE 27 CBR SUPPLEMENTARY PRINT 120 CDR SUPPLEMENTARY 81 CFOOD CONSTANT 270 CIO CONSTANT 248 CMI AUXILIARY 114 CMIT TABLE 115 CMPLE AUXILIARY 128 CMPLET TABLE 177 CSO CONSTANT 264 DCFS AUXILIARY 127 DCFSN CONSTANT 129 DIOPC1 AUXILIARY 148 DT PARAMETER 11 DTF AUXILIARY 123 D1 RATE 28 D2 RATE 42 D3 RATE 55 D4 RATE 68 EHSPC AUXILIARY 94 F AUXILIARY 265 FCAPC AUXILIARY 211 FCE AUXILIARY 122 FCEST CONSTANT 209 FCET TABLE 207 FCFPC AUXILIARY 208 FIE AUXILIARY 133 FM AUXILIARY 125 FMT TABLE 126 FPC AUXILIARY PRINT 88 FPU AUXILIARY 112 FPUT TABLE 113 FRSN AUXILIARY 130 FRSNT TABLE 132 FSAFC AUXILIARY 238 FSAFCT TABLE 239 F1 AUXILIARY 267 F11 AUXILIARY 269 F12 AUXILIARY 268 F2 CONSTANT 266 HSAPC RATE 97 HSAPCT TABLE 108 HSID CONSTANT 106 HSPC AUXILIARY 107 IEAT CONSTANT 145 IO AUXILIARY 241 IOPC RATE PRINT 136 IOPC1 AUXILIARY 116 IO1 AUXILIARY 243 IO11 AUXILIARY 246 IO12 AUXILIARY 245 IO2 CONSTANT 242 IPHST CONSTANT 92 LE AUXILIARY 41 LEN CONSTANT 82 LMC AUXILIARY 86 LMF AUXILIARY 83 LMFT TABLE 87 LMHS AUXILIARY 84 LMHS1 AUXILIARY 91 LMHS1T TABLE 93 LMHS2 AUXILIARY 90 LMHS2T TABLE 95 LMP AUXILIARY 85 LMPT TABLE 110 LPD CONSTANT 181 LT CONSTANT 244 LT2 CONSTANT 247 MAT1 RATE 29 MAT2 RATE 43 MAT3 RATE 56 MTF AUXILIARY 121 MTFN CONSTANT 124 M1 AUXILIARY 39 M1T TABLE 40 M2 AUXILIARY 53 M2T TABLE 54 M3 AUXILIARY 66 M3T TABLE 67 M4 AUXILIARY 78 M4T TABLE 79 NFC AUXILIARY 240 PET CONSTANT 119 PLE AUXILIARY 178 PLTPER PARAMETER 16 POP AUXILIARY PRINT 22 POP1 AUXILIARY 271 POP2 AUXILIARY 272 POP3 AUXILIARY 273 POP4 AUXILIARY 274 PPOLX AUXILIARY 111 PRTPER PARAMETER 15 PS CONSTANT 249 PT CONSTANT 250 P1 LEVEL 23 P1I CONSTANT 37 P2 LEVEL 24 P2I CONSTANT 51 P3 LEVEL 25 P3I CONSTANT 64 P4 LEVEL 26 P4I CONSTANT 76 RLT CONSTANT 118 SAD CONSTANT 151 SFPC CONSTANT 89 SFSN AUXILIARY 146 SFSNT TABLE 147 SO AUXILIARY 259 SOPC AUXILIARY PRINT 109 SO1 AUXILIARY 261 SO11 AUXILIARY 263 SO12 AUXILIARY 262 SO2 CONSTANT 260 START PARAMETER 14 STOP PARAMETER 13 TF AUXILIARY 117 TIME LEVEL 12 ZPGT CONSTANT 131 1PAGE 6 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND ERRORS DURING RUN 1: OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1900 ERRORS DURING RUN 1: OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1900 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1901 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1902 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1903 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1904 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1905 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1906 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1907 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1908 OVERFLOW HAS OCCURRED IN STATEMENT 0159 AT TIME = 1909 RUN 1PAGE 7 WORLD POPULATION SECTOR IN NDTRAN VERSION TWO (C) 1978 UND TIME SOPC IOPC FPC POP CBR E+03 E+03 E+03 E+03 E+09 E+00 1.9000 2.0000 1.8750 1.5000 1.6000 29.777 1.9050 1.7743 1.6634 1.3307 1.8036 29.824