| 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 |
Parent Directory
| 
Revision Log
| 
Revision Graph