Google
 

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