TITLE WORLD POPULATION SECTOR IN NDTRAN VERSION TWO 00000010 * EULER 00000050 * SYMBOL 00000060 * NOSTATS 00000070 * CHECK 00000080 NOTE LIMITS TO GROWTH BY DONELLA MEADOWS ET.AL. 00000020 NOTE THIS IS THE VERSION USED IN THE TECHNICAL REPORT 00000030 NOTE LAST REVISED DECEMBER 14, L972 00000040 NOTE 00000090 NOTE FOUR POPULATION LEVELS 00000100 NOTE 00000110 A POP.K=P1.K+P2.K+P3.K+P4.K 00000120 L P1.K=INTGRL(B.JK-D1.JK-MAT1.JK) 00000130 N P1=P1I 00000140 C P1I=65E7 00000150 R D1.KL=AD1.K 00000160 A AD1.K=P1.K*M1.K 00000170 A M1.K=TABHL(M1T,LE.K,20,70,10) 00000180 T M1T=.0567,.0366,.0243,.0155,.0082,.0023 00000190 R MAT1.KL=(P1.K)(1-M1.K)/15 00000200 L P2.K=INTGRL(MAT1.JK-D2.JK-MAT2.JK) 00000210 N P2=P2I 00000220 C P2I=70E7 00000230 R D2.KL=AD2.K 00000240 A AD2.K=P2.K*M2.K 00000250 A M2.K=TABHL(M2T,LE.K,20,70,10) 00000260 T M2T=.0266,.0171,.0110,.0065,.0040,.0016 00000270 R MAT2.KL=(P2.K)(1-M2.K)/30 00000280 L P3.K=INTGRL(MAT2.JK-D3.JK-MAT3.JK) 00000290 N P3=P3I 00000300 C P3I=19E7 00000310 R D3.KL=AD3.K 00000320 A AD3.K=P3.K*M3.K 00000330 A M3.K=TABHL(M3T,LE.K,20,70,10) 00000340 T M3T=.0562,.0373,.0252,.0171,.0118,.0083 00000350 R MAT3.KL=(P3.K)(1-M3.K)/20 00000360 L P4.K=INTGRL(MAT3.JK-D4.JK) 00000370 N P4=P4I 00000380 C P4I=6E7 00000390 R D4.KL=AD4.K 00000400 A AD4.K=P4.K*M4.K 00000410 A M4.K=TABHL(M4T,LE.K,20,70,10) 00000420 T M4T=.13,.11,.09,.07,.07,.07 00000430 NOTE 00000440 NOTE DEATH RATE EQUATIONS 00000450 NOTE 00000460 A AD.K=AD1.K+AD2.K+AD3.K+AD4.K 00000470 S CDR.K=1000*AD.K/POP.K 00000480 A LE.K=LEN*LMF.K*LMHS.K*LMP.K*LMC.K 00000490 C LEN=28 00000500 A LMF.K=TABHL(LMFT,FPC.K/SFPC,0,5,1) 00000510 T LMFT=0,1,1.2,1.3,1.35,1.4 00000520 C SFPC=230 00000530 A LMHS.K=CLIP(LMHS2.K,LMHS1.K,TIME.K,IPHST) 00000540 C IPHST=40 00000550 A LMHS1.K=TABHL(LMHS1T,EHSPC.K,0,100,20) 00000560 T LMHS1T=1,1.1,1.4,1.6,1.7,1.8 00000570 A LMHS2.K=TABHL(LMHS2T,EHSPC.K,0,100,20) 00000580 T LMHS2T=1,1.4,1.6,1.8,1.95,2.0 00000590 EXPND SMOOTH(EHSPC,HSAPC,HSID) 00000600 C HSID=20 00000610 R HSAPC.KL=HSPC.K 00000620 A HSPC.K=TABHL(HSAPCT,SOPC.K,0,2000,250) 00000630 T HSAPCT=0,20,50,95,140,175,200,220,230 00000640 A LMP.K=TABHL(LMPT,PPOLX.K,0,100,10) 00000650 T LMPT=1.0,.99,.97,.95,.90,.85,.75,.65,.55,.40,.20 00000660 A FPU.K=TABHL(FPUT,POP.K,0,16E9,2E9) 00000670 T FPUT=0,.2,.4,.5,.58,.65,.72,.78,.80 00000680 A LMC.K=1-(CMI.K*FPU.K) 00000690 A CMI.K=TABHL(CMIT,IOPC1.K,0,1600,200) 00000700 T CMIT=.5,.05,-.1,-.08,-.02,.05,.1,.15,.2 00000710 NOTE 00000720 NOTE BIRTH RATE EQUATIONS 00000730 NOTE 00000740 R B.KL=CLIP(AD.K,(TF.K*P2.K*0.5/RLT),TIME.K,PET) 00000750 C RLT=30 00000760 C PET=4000 00000770 S CBR.K=1000*B.JK/POP.K 00000780 A TF.K=MIN(MTF.K,(MTF.K*(1-FCE.K)+DTF.K*FCE.K)) 00000790 A MTF.K=MTFN*FM.K 00000800 C MTFN=12 00000810 A FM.K=TABHL(FMT,LE.K,0,80,10) 00000820 T FMT=0,.2,.4,.6,.8,.9,1,1.05,1.1 00000830 A DTF.K=DCFS.K*CMPLE.K 00000840 A DCFS.K=CLIP(2,DCFSN*FRSN.K,TIME.K,ZPGT) 00000850 C ZPGT=4000 00000860 C DCFSN=4 00000870 A FRSN.K=TABHL(FRSNT,FIE.K,-.2,.2,.1) 00000880 T FRSNT=.5,.6,.7,.85,1 00000890 A FIE.K=(IOPC1.K/AIOPC.K) 00000900 EXPND SMOOTH(AIOPC,IOPC,IEAT) 00000910 C IEAT=3 00000920 A SFSN.K=TABHL(SFSNT,DIOPC1.K,0,800,200) 00000930 T SFSNT=1.25,1,.9,.8,.75 00000940 EXPND DLINF3(DIOPC1,IOPC1,SAD) 00000950 C SAD=20 00000960 A CMPLE.K=TABHL(CMPLET,PLE.K,0,80,10) 00000970 T CMPLET=3,2.1,1.6,1.4,1.3,1.2,1.1,1.05,1 00000980 EXPND DLINF3(PLE,LE,LPD) 00000990 C LPD=20 00001000 A FCE.K=CLIP(1.0,(TABHL(FCET,FCFPC.K,0,3,.5)),TIME.K,FCEST) 00001010 C FCEST=4000 00001020 T FCET=.75,.85,.9,.95,.98,.99,1 00001030 EXPND DLINF3(FCFPC,FCAPC,HSID) 00001040 A FCAPC.K=FSAFC.K*SOPC.K 00001050 A FSAFC.K=TABHL(FSAFCT,NFC.K,0,10,2) 00001060 T FSAFCT=0,.005,.015,.025,.03,.035 00001070 A NFC.K=(MTF.K/DTF.K)-1 00001080 NOTE 00001090 NOTE INTERFACE AND PROGRAM VARIABLES 00001100 NOTE 00001110 A IO.K=CLIP(IO2,IO1.K,TIME.K,LT) 00001120 C LT=500 00001130 A IO1.K=CLIP(IO12.K,IO11.K,TIME.K,LT2) 00001140 C LT2=500 00001150 A IO11.K=.7E11*EXP(TIME.K*.037) 00001160 A IO12.K=POP.K*CIO 00001170 C CIO=100 00001180 C IO2=30E11 00001190 R IOPC.KL=IOPC1.K 00001200 A IOPC1.K=IO.K/POP.K 00001210 A PPOLX.K=1+RAMP(PS,PT) 00001220 C PS=0 00001230 C PT=10 00001240 A SO.K=CLIP(SO2,SO1.K,TIME.K,LT) 00001250 A SO1.K=CLIP(SO12.K,SO11.K,TIME.K,LT2) 00001260 A SO11.K=1.5E11*EXP(TIME.K*.030) 00001270 A SO12.K=POP.K*CSO 00001280 C CSO=150 00001290 C SO2=32E11 00001300 A SOPC.K=SO.K/POP.K 00001310 A F.K=CLIP(F2,F1.K,TIME.K,LT) 00001320 A F1.K=CLIP(F12.K,F11.K,TIME.K,LT2) 00001330 A F11.K=4.0E11*EXP(TIME.K*.020) 00001340 A F12.K=POP.K*CFOOD 00001350 C CFOOD=250 00001360 C F2=24E11 00001370 A FPC.K=F.K/POP.K 00001380 A POP1.K=P1.K/POP.K 00001390 A POP2.K=P2.K/POP.K 00001400 A POP3.K=P3.K/POP.K 00001410 A POP4.K=P4.K/POP.K 00001420 PARM DT=1 00001430 PARM START=1900 00001440 PARM STOP=2100 00001450 PARM PRTPER=5 00001460 PRINT SOPC,IOPC,FPC,POP,CBR 00001470