| 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