| 1 |
jscott |
1.1 |
|
| 2 |
|
|
#include "ctrparam.h" |
| 3 |
|
|
|
| 4 |
|
|
! ========================================================== |
| 5 |
|
|
! |
| 6 |
|
|
! TRENDS.F: This file contains subroitines a, b, & ctrend |
| 7 |
|
|
! in the old MD2G04.f. |
| 8 |
|
|
! |
| 9 |
|
|
! ---------------------------------------------------------- |
| 10 |
|
|
! |
| 11 |
|
|
! Revision History: |
| 12 |
|
|
! |
| 13 |
|
|
! When Who What |
| 14 |
|
|
! ---- ---------- ------- |
| 15 |
|
|
! 073100 Chien Wang repack based on CliChem3 & M24x11, |
| 16 |
|
|
! and add cpp. |
| 17 |
|
|
! |
| 18 |
|
|
! ========================================================== |
| 19 |
|
|
|
| 20 |
|
|
|
| 21 |
|
|
SUBROUTINE ATREND(XGAS,YEAR,NGAS) 9139.1 |
| 22 |
|
|
C 9139.2 |
| 23 |
|
|
C-----------------------------------------------------------------------9139.3 |
| 24 |
|
|
C T-GAS SCENARIO A: PRESENT TREND FOR TRACE GAS INCREASES 9139.4 |
| 25 |
|
|
C (EXPONENTIAL FORCING PROJECTED BEYOND 2000) 9139.5 |
| 26 |
|
|
C 9139.6 |
| 27 |
|
|
C CO2: OBSERVED TREND 1958-1984, 9139.7 |
| 28 |
|
|
C 0.015/YR RATE OF INCREASE 9139.8 |
| 29 |
|
|
C ON 1.5 PPM INCREMENT BEYOND 1984 9139.9 |
| 30 |
|
|
C 9140. |
| 31 |
|
|
C N2O: WEISS FORMULA THRU 1978, 9140.1 |
| 32 |
|
|
C 0.035/YR INCREASE 9140.2 |
| 33 |
|
|
C ON 1978 INCREMENT BEYOND 1978 9140.3 |
| 34 |
|
|
C 9140.4 |
| 35 |
|
|
C CH4: OBSERVED RATES THRU 1980, 9140.5 |
| 36 |
|
|
C 0.015/YR INCREASE BEYOND 1980 9140.6 |
| 37 |
|
|
C 9140.7 |
| 38 |
|
|
C CFC: OBSERVED CONCENTRATION THRU 1984, 9140.8 |
| 39 |
|
|
C 0.03 EMISSION RATE INCREASE BEYOND 1984, 9140.9 |
| 40 |
|
|
C DOUBLED FORCING TO ACCOUNT FOR OTHER CFCS 9141. |
| 41 |
|
|
C SLOW 75YR, 150YR DECAY RATE BEYOND 2000 9141.1 |
| 42 |
|
|
C 9141.2 |
| 43 |
|
|
C-----------------------------------------------------------------------9141.3 |
| 44 |
|
|
C 9141.4 |
| 45 |
|
|
DIMENSION XGAS(1) 9141.5 |
| 46 |
|
|
DIMENSION XXF11(45),XXF12(45),XXCO2(27) 9141.6 |
| 47 |
|
|
DATA XXF11/ 9141.7 |
| 48 |
|
|
10.00001,0.00001,0.00001,0.00001,0.00001,0.00004,0.00010,0.00023, 9141.8 |
| 49 |
|
|
20.00041,0.00066,0.00100,0.00148,0.00215,0.00297,0.00397,0.00523, 9141.9 |
| 50 |
|
|
30.00662,0.00791,0.00921,0.01093,0.01316,0.01596,0.01939,0.02345, 9142. |
| 51 |
|
|
40.02807,0.03321,0.03904,0.04569,0.05343,0.06224,0.07186,0.08267, 9142.1 |
| 52 |
|
|
50.09505,0.10867,0.12174,0.13422,0.14665,0.15820,0.16871,0.17841, 9142.2 |
| 53 |
|
|
60.18919,0.19983,0.21032,0.22068,0.23090/ 9142.3 |
| 54 |
|
|
DATA XXF12/ 9142.4 |
| 55 |
|
|
10.00001,0.00001,0.00001,0.00040,0.00112,0.00198,0.00298,0.00411, 9142.5 |
| 56 |
|
|
20.00539,0.00679,0.00836,0.00998,0.01181,0.01387,0.01619,0.01888, 9142.6 |
| 57 |
|
|
30.02195,0.02513,0.02868,0.03292,0.03766,0.04312,0.04954,0.05702, 9142.7 |
| 58 |
|
|
40.06546,0.07485,0.08545,0.09756,0.11130,0.12634,0.14250,0.16011, 9142.8 |
| 59 |
|
|
50.17962,0.20082,0.22144,0.24105,0.25951,0.27630,0.29247,0.30831, 9142.9 |
| 60 |
|
|
60.32585,0.34328,0.36059,0.37778,0.39486/ 9143. |
| 61 |
|
|
DATA XXCO2/ 9143.1 |
| 62 |
|
|
1 315.00, 315.69, 316.66, 317.27, 318.24, 318.93, 319.33, 319.95, 9143.2 |
| 63 |
|
|
2 320.80, 321.44, 322.25, 323.44, 324.82, 326.01, 326.97, 328.86, 9143.3 |
| 64 |
|
|
3 330.35, 330.71, 331.69, 332.83, 334.58, 335.95, 337.65, 339.04, 9143.4 |
| 65 |
|
|
4 340.01, 341.56, 343.81/ 9143.5 |
| 66 |
|
|
C 9143.6 |
| 67 |
|
|
DO 9000 N=1,NGAS 9143.7 |
| 68 |
|
|
GO TO (100,200,300,400,500),N 9143.8 |
| 69 |
|
|
C 9143.9 |
| 70 |
|
|
C CO2 SCENARIOS 9144. |
| 71 |
|
|
C------------------------------------------------------------------- 9144.1 |
| 72 |
|
|
100 CONTINUE 9144.2 |
| 73 |
|
|
IF(YEAR.GT.1958.0) GO TO 110 9144.3 |
| 74 |
|
|
DT=(YEAR-1850.0)/(1958.0-1850.0) 9144.4 |
| 75 |
|
|
XX=270.0*(315.0/270.0)**DT 9144.5 |
| 76 |
|
|
GO TO 120 9144.6 |
| 77 |
|
|
110 IF(YEAR.GT.1984.0) GO TO 115 9144.7 |
| 78 |
|
|
I=YEAR 9144.8 |
| 79 |
|
|
DELTA=YEAR-I 9144.9 |
| 80 |
|
|
I=I-1957 9145. |
| 81 |
|
|
XX=XXCO2(I)+DELTA*(XXCO2(I+1)-XXCO2(I)) 9145.1 |
| 82 |
|
|
GO TO 120 9145.2 |
| 83 |
|
|
115 DT=YEAR-1984.0 9145.3 |
| 84 |
|
|
XX0=343.81 9145.4 |
| 85 |
|
|
ETA=1.50 9145.5 |
| 86 |
|
|
XMU=0.015 9145.6 |
| 87 |
|
|
XX=XX0+ETA*(EXP(XMU*DT)-1.0)/XMU 9145.7 |
| 88 |
|
|
120 CONTINUE 9145.8 |
| 89 |
|
|
GO TO 8000 9145.9 |
| 90 |
|
|
C 9146. |
| 91 |
|
|
C N2O SCENARIOS 9146.1 |
| 92 |
|
|
C------------------------------------------------------------------- 9146.2 |
| 93 |
|
|
200 XSS=282.0 9146.3 |
| 94 |
|
|
TAU=150.0 9146.4 |
| 95 |
|
|
BNA=16.7 9146.5 |
| 96 |
|
|
IF(YEAR.GT.1912.0) GO TO 205 9146.6 |
| 97 |
|
|
X0=285.355 9146.7 |
| 98 |
|
|
XMU=0.04 9146.8 |
| 99 |
|
|
T0=1912.0 9146.9 |
| 100 |
|
|
XJ0=11.0*EXP(-0.01*33.0-0.035*33.0) 9147. |
| 101 |
|
|
GO TO 219 9147.1 |
| 102 |
|
|
205 IF(YEAR.GT.1945.0) GO TO 210 9147.2 |
| 103 |
|
|
X0=289.96 9147.3 |
| 104 |
|
|
XMU=0.01 9147.4 |
| 105 |
|
|
T0=1945.0 9147.5 |
| 106 |
|
|
XJ0=11.0*EXP(-0.035*33.0) 9147.6 |
| 107 |
|
|
GO TO 219 9147.7 |
| 108 |
|
|
210 IF(YEAR.GT.1978.0) GO TO 215 9147.8 |
| 109 |
|
|
X0=300.2 9147.9 |
| 110 |
|
|
XMU=0.035 9148. |
| 111 |
|
|
T0=1978.0 9148.1 |
| 112 |
|
|
XJ0=11.0 9148.2 |
| 113 |
|
|
GO TO 219 9148.3 |
| 114 |
|
|
215 CONTINUE 9148.4 |
| 115 |
|
|
X0=300.2 9148.5 |
| 116 |
|
|
XMU=0.035 9148.6 |
| 117 |
|
|
T0=1978.0 9148.7 |
| 118 |
|
|
XJ0=11.0 9148.8 |
| 119 |
|
|
219 DT=YEAR-T0 9148.9 |
| 120 |
|
|
XX=XSS+(X0-XSS)*EXP(-DT/TAU) 9149. |
| 121 |
|
|
+ +XJ0*TAU/BNA/(1.0+XMU*TAU)*(EXP(XMU*DT)-EXP(-DT/TAU)) 9149.1 |
| 122 |
|
|
220 XX=XX*1.E-03 9149.2 |
| 123 |
|
|
GO TO 8000 9149.3 |
| 124 |
|
|
C 9149.4 |
| 125 |
|
|
C CH4 SCENARIOS 9149.5 |
| 126 |
|
|
C------------------------------------------------------------------- 9149.6 |
| 127 |
|
|
300 CONTINUE 9149.7 |
| 128 |
|
|
IF(YEAR.GT.1958.0) GO TO 305 9149.8 |
| 129 |
|
|
DT=(YEAR-1850.0)/108.0 9149.9 |
| 130 |
|
|
XX=1.00*(1.40/1.00)**DT 9150. |
| 131 |
|
|
GO TO 320 9150.1 |
| 132 |
|
|
305 IF(YEAR.GT.1970.0) GO TO 310 9150.2 |
| 133 |
|
|
DT=(YEAR-1958.0)/12.0 9150.3 |
| 134 |
|
|
XX=1.40*(1.50/1.40)**DT 9150.4 |
| 135 |
|
|
GO TO 320 9150.5 |
| 136 |
|
|
310 IF(YEAR.GT.1980.0) GO TO 315 9150.6 |
| 137 |
|
|
DT=(YEAR-1970.0)/10.0 9150.7 |
| 138 |
|
|
XX=1.50*(1.65/1.50)**DT 9150.8 |
| 139 |
|
|
GO TO 320 9150.9 |
| 140 |
|
|
315 DT=YEAR-1980.0 9151. |
| 141 |
|
|
XX=1.65*1.015**DT 9151.1 |
| 142 |
|
|
320 CONTINUE 9151.2 |
| 143 |
|
|
GO TO 8000 9151.3 |
| 144 |
|
|
C 9151.4 |
| 145 |
|
|
C F11 SCENARIOS 9151.5 |
| 146 |
|
|
C------------------------------------------------------------------- 9151.6 |
| 147 |
|
|
400 CONTINUE 9151.7 |
| 148 |
|
|
IF(YEAR.LT.1985.0) GO TO 410 9151.8 |
| 149 |
|
|
UAT=1.0/75.0 9151.9 |
| 150 |
|
|
XMU=0.03 9152. |
| 151 |
|
|
XX0=XXF11(45) 9152.1 |
| 152 |
|
|
ETA=XXF11(45)-XXF11(44)*EXP(-UAT) 9152.2 |
| 153 |
|
|
DT=YEAR-1985.0 9152.3 |
| 154 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9152.4 |
| 155 |
|
|
GO TO 420 9152.5 |
| 156 |
|
|
410 XX=0.0 9152.6 |
| 157 |
|
|
IF(YEAR.LT.1945.0) GO TO 420 9152.7 |
| 158 |
|
|
I=YEAR 9152.8 |
| 159 |
|
|
DELTA=YEAR-I 9152.9 |
| 160 |
|
|
I=I-1940 9153. |
| 161 |
|
|
XX=XXF11(I)+DELTA*(XXF11(I+1)-XXF11(I)) 9153.1 |
| 162 |
|
|
420 XX=XX*2.0 9153.2 |
| 163 |
|
|
GO TO 8000 9153.3 |
| 164 |
|
|
C 9153.4 |
| 165 |
|
|
C F12 SCENARIOS 9153.5 |
| 166 |
|
|
C------------------------------------------------------------------- 9153.6 |
| 167 |
|
|
500 CONTINUE 9153.7 |
| 168 |
|
|
IF(YEAR.LT.1985.0) GO TO 510 9153.8 |
| 169 |
|
|
UAT=1.0/150.0 9153.9 |
| 170 |
|
|
XMU=0.03 9154. |
| 171 |
|
|
XX0=XXF12(45) 9154.1 |
| 172 |
|
|
ETA=XXF12(45)-XXF12(44)*EXP(-UAT) 9154.2 |
| 173 |
|
|
DT=YEAR-1985.0 9154.3 |
| 174 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9154.4 |
| 175 |
|
|
GO TO 520 9154.5 |
| 176 |
|
|
510 XX=0.0 9154.6 |
| 177 |
|
|
IF(YEAR.LT.1943.0) GO TO 520 9154.7 |
| 178 |
|
|
I=YEAR 9154.8 |
| 179 |
|
|
DELTA=YEAR-I 9154.9 |
| 180 |
|
|
I=I-1940 9155. |
| 181 |
|
|
XX=XXF12(I)+DELTA*(XXF12(I+1)-XXF12(I)) 9155.1 |
| 182 |
|
|
520 XX=XX*2.0 9155.2 |
| 183 |
|
|
GO TO 8000 9155.3 |
| 184 |
|
|
C 9155.4 |
| 185 |
|
|
8000 XGAS(N)=XX 9155.5 |
| 186 |
|
|
9000 CONTINUE 9155.6 |
| 187 |
|
|
RETURN 9155.7 |
| 188 |
|
|
END 9155.8 |
| 189 |
|
|
SUBROUTINE BTREND(XGAS,YEAR,NGAS) 9155.9 |
| 190 |
|
|
C 9156. |
| 191 |
|
|
C-----------------------------------------------------------------------9156.1 |
| 192 |
|
|
C T-GAS SCENARIO B: REDUCED TREND FOR TRACE GAS INCREASES 9156.2 |
| 193 |
|
|
C (REDUCED FORCING BEYOND 2000) 9156.3 |
| 194 |
|
|
C 9156.4 |
| 195 |
|
|
C CO2: OBSERVED TREND 1958-1984, 9156.5 |
| 196 |
|
|
C 0.015/YR INCREASE IN 1.5 PPM INCREMENT 9156.6 |
| 197 |
|
|
C THRU 1990, 0.010/YR THRU 2000, 0.005/YR 9156.7 |
| 198 |
|
|
C 2010, CONSTANT INCREMENT BEYOND 2010 9156.8 |
| 199 |
|
|
C 9156.9 |
| 200 |
|
|
C N2O: WEISS FORMULA THRU 1978, 0.035/YR INCREASE 9157. |
| 201 |
|
|
C IN 1978 INCREMENT THRU 1990, 0.025/YR 9157.1 |
| 202 |
|
|
C THRU 2000, 0.015/YR THRU 2010, 0.005/YR 9157.2 |
| 203 |
|
|
C INCREASE IN INCREMENT BEYOND 2010 9157.3 |
| 204 |
|
|
C 9157.4 |
| 205 |
|
|
C CH4: OBSERVED RATES THRU 1980, 9157.5 |
| 206 |
|
|
C 0.015/YR INCREASE THRU 1990, 0.010/YR 9157.6 |
| 207 |
|
|
C INCREASE THRU 2000, 0.005/YR BEYOND 2000 9157.7 |
| 208 |
|
|
C 9157.8 |
| 209 |
|
|
C CFC: OBSERVED CONCENTRATION THRU 1984, 9157.9 |
| 210 |
|
|
C 0.03 EMISSION RATE INCREASE TO 1990, 9158. |
| 211 |
|
|
C 0.02 EMISSION RATE INCREASE TO 2000, 9158.1 |
| 212 |
|
|
C 0.01 EMISSION RATE INCREASE TO 2010, 9158.2 |
| 213 |
|
|
C CONSTANT EMISSION RATE BEYOND 2010 9158.3 |
| 214 |
|
|
C 9158.4 |
| 215 |
|
|
C-----------------------------------------------------------------------9158.5 |
| 216 |
|
|
C 9158.6 |
| 217 |
|
|
DIMENSION XGAS(1) 9158.7 |
| 218 |
|
|
DIMENSION XXF11(45),XXF12(45),XXCO2(27) 9158.8 |
| 219 |
|
|
DATA XXF11/ 9158.9 |
| 220 |
|
|
10.00001,0.00001,0.00001,0.00001,0.00001,0.00004,0.00010,0.00023, 9159. |
| 221 |
|
|
20.00041,0.00066,0.00100,0.00148,0.00215,0.00297,0.00397,0.00523, 9159.1 |
| 222 |
|
|
30.00662,0.00791,0.00921,0.01093,0.01316,0.01596,0.01939,0.02345, 9159.2 |
| 223 |
|
|
40.02807,0.03321,0.03904,0.04569,0.05343,0.06224,0.07186,0.08267, 9159.3 |
| 224 |
|
|
50.09505,0.10867,0.12174,0.13422,0.14665,0.15820,0.16871,0.17841, 9159.4 |
| 225 |
|
|
60.18919,0.19983,0.21032,0.22068,0.23090/ 9159.5 |
| 226 |
|
|
DATA XXF12/ 9159.6 |
| 227 |
|
|
10.00001,0.00001,0.00001,0.00040,0.00112,0.00198,0.00298,0.00411, 9159.7 |
| 228 |
|
|
20.00539,0.00679,0.00836,0.00998,0.01181,0.01387,0.01619,0.01888, 9159.8 |
| 229 |
|
|
30.02195,0.02513,0.02868,0.03292,0.03766,0.04312,0.04954,0.05702, 9159.9 |
| 230 |
|
|
40.06546,0.07485,0.08545,0.09756,0.11130,0.12634,0.14250,0.16011, 9160. |
| 231 |
|
|
50.17962,0.20082,0.22144,0.24105,0.25951,0.27630,0.29247,0.30831, 9160.1 |
| 232 |
|
|
60.32585,0.34328,0.36059,0.37778,0.39486/ 9160.2 |
| 233 |
|
|
DATA XXCO2/ 9160.3 |
| 234 |
|
|
1 315.00, 315.69, 316.66, 317.27, 318.24, 318.93, 319.33, 319.95, 9160.4 |
| 235 |
|
|
2 320.80, 321.44, 322.25, 323.44, 324.82, 326.01, 326.97, 328.86, 9160.5 |
| 236 |
|
|
3 330.35, 330.71, 331.69, 332.83, 334.58, 335.95, 337.65, 339.04, 9160.6 |
| 237 |
|
|
4 340.01, 341.56, 343.81/ 9160.7 |
| 238 |
|
|
C 9160.8 |
| 239 |
|
|
DO 9000 N=1,NGAS 9160.9 |
| 240 |
|
|
GO TO (100,200,300,400,500),N 9161. |
| 241 |
|
|
C 9161.1 |
| 242 |
|
|
C CO2 SCENARIOS 9161.2 |
| 243 |
|
|
C------------------------------------------------------------------- 9161.3 |
| 244 |
|
|
100 CONTINUE 9161.4 |
| 245 |
|
|
IF(YEAR.GT.1958.0) GO TO 110 9161.5 |
| 246 |
|
|
DT=(YEAR-1850.0)/(1958.0-1850.0) 9161.6 |
| 247 |
|
|
XX=270.0*(315.0/270.0)**DT 9161.7 |
| 248 |
|
|
GO TO 120 9161.8 |
| 249 |
|
|
110 IF(YEAR.GT.1984.0) GO TO 115 9161.9 |
| 250 |
|
|
I=YEAR 9162. |
| 251 |
|
|
DELTA=YEAR-I 9162.1 |
| 252 |
|
|
I=I-1957 9162.2 |
| 253 |
|
|
XX=XXCO2(I)+DELTA*(XXCO2(I+1)-XXCO2(I)) 9162.3 |
| 254 |
|
|
GO TO 120 9162.4 |
| 255 |
|
|
115 XX0=343.81 9162.5 |
| 256 |
|
|
ETA=1.50 9162.6 |
| 257 |
|
|
XMU=0.015 9162.7 |
| 258 |
|
|
IF(YEAR.GT.1990.0) GO TO 116 9162.8 |
| 259 |
|
|
DT=YEAR-1984.0 9162.9 |
| 260 |
|
|
GO TO 119 9163. |
| 261 |
|
|
116 DT=1990.0-1984.0 9163.1 |
| 262 |
|
|
XX0=XX0+ETA*(EXP(XMU*DT)-1.0)/XMU 9163.2 |
| 263 |
|
|
ETA=ETA*EXP(XMU*DT) 9163.3 |
| 264 |
|
|
XMU=0.010 9163.4 |
| 265 |
|
|
IF(YEAR.GT.2000.0) GO TO 117 9163.5 |
| 266 |
|
|
DT=YEAR-1990.0 9163.6 |
| 267 |
|
|
GO TO 119 9163.7 |
| 268 |
|
|
117 DT=2000.0-1990.0 9163.8 |
| 269 |
|
|
XX0=XX0+ETA*(EXP(XMU*DT)-1.0)/XMU 9163.9 |
| 270 |
|
|
ETA=ETA*EXP(XMU*DT) 9164. |
| 271 |
|
|
XMU=0.005 9164.1 |
| 272 |
|
|
IF(YEAR.GT.2010.0) GO TO 118 9164.2 |
| 273 |
|
|
DT=YEAR-2000.0 9164.3 |
| 274 |
|
|
GO TO 119 9164.4 |
| 275 |
|
|
118 DT=2010.0-2000.0 9164.5 |
| 276 |
|
|
XX0=XX0+ETA*(EXP(XMU*DT)-1.0)/XMU 9164.6 |
| 277 |
|
|
ETA=ETA*EXP(XMU*DT) 9164.7 |
| 278 |
|
|
DT=YEAR-2010.0 9164.8 |
| 279 |
|
|
XX=XX0+ETA*DT 9164.9 |
| 280 |
|
|
GO TO 120 9165. |
| 281 |
|
|
119 XX=XX0+ETA*(EXP(XMU*DT)-1.0)/XMU 9165.1 |
| 282 |
|
|
120 CONTINUE 9165.2 |
| 283 |
|
|
GO TO 8000 9165.3 |
| 284 |
|
|
C 9165.4 |
| 285 |
|
|
C N2O SCENARIOS 9165.5 |
| 286 |
|
|
C------------------------------------------------------------------- 9165.6 |
| 287 |
|
|
200 XSS=282.0 9165.7 |
| 288 |
|
|
UAT=1.0/150.0 9165.8 |
| 289 |
|
|
IF(YEAR.GT.1912.0) GO TO 205 9165.9 |
| 290 |
|
|
XX0=285.355 9166. |
| 291 |
|
|
XMU=0.04 9166.1 |
| 292 |
|
|
T0=1912.0 9166.2 |
| 293 |
|
|
ETA=11.0/16.7*EXP(-0.01*33.0-0.035*33.0) 9166.3 |
| 294 |
|
|
GO TO 219 9166.4 |
| 295 |
|
|
205 IF(YEAR.GT.1945.0) GO TO 210 9166.5 |
| 296 |
|
|
XX0=289.96 9166.6 |
| 297 |
|
|
XMU=0.01 9166.7 |
| 298 |
|
|
T0=1945.0 9166.8 |
| 299 |
|
|
ETA=11.0/16.7*EXP(-0.035*33.0) 9166.9 |
| 300 |
|
|
GO TO 219 9167. |
| 301 |
|
|
210 XX0=300.2 9167.1 |
| 302 |
|
|
XMU=0.035 9167.2 |
| 303 |
|
|
T0=1978.0 9167.3 |
| 304 |
|
|
ETA=11.0/16.7 9167.4 |
| 305 |
|
|
IF(YEAR.GT.1990.0) GO TO 215 9167.5 |
| 306 |
|
|
GO TO 219 9167.6 |
| 307 |
|
|
215 DT=1990.0-T0 9167.7 |
| 308 |
|
|
XX0=XSS+(XX0-XSS)*EXP(-DT*UAT) 9167.8 |
| 309 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9167.9 |
| 310 |
|
|
ETA=ETA*EXP(XMU*DT) 9168. |
| 311 |
|
|
XMU=0.025 9168.1 |
| 312 |
|
|
T0=1990.0 9168.2 |
| 313 |
|
|
IF(YEAR.GT.2000.0) GO TO 216 9168.3 |
| 314 |
|
|
GO TO 219 9168.4 |
| 315 |
|
|
216 DT=2000.0-T0 9168.5 |
| 316 |
|
|
XX0=XSS+(XX0-XSS)*EXP(-DT*UAT) 9168.6 |
| 317 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9168.7 |
| 318 |
|
|
ETA=ETA*EXP(XMU*DT) 9168.8 |
| 319 |
|
|
XMU=0.015 9168.9 |
| 320 |
|
|
T0=2000.0 9169. |
| 321 |
|
|
IF(YEAR.GT.2010.0) GO TO 217 9169.1 |
| 322 |
|
|
GO TO 219 9169.2 |
| 323 |
|
|
217 DT=2010.0-T0 9169.3 |
| 324 |
|
|
XX0=XSS+(XX0-XSS)*EXP(-DT*UAT) 9169.4 |
| 325 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9169.5 |
| 326 |
|
|
ETA=ETA*EXP(XMU*DT) 9169.6 |
| 327 |
|
|
XMU=0.005 9169.7 |
| 328 |
|
|
T0=2010.0 9169.8 |
| 329 |
|
|
219 DT=YEAR-T0 9169.9 |
| 330 |
|
|
XX=XSS+(XX0-XSS)*EXP(-DT*UAT) 9170. |
| 331 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9170.1 |
| 332 |
|
|
220 XX=XX*1.E-03 9170.2 |
| 333 |
|
|
GO TO 8000 9170.3 |
| 334 |
|
|
C 9170.4 |
| 335 |
|
|
C CH4 SCENARIOS 9170.5 |
| 336 |
|
|
C------------------------------------------------------------------- 9170.6 |
| 337 |
|
|
300 CONTINUE 9170.7 |
| 338 |
|
|
IF(YEAR.GT.1958.0) GO TO 305 9170.8 |
| 339 |
|
|
DT=(YEAR-1850.0)/(1958.0-1850.0) 9170.9 |
| 340 |
|
|
XX=1.00*(1.40/1.00)**DT 9171. |
| 341 |
|
|
GO TO 320 9171.1 |
| 342 |
|
|
305 IF(YEAR.GT.1970.0) GO TO 310 9171.2 |
| 343 |
|
|
DT=(YEAR-1958.0)/12.0 9171.3 |
| 344 |
|
|
XX=1.40*(1.50/1.40)**DT 9171.4 |
| 345 |
|
|
GO TO 320 9171.5 |
| 346 |
|
|
310 IF(YEAR.GT.1980.0) GO TO 315 9171.6 |
| 347 |
|
|
DT=(YEAR-1970.0)/10.0 9171.7 |
| 348 |
|
|
XX=1.50*(1.65/1.50)**DT 9171.8 |
| 349 |
|
|
GO TO 320 9171.9 |
| 350 |
|
|
315 IF(YEAR.GT.1990.0) GO TO 316 9172. |
| 351 |
|
|
DT=YEAR-1980.0 9172.1 |
| 352 |
|
|
XX=1.65*1.015**DT 9172.2 |
| 353 |
|
|
GO TO 320 9172.3 |
| 354 |
|
|
316 IF(YEAR.GT.2000.0) GO TO 317 9172.4 |
| 355 |
|
|
DT=YEAR-1990.0 9172.5 |
| 356 |
|
|
XX=1.65*1.015**10.0*1.010**DT 9172.6 |
| 357 |
|
|
GO TO 320 9172.7 |
| 358 |
|
|
317 IF(YEAR.GT.2010.0) GO TO 318 9172.8 |
| 359 |
|
|
DT=YEAR-2000.0 9172.9 |
| 360 |
|
|
XX=1.65*1.015**10.0*1.010**10.0*1.005**DT 9173. |
| 361 |
|
|
GO TO 320 9173.1 |
| 362 |
|
|
318 DT=YEAR-2010.0 9173.2 |
| 363 |
|
|
XX=1.65*1.015**10.0*1.010**10.0*1.005**10.0*1.005**DT 9173.3 |
| 364 |
|
|
320 CONTINUE 9173.4 |
| 365 |
|
|
GO TO 8000 9173.5 |
| 366 |
|
|
C 9173.6 |
| 367 |
|
|
C F11 SCENARIOS 9173.7 |
| 368 |
|
|
C------------------------------------------------------------------- 9173.8 |
| 369 |
|
|
400 XX=0.0 9173.9 |
| 370 |
|
|
IF(YEAR.LT.1945.0) GO TO 420 9174. |
| 371 |
|
|
IF(YEAR.GT.1985.0) GO TO 410 9174.1 |
| 372 |
|
|
I=YEAR 9174.2 |
| 373 |
|
|
DELTA=YEAR-I 9174.3 |
| 374 |
|
|
I=I-1940 9174.4 |
| 375 |
|
|
XX=XXF11(I)+DELTA*(XXF11(I+1)-XXF11(I)) 9174.5 |
| 376 |
|
|
GO TO 420 9174.6 |
| 377 |
|
|
410 UAT=1.0/75.0 9174.7 |
| 378 |
|
|
XMU=0.03 9174.8 |
| 379 |
|
|
XX0=XXF11(45) 9174.9 |
| 380 |
|
|
ETA=XXF11(45)-XXF11(44)*EXP(-UAT) 9175. |
| 381 |
|
|
IF(YEAR.GT.1990.0) GO TO 415 9175.1 |
| 382 |
|
|
DT=YEAR-1985.0 9175.2 |
| 383 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9175.3 |
| 384 |
|
|
GO TO 420 9175.4 |
| 385 |
|
|
415 DT=1990.0-1985.0 9175.5 |
| 386 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9175.6 |
| 387 |
|
|
ETA=ETA*EXP(XMU*DT) 9175.7 |
| 388 |
|
|
XMU=0.02 9175.8 |
| 389 |
|
|
IF(YEAR.GT.2000.0) GO TO 416 9175.9 |
| 390 |
|
|
DT=YEAR-1990.0 9176. |
| 391 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9176.1 |
| 392 |
|
|
GO TO 420 9176.2 |
| 393 |
|
|
416 DT=2000.0-1990.0 9176.3 |
| 394 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9176.4 |
| 395 |
|
|
ETA=ETA*EXP(XMU*DT) 9176.5 |
| 396 |
|
|
XMU=0.01 9176.6 |
| 397 |
|
|
IF(YEAR.GT.2010.0) GO TO 417 9176.7 |
| 398 |
|
|
DT=YEAR-2000.0 9176.8 |
| 399 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9176.9 |
| 400 |
|
|
GO TO 420 9177. |
| 401 |
|
|
417 DT=2010.0-2000.0 9177.1 |
| 402 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9177.2 |
| 403 |
|
|
ETA=ETA*EXP(XMU*DT) 9177.3 |
| 404 |
|
|
XMU=0.0 9177.4 |
| 405 |
|
|
DT=YEAR-2010.0 9177.5 |
| 406 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9177.6 |
| 407 |
|
|
420 XX=XX*1.0 9177.7 |
| 408 |
|
|
GO TO 8000 9177.8 |
| 409 |
|
|
C 9177.9 |
| 410 |
|
|
C F12 SCENARIOS 9178. |
| 411 |
|
|
C------------------------------------------------------------------- 9178.1 |
| 412 |
|
|
500 XX=0.0 9178.2 |
| 413 |
|
|
IF(YEAR.LT.1943.0) GO TO 520 9178.3 |
| 414 |
|
|
IF(YEAR.GT.1985.0) GO TO 510 9178.4 |
| 415 |
|
|
I=YEAR 9178.5 |
| 416 |
|
|
DELTA=YEAR-I 9178.6 |
| 417 |
|
|
I=I-1940 9178.7 |
| 418 |
|
|
XX=XXF12(I)+DELTA*(XXF12(I+1)-XXF12(I)) 9178.8 |
| 419 |
|
|
GO TO 520 9178.9 |
| 420 |
|
|
510 UAT=1.0/150.0 9179. |
| 421 |
|
|
XMU=0.03 9179.1 |
| 422 |
|
|
XX0=XXF12(45) 9179.2 |
| 423 |
|
|
ETA=XXF12(45)-XXF12(44)*EXP(-UAT) 9179.3 |
| 424 |
|
|
IF(YEAR.GT.1990.0) GO TO 515 9179.4 |
| 425 |
|
|
DT=YEAR-1985.0 9179.5 |
| 426 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9179.6 |
| 427 |
|
|
GO TO 520 9179.7 |
| 428 |
|
|
515 DT=1990.0-1985.0 9179.8 |
| 429 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9179.9 |
| 430 |
|
|
ETA=ETA*EXP(XMU*DT) 9180. |
| 431 |
|
|
XMU=0.02 9180.1 |
| 432 |
|
|
IF(YEAR.GT.2000.0) GO TO 516 9180.2 |
| 433 |
|
|
DT=YEAR-1990.0 9180.3 |
| 434 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9180.4 |
| 435 |
|
|
GO TO 520 9180.5 |
| 436 |
|
|
516 DT=2000.0-1990.0 9180.6 |
| 437 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9180.7 |
| 438 |
|
|
ETA=ETA*EXP(XMU*DT) 9180.8 |
| 439 |
|
|
XMU=0.01 9180.9 |
| 440 |
|
|
IF(YEAR.GT.2010.0) GO TO 517 9181. |
| 441 |
|
|
DT=YEAR-2000.0 9181.1 |
| 442 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9181.2 |
| 443 |
|
|
GO TO 520 9181.3 |
| 444 |
|
|
517 DT=2010.0-2000.0 9181.4 |
| 445 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9181.5 |
| 446 |
|
|
ETA=ETA*EXP(XMU*DT) 9181.6 |
| 447 |
|
|
XMU=0.0 9181.7 |
| 448 |
|
|
DT=YEAR-2010.0 9181.8 |
| 449 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9181.9 |
| 450 |
|
|
520 XX=XX*1.0 9182. |
| 451 |
|
|
GO TO 8000 9182.1 |
| 452 |
|
|
C 9182.2 |
| 453 |
|
|
8000 XGAS(N)=XX 9182.3 |
| 454 |
|
|
9000 CONTINUE 9182.4 |
| 455 |
|
|
RETURN 9182.5 |
| 456 |
|
|
END 9182.6 |
| 457 |
|
|
|
| 458 |
|
|
SUBROUTINE CTREND(XGAS,YEAR,NGAS) 9182.7 |
| 459 |
|
|
C 9182.8 |
| 460 |
|
|
C-----------------------------------------------------------------------9182.9 |
| 461 |
|
|
C T-GAS SCENARIO C: CURTAILED TREND FOR TRACE GAS INCREASES 9183. |
| 462 |
|
|
C (EFFECTIVE FORCING CUTOFF BEYOND 2000) 9183.1 |
| 463 |
|
|
C 9183.2 |
| 464 |
|
|
C CO2: OBSERVED TREND 1958-1984, 9183.3 |
| 465 |
|
|
C 1.5PPM FIXED INCREMENTS 1984-2000 9183.4 |
| 466 |
|
|
C FIXED CONCENTRATION BEYOND 2000 9183.5 |
| 467 |
|
|
C 9183.6 |
| 468 |
|
|
C N2O: WEISS FORMULA THRU 1984, 9183.7 |
| 469 |
|
|
C FIXED INCREMENT 1984-2000, 9183.8 |
| 470 |
|
|
C FIXED CONCENTRATION BEYOND 2000 9183.9 |
| 471 |
|
|
C 9184. |
| 472 |
|
|
C CH4: OBSERVED RATES THRU 1980, 9184.1 |
| 473 |
|
|
C 0.01 INCREMENTS TO 1990, 0.005 TO 2000 9184.2 |
| 474 |
|
|
C FIXED CONCENTRATION BEYOND 2000 9184.3 |
| 475 |
|
|
C 9184.4 |
| 476 |
|
|
C CFC: OBSERVED CONCENTRATION THRU 1984, 9184.5 |
| 477 |
|
|
C 0.03 EMISSION RATE INCREASE TO 1990, 9184.6 |
| 478 |
|
|
C EMISSIONS TO ZERO DURING 1990-2000, 9184.7 |
| 479 |
|
|
C SLOW 75YR, 150YR DECAY RATE BEYOND 2000 9184.8 |
| 480 |
|
|
C 9184.9 |
| 481 |
|
|
C-----------------------------------------------------------------------9185. |
| 482 |
|
|
C 9185.1 |
| 483 |
|
|
DIMENSION XGAS(1) 9185.2 |
| 484 |
|
|
DIMENSION XXF11(45),XXF12(45),XXCO2(27) 9185.3 |
| 485 |
|
|
DATA XXF11/ 9185.4 |
| 486 |
|
|
10.00001,0.00001,0.00001,0.00001,0.00001,0.00004,0.00010,0.00023, 9185.5 |
| 487 |
|
|
20.00041,0.00066,0.00100,0.00148,0.00215,0.00297,0.00397,0.00523, 9185.6 |
| 488 |
|
|
30.00662,0.00791,0.00921,0.01093,0.01316,0.01596,0.01939,0.02345, 9185.7 |
| 489 |
|
|
40.02807,0.03321,0.03904,0.04569,0.05343,0.06224,0.07186,0.08267, 9185.8 |
| 490 |
|
|
50.09505,0.10867,0.12174,0.13422,0.14665,0.15820,0.16871,0.17841, 9185.9 |
| 491 |
|
|
60.18919,0.19983,0.21032,0.22068,0.23090/ 9186. |
| 492 |
|
|
DATA XXF12/ 9186.1 |
| 493 |
|
|
10.00001,0.00001,0.00001,0.00040,0.00112,0.00198,0.00298,0.00411, 9186.2 |
| 494 |
|
|
20.00539,0.00679,0.00836,0.00998,0.01181,0.01387,0.01619,0.01888, 9186.3 |
| 495 |
|
|
30.02195,0.02513,0.02868,0.03292,0.03766,0.04312,0.04954,0.05702, 9186.4 |
| 496 |
|
|
40.06546,0.07485,0.08545,0.09756,0.11130,0.12634,0.14250,0.16011, 9186.5 |
| 497 |
|
|
50.17962,0.20082,0.22144,0.24105,0.25951,0.27630,0.29247,0.30831, 9186.6 |
| 498 |
|
|
60.32585,0.34328,0.36059,0.37778,0.39486/ 9186.7 |
| 499 |
|
|
DATA XXCO2/ 9186.8 |
| 500 |
|
|
1 315.00, 315.69, 316.66, 317.27, 318.24, 318.93, 319.33, 319.95, 9186.9 |
| 501 |
|
|
2 320.80, 321.44, 322.25, 323.44, 324.82, 326.01, 326.97, 328.86, 9187. |
| 502 |
|
|
3 330.35, 330.71, 331.69, 332.83, 334.58, 335.95, 337.65, 339.04, 9187.1 |
| 503 |
|
|
4 340.01, 341.56, 343.81/ 9187.2 |
| 504 |
|
|
C 9187.3 |
| 505 |
|
|
DO 9000 N=1,NGAS 9187.4 |
| 506 |
|
|
GO TO (100,200,300,400,500),N 9187.5 |
| 507 |
|
|
C 9187.6 |
| 508 |
|
|
C CO2 SCENARIOS 9187.7 |
| 509 |
|
|
C------------------------------------------------------------------- 9187.8 |
| 510 |
|
|
100 CONTINUE 9187.9 |
| 511 |
|
|
IF(YEAR.GT.1958.0) GO TO 110 9188. |
| 512 |
|
|
DT=(YEAR-1850.0)/(1958.0-1850.0) 9188.1 |
| 513 |
|
|
XX=270.0*(315.0/270.0)**DT 9188.2 |
| 514 |
|
|
GO TO 120 9188.3 |
| 515 |
|
|
110 IF(YEAR.GT.1984.0) GO TO 115 9188.4 |
| 516 |
|
|
I=YEAR 9188.5 |
| 517 |
|
|
DELTA=YEAR-I 9188.6 |
| 518 |
|
|
I=I-1957 9188.7 |
| 519 |
|
|
XX=XXCO2(I)+DELTA*(XXCO2(I+1)-XXCO2(I)) 9188.8 |
| 520 |
|
|
GO TO 120 9188.9 |
| 521 |
|
|
115 XX0=343.81 9189. |
| 522 |
|
|
ETA=1.50 9189.1 |
| 523 |
|
|
XMU=0.0 9189.2 |
| 524 |
|
|
IF(YEAR.GT.2000.0) GO TO 116 9189.3 |
| 525 |
|
|
DT=YEAR-1984.0 9189.4 |
| 526 |
|
|
GO TO 119 9189.5 |
| 527 |
|
|
116 DT=2000.0-1984.0 9189.6 |
| 528 |
|
|
XX0=XX0+ETA*DT 9189.7 |
| 529 |
|
|
DT=0.0 9189.8 |
| 530 |
|
|
119 XX=XX0+ETA*DT 9189.9 |
| 531 |
|
|
120 CONTINUE 9190. |
| 532 |
|
|
GO TO 8000 9190.1 |
| 533 |
|
|
C 9190.2 |
| 534 |
|
|
C N2O SCENARIOS 9190.3 |
| 535 |
|
|
C------------------------------------------------------------------- 9190.4 |
| 536 |
|
|
200 XSS=282.0 9190.5 |
| 537 |
|
|
UAT=1.0/150.0 9190.6 |
| 538 |
|
|
IF(YEAR.GT.1912.0) GO TO 205 9190.7 |
| 539 |
|
|
XX0=285.355 9190.8 |
| 540 |
|
|
XMU=0.04 9190.9 |
| 541 |
|
|
T0=1912.0 9191. |
| 542 |
|
|
ETA=11.0/16.7*EXP(-0.01*33.0-0.035*33.0) 9191.1 |
| 543 |
|
|
GO TO 219 9191.2 |
| 544 |
|
|
205 IF(YEAR.GT.1945.0) GO TO 210 9191.3 |
| 545 |
|
|
XX0=289.96 9191.4 |
| 546 |
|
|
XMU=0.01 9191.5 |
| 547 |
|
|
T0=1945.0 9191.6 |
| 548 |
|
|
ETA=11.0/16.7*EXP(-0.035*33.0) 9191.7 |
| 549 |
|
|
GO TO 219 9191.8 |
| 550 |
|
|
210 XX0=300.2 9191.9 |
| 551 |
|
|
XMU=0.035 9192. |
| 552 |
|
|
T0=1978.0 9192.1 |
| 553 |
|
|
ETA=11.0/16.7 9192.2 |
| 554 |
|
|
IF(YEAR.GT.1984.0) GO TO 215 9192.3 |
| 555 |
|
|
GO TO 219 9192.4 |
| 556 |
|
|
215 DT=1984.0-T0 9192.5 |
| 557 |
|
|
XX0=XSS+(XX0-XSS)*EXP(-DT*UAT) 9192.6 |
| 558 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9192.7 |
| 559 |
|
|
ETA=ETA*EXP(XMU*DT) 9192.8 |
| 560 |
|
|
XMU=0.0 9192.9 |
| 561 |
|
|
T0=1984.0 9193. |
| 562 |
|
|
IF(YEAR.GT.2000.0) GO TO 216 9193.1 |
| 563 |
|
|
GO TO 219 9193.2 |
| 564 |
|
|
216 DT=2000.0-T0 9193.3 |
| 565 |
|
|
XX0=XSS+(XX0-XSS)*EXP(-DT*UAT) 9193.4 |
| 566 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9193.5 |
| 567 |
|
|
ETA=ETA*EXP(XMU*DT) 9193.6 |
| 568 |
|
|
T0=YEAR 9193.7 |
| 569 |
|
|
219 DT=YEAR-T0 9193.8 |
| 570 |
|
|
XX=XSS+(XX0-XSS)*EXP(-DT*UAT) 9193.9 |
| 571 |
|
|
+ +ETA/(XMU+UAT)*(EXP(XMU*DT)-EXP(-DT*UAT)) 9194. |
| 572 |
|
|
220 XX=XX*1.E-03 9194.1 |
| 573 |
|
|
GO TO 8000 9194.2 |
| 574 |
|
|
C 9194.3 |
| 575 |
|
|
C CH4 SCENARIOS 9194.4 |
| 576 |
|
|
C------------------------------------------------------------------- 9194.5 |
| 577 |
|
|
300 CONTINUE 9194.6 |
| 578 |
|
|
IF(YEAR.GT.1958.0) GO TO 305 9194.7 |
| 579 |
|
|
DT=(YEAR-1850.0)/(1958.0-1850.0) 9194.8 |
| 580 |
|
|
XX=1.00*(1.40/1.00)**DT 9194.9 |
| 581 |
|
|
GO TO 320 9195. |
| 582 |
|
|
305 IF(YEAR.GT.1970.0) GO TO 310 9195.1 |
| 583 |
|
|
DT=(YEAR-1958.0)/12.0 9195.2 |
| 584 |
|
|
XX=1.40*(1.50/1.40)**DT 9195.3 |
| 585 |
|
|
GO TO 320 9195.4 |
| 586 |
|
|
310 IF(YEAR.GT.1980.0) GO TO 315 9195.5 |
| 587 |
|
|
DT=(YEAR-1970.0)/10.0 9195.6 |
| 588 |
|
|
XX=1.50*(1.65/1.50)**DT 9195.7 |
| 589 |
|
|
GO TO 320 9195.8 |
| 590 |
|
|
315 IF(YEAR.GT.1990.0) GO TO 316 9195.9 |
| 591 |
|
|
DT=YEAR-1980.0 9196. |
| 592 |
|
|
XX=1.65*1.010**DT 9196.1 |
| 593 |
|
|
GO TO 320 9196.2 |
| 594 |
|
|
316 IF(YEAR.GT.2000.0) GO TO 317 9196.3 |
| 595 |
|
|
DT=YEAR-1990.0 9196.4 |
| 596 |
|
|
XX=1.65*1.010**10.0*1.005**DT 9196.5 |
| 597 |
|
|
GO TO 320 9196.6 |
| 598 |
|
|
317 DT=YEAR-2000.0 9196.7 |
| 599 |
|
|
XX=1.65*1.010**10.0*1.005**10.0 9196.8 |
| 600 |
|
|
320 CONTINUE 9196.9 |
| 601 |
|
|
GO TO 8000 9197. |
| 602 |
|
|
C 9197.1 |
| 603 |
|
|
C F11 SCENARIOS 9197.2 |
| 604 |
|
|
C------------------------------------------------------------------- 9197.3 |
| 605 |
|
|
400 XX=0.0 9197.4 |
| 606 |
|
|
IF(YEAR.LT.1945.0) GO TO 420 9197.5 |
| 607 |
|
|
IF(YEAR.GE.1985.0) GO TO 410 9197.6 |
| 608 |
|
|
I=YEAR 9197.7 |
| 609 |
|
|
DELTA=YEAR-I 9197.8 |
| 610 |
|
|
I=I-1940 9197.9 |
| 611 |
|
|
XX=XXF11(I)+DELTA*(XXF11(I+1)-XXF11(I)) 9198. |
| 612 |
|
|
GO TO 420 9198.1 |
| 613 |
|
|
410 UAT=1.0/75.0 9198.2 |
| 614 |
|
|
TAU=75.0 9198.3 |
| 615 |
|
|
XMU=0.03 9198.4 |
| 616 |
|
|
XX0=XXF11(45) 9198.5 |
| 617 |
|
|
ETA=XXF11(45)-XXF11(44)*EXP(-UAT) 9198.6 |
| 618 |
|
|
IF(YEAR.GT.1990.0) GO TO 415 9198.7 |
| 619 |
|
|
DT=YEAR-1985.0 9198.8 |
| 620 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9198.9 |
| 621 |
|
|
GO TO 420 9199. |
| 622 |
|
|
415 DT=1990.0-1985.0 9199.1 |
| 623 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9199.2 |
| 624 |
|
|
ETA=ETA*EXP(XMU*DT) 9199.3 |
| 625 |
|
|
IF(YEAR.GT.2000.0) GO TO 416 9199.4 |
| 626 |
|
|
DT=YEAR-1990.0 9199.5 |
| 627 |
|
|
TAUETA=TAU*ETA+TAU*TAU*ETA/10.0 9199.6 |
| 628 |
|
|
XX=(XX0-TAUETA)*EXP(-UAT*DT)+TAUETA-TAU*ETA/10.0*DT 9199.7 |
| 629 |
|
|
GO TO 420 9199.8 |
| 630 |
|
|
416 DT=2000.0-1990.0 9199.9 |
| 631 |
|
|
TAUETA=TAU*ETA+TAU*TAU*ETA/10.0 9200. |
| 632 |
|
|
XX0=(XX0-TAUETA)*EXP(-UAT*DT)+TAUETA-TAU*ETA/10.0*DT 9200.1 |
| 633 |
|
|
DT=YEAR-2000.0 9200.2 |
| 634 |
|
|
XX=XX0*EXP(-UAT*DT) 9200.3 |
| 635 |
|
|
420 XX=XX*1.0 9200.4 |
| 636 |
|
|
GO TO 8000 9200.5 |
| 637 |
|
|
C 9200.6 |
| 638 |
|
|
C F12 SCENARIOS 9200.7 |
| 639 |
|
|
C------------------------------------------------------------------- 9200.8 |
| 640 |
|
|
500 XX=0.0 9200.9 |
| 641 |
|
|
IF(YEAR.LT.1943.0) GO TO 520 9201. |
| 642 |
|
|
IF(YEAR.GE.1985.0) GO TO 510 9201.1 |
| 643 |
|
|
I=YEAR 9201.2 |
| 644 |
|
|
DELTA=YEAR-I 9201.3 |
| 645 |
|
|
I=I-1940 9201.4 |
| 646 |
|
|
XX=XXF12(I)+DELTA*(XXF12(I+1)-XXF12(I)) 9201.5 |
| 647 |
|
|
GO TO 520 9201.6 |
| 648 |
|
|
510 UAT=1.0/150.0 9201.7 |
| 649 |
|
|
TAU=150.0 9201.8 |
| 650 |
|
|
XMU=0.03 9201.9 |
| 651 |
|
|
XX0=XXF12(45) 9202. |
| 652 |
|
|
ETA=XXF12(45)-XXF12(44)*EXP(-UAT) 9202.1 |
| 653 |
|
|
IF(YEAR.GT.1990.0) GO TO 515 9202.2 |
| 654 |
|
|
DT=YEAR-1985.0 9202.3 |
| 655 |
|
|
XX=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9202.4 |
| 656 |
|
|
GO TO 520 9202.5 |
| 657 |
|
|
515 DT=1990.0-1985.0 9202.6 |
| 658 |
|
|
XX0=XX0*EXP(-UAT*DT)+ETA/(UAT+XMU)*(EXP(XMU*DT)-EXP(-UAT*DT)) 9202.7 |
| 659 |
|
|
ETA=ETA*EXP(XMU*DT) 9202.8 |
| 660 |
|
|
IF(YEAR.GT.2000.0) GO TO 516 9202.9 |
| 661 |
|
|
DT=YEAR-1990.0 9203. |
| 662 |
|
|
TAUETA=TAU*ETA+TAU*TAU*ETA/10.0 9203.1 |
| 663 |
|
|
XX=(XX0-TAUETA)*EXP(-UAT*DT)+TAUETA-TAU*ETA/10.0*DT 9203.2 |
| 664 |
|
|
GO TO 520 9203.3 |
| 665 |
|
|
516 DT=2000.0-1990.0 9203.4 |
| 666 |
|
|
TAUETA=TAU*ETA+TAU*TAU*ETA/10.0 9203.5 |
| 667 |
|
|
XX0=(XX0-TAUETA)*EXP(-UAT*DT)+TAUETA-TAU*ETA/10.0*DT 9203.6 |
| 668 |
|
|
DT=YEAR-2000.0 9203.7 |
| 669 |
|
|
XX=XX0*EXP(-UAT*DT) 9203.8 |
| 670 |
|
|
520 XX=XX*1.0 9203.9 |
| 671 |
|
|
GO TO 8000 9204. |
| 672 |
|
|
C 9204.1 |
| 673 |
|
|
8000 XGAS(N)=XX 9204.2 |
| 674 |
|
|
9000 CONTINUE 9204.3 |
| 675 |
|
|
RETURN 9204.4 |
| 676 |
|
|
END 9204.5 |
Parent Directory
| 
Revision Log
| 
Revision Graph