1 |
jscott |
1.1 |
C $Header$ |
2 |
|
|
C $Name$ |
3 |
|
|
|
4 |
|
|
#include "ctrparam.h" |
5 |
|
|
|
6 |
|
|
subroutine eppaemission (iyr) |
7 |
|
|
! =================================== |
8 |
|
|
|
9 |
|
|
#include "chem_para" |
10 |
|
|
#include "chem_com" |
11 |
|
|
#include "BD2G04.COM" |
12 |
|
|
#include "chem_meta" |
13 |
|
|
parameter (neppa=124) |
14 |
|
|
character * 120 emiss_data,SO2ERATIO |
15 |
|
|
namelist /EPPA/ LYEAREM,emiss_data,SO2ERATIO |
16 |
|
|
common /SO2EMIN/SO2EM |
17 |
|
|
data ifirst / 0 / |
18 |
|
|
common /EPPAEMISS/ edailyf11eppa (nlon,nlat,neppa), |
19 |
|
|
& edailyf12eppa (nlon,nlat,neppa), |
20 |
|
|
& edailyn2oeppa (nlon,nlat,neppa), |
21 |
|
|
& edailycoeppa (nlon,nlat,neppa), |
22 |
|
|
& edailynoxeppa (nlon,nlat,neppa), |
23 |
|
|
& edailych4eppa (nlon,nlat,neppa), |
24 |
|
|
& edailyso2eppa (nlon,nlat,neppa), |
25 |
|
|
& edailyco2eppa (nlon,nlat,neppa), |
26 |
|
|
& edailyhfc134aeppa (nlon,nlat,neppa), |
27 |
|
|
& edailypfceppa (nlon,nlat,neppa), |
28 |
|
|
& edailysf6eppa (nlon,nlat,neppa), |
29 |
|
|
& edailyucoeppa (nlon,nlat,neppa), |
30 |
|
|
& edailyunmveppa (nlon,nlat,neppa), |
31 |
|
|
& edailyunoxeppa (nlon,nlat,neppa), |
32 |
|
|
& edailyusoxeppa (nlon,nlat,neppa), |
33 |
|
|
& edailybceppa (nlon,nlat,neppa), |
34 |
|
|
& edailynh3eppa (nlon,nlat,neppa), |
35 |
|
|
& edailyoceppa (nlon,nlat,neppa), |
36 |
|
|
& edailyubceppa (nlon,nlat,neppa), |
37 |
|
|
& edailyunh3eppa (nlon,nlat,neppa), |
38 |
|
|
& edailyuoceppa (nlon,nlat,neppa), |
39 |
|
|
& n_total_urbaneppa (nlat,neppa), |
40 |
|
|
& n_urbaneppa (3,nlat,neppa) |
41 |
|
|
|
42 |
|
|
if (ifirst.eq.0) then |
43 |
|
|
open(535,file='eppaemis.dat') |
44 |
|
|
read (535,NML=EPPA) |
45 |
|
|
close (535) |
46 |
|
|
open(164,file=emiss_data, |
47 |
|
|
& form='unformatted', |
48 |
|
|
& status='old') |
49 |
|
|
open(664,file=SO2ERATIO, |
50 |
|
|
& form='formatted', |
51 |
|
|
& status='old') |
52 |
|
|
|
53 |
|
|
read(664,'(f10.6)')SO2EM |
54 |
|
|
print *,' SO2EM=',SO2EM |
55 |
|
|
|
56 |
|
|
read(164)edailyf11eppa, |
57 |
|
|
& edailyf12eppa, |
58 |
|
|
& edailyn2oeppa, |
59 |
|
|
& edailycoeppa, |
60 |
|
|
& edailynoxeppa, |
61 |
|
|
& edailych4eppa, |
62 |
|
|
& edailyso2eppa, |
63 |
|
|
& edailyco2eppa, |
64 |
|
|
& edailyhfc134aeppa, |
65 |
|
|
& edailypfceppa, |
66 |
|
|
& edailysf6eppa, |
67 |
|
|
& edailyucoeppa, |
68 |
|
|
& edailyunmveppa, |
69 |
|
|
& edailyunoxeppa, |
70 |
|
|
& edailyusoxeppa, |
71 |
|
|
& edailybceppa, |
72 |
|
|
& edailynh3eppa, |
73 |
|
|
& edailyoceppa, |
74 |
|
|
& edailyubceppa, |
75 |
|
|
& edailyunh3eppa, |
76 |
|
|
& edailyuoceppa, |
77 |
|
|
& n_total_urbaneppa, |
78 |
|
|
& n_urbaneppa |
79 |
|
|
|
80 |
|
|
close (164) |
81 |
|
|
close (664) |
82 |
|
|
print *,'Emissioms for ',neppa,' years' |
83 |
|
|
ifirst = 1 |
84 |
|
|
return |
85 |
|
|
endif ! first |
86 |
|
|
iyremi = iyr-1976 !year from starting point |
87 |
|
|
myyearlast = min(LYEAREM-1976,neppa) !last year of emission |
88 |
|
|
print *,' EPPAEMISION myyearlast=',myyearlast |
89 |
|
|
iyremi = min(iyremi,myyearlast) |
90 |
|
|
print *,' EPPAEMISION iyr=',iyr,' iyremi=',iyremi |
91 |
|
|
print *,' EPPAEMISION nlat=',nlat |
92 |
|
|
print *,'CO2 emissions for ',iyremi |
93 |
|
|
print *,(edailyco2eppa(1,j,iyremi),j=1,nlat) |
94 |
|
|
do j=1,nlat |
95 |
|
|
edailyf11(1,j,1)=edailyf11eppa(1,j,iyremi) |
96 |
|
|
edailyf12(1,j,1)=edailyf12eppa(1,j,iyremi) |
97 |
|
|
edailyn2o(1,j,1)=edailyn2oeppa(1,j,iyremi) |
98 |
|
|
edailyco(1,j,1)=edailycoeppa(1,j,iyremi) |
99 |
|
|
edailynox(1,j,1)=edailynoxeppa(1,j,iyremi) |
100 |
|
|
edailych4(1,j,1)=edailych4eppa(1,j,iyremi) |
101 |
|
|
edailyso2(1,j,1)=edailyso2eppa(1,j,iyremi) |
102 |
|
|
edailyco2(1,j,1)=edailyco2eppa(1,j,iyremi) |
103 |
|
|
edailyhfc134a(1,j,1)=edailyhfc134aeppa(1,j,iyremi) |
104 |
|
|
edailypfc(1,j,1)=edailypfceppa(1,j,iyremi) |
105 |
|
|
edailysf6(1,j,1)=edailysf6eppa(1,j,iyremi) |
106 |
|
|
edailyuco(1,j,1)=edailyucoeppa(1,j,iyremi) |
107 |
|
|
edailyunmv(1,j,1)=edailyunmveppa(1,j,iyremi) |
108 |
|
|
edailyunox(1,j,1)=edailyunoxeppa(1,j,iyremi) |
109 |
|
|
edailyusox(1,j,1)=edailyusoxeppa(1,j,iyremi) |
110 |
|
|
edailybc(1,j,1)=edailybceppa(1,j,iyremi) |
111 |
|
|
edailynh3(1,j,1)=edailynh3eppa(1,j,iyremi) |
112 |
|
|
edailyoc(1,j,1)=edailyoceppa(1,j,iyremi) |
113 |
|
|
edailyubc(1,j,1)=edailyubceppa(1,j,iyremi) |
114 |
|
|
edailyunh3(1,j,1)=edailyunh3eppa(1,j,iyremi) |
115 |
|
|
edailyuoc(1,j,1)=edailyuoceppa(1,j,iyremi) |
116 |
|
|
n_total_urban (j,1)=n_total_urbaneppa(j,iyremi) |
117 |
|
|
do k=1,3 |
118 |
|
|
n_urban (k,j,1)= n_urbaneppa(k,j,iyremi) |
119 |
|
|
enddo |
120 |
|
|
enddo |
121 |
|
|
|
122 |
|
|
|
123 |
|
|
return |
124 |
|
|
end |