Google
 

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