Trailing-Edge
-
PDP-10 Archives
-
decuslib20-06
-
decus/20-160/pop.dat
There are no other files named pop.dat in the archive.
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