From a3b6dc3c2b84db591ca0a88e4c9bc35ef7374ac2 Mon Sep 17 00:00:00 2001 From: chen-0126 <252941087@qq.com> Date: Thu, 23 Feb 2023 09:37:30 +0800 Subject: [PATCH 1/4] fix a bug in the read data file --- .../run_meic_cams_htap/meic/emission$py.class | Bin 8272 -> 0 bytes .../run_meic_cams_htap/meic/emission.py | 56 ------------ .../meic/emission_meic_2017.py | 4 +- .../meic/grid_spec$py.class | Bin 9112 -> 0 bytes .../run_meic_cams_htap/meic/grid_spec.py | 80 ------------------ .../run_meic_cams_htap/meic/run_VOC.py | 2 +- .../run_meic_cams_htap/meic/run_pollutants.py | 8 +- 7 files changed, 7 insertions(+), 143 deletions(-) delete mode 100644 emips/run/run_with_script/run_meic_cams_htap/meic/emission$py.class delete mode 100644 emips/run/run_with_script/run_meic_cams_htap/meic/emission.py delete mode 100644 emips/run/run_with_script/run_meic_cams_htap/meic/grid_spec$py.class delete mode 100644 emips/run/run_with_script/run_meic_cams_htap/meic/grid_spec.py diff --git a/emips/run/run_with_script/run_meic_cams_htap/meic/emission$py.class b/emips/run/run_with_script/run_meic_cams_htap/meic/emission$py.class deleted file mode 100644 index b13e0550ee0712bcb491eb3755a8d34b40283659..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 8272 zcmds6cX(UX6+g$4t!G6EFV7h`dpO0(v!cW#5bQt_Vh04VVl3B?5O{sDpB)A1Jt4^u z)6zxh-aB11Z7JP_0E$aW_oCzL`ucV6N%vl~(BHXFQWQxx`P%>LFIVsC-E+@9<9EhY zKKtkcA0nb|&M^hXJ2o*5b?)_Y6S3@6;jrh%M!lR98=TscvolU_-NNwUoMVqA3wFWj zU0=64IJNK45off}%d~g~mLUo-h4M}zopK$fU}t>&NQ9cHB}_pIG0kx@sce4ZWFeK# zYt+i5AF*$;V`IF;F|xIlh5L4?7tvpLIqxTFNwc&`YN$3$|PEwO5dkF{1ElQZR6HN+B;;Y!E?!Mj@Uz&wL;a) zffd-x0yy&4xZBYU`SHDk%jJm1*_d2Od&OP@h+9!vqkRZi%K`(tP{>(Ugsu{ccGD13 zbK1LE0&`ed2(V9E7Sl!bx#lcWN`yw}YTrMO#ga;^ieXA5edLvMs!^1_s%EOZ|p+LcHS?i(7`2-$vK zA!oaJ$$Ur&o~uy-QcY&F;-#A~9USPkwrX?~-o+HNTktY5FtTsAMo)o4YVuU3g&8L` z8tvY+`9f=5%37C?uFG$9kxic7FzXLbX46iHo(ZFGPr0cAVqrz+Ea(^aJ&SG+({1!@ zNFuQoq38NWejZ*RpB>id1(+T2Qc&VWu-I5Crw}a#L8F&2W$L|I4YQfSGxN%Qw6IdI~AhW0t5|?+U<=w@Zc5nx>hVLKF%#zK#@;DxO zt9ayXOv(RN+Q;GlcHw^){CtZ>?^MwbtECDk%gBT8_5u4IfK%nid^$CX5;40+XI6FQ zm05*+AH83Myc{=QkU&ZDMR!E@w~WL-Y|yh@Jr{ zv=a&DbW>M?c>NpU_%Zr;m_ACMh|pX7sQVQ1W(DmUeFhCmB^oBsIdyD7LCfa6BPn;> z)97<8JFmL!eEGv`7L#By(>sQC*5J3>Deko_b=pY<*Hv+5p0Os3H<{G5W^C=D=rm1Ui?^86+JO6vEoBXp6^^SAIy z#=gPXjpT~$c7slD7L^)+N~PLaW%xb)L45Q_rd54C2d{Ld`p};qEO&_q2m0f?4@wa{ zD7Ej;66mU_66{Dt{cE)q&}nUr(BFMs{&7Y=y#T?PR~MI_RZqTEqkqAr<6bUfqY-`- z0hq;|r;v7;me+f<5>G%f5)ngQZi>(&>KqQ%LHEK4xA<$gX@$i#mJNCH8n>dXH?ZmY z5vG+hx1huJA=Gxp+30qpZu7h_&t=@bf_W!B4h!rRn@Cwk+3St)d^S+oxDBzg-EX2V z(HM855F$rB!gNPxya52KVz&gz`Zw_d2O9Xf+$}3ndmHg?a9q4$e)Ov|FJ@fDc5r)y z*tg|U4SG?F??jDP0IJ#Ybxh+`FvgT)=QLh}0L^$V+86Es1bljRchL9=XwNE}O7A*x zI~#X3-T>E*Yup7!az)a3BNp7F8sqNXz2lrOF9_{vA*kMwBlNfTw(};1zL4};DyE@=hVry{I zI8%}cU&(vJd2ukOdT>=7&8~!H5eapa*{{jjOLr{^rREO zk68FwVQEA9f#DQ>y=o@n+aKBX)CbV=G$bo~3QQ05^~46I26Ns73>H((!+=w8yf~_S zpL0{J;Cbmo-qBe9KzuM6LsrY?W5A^hOYTF~sGZ4MhYNN#CSQd@<-LvBDF|5$pQ=(E z;u{c`;%21lz4%cqF_}5!6v7Y+49>!Y41;OdD&rJg7x@|2>1C-?~* zz%@hEbSbsqXFK9Y&0vu*eoevus`#ENxyh(5l1#o&E z+okpyf0?Em)iyb;C)geg7%xpmW!FzR0k(fD$*!B&d4hH)PEr4U^FE5>kbm+?x&t)F z%QUe%8g4GpJXEIHQct6mXdW)pKt4$MH%MtbT%=@Ik@mN7xHL- zE|^_>c2h*=7)=-wok0vB1UMm}kf;;G#E4;R!Po|BBo!Y|^hVIhj$2Km*EBAM1aO_H zod-cYN2dp*skYn5MHPxJE{<%{&2)9TDqD$IdE?HFqzDzct0D*u%<5tWaddIW4G z020O^jAT5qVkr%FItb+HbITNPmtG{xut5{GisW#jNU39@BoL$-HvkVns=x0py2RAa zsok2d-wKwRF?r_JfT_=yt%$9XTsK$b(p$=c8D$a(wn7A(s^PxW)Z2vnlNI;!#=B^T zsV}JAdu#pPg|hc)HG3i2B6exYkvvJyDbn*{ z!%$h~#d5CEA=17OCmJsnz082<`BF?liD{8ukuU|HQE_Jzd-$9~!GFh6#0|TgcB@En1A(j3vHuWV+sN)c-OF==|-4qR! zUPZdI;#N~%s)#o!V$rrp?|_08r!13YrM>Ykn5anaW#%H?Q*DIhpqU?7HBI;9Z3Sm= zp*ElY?7wX-8*HWGCsr{|DqHDDT(V6a4LhB-?WzhXQh9d*!Hs{&qxbg&PtgbWUw=(| z@Faa$Wh{RUpEXNCD=AKm!f`6=LHz0yR##4_@{gZA)*z&amS8N!Scb72V+F=aj8zz` zVSY3jrEKvaj4!$6PE!XOfQdC?Q9a~mWTY+2_`K+g{G}vbKwp@quS6|KxPZQ@rp;E; zw_?<4k+)w{Z$qaEG!Q6*MlBr_II_SQc0CKkR*w3v#CQ|Kl@q+;8a3(4W<;F8nW#ITeMsywy6@^0aFLs zG~?fZlkr%ED{wc>h{nyQk9?@Icos*7aS+iaksAl1zbhTP6*JLhxv6hG&AMDm>iq5- zWGf`cy~Ah0_aG!uKlsiKCt2TcFE0oLVEKVARG&qGK;Q&kFbC%AG|fxy35p0sm*~lUY_gEW|u2o}OH>NAYtfm1f)(b~r63;T!kcE!G>iv6d9(hUrB`QnLf-+DE%HTLCBZ8Z1msYjwpO3riAw8yc=~V)8bH%(5>s$Qg zDp^>&2fmgT0GGuQ0(3Z2P#!$N%OnA#`11vPVQFjO2FfPZ~TjZg;i#&`mf^iMT0gQtf*JC{CZkF>{)sJ6O>tkqIf+6(P2-gs&DS){@ qeGEGNKdqoF4RL!KHjh%~y5;1!^6wPJg_Oy?#W2za*#3!S` z839XHMx!-Mwabr@ah8wZzD?7p*ns)BruTEb0tR#+>@9-ws; z^pj4tOqvLqnusNsLKV(jnTRN518o$gY+}03QOayCp2%uc2Wg_vaw^kr&n7jx9%4HP z*>o%yk0s1l0*7j}4Gf5m`m8M8wlh_wuzd$pwWEw7vB1&S##d4Bm01qZE@}|Dcf-Ii z;#?wrMx#A=nT+LPP=rQJcmc~+CJufK@MEVnYGn#8x7zTWk3vj={$x5=@Tn}*^{Z@#tSxuk=baOQ<&8^AP0aRG&s&8y|4A$IKKRRO-rVE%ED>s|51qTFx z;9xGHrqMyB&E--RgW*=FY6_kmwsTg-wkDa{Dp<8rygC%R19X^r{B($JgUdkFg3oS4 zbgrhKKI#__Jpz?tlcUque)dAV9Hp?3xgBSRBYh{Lqa))Q4OW!c+2f-Kg2FU$Xe^h@ zm}Y>6g~kxQ9!br#-o&)OQc^3NxT;K18Wln()B$$NzSqj6OcP`qn@$DjB;6sa<9LtT zaN;}R(Nkt3md(YPI2@oi&>O{Nu;IF%fst_UiNVnmBmEl1up>5&BsIb0?T_@23_wPN zpDUfZN26(=kMu{wef^QqaL=GdpytY0h}OHAHkKQHFrLl%2$`&UbfhN|9UdAP)hLBM zJ>e0JAigJhG&~&D2p;Y|7`Nj&xX2BStHh5aoLM^OC!6j89W!nRXwFf^y?D=MVjA6# zd>u*ImPT)e!%ZWzs1)6hMzlTzb5(Ga-qfT8@9L(vRMQ1|E23r20VLaFvqy)#&Xo9Oi6Pxd1)pIN&>wNJp$(Cbhh%@ZpA66wj!S+DhA1(TMiWwH%>;VQAo$*Fh(p8p)ed`hFwBJ{8r zkDWsSZh5c~HW&4>?c^ete_BGX|04`Af!P_iS0@gHIwOQ$UxUh>CC8>T60vaRzuN!x$o z=;BWUblI`TpO*%Ws9&SMfSA&CI>`q0+hy7a+P`wL(_gD>l?~AMocF&~kp&Cln4}!> z_s;x3sIxJjO~fY$m3?UCl~r~7PxQ~C+kXkrUpj*RD?DgHDvkcF)a@ju8uoh+ zef%l;dm7FK_c$Wj7w77$bdXwwYy1cOr?|!|-b0^rHhaG~cK_G$_YK(LyqU_xr#1Sq z8c$e>ghsEzZ0AjBB;ncQopVsQN{tl?yMZ$2pm>?4SNE`0l6Zv{yOJo2{>U11#?=A( zccnAFuGBeH;^8$;{GtX+o1G2tTJb11>j8SVvndEgCu4V8y{OVChp`Hs(JDMwFt@5{ zt1-^&c>|2V8;ds=Y1~6SnM`ah>*LJ`ALJ5T2^vw57w;hjvCC(Wh$+%G7?U z>ojgh_URETYsE4XGpGBHqtu_4nx!HnE~(cxrRi=Q8|rP+_y(qRDo?9*Q@jq>C9W^k zChF6H0Pl2EwyV?`RSvv4z`LEb8<#|5XMh`>dGz*9lPuLvIq*&G0XCe)mH@q3dg*p< zLn*Cbfu=F0?d5$=v27ofJnLh0-Dtz4u=DHu+`-7q?yQxVLVX{UBa7`wcYycv0dyRE zGh*kUQ=o6vxC{9Tkt7CW`cPxIRQ-icxJ3Oc&V)~$s^At&c7H<%KK^D%6^VFKQl zv!b|wCal|1sT_oer7M=L;W@=hd7#Pk!d9624j%XONj?SHGPAZL+rIK7wRBqrTS%Op z<~NAS-{?5=k}XU-Dky9TN;kToawi^{wQcd=1)hC;hN%gq=m=U0skFT*8+w?PEB%Ny zmB2L<5xjuT)KMBSDvwq8__CCf08g^z=Ltq|tWo{sRN)FNZDqKU|Ng(e_=t~Ii4&Za zZeNWOl7iI;cZnyQQCbSUn#)KfRHK~acQraE8Ay%JOTJN~dnLoD(fyK6)Ckv)3U~aX zo_+ifQ~gp2lWC|r6Y3cbk6D>4vMf~HVaeM{y5b-0geEVD9 zM{8CQt!PIO9X{0+lJp#!j_1sDCMBQ4Lb+5baVB*>q%w6#eX9wu(sT#O6~wYehjfHrbo&i#3!}LZ1RkiQTY;v$A`L!yN}1xg?E26 zHJh2Rgk7dO=qm^D#^TvHytBu)Q@D*lHCfsLJ9ZHB1qzU)oCDC&W=TToiOZaNqQ0sp zI%)Mp2dAEBaMTk;Tb@!>)Ds0nJy9^!6SVCj|ox#33x7yeRFLibXW#{CQowa2?vJcYHKyj`L5O0N45cSX52%NY+W zIb(~Q5quCTENTfRsN%=*yOumZrd~j&r?c}%_!cHM?Ld|xs!_|54+tO7iGqMyzun~N(5&T&=K4R9;w}Xvk0ru(TzP!5}=ex^&d2cz+ z50?9KcNq=`miuy^*Dd;SP;tB+T<-?#1$09u6v={2pMsByOph6QJ!nfdgkOqx2A@Q1 zJcX2`?~sMyqtJ7UTF=u5BZdd7&&cXdwTjIz* zQh$UW7zlZ)UDkcF={vIQph9COpCy-}i#ou@pGlhk%B$B5y`zX;Vn-(!YpCic^q`N1 z@b8E8a!G!@5jNceXaqC?ngIr&1<(p;1B3wW09i+`AkAFCcVC75#doh7`aU^a8S+W+ z8c}X4DfMqSG7iN*O^qA50^g}i4m7#;+ z1A`DaQk(L@b}P@pcCQm zG+*I;&ChYC%LTI>Z-L6XT`t#UOfb9F=`_!Wp5q?3bYaI^T)`)!ElA&bwL)<#6m)Qh zo}-701cHz2+tjM3gg~vw=Mbpn+oCOQAp?{5MO$F;BSP|$zFlqdmXQ3kS1;cc{AO_z zyyR%K1x|8%6wX=8gF@#;eYc|HE1~mUZ#fx-UOfCIeXrW2m2CP!QDnt2)Ka5EG0=-% z!XgKXSdL;Y>Aen4$?@Fa@hmFOE&L0c0}f)zk-TRKVnLvn^kWWY^&)fJTvWEOHRP=O zOV(R2!9S$*;Uh=EYn&?!TF=YeBPAV;Cl4{4Las0>teWrQSOGn#kkP*eFU<3_Lymcda2jwIU^35lU(9n7kOI63kOj;F?g7jJ?gKmkc<>?%A1o~!cTwOP4v++h z#3vb$TnowY8f5+H%1Syw)BA4dv5Ft%_u#n=|Mlf9(CS;m{4hTvgWLIS{B{{+`O*IY DRreMS diff --git a/emips/run/run_with_script/run_meic_cams_htap/meic/grid_spec.py b/emips/run/run_with_script/run_meic_cams_htap/meic/grid_spec.py deleted file mode 100644 index a861c3d..0000000 --- a/emips/run/run_with_script/run_meic_cams_htap/meic/grid_spec.py +++ /dev/null @@ -1,80 +0,0 @@ -from emips.chem_spec import GridSpecReader -from emips.utils import SectorEnum -from emips.spatial_alloc import GridDesc -import os -from mipylib import dataset -from mipylib import numeric as np - - -def get_sector_str(sector): - """ - Get sector string. - :param sector: (*Sector*) The sector. - :return: (*str*) Sector string. - """ - if sector == SectorEnum.INDUSTRY: - return "inc" - elif sector == SectorEnum.AGRICULTURE: - return "agr" - elif sector == SectorEnum.ENERGY: - return "pow" - elif sector == SectorEnum.RESIDENTIAL: - return "res" - elif sector in [SectorEnum.TRANSPORT, SectorEnum.AIR, SectorEnum.SHIPS]: - return "tra" - else: - return None - - -class RetroGridSpecReader(GridSpecReader): - - def get_spec_fn(self, sector): - self.sector = sector - sector_str = get_sector_str(sector) - fn = os.path.join(self.dir_grid, 'retro_nmvoc_ratio_{}_2000_0.1deg.nc'.format(sector_str)) - return fn - - def get_spec_vars(self, sector, dims): - fn = self.get_spec_fn(sector) - self.f = dataset.addfile(fn) - spec_vars = [] - for var in self.f.variables: - if var.ndim == 2: - spec_var = dataset.DimVariable() - spec_var.name = var.name - spec_var.dtype = np.dtype.float - spec_var.dims = dims - spec_var.addattr('units', 'g/m2/s') - spec_vars.append(spec_var) - - return spec_vars - - def read_spec(self, sector, spec_var): - if self.sector != sector or self.f is None: - fn = self.get_spec_fn(sector) - self.f = dataset.addfile(fn) - - return self.f[spec_var.name][:] - - def get_spec_grid(self): - return GridDesc(x_orig=0.05, x_cell=0.1, x_num=3600, - y_orig=-89.95, y_cell=0.1, y_num=1800) - - -_grid_spec_reader = RetroGridSpecReader(dir_grid="Z:\chen\Research\EMIPS\Grid_speciation_data(VOC)") - - -def get_spec_fn(sector): - return _grid_spec_reader.get_spec_fn(sector) - - -def get_spec_vars(sector, dims): - return _grid_spec_reader.get_spec_vars(sector, dims) - - -def read_spec(sector, spec_var): - return _grid_spec_reader.read_spec(sector, spec_var) - - -def get_spec_grid(): - return _grid_spec_reader.get_spec_grid() diff --git a/emips/run/run_with_script/run_meic_cams_htap/meic/run_VOC.py b/emips/run/run_with_script/run_meic_cams_htap/meic/run_VOC.py index 52b2dc8..24056ca 100644 --- a/emips/run/run_with_script/run_meic_cams_htap/meic/run_VOC.py +++ b/emips/run/run_with_script/run_meic_cams_htap/meic/run_VOC.py @@ -50,7 +50,7 @@ def run(year, month, dir_inter, emission, model_grid): scc = emis_util.get_scc(sector) print('Read emission data...') - emis_data = emission.read_emis(sector, pollutant, month) + emis_data = emission.read_emis(sector, pollutant, year, month) #### Spatial allocation print('Spatial allocation...') diff --git a/emips/run/run_with_script/run_meic_cams_htap/meic/run_pollutants.py b/emips/run/run_with_script/run_meic_cams_htap/meic/run_pollutants.py index 61faa64..1675fce 100644 --- a/emips/run/run_with_script/run_meic_cams_htap/meic/run_pollutants.py +++ b/emips/run/run_with_script/run_meic_cams_htap/meic/run_pollutants.py @@ -61,15 +61,15 @@ def run(year, month, dir_inter, emission, model_grid): print(pollutant) print('Read emission data...') - emis_data = emission.read_emis(sector, pollutant, month) + emis_data = emission.read_emis(sector, pollutant, year, month) # Remove PM2.5 included in PM10, Remove BC and OC included in PM2.5 if pollutant == PollutantEnum.PM10: - emis_data_pm25 = emission.read_emis(sector, PollutantEnum.PM2_5, month) + emis_data_pm25 = emission.read_emis(sector, PollutantEnum.PM2_5, year, month) emis_data = emis_data - emis_data_pm25 if pollutant == PollutantEnum.PM2_5: - emis_data_bc = emission.read_emis(sector, PollutantEnum.BC, month) - emis_data_oc = emission.read_emis(sector, PollutantEnum.OC, month) + emis_data_bc = emission.read_emis(sector, PollutantEnum.BC, year, month) + emis_data_oc = emission.read_emis(sector, PollutantEnum.OC, year, month) emis_data = emis_data - emis_data_bc - emis_data_oc # Spatial allocation From 4f1e1c82bdff98c32b7bd6e4d07da061e4df90fc Mon Sep 17 00:00:00 2001 From: chen-0126 <252941087@qq.com> Date: Mon, 27 Feb 2023 16:46:38 +0800 Subject: [PATCH 2/4] fixed lack of molar mass for NOx --- emips/chem_spec/species_profile$py.class | Bin 6834 -> 7237 bytes emips/chem_spec/species_profile.py | 7 ++++++- 2 files changed, 6 insertions(+), 1 deletion(-) diff --git a/emips/chem_spec/species_profile$py.class b/emips/chem_spec/species_profile$py.class index e34cc0816a1a2408a56f8b36681e57907693477f..b666a1eeb4d4d5add75a84e5d2ee997c80787bb3 100644 GIT binary patch delta 1808 zcmbtUYitx%6#mZ6?C$Jz+Rjc(Tei5ZKuPKDQXm${rf5irrbGtJ8i<%MuI)k#yIs4x z#n+Dpe`pK|MsNHf7((!YPbe*H7Z4B?m4||FMG>DUz6GBYsAqNyH1UTf{5W&Zo%4O? zJNKMB)7P|bhI-)qj@;l^P`cxtbW+M?kB zq&(YlKl^C6$SLvl z1sb|>NX0?CA+~FiD-Pp`UD8_wr8$vo>yF2+R`8A((l*`l9^SWgKF~k~`HvJFCOFbd z6nw&K+X@Ar@|xiNX9P+q_<}WBX2y(omGp67K!IOx!wnsZV(g_9Ma-c1mi7NO` z6qlG~`&Cx|Rl^VX(e_Rgd}41&x##+1v@Mo$<9D&DWP$Q0{!;M={uU=oW-Eg@Yu9y7 zG?gxP5Q&FNZ}pHwvPuqeKBRYs$YXcR!TG32Lvz6^)^K;hE42rak9uvjWT{VbbB}WR z9Ou2OjlBzn+dHm&_M&iOYkRD-acM`!x-FGl7HiFj7Vix-My&NNEItULa=5@9!4Wa! zZ6d#DGS=7a&f%r7955P2%USp+EQctE*8`4ii1#7c%lD3`Mw~IK8DCMhHn`PzFo)R# zxWfz>*(~aVp)3}eA>QTK^N|_I;)M&eV`Vku?>|?zjeND7sSK_`Ll9M3j^$H?dA0{npEb}rt@Yn|9wvAl;$1A0>yO2C!HTxSkz zvRG#tk8x=sD~l(|%wpjqm>JX)5F zry?P{9nX?!1{6NlePnLLnt(q2f@04_$>;wqIpaSiC;d;!`hXt3lx#%lB3Nlqv$&?R zo@R>n$~xLCR#*BKZ?=z_hq((a;n{>NmlX*ZmqN;jklTvz@%qgCyhHZb98}9s4bHPT z&NJ#7@2DqxER#?qsH(DO8vE1PX0V0Xn%QRZpt`ZYFesctF(iI_qG14UHuU3dN#e%N z4F-&NB}wYzJ(1+qk;8}mIC?SIpyY${h~RGBrMvQ6$@n-Dlnph~6LcD5B0-mNEE04Z z$0I?-_&ma?@6z3Qx-(DbE6lkwGEba^hsj^sxK;YNJkB)&cmE}ouWglYA}n&!7Wv*5 zQEd^{LIl^OB!1Vo6ctR zQsE1%W%N3@hENaL%|916FLVVY_6zP1*K)a$<(HR^Z`_q!rh}?!EU#Q0jp8#eZlMaQ Ow4FIrMLyd}Q}sXNKAqqI delta 1491 zcmbtTYfPJE6n;+owf*|>LBD>K%fNy0lGR?qHb)?cOGr=yNl{ZXqlsXJ!pM3FZ9%7r z{PKq}(U>=jQIo~EO}D97ij=y^P}#g}x_O(rd3TEUsq?OLc;1#c{_n^6-sgPJd(Ly- z^PcxWYpz*4cxKN&0BfmcQk&Cs+h8mi$*6dgaGh9VYY;|Rz1XZauHCNT0_@OnA5u+qKXC&d+dN}a+p z8rEQ1%;{A$BaZ6p)Mv3rLjik*yS&}Kk6`s)vO&ivp4Tvk7esIQ>gsvym)ZviYI{66 zFdB=lSMjpgUOtg{6$fR)LprFW^SX+8f+anu;tgI0##Fq;Yn=Ue2$WLsE^GA1(kk93 zD6!#0RK;Z^>k>LcP&MV075q&P3I99RBnANzZ zFLiv1&m{Uv7_LV8Myz&iP`|_X8otF3;$By$`XhdlCH^dqyLv77O|%)eIPnMm)bKn0 z5<8595C6(GS!9+)9qmO#q~nefvdZpDmTcm3cOyB(P3~^>uIzl~5TQZb?5^@%IUI{7 zBJt>DgsVH-nnIheNhOC;k$6*gWGvDYie{q8uEb!nGddmQ+of#l*qYf|*jl**qnMQ)7-#ns%8LzJ@Xg>w-xOZ5T4$k$0)F!qtJRugpUCPq zc=lhL!r}i(pIRi9Wn|8<8}?$X)jSdoD5e$;*v#|70lWETIN&hf4hK~8XqYqKZa9k3 ze__sq%qT`ADjKB2`#am7;ushOZdxSQfl< zl8-E>fu93tocOw6uq($y$>FfNbEx`30u0#22Z zIq?gAUD|CH_X2ivq=3Im*$DsO^b&uP2ZKig>wSS>{ch1xf6+M3Wb*_z3!B1b&r)gN z Date: Sat, 24 Jun 2023 10:38:00 +0800 Subject: [PATCH 3/4] Add support for CMAQ and bug fixes --- README.md | 5 +- emips/__init__$py.class | Bin 3340 -> 3340 bytes emips/__init__.py | 4 +- emips/chem_spec/chemical_speciation$py.class | Bin 13375 -> 15082 bytes emips/chem_spec/chemical_speciation.py | 36 ++++- emips/loadApp.py | 2 +- emips/run/run_with_script/for_CMAQ/README.txt | 8 + .../run/run_with_script/for_CMAQ/__init__.py | 0 .../run/run_with_script/for_CMAQ/for_CMAQ.py | 137 ++++++++++++++++++ emips/run/run_with_script/for_CMAQ/height.py | 93 ++++++++++++ emips/run/run_with_script/for_CMAQ/merge.py | 92 ++++++++++++ emips/run/run_with_script/for_CMAQ/proj.py | 117 +++++++++++++++ .../run/run_with_script/for_CMAQ/transform.py | 90 ++++++++++++ .../for_CUACE/convert_grads$py.class | Bin 10268 -> 0 bytes .../for_CUACE/for_CUACE$py.class | Bin 5372 -> 0 bytes .../run_with_script/for_CUACE/for_CUACE.py | 63 ++++---- .../for_CUACE/for_CUACE_bak.py | 42 ------ .../for_CUACE/write_ctl$py.class | Bin 9824 -> 0 bytes .../for_WRFChem/for_WRFChem.py | 106 ++++++++------ .../run/run_with_script/for_WRFChem/height.py | 8 +- .../run/run_with_script/for_WRFChem/merge.py | 4 +- .../for_WRFChem/proj_2_file.py | 2 +- .../run_with_script/for_WRFChem/transform.py | 2 +- 23 files changed, 675 insertions(+), 136 deletions(-) create mode 100644 emips/run/run_with_script/for_CMAQ/README.txt create mode 100644 emips/run/run_with_script/for_CMAQ/__init__.py create mode 100644 emips/run/run_with_script/for_CMAQ/for_CMAQ.py create mode 100644 emips/run/run_with_script/for_CMAQ/height.py create mode 100644 emips/run/run_with_script/for_CMAQ/merge.py create mode 100644 emips/run/run_with_script/for_CMAQ/proj.py create mode 100644 emips/run/run_with_script/for_CMAQ/transform.py delete mode 100644 emips/run/run_with_script/for_CUACE/convert_grads$py.class delete mode 100644 emips/run/run_with_script/for_CUACE/for_CUACE$py.class delete mode 100644 emips/run/run_with_script/for_CUACE/for_CUACE_bak.py delete mode 100644 emips/run/run_with_script/for_CUACE/write_ctl$py.class diff --git a/README.md b/README.md index c4b5d35..de2d0b7 100644 --- a/README.md +++ b/README.md @@ -1,8 +1,6 @@ EMIPS：EMission Inventory Processing System ============================================ -[![Join the chat at https://gitter.im/meteoinfo/community](https://badges.gitter.im/meteoinfo/community/meteoinfo.svg)](https://gitter.im/meteoinfo/community) - Installation ------------ @@ -34,6 +32,7 @@ Learn more about MeteoInfo and EMIPS in its official documentation at http://met Version ------- +EMIPS-1.1 was released (2023-06-24). Added support for CMAQ and bug fixes. EMIPS-1.0 was released (2023-01-26). Author @@ -61,7 +60,7 @@ Sci. Total Environ. 870, 161909. https://doi.org/10.1016/j.scitotenv.2023.161909 License ------- -Copyright 2019-2022, EMIPS Developers +Copyright 2019-2023, EMIPS Developers Licensed under the LGPL License, Version 3.0 (the "License"); you may not use this file except in compliance with the License. diff --git a/emips/__init__$py.class b/emips/__init__$py.class index 65e9f72f7ccfff3c620ae19a8ac209c8d96be58e..8e12bef7b77e11c7275af64ae196fcff9fe735f9 100644 GIT binary patch delta 31 mcmeB?>XF*O!oqH>XJDvjWVV@`MUs`Z;~|@L{$@?i$IJkFGzbj< delta 31 ncmeB?>XF*O!oqH(XJDacV6vH;MUs`Z_1onOW}7uRA2S00g$D_0 diff --git a/emips/__init__.py b/emips/__init__.py index 7685688..0c2c5ad 100644 --- a/emips/__init__.py +++ b/emips/__init__.py @@ -4,10 +4,10 @@ import chem_spec import vertical_alloc -__version__ = '1.0' +__version__ = '1.1' import os ge_data_dir = os.path.join(os.path.dirname(os.path.abspath(__file__)), 'ge_data') -version = 'EMIPS-1.0 (2022.08.04)' +version = 'EMIPS-1.1 (2023.06.24)' __all__ = ['ge_data_dir', 'utils', 'temp_alloc', 'spatial_alloc', 'chem_spec', 'vertical_alloc', 'version'] diff --git a/emips/chem_spec/chemical_speciation$py.class b/emips/chem_spec/chemical_speciation$py.class index a43bf84d2984be3a2068e456e4b3aecf4c64d95b..bf8b988af087facc93fe66088957a9e2745a740c 100644 GIT binary patch literal 15082 zcmdT~33yw@wVoL-8aqmIeQy%7Ku{8s*pVG4&H^|g#8I5s!5fyd0##*OP6V=KBsmaD zODSC_G_?|biiFCXLkYv$aU zIp_T6%yRk6zux~rAw;u!vlL~+wX3D*tLw~U4>XM)&L7OAn+7x4L{snK?b-N9Vq;N5 z|G{h`J`~Nz^NEf1MXkMudj>8|4CXgVaoQ=m42Uu*0=Yy!l}sn3D6i|N-yIYcqOw|) zi+~h0xzWU6GLegA6T^vYB0ZSUqDl(q()gixQ!1W5&=k#Qlj#EiF+++uQ5Np=WjpeT z5iMpy{^S-dW=r9nJgMZ!Xh76RQ4NsE=C@}vBT`&k*Lm`28!sxZDxp0at)u>;5^dHp zLE(ux)xs5XrO;%C%wRkvMN^3}Pu&Tb>@+c-OtwIZ*}lm}vzg&!3I3pmpdM0eu>sJ?qmq<}n5Y}QDeinqZI16KrBu5XYk^`&KVFLM>XQsN9rH zC-a-7SX?&^2Z=rrv8`IH7HtT|SS+4O#bR3QkisptT|jg|`L>bKyp3-7cF{Dw2gL=V zvr6m~T~aJ8?r!^7dXQ8Nh#pvCICdbF8HlIQ=fde2aG~g{0^6fVEG7L{yGF|L>Ja|z zNaRBW_l7LbL*%;>)hTc$G;nxQOB!i!mFcKuG>#l-ntYbzT0(Uz>7AAJq5gx3Nw4|i zP-D^SCgkYVoJn1wOg2=ArlOuEO(2`dk7d*JVnN(SDBdbm`8x8jTI{BDmP@8{2z1Jx zdtqkSjP3mFd0WU=`^80+Xf8&+LN_*Foq@bPoxh$VV)W8;F_(dOF2RaxF@R2DF=RaO zO^T7?g=g~ma_at1-7@W@hD4%TG>c(ms#q+YNne%7W@sK|b@ls$;-E-U>!nge|J`~0 zC!?pRe=IPI>pG^D~XIlvpl_S+U~O=N8RU*a&S~UpT}!Kafsa4EUv&> zOVrZhuoM+zqoawe7SEI7-05|lcE};|eCqlIuwX24xfC0U`zq`RrtRp3;#%tHMMy@& zu~a;lPf8&>g5t&EB^329MgED!$PVemm3dO?dRF%3P^vqVhEg|5u{?n_Dc3ZJJZbW# zK}m3pv#80hwWEgv;s`Q*3HwKoagnBXiC0#Mo5ZW6xO5tIrpbS&>bgX|jWN#0udWuo z;x%wzVRfi&mg2i67xiJmD9$!HPHTmoQYfB=6fgK8T3&LK6B)!}!^kF-xF%j9vB^^Q z^@8N&?k3CS1_xV0(aczOFcBJr%W}zK*qRc1Xb^^EQ9pqdrsRWk2+u+H} zD52Mbo|G=I&fh`AO67y%P0WuwscAHo%xm#xB$jmI3R;d2B(f%?N-@l~!>w->cUNOg zy2lR>uxW&tbT4=PcHi!Prn(=CX&h0K9@65S75|K`G9jCQPS~xi5(?Q`-4TP<)X6@ipvhkUC3mzG~M=xx%<-2eue? zXBHYu4<)jpq2#bRt>i11vs7fPl>*uXWCSJ{?zZL!|2HLF`fx%g;y z@Z577LW2XV*N3P!M%CH)&`5JbXe4u0Je#lQSmNDCq1PBC^NBr$-9}JlC^KAQ7DZ09 z+nCmrWmzIuml@0ZJI;EOOXB3AOlnF?43+g(?GnoJI9hxYUQK1LpjGDEa3SfZ#dpDg zwpeSVtzV1pLjVf5_G$4$uI{{GtrkD_@j5Mj>SMh6xDoH~>D1!q2JYPl_)84R7K!xqv}^HujE^5ZTKv(7MEcsITKw7c($)`wznW^mRxSQ+s&{nm z)8e0|dV5b_mlh{Xhwa?q)8UTM*6hGF&FxEwOycP1y{FW zGBV)fW-Y4>ZtJJ|8C*{T%bDiY?)J4>;?Rm-+9RD>;+P6Gx_h+rxV;Bm;1mgqcJy?% zX?dEdHgJKd?rm@D)^d@VEln2kI0cdiI=AiB@=PNirEZq6OmtW0b}ez*gG_&0A2d3P z8@ig;cWAky@MEWz_2%ui2)kH>lk%1S0 z=<*UfMZQZKnW&Y2d zG5R^e^Xfz=84J<(vp zIQD#B3U6FvF!;7c%WHj%H?A|dwXX;DFYz(zUuJN`nCj&Q_ZfU6)^czd5!7eKIAYqO zZP5$0e5Ju#TWRVy8)c0g0#Cliyw*iw8|j1fUu!h(+KqsFotZPjKSrT-yZO~;f4u=c za$?o;jizO5$Bx!cymhDl%go} zzG2o!`5Ek-i`PfGeW+b8#UH9e`1fxtj?{r&pZ)(bn8UMRfK?(T#0K2;TwEQ7}N#Q>Rpw&_Slt3Gj+~Wkg z0Dcbl^>+32x1;r!=)4#CQ_Ejda}Sz-OI1jB|7nELasq7X?T++xY56qOZAaWcL!h-&D@h>Qp;Z}yet-%BUENv*Xb_4DRfQ^A zSJVtt8Oo`dKxRo%v(ZTFxT*o5Ra|*!>!Q_B%>^>+p*qbVtjKBskXf+QA|SKWsWb2@ z<#Kf0 zfiAv)?%2P-#QxU5Uo%>4EY)|3;W|p*=J1JFtySxAbW-bsLb8F+#x|W!6oqRst(&iY=(MdB(4r#TE#Qk=y_F#ge zqe-py;czxtORX+~Ubz@g(pFllLEM*5s-PXTR&ad)xwH#s zMBIIryxTGMhw>vvW*)9|i&#;=TZ<#O%~dO^L9oVx4G#Bp8;C|oVg z9mGXLJ`u&GeJZgvlgZ=g8%KSK8=>CAW-7%eY^pB23c6=iu3BBLa0wO2j-^?*^~LA2 zsjrn(E2s{sD~NVi;sXWjH02d3mXt8q6cKQ6*~tdzJ=`~zPE&Lj<{VI0OHq&hx8oj! z2olNQ;xs!pn9pSU62mEc%)*sML0pRF$@wHg7at|QgH4@DP`yB1Q>~t_UWm9ccb~)d z4xFwUJ4(gyf9vBL#i=3^yqIpY40R)qy{oqHnnwi^hTSWMOI!ldKgnWnHZC9odXINQ03%6 z1Z^ksPxsi!Kq70hDcpxIpM#fm#YgS$YQ6!ZY2v#L*p-J`yK(ytj9xD_%_+(wuuc3t$Obq&uMkcVT@%N;_fC882vO@;h% zxO`k_x4Pm{5$tX(b7w^x%Uvhh7%o39=Aojwr9AA;S^%MivZcJS-28U@-{A^+`Sf!8 z^5T=;TU?Nkd2g`XgfX zKK-a@-dEH3G{E0^=gx z9WKWVFPGi-in6fdkw6(<9T#J$iH9r5$CrnlIn)4p#ogL*aix0|{ka-{8pg#nEfr?8 zsg(;VV&mdE*@Aa4_WZCj*Y5%1WQ+wB5Ab*x{)U$F1?8W7=9aKC&&==@nBjyN6jka? z_*)~&PKfTOMNL37L2?v0A`yg*{;OH<(to9Biun*0ZP;;b9H> z+ptq`%2QWN!;QEd^Qk z6#m(0B&Pbu2~L=MQ=uKpf%Ts=bE5q`YBy>eXLkglZP=MgcN&R$4?EjU8_nKtv;8wii4~1LN1y(E*x6y^ zi6)at+CEV`iKQK;md5QiDXo#RJEyAUqUvdp<2?10*jw>W@v5i9UHJW!XsrmJkPA+fdq2wxiJdJ5btDI*>E5x3nwHF3y>2%#`%&XKN`CHR@U98T3Ro4ZO4#W)BAh)70b9teyQr0NkoyVZ1x1lzFcvLzzuUZv z3~`D3DMK5QG~E}{?;a!a88ZW-K3aM0+8ow4)9xh^njJ5SoTKh@KWC=;MevLw&Ln9f z>{m_hrEUwnc;;Bx!K$i_V3Dh53H7wN-(Qkij7#Bun6M8eiqem=3uQOT9+bWIxeq4B zrkzb3yao$DwgAC}VP^n?ncW6#FzC^VIGy%+slKf5jl+{)T~ef5)K8{R4w) z_fHIFxPM^~bpOU+ruz>Dv)m^c%yyq*;JE){P=m`NMy@BZeZ=$}&yz@Rgywj%lF?jG zYDV)sSGFsuu21a$BT*au~ldBo6^kg%muqSzE z>@;|C1J^Zr@*GC1JlVpi$&1>^^FC*-C$~|ZSvc+HF6?A{ zF3dJu@aJ3ZPq)IKZKXfifInAVm@3-rRMBXsidH+-s=`#O3sap@m};;v)tQB<&YC<` zn=$Uye+oO<02t4zBw0+eANK^b;ACd?iWSe|0GJK12BK2M=a zCs(1dn-Ns$WmN440W|63DrnNr2%79>G|LZ~*`D0XRgNe3GXl#mW(1Z$hY?u*Tt;BI z2q+u>4X+V0nrWSe|ACmfPq8%SXAc&XZY2V0oUz}d()Q6X4M1(t(tJ0RTGZ0YQk|=IXKQJXH@c; zTG<|Uu9?KtLr{7$Qy&l4i0gc<&$e1KLqXLmjj=%0s|g$O z3m8G=7czp%FJc6hU(5(9zmySFzMc_OwrqsTH*g(Pegz}P{Q(1LT*(L;-^2)1UCszp zeHEj*e%=RD=WrF6I*$>UYBT8qe^mfcH8(EuR|D{L7FU6%uVw_E-ogkxy_FGodK)9~ zbTK3FbSI-_{@PUQuPI>Z>$wdqy@L^0>M#OJ>llHhD;XgdGvS=D&p67(mU4{+${`n9 z${`n9${`n9${`n9${`n9dLb8EdaWzaYkh%U8z!eeKU3F;8+^+AFQsZbSFpy;6|AwL z0M^)0D4D9uc7WkWCPl!TV9-esP__c?H~Z?WHxXcog8*ROAOM&)2!OkcKS1-h7yy>v z%^pJA+Zdsd1IL*X0~|vRV&Ls82W{?W1a02Q2->`x5wx*1M+|U8Itcp*xD`Qglo5j9 zAw~!S%OC`S;yUE(aYl#&%PK^HW!3*czCO;QAr9Wp2yyTMM#$G6WQ2S@$_SBQlGZ*4 zk?>)zLnK&#BND8?5eXmVRz!mJJR-q*9+5D>t%w9pQO+tqh!F^naU%kuoDl-S1~vl0 zhVEKFIuQq-Ah9WHLyfr2C&_82TpTz=!Gu+27!ZL26g+6=fzf8;6OA@L(Q4xpnbO85 zSYj;+mRL)ICDxK)iM1qHVl4@lSWAK>){Qj8XQ%0|p0PuX;Ix~DX^p5du# zMvFZaWOSydY-EKz<#62+Pq~bidTI`%WuBVHsMb^S8J*>+g^ZSa>U2gcJhhlnou@*K z>OHlT&`AM`8@s8yH#R6d$M%Hh*qrbVTNBP_mKc1(?otITmVw(__-%$-%KtJL$XxNC%ksB~{S zuCC?oqm6M1i0vyI8W#~&@NbLx%in$Kfa!K`H45Ft(GA|FarKf1#?{ME zUXF4j$`O=TJ|L-oX+Db_7lCK+VtGK^ByoilP;U@rzzm;(4*yqMaQ}HHv{34Hbr+sz bh%)tR==7Qnb&Gl}ft~7BbsK@4y8S-@I7H9g literal 13375 zcmc&(33yY-wLWuYTVo;M`U(h!ga~lJ_O&sXC5a&fBMdh2R?85a(5kYz2Bos3N-~Ml zv}v2JX}Y9q+NS$LlWr+#4Q-mP>Avs#-gHa%y&-A(pP9L`jD=sm_rCWYALD!O%$YN1 z`OlfTqtE{P10N!yc59EIz<9$(K|_tb+5FMA+*I*cHq$ni&3kQwQ+xBt32*zHf#G9$ zFL@+hOcuTEO>;&Erv^r^@y3eV1)X;ekwFRw3KqO#I+gJRRWuf`aM_ zId3fG6%u)G+{=5JF;7yhAp4r+@nl;%nK{}PFXmI3qd{6AXh~c%5Ba=3MQ=jVLhzqG zB5ARp$k|GzCUQZl6BGiG%@_CPvlD_2H};;r+V-pGc2%%mJFTbb>I!4VGhvF*k`Otx zRFGshWXFb*?F{_T~;q>vG1~6K0BUD!(l50Eiw)m^!qM=Z*$pnda{_z z6cZ)Z3*oClzq`b|T2O6?S<+g(mY5`6gf%Bpxv6w&bVFux!po<|B&`>;Kznf@<8>AU zbw4_je4m}$+AuXyGrPXhQ|_DrJP0usuu{(8XvJJ{dnYqvY)X*&VCi_`XgWKZOhe|1 z`64_(gS9~OV^BCM#;0ME5LKZl{_XRM(Ng%L2FFp3O@;0>FcBS{DjUh5tZh(&I&w+W zI5p&>pwJO*Z00%5P0`_FUfDE%9Uh!BU8NLHEvV8I&E}(}Y?&jeOo6;toXlri$$l7_il6;hc35}B*xLAV({(|&&RxGdG9Lo~vL%DxI z^$0XJT@|KrI?AKR1jYZmdi@uxU&HdBz{%M&j|HlBPq6ng%3CQ^6xQV8VN%plx$y zm3I7j1wDmsWI<0wF&a;#lZ9eRkmw21)9EJ8`e&g2Boge0jCWm;je3?g_SrD1KbwJ3 z&lR-J!=6)U8$0GrsJf`A6tJeNq^h3{xv3yM50$)v|6{1AD8q4jel6WXFA#LiJm$<( z|IQ`4Lchr|U5{TBqJFv+;VV`DhIT>UQ?+Oa8^YXbQ_y9tbP|dtGvLx2dX%@1LU|zr zTVfp5gbUZq6lN>uGSf?(<>Jcu3S(m((Rg+;KjuZp5VAsQ9KPnl9vy=tHLG91=F}>! zn535iI%yBv+RG8iY>uhh08cI#*vDVRz)Baxbi2mKYj`M^P8B7+4uvJ-UCZ0;Q4g3$ zN>xIb8An*(KyM6TPkNJ|9zfFsBk3(#^jm$u`<3c#*i4hilFSiFcVgY~^km_fq<4TE zOJ@s!+`H7)naL_Z-lYNg9^iJZq`Q&0W7$kGi7G7VeQZW56Q)yotIx{$HQ#q0>CV5*(=JUy^S$hGGL{;~!3AwY+Gw0Y;DsE~oB);* zrjMk?)h9_Y$~`K66W-XdWF}RZ&^~PJN(}aG+Sm~- zBy;((OD}DSj*V{I8s**ucjuEwCfZw~6WQyN`C?Ofx|IYaxQNf3CCu@GJEPh0b1vL5 z6FhT#=}a4Rd4xvqlo=Y>G~1j(GVB%M`SMlD>gxrnA@b?(!{! zj_s559Ynq>(HZL+mh?R^z~IgyNk7oty;p3O^dq0%BIzeSjj2CV?85`SlI~aJp%IY3 zR2;iId-q8Cwc*~=F6jZ}WLF{<8yM)8^gFDNmjOwCP%N>buDGN>DJfmUVEBvb2H7d; zZ>oD=?}(&-sP4T3Lw%C|r37|sfoIf9&*0DimVH=p^~47TB|W0(zRuW?6haB=>$*aU zfbQ;sGNRh2+ocF9x@(yGYjr;_EEcG#{oR|TSg6SESg#cLE`=HW15(uK@d1d4sA)X| zyt7nNT2OR-+f?uwzm$&l!E)!c!hedxyL7uEueetd3$jt^`>|0bVCKS~|? z`v*3of2&WU|6-pfCdH+GKWKagX_=reGs}X!Ty=9!NzsXJl|WL&6bE7>MVBID-CLyC zr^xOY+}mTwtx{ZJ$V;Ta2MM4#&57;;#O#~2ct#uXWhBPs{f%H6SUtZ}94W-y8o zMP31*i>r;woe=j}<3VJ(6bUuZ-`NR!lRgdUW7sv6-+B_AQh5H`ettWOL%bQ=y<3Wu z|GGcHuTRi=^+L#mAt5`f7L9dpW)wXUFQ)29QBbSK+Pk+&F{$az-EhFQem_*6QtKl* zy83~!C;42M@nl5MQ!ruA{!#)jbimne(-4kF=RrgTR7CpfiUq27kf z+bbba7gi}=twwhD?Cb2stk?K&N{iR2mEjv8<_%h2d-qnB_a;5k-Ht`yqDSCc2EYlW z5_@_8f26p>_dHK~yHXz&?0}O0PW8IGtN$wX`feo%MNPfEM=jAg*n|ANThaa9oI3CG zX@vEZvbVnz-S0OfdhgN74KP^l)yoekYklVKS?U*X5NcG1UGb%m_@ekSDqF!z zj|*z-eRPNS?P2j%@im;!#n<7t9oVNQ)80-gz6pzQ9ONUspw~6_RJx9vb7wdA%C6_C zt16rr{2L2p#r8_FD=gPj@jR-(F~qmUcW@XL-wl({4*ot4K72+iT{KAXLzrX^-%|V- zfz#)ADSisD&79Ju_&L0a)3_ABfJ^vj!dIhG{0c|Wj1<2CMAS?vehYW=t*R8i2gQdr zDgKB#Jj{nFDgKNLpgBjyf)sxRW=l&-@psmMwdG?6hW}ZnQR;{k4>G$yF2zGoFqcb7 z@oyY=%dM0a?5`;#lq!}5g;Ub1!rkQT(~h)iaONvJq)1ERhOw-J&nD6e!DIOF?Zc@X zw@DShmX-fu{Dk6JhD)+JsAxJU>3O&x6Q2*!EmmFF3R?@Aid4AGVsX9MB#uNuYbi2g z44-aAFOF-@w6{B(E#hM$nXCAXJvgOQnfuMc-KDA9!`6A$@{qO6s)ud)$&9vb>)aFX z>}eIv3R^3!^8q{S0$l0BPqU^7T3x|mvsktZT01wzhbA)_&hC=VL2I?3Cdl85FG~hU zEQ_nC{Nz|Mn;-JV)42Y}M`DRxPZw!~dh7D>zRr20;_w4=Bpzk2poO^5=u1+PGjUYEuQU3{LsJ3R;^5MQ3X{k%PPP zZ7BT*y?g;4g1f7aC)1N&0Kdkp$CYbE{=VT9euk+g63@K$jt3xiRYf*G77!ZW)6v#9 zHJHyHg-hD@^Uo%IUeU|;WX7|5yyI=fY&Jcby{@gRuV*mc_9)8oh8&IxE6+wo*+pK^ zx(q4D5k%JZq5$ou-z4*V)j@>h-%ePrrvYjs)+e z_K|?IRIjVRL>e`7e7zY9!%ZqPYq$}N$Tythx$j{_tH@7NqV7k0F0;T z#o*iG^VQGgyT^H1iEoF`*E^T*YtDyCd^>%wWahlj3uCA&|n_ zg}==RC}=bbjq{}n4I!l<8jsbZ^`adH7ZTV6ek1>`h~E2M+s4w$#4BC94z%lBIZao$ zFd`DJlLS3Y-ZZ6NNe!;;D28<-f+0_|q~?PwuEztIh7LuZE1Li4M%c+RI&gsV@QS?1=j6yYzZk{y800 z^;MsG8ocPiD4{}h2LCzec8H-5euUq`H@iG zeXZ}^niB-qHO&ZoeVM@EjKF5a!Se&B>CMx00@B!^*|F+zZC9DihwW*|eqUD;1d7_L<&y!q^j3;dBSJkdsl(4nwuWNu-3H~Dw{o4 zTiC9>lewH#m4mxzM%v|SHJ;31T!5#t>is_Iu5j&LN)$Ycf@9Z~FUC_A&&cXjvKE8G zF-QTpICenF^{yRLa*%*5N2|+7FPE_AoDy8SOO5gRzI0xRrMvcCHOAWeV@AHlD5Ij) zFF+6f!L|1(euhaoNV{ix%TU^*dU@S`Rnl6yvg=Z{j%pvKg052!(V^;p(Tg6U*W>*m z>a2Fph!tn(_=Dixg&8lWi|{v!wgzo2S_2yQtwURn)`ZrK=AyNr@mMR`2DD9R?P#0P zwxDfA+lF=tS_j&8v>j-dq3uM&!GU(6aY2ot?Lq58+l$8Y_n~#8^`K_rcxZN<1G;94 zTGt+e)oK%jcjN2<=OZd`T>PAm0eEnHzgrEjVpnnPVa1}WXDMMzweQ#u=xIGz0c4s`j$j5-uJ&<=P`Ggc;|AdJ$O<)1Cg7)M4#{^ALXg9YTwv4Wk`E zJBW5A+9C5Cfs1kO)-Db|1)Dy$|8Tu)k76-(et-|geC?fYvkTPL6!4vU#2@T;zKgQc z3iF>P`Z?3B54!d-jB*vVx9FqP_rbzPsUNCmwew^3tZ{y-oM5OHsHa`y{7Mr+ z=Qo;=&TlnQ>-=65A?J^pSm69w6Jh7Cnpo)kT@#C(e`;c}^Pnbd=OInhIsev#gRfLg zMIr(xLg;BPi3rp=rj|xTji#1Ggw)h|5fReV@`wm)sy-qXX=+78*qT}y5ss$LkBB9j zx*#H!X=+tOEZ5Y95wSv3(TF%-Q>!Clm8RB2L{w92BVvuF8X}@WQx`?VI!&#Mh(=AV zkBDYXHAX~>rkWyRgQl7zVxy+qh-lYTOGN0yuH70D+jQTCh`2;kZ4uF-sf`h#55xAR zh}fz7+9N_8h)3+r5wV;5)V^uXcCMZEp|Ds%!PmFSSGU^Nw#HXB=<6y=s$$Hjib10) zMvbbqB~?QuRTq?04VP41SW-1T1xtoB1xtoCwaCwz z#Sw8(ci9nfNKuF%yHvOpF>c5jJBc5NOOq+>Dus zn=un{GiCyT#!Ljxn2Ep{GZ8pr90I3|QzrRHt?72{8_SSdgwf?lz2B{)n|!M;Hdbqd z0;u*gB{+0ri-loTZ z(pPB;C|#r}pmd$4fYL@yp%!bv*{+W`yQKuV))M4Uiw)#Ziw)#Ziw)#Ziw)#Ziw(R` ziw(TCl;E|s1g~x7<n{+U@uj!B0o#bMG)M_hISug^r1ekYEx2;ELyEmhu?uFW7uMLe&;K*< zN5uVV!Y|`32;r~eEePKO@fL*cckvb*U4Os>eSg9Ooqxdt!TcK@h~z);Kp6jp2V!^z z4+QXGJP^G{@K_wN1S_(`u2+gwmnm{jwOvvi!mFGp^l(0=)PMrPh!tQptZEjp)9zIQ zYRUwBt;(No1=$d*mJP8MF#j(5u;LHS@Qcs92cG83U;H$uJ&kKizR}0U?@93kzCIwT zwH8^6&FE@Gfk&6#L${1XEc=w@Oj}D%TPvDRTdTCNTI)h9YS8Sgx@6 z{@}G#^9UwX1@V&xK3;>?P6~k5r!hw69sF>e*pFrJnaAg diff --git a/emips/chem_spec/chemical_speciation.py b/emips/chem_spec/chemical_speciation.py index 63d5cb0..4afdc68 100644 --- a/emips/chem_spec/chemical_speciation.py +++ b/emips/chem_spec/chemical_speciation.py @@ -4,7 +4,7 @@ from .pollutant import Pollutant import mipylib.numeric as np -__all__ = ['get_pollutant_profile', 'read_file', 'speciation', 'get_model_species_wrf'] +__all__ = ['get_pollutant_profile', 'read_file', 'speciation', 'get_model_species_WRFChem', 'get_model_species_CMAQ'] def get_pollutant_profile(poll_profiles, pollutant): @@ -89,9 +89,9 @@ def speciation(data, pollutant_profile): return s_data -def get_model_species_wrf(mechanism_name): +def get_model_species_WRFChem(mechanism_name): """ - Get species in wrf-chem in different chemical mechanism. + Get species under different chemical mechanisms in WRF-Chem. (E_PM10: saprc99, cb05 E_PM_10: radm2, mozart) :param mechanism_name: (*string*) The name of chemical mechanism. :return: (*list*) All species(out_species) and aerosols(out_species_aer) under chemical mechanisms. @@ -99,7 +99,7 @@ def get_model_species_wrf(mechanism_name): mechanism_name = mechanism_name.lower() if mechanism_name == 'cb05': ########################################## - # ------CB05, emiss_opt=15, (52, 15)------# + #------CB05, emiss_opt=15, (52, 15)------# ########################################## out_species = ['E_ACET', 'E_PAR', 'E_ALK3', 'E_ALK4', 'E_ALK5', 'E_TOL', 'E_XYL', 'E_BALD', 'E_ALD2', 'E_CCOOH', 'E_CO', 'E_CRES', 'E_ETH', 'E_ETHA', 'E_GLY', 'E_FORM', @@ -113,7 +113,7 @@ def get_model_species_wrf(mechanism_name): elif mechanism_name == 'radm2': ########################################## - # ------RADM2, emiss_opt=3, (42, 19)------# + #------RADM2, emiss_opt=3, (42, 19)------# ########################################## out_species = ['E_ISO', 'E_SO2', 'E_NO', 'E_NO2', 'E_CO', 'E_CH4', 'E_ETH', 'E_HC3', 'E_HC5', 'E_HC8', 'E_XYL', 'E_OL2', 'E_OLT', 'E_OLI', 'E_TOL', 'E_CSL', @@ -126,7 +126,7 @@ def get_model_species_wrf(mechanism_name): 'E_ORGJ_BB', 'E_CLI', 'E_CLJ'] elif mechanism_name == 'saprc99': ########################################### - # -----SAPRC99, emiss_opt=13, (55, 15)-----# + #-----SAPRC99, emiss_opt=13, (55, 15)-----# ########################################### out_species = ['E_SO2', 'E_C2H6', 'E_C3H8', 'E_C2H2', 'E_ALK3', 'E_ALK4', 'E_ALK5', 'E_ETHENE', 'E_C3H6', 'E_OLE1', 'E_OLE2', 'E_ARO1', 'E_ARO2', 'E_HCHO', 'E_CCHO', 'E_RCHO', @@ -140,7 +140,7 @@ def get_model_species_wrf(mechanism_name): elif mechanism_name == 'mozart': ########################################## - # -----MOZART, emiss_opt=10, (53, 23)-----# + #-----MOZART, emiss_opt=10, (53, 23)-----# ########################################## out_species = ['E_CO', 'E_NO', 'E_NO2', 'E_BIGALK', 'E_BIGENE', 'E_C2H4', 'E_C2H5OH', 'E_C2H6', 'E_C3H6', 'E_C3H8', 'E_CH2O', 'E_CH3CHO', 'E_CH3COCH3', 'E_CH3OH', 'E_MEK', 'E_SO2', @@ -154,3 +154,25 @@ def get_model_species_wrf(mechanism_name): 'E_CO_A', 'E_ORGI_A', 'E_ORGJ_A', 'E_CO_BB', 'E_ORGI_BB', 'E_ORGJ_BB', 'E_PM_10'] return out_species, out_species_aer + + +def get_model_species_CMAQ(mechanism_name): + """ + Get species under different chemical mechanisms in CMAQ. + :param mechanism_name: (*string*) The name of chemical mechanism. + :return: (*list*) All species(out_species) and aerosols(out_species_aer) under chemical mechanisms. + """ + mechanism_name = mechanism_name.lower() + if mechanism_name == 'cb05': + ######################################## + #-----CB05, emiss_opt=10, (42, 19)-----# + ######################################## + out_species = ['ALD2', 'CO', 'ETH', 'FORM', 'ISOP', 'NH3', 'NO', 'NO2', 'UNR', + 'OLE', 'PAR', 'PEC', 'PMC', 'PMOTHR', 'PNO3', 'POC', 'PSO4', 'PCL', + 'PNH4', 'PNA', 'PMG', 'PK', 'PCA', 'PNCOM', 'PFE', 'PAL', 'PSI', 'PTI', + 'PMN', 'PH2O', 'SO2', 'SULF', 'TERP', 'TOL', 'XYL', 'MEOH', 'ETOH', + 'ETHA', 'ALDX', 'IOLE', 'CH4', 'AACD'] + out_species_aer = ['PEC', 'PMC', 'PMOTHR', 'PNO3', 'POC', 'PSO4', 'PCL', 'PNH4', 'PNA', + 'PMG', 'PK', 'PCA', 'PNCOM', 'PFE', 'PAL', 'PSI', 'PTI', 'PMN', 'PH2O'] + + return out_species, out_species_aer diff --git a/emips/loadApp.py b/emips/loadApp.py index 4edf49a..ec79639 100644 --- a/emips/loadApp.py +++ b/emips/loadApp.py @@ -10,7 +10,7 @@ class LoadApp(PluginBase): def __init__(self): self.setName("EMIPS") - self.setAuthor("Yaqiang Wang & Wenchong Chen") + self.setAuthor("Yaqiang Wang & Wencong Chen") self.setVersion("1.0") self.setDescription("EMission Inventory Processing System") self.app_menu_item = None diff --git a/emips/run/run_with_script/for_CMAQ/README.txt b/emips/run/run_with_script/for_CMAQ/README.txt new file mode 100644 index 0000000..325214f --- /dev/null +++ b/emips/run/run_with_script/for_CMAQ/README.txt @@ -0,0 +1,8 @@ +Read before running +---------------------- + +1, Set parameters(time, directory and dimension etc.) in "for_CMAQ.py". + + + + diff --git a/emips/run/run_with_script/for_CMAQ/__init__.py b/emips/run/run_with_script/for_CMAQ/__init__.py new file mode 100644 index 0000000..e69de29 diff --git a/emips/run/run_with_script/for_CMAQ/for_CMAQ.py b/emips/run/run_with_script/for_CMAQ/for_CMAQ.py new file mode 100644 index 0000000..88750c4 --- /dev/null +++ b/emips/run/run_with_script/for_CMAQ/for_CMAQ.py @@ -0,0 +1,137 @@ +""" +# Author: Wencong Chen +# Date: 2023-06-21 +# Purpose: Convert model-ready emission file to meet CMAQ's input requirements. + (height allocation, merge files of different sectors, convert pollutants and projection) +""" + +import time +start_time = time.clock() +#Set current working directory +from inspect import getsourcefile +dir_run = os.path.dirname(os.path.abspath(getsourcefile(lambda:0))) +if not dir_run in sys.path: + sys.path.append(dir_run) + +import os +from emips.spatial_alloc import GridDesc +from emips import ge_data_dir +from emips.utils import SectorEnum +from emips.chem_spec import get_model_species_CMAQ +from collections import OrderedDict +import datetime +import height +import merge +import transform +import proj + +#Choose Chemical mechanism +mechanism_name = 'cb05' +print('--------CHEMICAL MECHANISM: {}--------'.format(mechanism_name.upper())) + +#Set year and month +year = 2019 +month = 1 + +#Set input&output file path +dir_inter = os.path.join(r'G:\CMAQ', mechanism_name, r'merge\{0:}\{0:}{1:>02d}'.format(year, month)) + +#Set original grid +original_proj = geolib.projinfo() +original_grid = GridDesc(original_proj, x_orig=64., x_cell=0.25, x_num=324, + y_orig=15., y_cell=0.25, y_num=180) + +#Set target grid +ncols = 164 #the number of x dimension +nrows = 97 #the number of y dimension +nlays = 7 #the number of z dimension +xcell = 36000.0 #Grid spacing in x dimensions +ycell = 36000.0 #Grid spacing in y dimensions +proj_name = 'lcc' #Name of the projection ('lcc' is Lambert Conformal Conic projection) +xcent = 110.0 #Longitude of projection center +ycent = 34.0 #Latitude of projection center +p_alp = 25.0 #First standard parallel +p_bet = 40.0 #Second standard parallel +p_gam = 110.0 +radius = 6370000.0 #Radius of the sphere(unit: meter) + +if ncols % 2 == 0: + xorig = (ncols-2) / 2.0 * xcell + xcell / 2.0 +else: + xorig = (ncols-1) / 2.0 * xcell +if nrows % 2 == 0: + yorig = (nrows-2) / 2.0 * ycell + ycell / 2.0 +else: + yorig = (nrows-1) / 2.0 * ycell +target_proj = geolib.projinfo(proj=proj_name, lon_0=xcent, lat_0=ycent, lat_1=p_alp, lat_2=p_bet, a=radius, b=radius) +target_grid = GridDesc(target_proj, x_orig=-xorig, x_cell=xcell, x_num=ncols, + y_orig=-yorig, y_cell=ycell, y_num=nrows) + +#Set path of the vertical allocate file +vert_prof_file = os.path.join(ge_data_dir, 'height.txt') + +#Set sectors that need to be processed +sectors = [SectorEnum.INDUSTRY, SectorEnum.AGRICULTURE, SectorEnum.ENERGY, + SectorEnum.RESIDENTIAL, SectorEnum.TRANSPORT] + +#Acquisition of model species (gases and aerosol) and aerosol names based on chemical mechanisms +out_species, out_species_aer = get_model_species_CMAQ(mechanism_name) +varlist = '' +for i in out_species: + varlist = varlist + '{:<16s}'.format(i) + +#Set global attributes +gattrs = OrderedDict() +#gattrs['Conventions'] = 'CF-1.6' +gattrs['TITLE'] = "Created using MeteoInfo, mechanism: {}".format(mechanism_name.upper()) +gattrs['IOAPI_VERSION'] = "{:<80s}".format('$Id: @(#) ioapi library version 3.1 $') +gattrs['EXEC_ID'] = "{:<80s}".format('????????????????') +gattrs['FTYPE'] = 1 +gattrs['CDATE'] = datetime.date(year=year, month=month, day=1).strftime("%Y%j") #test +gattrs['CTIME'] = 81255 #test +gattrs['WDATE'] = datetime.date(year=year, month=month, day=1).strftime("%Y%j") #test +gattrs['WTIME'] = 81255 #test +gattrs['SDATE'] = datetime.date(year=year, month=month, day=1).strftime("%Y%j") +gattrs['STIME'] = 0 +gattrs['TSTEP'] = 10000 +gattrs['NTHIK'] = 1 +gattrs['NCOLS'] = ncols +gattrs['NROWS'] = nrows +gattrs['NLAYS'] = nlays +gattrs['NVARS'] = len(out_species) +gattrs['GDTYP'] = 2 +gattrs['P_ALP'] = p_alp +gattrs['P_BET'] = p_bet +gattrs['P_GAM'] = p_gam +gattrs['XCENT'] = xcent +gattrs['YCENT'] = ycent +gattrs['XORIG'] = -xorig +gattrs['YORIG'] = -yorig +gattrs['XCELL'] = xcell +gattrs['YCELL'] = ycell +gattrs['VGTYP'] = 2 #test +gattrs['VGTOP'] = 10000.0 +gattrs['VGLVLS'] = 1.0, 0.995, 0.988, 0.98, 0.97, 0.956, 0.938, 0.893, 0.839, 0.777, 0.702, 0.582, 0.4, 0.2, 0.0 +gattrs['GDNAM'] = "{:<16s}".format('M_36_08CHINA') +gattrs['UPNAM'] = "" +gattrs['VAR-LIST'] = varlist +gattrs['FILEDESC'] = "" +gattrs['HISTORY'] = "" + +#Run all scripts +print('Allocate according to height...') +height.run(year, month, dir_inter, original_grid, sectors, nlays, vert_prof_file) +print('Merge data from different sectors...') +merge.run(year, month, dir_inter, original_grid, sectors, nlays) +print('Distribution of particulate matter and change unit...') +transform.run(year, month, dir_inter, original_grid, out_species, out_species_aer, nlays) +print('Conversion projection...') +proj.run(year, month, dir_inter, original_grid, target_grid, out_species, out_species_aer, gattrs, mechanism_name, nlays) +print('-----------------------------') +print('--------All finished!--------') +print('-----------------------------') + +#Calculate running time +end_time = time.clock() +time = (end_time - start_time) / 60.0 +print('Time: {:.2f}min'.format(time)) \ No newline at end of file diff --git a/emips/run/run_with_script/for_CMAQ/height.py b/emips/run/run_with_script/for_CMAQ/height.py new file mode 100644 index 0000000..1331145 --- /dev/null +++ b/emips/run/run_with_script/for_CMAQ/height.py @@ -0,0 +1,93 @@ +from emips.spatial_alloc import GridDesc +from collections import OrderedDict +from mipylib import dataset +import mipylib.numeric as np +import os +from emips.utils import emis_util +import emips + +out_species_unit = ['PEC', 'POA', 'PMFINE', 'PNO3', 'PSO4', 'PMC'] + + +def run(year, month, dir_inter, model_grid, sectors, nlays, z_file): + """ + Allocate data to different heights. + + :param year: (*int*) Year. + :param month: (*int*) Month. + :param dir_inter: (*string*) The directory where data is stored. + :param model_grid: (*GridDesc*) Model data grid describe. + :param sectors: (*GridDesc*) The sectors need to be processed. + :param nlays: (*int*) The z-dimension of the output data. + :param z_file: (*string*) The path of the vertical allocate file. + """ + print('Define dimension and global attributes...') + tdim = np.dimension(np.arange(24), 'HOUR') + ydim = np.dimension(model_grid.y_coord, 'ROW', 'Y') + xdim = np.dimension(model_grid.x_coord, 'COL', 'X') + zdim = np.dimension(np.arange(nlays), 'LAY') + dims = [tdim, zdim, ydim, xdim] + + gattrs = OrderedDict() + gattrs['Conventions'] = 'CF-1.6' + gattrs['Tools'] = 'Created using MeteoInfo' + + for sector in sectors: + fn = dir_inter + '\emis_{}_{}_{}_hour.nc'.format(sector.name, year, month) + print('File input: {}'.format(fn)) + dimvars = [] + if os.path.exists(fn): + f = dataset.addfile(fn) + for var in f.varnames: + if var == 'lat' or var == 'lon': + continue + else: + dimvar = dataset.DimVariable() + dimvar.name = var + dimvar.dtype = np.dtype.float + dimvar.dims = dims + dimvar.addattr('description', "EMISSION_{}".format(var)) + if var in out_species_unit: + dimvar.addattr('units', 'g/m2/s') + else: + dimvar.addattr('units', 'mole/m2/s') + dimvars.append(dimvar) + + out_fn = dir_inter + '\emis_{}_{}_{}_hour_height.nc'.format(sector.name, year, month) + print('Create output file:{}'.format(out_fn)) + ncfile = dataset.addfile(out_fn, 'c', largefile=True) + ncfile.nc_define(dims, gattrs, dimvars) + + data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + dd = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + # get vertical profiles + scc = emis_util.get_scc(sector) + vertical_pro = emips.vertical_alloc.read_file(z_file, scc) + if nlays != len(vertical_pro.get_ratios()): + print('----------Warning-----------\nThe z-dimensions in the vertical profile file and the output file do not match.\n----------Warning-----------') + # read, merge and output + if round(vertical_pro.get_ratios()[0], 2) != 1.0: + print('Allocating: {}'.format(sector.name)) + else: + print('Do not need to be allocated: {}'.format(sector.name)) + print('Write data to file...') + for var in f.varnames: + if var == 'lat' or var == 'lon': + continue + else: + print(var) + dd[:, 0, :, :] = f[var][:] + if round(vertical_pro.get_ratios()[0], 2) == 1.0: + data[:, 0, :, :] = dd[:, 0, :, :] + else: + for lay in np.arange(len(vertical_pro.get_ratios())): + data[:, lay, :, :] = dd[:, 0, :, :] * vertical_pro.get_ratios()[lay] + # Turn nan to zero + data[data == np.nan] = 0 + ncfile.write(var, data) + ncfile.close() + f.close() + else: + print('File not exist: {}'.format(fn)) + continue + print('Allocate of height finished!') diff --git a/emips/run/run_with_script/for_CMAQ/merge.py b/emips/run/run_with_script/for_CMAQ/merge.py new file mode 100644 index 0000000..19cead3 --- /dev/null +++ b/emips/run/run_with_script/for_CMAQ/merge.py @@ -0,0 +1,92 @@ +from emips.spatial_alloc import GridDesc +from collections import OrderedDict +from mipylib import dataset +import mipylib.numeric as np +import os + +out_species_aer = ['PEC', 'POA', 'PMFINE', 'PNO3', 'PSO4', 'PMC'] + +def run(year, month, dir_inter, model_grid, sectors, nlays): + """ + Combine all sectors into one file + + :param year: (*int*) Year. + :param month: (*int*) Month. + :param dir_inter: (*string*) The directory where data is stored. + :param model_grid: (*GridDesc*) Model data grid describe. + :param sectors: (*list*) Sectors that needs to be merged. + :param nlays: (*int*) The z-dimension of the output data. + """ + print('Define dimension and global attributes...') + tdim = np.dimension(np.arange(24), 'HOUR') + zdim = np.dimension(np.arange(nlays), 'LAY') + ydim = np.dimension(model_grid.y_coord, 'ROW', 'Y') + xdim = np.dimension(model_grid.x_coord, 'COL', 'X') + dims = [tdim, zdim, ydim, xdim] + #set the definition of the output variable + print('Define variables...') + dimvars = [] + count = [] + dict_spec = {} + print('Add files...') + for sector in sectors: + fn = dir_inter + '\emis_{}_{}_{}_hour_height.nc'.format(sector.name, year, month) + if os.path.exists(fn): + print(fn) + f = dataset.addfile(fn) + for var in f.variables: + if var.ndim == 4: + if dict_spec.has_key(var.name): + dict_spec[var.name].append(fn) + else: + dict_spec[var.name] = [fn] + for var in f.varnames: + if var == 'ROW' or var == 'COL': + continue + if var in count: + continue + else: + dimvar = dataset.DimVariable() + dimvar.name = var + dimvar.dtype = np.dtype.float + dimvar.dims = dims + dimvar.addattr('description', "EMISSION_{}".format(var)) + if var in out_species_aer: + dimvar.addattr('units', 'g/m2/s') + else: + dimvar.addattr('units', 'mole/m2/s') + count.append(var) + dimvars.append(dimvar) + f.close() + else: + print('File not exist: {}'.format(fn)) + continue + + #set dimension and define ncfile + out_fn = dir_inter + '\emis_{}_{}_hour.nc'.format(year, month) + gattrs = OrderedDict() + gattrs['Conventions'] = 'CF-1.6' + gattrs['Tools'] = 'Created using MeteoInfo' + print('Create output data file:{}'.format(out_fn)) + ncfile = dataset.addfile(out_fn, 'c', largefile=True) + ncfile.nc_define(dims, gattrs, dimvars) + #read, merge and output + print('Write variables data...') + for sname, fns in dict_spec.iteritems(): + print(sname) + spec_data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + dd = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + for fn in fns: + f = dataset.addfile(fn) + dd = f[sname][:] + #turn nan to zero + dd[dd==np.nan] = 0.0 + if spec_data.sum() == 0.0: + spec_data = dd + else: + spec_data = spec_data + dd + f.close() + ncfile.write(sname, spec_data) + ncfile.close() + print('Merge data finished!') + \ No newline at end of file diff --git a/emips/run/run_with_script/for_CMAQ/proj.py b/emips/run/run_with_script/for_CMAQ/proj.py new file mode 100644 index 0000000..db48243 --- /dev/null +++ b/emips/run/run_with_script/for_CMAQ/proj.py @@ -0,0 +1,117 @@ +from emips.spatial_alloc import GridDesc, transform +from collections import OrderedDict +from mipylib import dataset +import mipylib.numeric as np +import datetime + +def run(year, month, dir_inter, model_grid, target_grid, out_species, out_species_aer, global_attributes, mechanism_name, nlays): + """ + Write Times variable, add global attributes, convert data's projection. + + :param year: (*int*) Year. + :param month: (*int*) Month. + :param dir_inter: (*string*) The directory where data is stored. + :param model_grid: (*GridDesc*) Model data grid describe. + :param target_grid: (*GridDesc*) Target data grid describe. + :param out_species: (*list*) The name of the output species (gases and aerosol). + :param out_species_aer: (*list*) The name of the output species (aerosol). + :param global_attributes: (*OrderedDict*) The global attributes of the output file. + :param mechanism_name: (*string*) The name of the chemical mechanism. + :param nlays: (*int*) The z-dimension of the output data. + """ + print('Add input file...') + fn_in = dir_inter + '\emis_{}_{}_hour_transform.nc'.format(year, month) + print(fn_in) + f_in = dataset.addfile(fn_in) + + #set dimension + tdim = np.dimension(np.arange(25), 'TSTEP') + ddim = np.dimension(np.arange(2), 'DATE-TIME') + ydim = np.dimension(target_grid.y_coord, 'ROW', 'Y') + xdim = np.dimension(target_grid.x_coord, 'COL', 'X') + zdim = np.dimension(np.arange(nlays), 'LAY') + vdim = np.dimension(np.arange(len(out_species)), 'VAR') + vardims = [tdim, zdim, ydim, xdim] + all_dims = [tdim, ddim, zdim, vdim, ydim, xdim] + + #Set variables + dimvars = [] + dimvar = dataset.DimVariable() + dimvar.name = 'TFLAG' + dimvar.dtype = np.dtype.int + dimvar.dims = [tdim, vdim, ddim] + dimvar.addattr('units', '') + dimvar.addattr('long_name', '{:<16s}'.format('TFLAG')) + dimvar.addattr('var_desc', '{:<80s}'.format('Timestep-valid flags: (1) YYYYDDD or (2) HHMMSS')) + dimvars.append(dimvar) + for out_specie in out_species: + dimvar = dataset.DimVariable() + dimvar.name = out_specie + dimvar.dtype = np.dtype.float + dimvar.dims = vardims + dimvar.addattr('long_name', "{:<16s}".format(out_specie)) + if out_specie in out_species_aer: + dimvar.addattr('units', "{:<16s}".format('g/s')) + else: + dimvar.addattr('units', "{:<16s}".format('moles/s')) + dimvar.addattr('var_desc', "{:<80s}".format("Hourly emission")) + dimvars.append(dimvar) + + #Output file + fn_out = dir_inter + '\GR_EMIS_{}_{}{:0>2d}'.format(mechanism_name, year, month) + print('Create output data file...') + print(fn_out) + ncfile = dataset.addfile(fn_out, 'c', largefile=True) + print('Define dimensions, global attributes and variables...') + ncfile.nc_define(all_dims, global_attributes, dimvars, write_dimvars=False) + + #Times + print('Write Times variable...') + daytime = datetime.date(year=year, month=month, day=1) + daytime1 = daytime.strftime("%Y%j") + daytime2 = (daytime + datetime.timedelta(days=1)).strftime("%Y%j") + t_out = np.zeros((tdim.length, vdim.length, ddim.length)) + for k in range(0, ddim.length): + for i in range(0, tdim.length): + if k == 0: + if i != tdim.length - 1: + for j in range(0, vdim.length): + t_out[i, j, k] = int(daytime1) + else: + for j in range(0, vdim.length): + t_out[i, j, k] = int(daytime2) + else: + if i != tdim.length - 1: + for j in range(0, vdim.length): + t_out[i, j, k] = int(i * 10000) + else: + for j in range(0, vdim.length): + t_out[i, j, k] = int(0) + ncfile.write('TFLAG', t_out) + + # + print('Write variable data except times...') + for out_specie in out_species: + data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + dd = np.zeros((tdim.length, zdim.length, model_grid.y_num, model_grid.x_num)) + if out_specie in f_in.varnames: + print(out_specie) + for i in range(tdim.length): + j = i + if i == tdim.length - 1: + j = 0 + dd[i, :, :, :] = f_in[out_specie][j, :, :, :] + #Conversion + dd = transform(dd, model_grid, target_grid) + #Set the fourth dimension + #dd = dd.reshape(12, 1, ydim.length, xdim.length) + #Set default values + dd[dd==np.nan] = 0 + data[:, :, :, :] = dd + else: + print('{} no data!'.format(out_specie)) + ncfile.write(out_specie, data) + ncfile.close() + f_in.close() + print('Convert projection finished!') + \ No newline at end of file diff --git a/emips/run/run_with_script/for_CMAQ/transform.py b/emips/run/run_with_script/for_CMAQ/transform.py new file mode 100644 index 0000000..61fbb5c --- /dev/null +++ b/emips/run/run_with_script/for_CMAQ/transform.py @@ -0,0 +1,90 @@ +from emips.spatial_alloc import GridDesc +from collections import OrderedDict +from mipylib import dataset +import mipylib.numeric as np + +def run(year, month, dir_inter, model_grid, out_species, out_species_aer, nlays): + """ + Change unit and match variable name. + + :param year: (*int*) Year. + :param month: (*int*) Month. + :param dir_inter: (*string*) The directory where data is stored. + :param model_grid: (*GridDesc*) Model data grid describe. + :param out_species: (*list*) The name of the output species (gases and aerosol). + :param out_species_aer: (*list*) The name of the output species (aerosol). + :param nlays: (*int*) The z-dimension of the output data. + """ + print('Add input file...') + fn_in = dir_inter + '\emis_{}_{}_hour.nc'.format(year, month) + f_in = dataset.addfile(fn_in) + + #Set dimension + print('Define dimensions and global attributes...') + tdim = np.dimension(np.arange(24), 'TSTEP') + zdim = np.dimension(np.arange(nlays), 'LAY') + ydim = np.dimension(model_grid.y_coord, 'ROW', 'Y') + xdim = np.dimension(model_grid.x_coord, 'COL', 'X') + dims = [tdim, zdim, ydim, xdim] + gattrs = OrderedDict() + gattrs['Conventions'] = 'CF-1.6' + gattrs['Tools'] = 'Created using MeteoInfo' + + #Set the definition of the output variable and ncfile + fn_out = dir_inter + '\emis_{}_{}_hour_transform.nc'.format(year, month) + print('Define variables and output file...') + dimvars = [] + for out_specie in out_species: + dimvar = dataset.DimVariable() + dimvar.name = out_specie + dimvar.dtype = np.dtype.float + dimvar.dims = dims + dimvar.addattr('long_name', "{:<16s}".format(out_specie)) + if out_specie in out_species_aer: + #g/m2/s to g/grid/s + dimvar.addattr('units', "{:<16s}".format('g/s')) + else: + #mole/m2/s to moles/grid/s + dimvar.addattr('units', "{:<16s}".format('moles/s')) + dimvar.addattr('var_desc', "{:<80s}".format("Hourly emission")) + dimvars.append(dimvar) + print('Create output data file:{}'.format(fn_out)) + ncfile = dataset.addfile(fn_out, 'c', largefile=True) + ncfile.nc_define(dims, gattrs, dimvars) + + #calculate grid areas (meter) + grid_areas = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + for i in range(tdim.length): + for j in range(zdim.length): + grid_areas[i, j, :, :] = model_grid.grid_areas() + + #add data to ncfile + print('Process data and write to file...') + for name in out_species: + data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + sname = name + print(sname) + if sname in f_in.varnames: + data = f_in[sname][:, :, :, :] +# data = data * 1e6 + data = data * grid_areas + ncfile.write(name, data) + elif sname == 'UNR': + data = f_in['NR'][:, :, :] + data = data * grid_areas + ncfile.write(name, data) + elif sname == 'PMOTHR': + data = f_in['PMFINE'][:, :, :] + data = data * grid_areas + ncfile.write(name, data) + #####test##### + elif sname == 'POC': + data = f_in['POA'][:, :, :] + data = data * grid_areas + ncfile.write(name, data) + #####test##### + else: + ncfile.write(name, data) + f_in.close() + ncfile.close() + print('Distribution of particulate matter and change unit finised!') diff --git a/emips/run/run_with_script/for_CUACE/convert_grads$py.class b/emips/run/run_with_script/for_CUACE/convert_grads$py.class deleted file mode 100644 index b2244c61255e4b130da824971d06672255809188..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 10268 zcmb_i34B}Sk)N?->sgkaXNf?RgD53{ZP_s~0fHO?#&Kd1TQ-)G;BW}pdbX9wPePKD zm^7sqg&xaxTee%+Zs|drLTMWU$Sx_A(j$e^QnsZn-R*9>2i>#Xt6OgRpYJ{SR1EE| ze;K`RzHer}nfYesoA=u zV%%t7v=AH18Hr@HkSG}Kjf+8N-EkrZY|<(nSs?SKcV3Q>MY>YzQ~*_E05N zxyeZ`rix5nqiQDKLjKfv)fg)Ot51D69$H1K#n?4a!JHU3a;agB)-o+KpmTo1@R%_k z&u0xxu4Aek=-AfVrqNkUHP-q>A(h47qs^?_vrU1(`fuO8>*?9TA8Z0q|8)dE7O?`OK`}+Y@{|fHPQJ{ z8dl`v!xM>NL!%3rPG782Pbyz<(M3!v7Oz={?V^j}yiucIB2Ab!lsT^@!>ygRNSs^( zr*`ZZ2zT!4iR~Kb(&$o5#N&{gC=_z>xQ9A~FhrMQexNHF-qsa~g*$pQ+J;Ps4Rl1J z{e1&5jke>UqTS*Cs7Bpx+AMZD(wT!A?JR4Al`Aglh37{QabdGy;A8pbn%Cl1LG=|dck3Te+ajZzqNc5l*X zzYL=ADH6}o73}s!QMB}0is(l)!zZO_hKs6H2{h4$5gBjD*X`# zR}EtD#~4Jqq0*mV5EVl{4C?5vp6wd_DW>|nIyL$;_@=*idpLp}@4#|@qz@c_F3Wu# z8vTV#Mf+g#T{6`Rj!)v)66v%`_mq*>=+j6L(>N#rF>2%_vC7s>d%fL5pP|pXk&H)S zdcoA_9;TIxlVCd~C4neTB2 z@Yk6!4;Z;Zz%&ZO$&tW#CTXNwjI*$BBA+xX2HRa>T zE(jiakRFnPSOil%4s}udaGx|RSwyGl%Tg#FR#AHDK(g-TuZZ|ZnChD5{b}@76!I)u zheFz9TC+se5;Lsu3!;nV&26EUQ2NwM z9(sbF6vMv*E#P4*nLPAe>rVL|vOJq8jA`_JjP_+x$jZM&z9aNB`XN{jA2}S_dR|)+ zivNgd{14Ue)Cw+nXiv~jMEB>IZu$S;`486VuSBQ6#xWX)kjSWme}kxq$IVPZF5#bB zf%Xe^$C(-PMkAMz({(g7?)A{$(?1CJKjPrW6Z;KhHTZwY)ZDV@Z#$Ao&Sud;T2%fIq{D=nO$_gk8aGTB<}hx_ z&5Om(wDsY}CAusR{TKbWB**`h?fzpE*|g!J|8>)Er8}utUD)9xarJM*XfCPIi?DGh zWe%g17tou%EOj59J%+EMIgO!61<+{@87R!+lGhR5@%+RP(}hduyp$i9xdMJ*2aYms z+~DC#3$KFYNX9fY)|i&h(=fBq+LQUV)Bm7?mZFGV^gWVQRb*OKPxM$;x~9Op9l& z*pGcZyv~|<3mi6_&g7xRdWf;=sA(LM80K?S4D-2AAe;|~b;^s}2quYVC>z8uKbFXX zzeV;k4G*_kqVbxz#k#WWn<6ilUCNZ3LfuyhGg-Ph6rRikkURn0wZmiU)~(aH4at?C z%)@MrH({e}VnP(X5bna6t1C+5%^*tT3zMQk2t*_4OroIi7N**Cs!&K9al=fe#4G6e zow3{m=q{6}nJgHYT++yCj855|SJ7NZCDIzBQ660!9&?NE5;0xooolIhq9=H}#M2HY z|5>N{(>RRfZPJn8pj9I?oEb-Q;Eq_mHnKJ9+qRn0Np&9rZ^J7zL`316H_beZ`%FT@-9DmUH7 zgNP^ea2lg!#|1Ch5szzpt-{6;?1G`9It$2hjT1OT8O8I3L=FKkj1}p~G)9H0vU)6y zN8zX)p=%089whKpl~hbzcf219cJecQ{)F0);@?2CsL&HP?|AFdjHzlYe0YrH5`s<_Y?uS(Ub zW?OXxFHHE}B;P4e;yw&A+^HjQIo_INC^ZA`O^gt_Ln!gxLtca!Z%t0AWE$hOsmc`7+-%=WF8ONBOI6 zeuN*xwWJ#T5&QL4`u%XZ`2RnBOht@Sg#^JjP{pm$x5|V=|76J+uzOcz54gtr2*&h;B?Ens^R6ZJX7D@BEgB8=HM}l5= zOOa0FNU%y~Pj@~-De z=pd9=SeL>oZET0as%$JQ*c4qMB89VOng$?xH{7Dv+M=P~tmvy2FzlVfP*;xOngtBk z&0$zpj^X+R3`27myf%X#uytB5N+*Rw%vIQk!eFk##uNr~6}C@dE*ncLOtUdlVS0%) zYoR6DoQ0N%3l>^!^B%BJw~ZdM&>9=P(Ly~o`c@0Awb8d(Xq}C|-9ney=sPV`x6yZ7 zXuXZT$3nd}daH%{Z1jB=dYX-Xz(SYX=xr8yx{ZF&LeH?#4_W968@=5^{Wkg$3tefW zAGOd`HhQOpuC~#SS?C%Y{kVmmX``R8(6u)DDM77(yt@?RD{<<&f0{mrD808dxdxL+ zhsX)KZ!q{E9cwwxd%O=IR;GAwB2KC67li%5>4?UyEsy$=kMT82XzAk zJ|42m0!AgyQHyg3oS21T%|-f8k)142AfB{(*DsW6J{~De@5l5B+B|4wjuXq4Nj@Gc zks$W4={{&90@J1G@p;o`X&SMq)mQx5-X~j{&d41Y3*j)VDy>BHt+L%s3gD1V_9Tj_`h4Tx`oN(o@s)jiclURy@T0$9dow>Usq~ zYgPAGD<7(?@xIkzrr00o9ue(T3d!K=WygyK068fO`SZtM*7ADBL8+qBvdHU-$`cBx}e@ROZQK6N%qLO=;w=? z|5CDOlVJZgs{uBflB*R~{Zv_n`vrF3_i691L>O`jP6|~Xt?czzRZi1yLN31xcdi!H zU4Bies`*_`HKCOzg8D@ebAo<5SZUn=I_}UZDgHWGrjTYYsLeuZvvN-JzBH@jte3E- zQuSxhHslKG)qdb`7KU4tDAb^*y{{-I=+%n!DOq&S(p{?PUQ;LHy%!JV%M&5jQDiDZw%?e>_9L<`bemN)|SOXnoV;ag< zt7F-w^tf9oZ~}9V2ZFlC4_6=v(F#c7g8b=F4Q_ZpHU}eyJJinNLphY1qx^hOuk(8p zSGSV5(3x_;KCYQ#Nu8Rhg_BfBmf}zOv^JZMLxrFDh48;epY^uVt_2j_lqKk z%22I1R>^kRvZ?AQJsZ>`bEsF@)apPk3F>`VQwJg&C94a+u}&1& z4cY3z`ju47OhGy8aMb33(&lS}`js;~ivhJI1CG+!K|MOdZ|z*-Sn#(HATt zTy9f)^LC5<5{n&cxLtjGDNu~9+ZD;SL{nM*SYAkGrD&8e|)H!!8&k zIgF&nvAF&2ldp6Hb*V0>HJCvyJo)ndK^+YZYK3Sr??k1{qSAw~%3pQzrK8rJgp0c& zmhv>9d=}Z@e@b49RuJU!-NiG-1T;z1| z53G6Xm*`dcd5JIuLLb$XC0xjKfz8Bv<;%3|(sdOt(&?|$?Ns}Q+?C)Y{5Ojzo||5v zLoZSN3)Ju;)vI3NaXJt8cmtpj&;$qqngKGu5zq!WAFv5bka1kx73rGU!-9e~RL z+W=jF9e{2?7;pt(C!hxq0rUa-0apT|fEYmN2Lh{CTButc4##m!FgxmXn&uBYz_;mD zx>v0u>_(!^h`O<8b0wP3Xmgd1&1kdB$JuBzLL(PN-{|8H;(@sj;epB9@u>FkNAPg_ z_@j8#`1npdJU;#y9<@IHI39IA{sbP&eEcaqbRXZ1N4<|f8x_5-(bp)&oJui#UxH$L z?uQNN;pFp~ZJ7=`rKfWq=i82{)oQ+n@3mLg!D?Aueu86zKK}d^f1$`MU1A<$FXVtnkA3_4 zCy%l4vuZ^9$?-a^mbj0x75?Oy3DMV~>H{=N+I{LJqFMESIL diff --git a/emips/run/run_with_script/for_CUACE/for_CUACE$py.class b/emips/run/run_with_script/for_CUACE/for_CUACE$py.class deleted file mode 100644 index 2c1354fcd512771338400e6a5eb8b9b63429313f..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 5372 zcmbtXd3YOF5&yl&w$`$eDx1W)nr3s`mTfhLUUo=voupN9TobtoNubMGTG^XQyJEFU zBT6Y(p&SX6tCXu;g;K6GZ5`8)a$n^t_kEWe`cZEFq5R&iB*&JVFMR0x-tOqlyqS43 zznS;;;r~4N5P&}XTOgF(c%i_aj$y~0NEFK6q+=)2j%y}H%e!17Zw}Nk#wK0U$fP{O zGY2~BxTEE3#t)fkZ$Kcd>lr7l>k2{wk3=?TgC(kW`?Yx8M{xJLK4hP1n(1yo zK&)W^E37D55s5%So51pVJAQV;SSZk1G`)$OGj8Mr5)Ew3pQFybQJjZZ1dFhk`J$ej zWt##`9ZAVrf~AtROkg4NCG((9xJIU^V7b7Oxh`4xLKrLQ4hhHgu91qaIyZL$q$SUS zrYKh9{0LTIjesKkbkasHinW3B8m;=d)YD#0!Fuk%v+|~bjRNXiKS#@9Y-VzH<_lg0 ztVMy<2j>q{6dmYn#TIl4tg2VHt7NAoaKh-OA6Y%f1$C}Eo624!S{YmxU$a}XH*OVh ze0YvIvmiSM>jAn01p4Cr1%qX(QH9PNWvz*Rt)pvC$!6j3PsugAI<+wZoB~=Pz`0(3 zse&ib--#LUoGU;p6&}S!xL5+@l7@v4NQH5kK%1`9jNy5%u17H-SM_5XGpXxLV@}rv zF01cA=h$eX*pADk!lw#c(2<-!<`U9tafS4Fh_>mGk0OCJg!}x! z%ckKfND8#q17g@JdSN_GU{Osx@7Ug?g53fw+brAiwo6~n0y*jHH5iRx7*D4jnckok z=^2Z5Q^a+#V0f01)BT>v!n869t|gj|=yoZuU?0hHK<=m2Bf8^S6ABKJE?4Mh`SWQr z$C_aZI3LgktDE@o8!!8+U-P}ay$aH75(Tz5FK5fB%^x%AzKLvvH8a68hVCn{1jg!@ zoSEda#bbV3GXK}h{ImTn+2?okUbX@MSSb2#)7d3G^*;rLfH9bp(zaEcG&Accr@&lBr%d%vL&kO-~`4+rY?(%H{vAMT^Fz#eO+QvdI<;?90-pNQ4Gcr5^cXlKX z*54b>WhQyxK!XU(+QC30-3`no8-liyzjxz360vthA%aZ4kC>MoD>!cyyr2K&d8`*b z!}Sz=fKdwi=ZA=2zcuns@e$H&Zo0(vy^MgQ+4B_K$K_`>xd#Lm9HG!jzRP@nlSQac@jYu(ylq3vKqYhTDqbJNT|l z_4f#3#%b0PfprZ6HcN!x1RLs0YEQ|wWi{1&4&z4xoz%aJ_Z*wwAUum0hf>;e+&yMC z$FB?Cwrb)6eP`D15Ia<{!M%F6CsF(aKaJqW_?f)i$~#E5a*LB)G})+b{;#hgSsM!^ z2p*OF&ZmCe01l(&YBa+5t$;S$k9>gvaT9~1$-Sm4KYbX`=EFv=WQO=Tw}qG3hVr}R zu`$pLh+7`Je4JpmZuk7M}E*HcJ5A&rV_l-6p9G~DN*Zh zGgNP(Oef~1_>3bcjK8pUWGZ-?4f9)Lq?8{wUB8<$5&Wo3U%%TZRQ^X&PRUK1(lLRR zjEt<@y;jj8S_W;~;f<0(C}JJsdOHLBJB0YT#74_#tdPyffAe_ZJIwz~0mt9??*X)+ ziSMIvwMlK6Msx~lJk~OWcK_eLcq}aG74ZerxL_o%uIQPZhZ zs!ngJPPcLT1mYv7(6c{&KeqOqz-{W2Y3?LmPpyY4atj{7%TM4{$LWy3-}!6@vqcxB z+u)YKPeGuKTw*8KLTHRUZehWr@mNe+mul(vV$yMS(-bb1kEe7`VQ?hgL_<5I`%5G@ zG=;0wtH~URH&+=;$vB1K{Y|IQzdsa+j8G(X5_?ACEmg5)we@5E`f!!eK8K-H8Ov*o zz5a?=P1g#^kukqYeMZ2lXA{R$G;| zW`;E$u%cDgTC&s(<2sGWcubRmD!sl`Jw(CoSPNst5OM=OdPq;jX_>k_=`QJw;1-CT+bm@70rb(L`)ma`IK9E&-Ya4h9m#?j6pc`G

2pV;Y&#*epr=mW{XRG=#J|pTwe732d;4`XziqD1WXZT#CelA6ehhlBM$tK^V z@J-%!FP)PugZGPy%Q3zM=bBF7ZrM&`|&0XC+tM8Gv;0IW~r4^Y7*4?j; zPvJMmMbO;-!MRa{o?up*!Z;?vSQN%z5h5LsPte8xmBW_!PwHKO$MG+|*HFdp>A@e8 S_#^%#iL3D#{w#?i9{(>k9MO&d diff --git a/emips/run/run_with_script/for_CUACE/for_CUACE.py b/emips/run/run_with_script/for_CUACE/for_CUACE.py index 76ac600..7c9954e 100644 --- a/emips/run/run_with_script/for_CUACE/for_CUACE.py +++ b/emips/run/run_with_script/for_CUACE/for_CUACE.py @@ -4,40 +4,39 @@ # Purpose: Convert netcdf model-ready emission file to GrADS data format for CUACE model and write the description file of the output binary data files. """ -import convert_grads -import write_ctl import time +time_start = time.time() +#Set current working directory +from inspect import getsourcefile +dir_run = os.path.dirname(os.path.abspath(getsourcefile(lambda:0))) +if not dir_run in sys.path: + sys.path.append(dir_run) -def run(run_config): - """ - To CUACE model ready emission data files. - - :param run_config: (*RunConfigure*) The run configure. - """ - time_start = time.time() - - year = run_config.emission_year - months = [run_config.emission_month] - dir_in = run_config.run_output_dir - dir_out = run_config.run_output_dir - model_grid = run_config.spatial_model_grid - xn = model_grid.x_num - yn = model_grid.y_num - xmin = model_grid.x_orig - ymin = model_grid.y_orig - xdelta = model_grid.x_cell - ydelta = model_grid.y_cell +#set year and months +year = 2017 +months = [1] +##Set directory for input and output files +dir_in = r'G:\test' +dir_out = r'G:\test' +#Set dimension length of output file. +xn = 324 +yn = 180 +xmin = 64.0 +ymin = 15.0 +xdelta = 0.25 +ydelta = 0.25 - print('Convert to grads...') - convert_grads.run(year, months, dir_in, dir_out, xn, yn) - - print('Write .ctl files...') - write_ctl.run(year, months, dir_out, xn, yn, xmin, ymin, xdelta, ydelta) +print('Convert to grads...') +import convert_grads +convert_grads.run(year, months, dir_in, dir_out, xn, yn) - print('-------------------') - print('---All finished!---') - print('-------------------') +print('Write .ctl files...') +import write_ctl +write_ctl.run(year, months, dir_out, xn, yn, xmin, ymin, xdelta, ydelta) - time_end = time.time() - time_count = (time_end - time_start) / 60 - print('Time: {:.2f}min'.format(time_count)) +print('-------------------') +print('---All finished!---') +print('-------------------') +time_end = time.time() +time = (time_end - time_start) / 60 +print('Time: {:.2f}min'.format(time)) diff --git a/emips/run/run_with_script/for_CUACE/for_CUACE_bak.py b/emips/run/run_with_script/for_CUACE/for_CUACE_bak.py deleted file mode 100644 index 7c9954e..0000000 --- a/emips/run/run_with_script/for_CUACE/for_CUACE_bak.py +++ /dev/null @@ -1,42 +0,0 @@ -""" -# Author: Wencong Chen -# Date: 2022-11-27 -# Purpose: Convert netcdf model-ready emission file to GrADS data format for CUACE model - and write the description file of the output binary data files. -""" -import time -time_start = time.time() -#Set current working directory -from inspect import getsourcefile -dir_run = os.path.dirname(os.path.abspath(getsourcefile(lambda:0))) -if not dir_run in sys.path: - sys.path.append(dir_run) - -#set year and months -year = 2017 -months = [1] -##Set directory for input and output files -dir_in = r'G:\test' -dir_out = r'G:\test' -#Set dimension length of output file. -xn = 324 -yn = 180 -xmin = 64.0 -ymin = 15.0 -xdelta = 0.25 -ydelta = 0.25 - -print('Convert to grads...') -import convert_grads -convert_grads.run(year, months, dir_in, dir_out, xn, yn) - -print('Write .ctl files...') -import write_ctl -write_ctl.run(year, months, dir_out, xn, yn, xmin, ymin, xdelta, ydelta) - -print('-------------------') -print('---All finished!---') -print('-------------------') -time_end = time.time() -time = (time_end - time_start) / 60 -print('Time: {:.2f}min'.format(time)) diff --git a/emips/run/run_with_script/for_CUACE/write_ctl$py.class b/emips/run/run_with_script/for_CUACE/write_ctl$py.class deleted file mode 100644 index 5e37e0d0737bda57a393dd3d9f48d4f9d2f49d33..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 9824 zcmb_i2Yg%Sala!#;t_%*2($(@p-mDbK~cjZDNCX#QiQ}6f*O=8BOv)mLB=xhD2Y-O zXE{#yN=LfV?bcn6oKS3Br|C$YrYqfP(>6`hrn^y_wyE7Z&i}p#Jc1uU+DXfwj_?#As zuI0@5RMd)FW>;ewIx;sraowD>x|lXDqtZhzCQsh9l8Lm*gkI(b3;9a6Nt_9!|#7GvTO}OQdJY4KI2WS;1-_ zt)?|z^3ht@Atq*E3{&`x=`Xl0?yRGLxKqbemB|~l9v-I>*|}t5!k~>ze%)LT)ic$k z60=q!Y0wrXukBu6E|X$#WnKtOLADJ6S z!tilNUC`8R( zYN8f&n9U>%3ac0Nd;l-;1TyF#qIxWmPFTm8wl}O01sSx%bi_-0=wYU6Yc6YI+&Y*7 z<&Js?>Fw)HWvxQC=b5&wpxZ}}(D7P2Mkkmym#gcWO;1{hOxi=;NbBj?OfoYOPeSMB zl?JVcPS+yCdSPsCHofNnR}FZB__>%%SZ2_=W(KFs{A4bX71hD%MA8gqrWHRkYh`Dx z;6x%F&&>rf;&BD%TfIs_XEu(RADlDexz1oiQzC6OH3rA=)mjRtGHL5tUIdaz2}nq0 zNDnr!M0RQ_p;YD7%!KTiSeMA5PGn%g&1+`Pv`q*m^1+)qD`8G8QE)50XvwXXDdcB5 zFJ7U$f-~YD9FPDM?JsFF_XgYCN+ndQ>0~A@K8c-@Jc)QRn9QVS64vaL3ghy(5~V9hdW4C+TYm~l$rHk9d=((x>vlY%#Z zMzI(&Fc*t4ZEaYw%@7THX&;Tq)Ql{z+^(miqCScV0-Nj^^2Q5Bp)mXCJY5hq7qRSQ zdVo({kEST)0Dh6xX8VOJ?cHxPmv! zLq3|IN$DN6u=dW~JDtbd+NbUsG>sxQoyn!*mON(Drh~ydcZ0G~isw+; z?sT@cPv328lQQRgK>A&fo_oX6?`B$mTS*5G1dGZK1P?&bd$944EkH5bz8>Fc^QJ7( zGiFZB(@LpQz`K0(jr2_t%x}h|85|rm=v(n(re_TLc0_YFn>8_!-pkZcJ|?BJr&1rz z9G;->lA-@@rd^ipMf>5R7j12~+YVsn+(mqjw58H;>3vwLZ^UzX6|;lw;Cvs`?qINI z7~Xeu1bdM_vY5ocXh>zSjD%C|;d~4IQrg70VA8w?H&w`(eeqlx~o673#1 zt%(Xp6a7&o+S7kX6BUjo`s21}bU+scM-zQUiS~~5Yofx@M1M+&_V*mrM1`Y?{&Yd~ zaGB^~UG!%QqDRX_kLsd7&$K@n?CR zy)>?i3P%(D4JA6%uV<=oG|}HwqS3OsDjZGpx7GZQo*(Gb#z3)3BJjI*pn7{Wxe;(m z#C}f&YGkl4GNccFfmtH_2deXtp<&%8#VQf}BfInAZoPAGO2q!e))yVt%gKnaN(BGZ z);Cz@3^*lXpT-)CW4e@9gZ={5eKMKJW3BuQ>aD#4Oq;hX*;0Mh-sb%^4m7)${W0io zfZK-`S;XxV%VZ{%O`4WDwbP)_qr>S~GM={*SfC?5`g{5ZIRX767Ve=;+Vs&s%aQ09 z{R=Em3wI29_GITg^aZTn6*d|@Sgx^VouPlLr3duyOxMeoC1IwqC|kQK>8{XS zq~ZzGp2Xb4rzK%}eMmEz&U^+XN#&%$5x8z)>|8ZuxkMO$j%$$&V_vqzRUf3ISUpf?V8Pp1`BC>9pjoW3z{^Gc`hg zA8vP~682=$h{LR08;zRNN!%Y{WZ4=`b=^jE0AC9t|$FBb7x!a^dcg$OSVehZh92 z`bIoCYr0Sn_v1FQqWqv#n67G6;-@~g_2XA?oU15Sj5d^ZUhWJ}UNh5STz$o|xr|)x zhkFMjBhfHU!r6Qni4?|RG1hawyQeqo)Rork9E{M1xTL;%_-&XwG870q{xf0n@xgdj zRRa{l0i1V-W>XVpuE6m|Gqbr#Qw(FOM`Wb5UP$B<$d>MOI)mFA1R+m55!Vi^YB$M6 zhLkl#W0MpF^~5@lrw88{0!Y1s$Tdeb_^$_uhurgI1iYcD1@aGt0^XK+dWeTY)xzCS zA~ zxGT`GK&{bncc4u@s{`%oSra%U&--)~=5!XD8fbctP7b|5r^Z7cq`sDW^u<5~g7@ja zAot%zxk1Q<0 zvM+%B5I?lQ_0rFlBKx9ZZ!NHQ6xo*)J6K@vF0#j$u=f_(m%*;%{fLX`ZY;8|fQ=Z! z_hx@B{7??!t0m@&MIpFWy}Q+Wqk6Aa@9Wfije4(B?;iEOMZFv9eTRDY+d{kT4_j)l z{Qo5c zz4K%)(Cjnh3RQiG{^=h5>vNbHRs2qSe!N=srRr*a7d|f22K;6LK~h}YAP>j@4T6S1 z!=O=66f_1p54r%l2)YCs2VDkT0bK<>28ux|>YqIgJ)fdl(OC1*kbgZaRvC>*Y!E*& znFiEQ`Jvqo|2uO zkl$C2S&x43x>iL3i`@$Yx6-mgPH=s*oR^g|AWI(FmG&1+*2k4YVD!0~7@91nmOt2JHdu1?>av2Q`2iK~10#s2S7(l53MNs10-g z)DAibIs`flIs$qabQB~ub%MG;$3TyOWD|V?bQ07JItA(hod)%S`aowu{UGr{{E2|h zg3e);$MM%*!Qur$!&=?`@!0Q2Y9>j0CO)MwRLv$B8HFZ1|ZN|3GZzvPo$^;e@=!+kKwb=j2!M)&vA0Ekm$bFCd zWz$6XBc+RcsetZ;e}b-mfgX(4@%VkdGS62Rc(Q4MXKZD)e2o(Y@qiM>z9;Ztfit!` z<}5!^z~6)ImTb88H8r02d}'.format(year, month)) -#set model grid -model_proj = geolib.projinfo() -model_grid = GridDesc(model_proj, x_orig=64., x_cell=0.25, x_num=324, +#Set original grid +original_proj = geolib.projinfo() +original_grid = GridDesc(original_proj, x_orig=64., x_cell=0.25, x_num=324, y_orig=15., y_cell=0.25, y_num=180) + #set target grid -target_proj = geolib.projinfo(proj='lcc', lon_0=103.5, lat_0=36.500008, lat_1=30.0, lat_2=60.0, a=6370000, b=6370000) -target_grid = GridDesc(target_proj, x_orig=-2497499.597352108, x_cell=15000.0, x_num=334, - y_orig=-2047499.8096037393, y_cell=15000.0, y_num=274) +e_we = 335 #the number of x dimension +e_sn = 275 #the number of y dimension +zdim = 7 #the number of z dimension +dx = 15000.0 #Grid spacing in x dimensions +dy = 15000.0 #Grid spacing in y dimensions +proj_name = 'lcc' #Name of the projection ('lcc' is Lambert Conformal Conic projection) +ref_lat = 36.500008 #Latitude of projection center +ref_lon = 103.5 #Longitude of projection center +truelat1 = 30.0 #First standard parallel +truelat2 = 60.0 #Second standard parallel +stand_lon = 103.5 +radius = 6370000.0 #Radius of the sphere(unit: meter) + +if (e_we-1) % 2 == 0: + xmin = (e_we-3) / 2.0 * dx + dx / 2.0 +else: + xmin = (e_we-2) / 2.0 * dx +if (e_sn-1) % 2 == 0: + ymin = (e_sn-3) / 2.0 * dy + dy / 2.0 +else: + ymin = (e_sn-2) / 2.0 * dy +target_proj = geolib.projinfo(proj=proj_name, lon_0=ref_lon, lat_0=ref_lat, lat_1=truelat1, lat_2=truelat2, a=radius, b=radius) +target_grid = GridDesc(target_proj, x_orig=-xmin, x_cell=dx, x_num=e_we-1, + y_orig=-ymin, y_cell=dy, y_num=e_sn-1) -#set the number of z dimension -zdim = 7 -#set path of the vertical allocate file -z_file = os.path.join(ge_data_dir, 'height.txt') +#Set path of the vertical allocate file +vert_prof_file = os.path.join(ge_data_dir, 'height.txt') -#set sectors that need to be processed +#Set sectors that need to be processed sectors = [SectorEnum.INDUSTRY, SectorEnum.AGRICULTURE, SectorEnum.ENERGY, SectorEnum.RESIDENTIAL, SectorEnum.TRANSPORT] @@ -57,29 +80,29 @@ #gattrs['Conventions'] = 'CF-1.6' gattrs['TITLE'] = 'Created using MeteoInfo, mechanism: {}'.format(mechanism_name.upper()) gattrs['START_DATE'] = "{}-{:0>2d}-01_00:00:00".format(year, month) -gattrs['WEST-EAST_GRID_DIMENSION'] = 335 -gattrs['SOUTH-NORTH_GRID_DIMENSION'] = 275 -gattrs['BOTTOM-TOP_GRID_DIMENSION'] = 8 -gattrs['DX'] = 15000.0 -gattrs['DY'] = 15000.0 -gattrs['WEST-EAST_PATCH_START_UNSTAG '] = 1 -gattrs['WEST-EAST_PATCH_END_UNSTAG'] = 334 +gattrs['WEST-EAST_GRID_DIMENSION'] = e_we +gattrs['SOUTH-NORTH_GRID_DIMENSION'] = e_sn +gattrs['BOTTOM-TOP_GRID_DIMENSION'] = zdim + 1 +gattrs['DX'] = dx +gattrs['DY'] = dy +gattrs['WEST-EAST_PATCH_START_UNSTAG'] = 1 +gattrs['WEST-EAST_PATCH_END_UNSTAG'] = e_we - 1 gattrs['WEST-EAST_PATCH_START_STAG'] = 1 -gattrs['WEST-EAST_PATCH_END_STAG'] = 335 +gattrs['WEST-EAST_PATCH_END_STAG'] = e_we gattrs['SOUTH-NORTH_PATCH_START_UNSTAG'] = 1 -gattrs['SOUTH-NORTH_PATCH_END_UNSTAG'] = 274 +gattrs['SOUTH-NORTH_PATCH_END_UNSTAG'] = e_sn - 1 gattrs['SOUTH-NORTH_PATCH_START_STAG'] = 1 -gattrs['SOUTH-NORTH_PATCH_END_STAG'] = 275 +gattrs['SOUTH-NORTH_PATCH_END_STAG'] = e_sn gattrs['BOTTOM-TOP_PATCH_START_UNSTAG'] = 1 -gattrs['BOTTOM-TOP_PATCH_END_UNSTAG'] = 7 +gattrs['BOTTOM-TOP_PATCH_END_UNSTAG'] = zdim gattrs['BOTTOM-TOP_PATCH_START_STAG'] = 1 -gattrs['BOTTOM-TOP_PATCH_END_STAG'] = 8 -gattrs['CEN_LAT'] = 36.500008 -gattrs['CEN_LON'] = 103.5 -gattrs['TRUELAT1'] = 30.0 -gattrs['TRUELAT2'] = 60.0 -gattrs['MOAD_CEN_LAT'] = 36.500008 -gattrs['STAND_LON'] = 103.5 +gattrs['BOTTOM-TOP_PATCH_END_STAG'] = zdim + 1 +gattrs['CEN_LAT'] = ref_lat +gattrs['CEN_LON'] = ref_lon +gattrs['TRUELAT1'] = truelat1 +gattrs['TRUELAT2'] = truelat2 +gattrs['MOAD_CEN_LAT'] = ref_lat +gattrs['STAND_LON'] = stand_lon gattrs['POLE_LAT'] = 90.0 gattrs['POLE_LON'] = 0.0 gattrs['GMT'] = 0.0 @@ -93,25 +116,24 @@ gattrs['ISURBAN'] = 13 gattrs['ISOILWATER'] = 14 -######set out species(gases and aerosol) and out species(aerosol) -################################################################################################# +##Acquisition of model species (gases and aerosol) and aerosol names based on chemical mechanisms ###------CB05,emiss_opt=15; RADM2,emiss_opt=3; SAPRC99,emiss_opt=13; MOZART,emiss_opt=10------### -################################################################################################# -out_species, out_species_aer = get_model_species_wrf(mechanism_name) +out_species, out_species_aer = get_model_species_WRFChem(mechanism_name) -#run all scripts +#Run all scripts print('Allocate according to height...') -height.run(year, month, dir_inter, model_grid, sectors, zdim, z_file) +height.run(year, month, dir_inter, original_grid, sectors, zdim, vert_prof_file) print('Merge data from different sectors...') -merge.run(year, month, dir_inter, model_grid, sectors, zdim) +merge.run(year, month, dir_inter, original_grid, sectors, zdim) print('Distribution of particulate matter and change unit...') -transform.run(year, month, dir_inter, model_grid, out_species, out_species_aer, zdim) +transform.run(year, month, dir_inter, original_grid, out_species, out_species_aer, zdim) print('Conversion projection...') -proj.run(year, month, dir_inter, model_grid, target_grid, out_species, out_species_aer, gattrs, mechanism_name, zdim) +proj.run(year, month, dir_inter, original_grid, target_grid, out_species, out_species_aer, gattrs, mechanism_name, zdim) print('-----------------------------') print('--------All finished!--------') print('-----------------------------') +#Calculate running time end_time = time.clock() -time = (end_time - start_time)/60 +time = (end_time - start_time) / 60.0 print('Time: {:.2f}min'.format(time)) \ No newline at end of file diff --git a/emips/run/run_with_script/for_WRFChem/height.py b/emips/run/run_with_script/for_WRFChem/height.py index 06fe943..80c017c 100644 --- a/emips/run/run_with_script/for_WRFChem/height.py +++ b/emips/run/run_with_script/for_WRFChem/height.py @@ -54,15 +54,17 @@ def run(year, month, dir_inter, model_grid, sectors, z, z_file): dimvars.append(dimvar) out_fn = dir_inter + '\emis_{}_{}_{}_hour_height.nc'.format(sector.name, year, month) - print('Create output data file:{}'.format(out_fn)) + print('Create output file:{}'.format(out_fn)) ncfile = dataset.addfile(out_fn, 'c', largefile=True) ncfile.nc_define(dims, gattrs, dimvars) - data = np.zeros((tdim.length, z, ydim.length, xdim.length)) - dd = np.zeros((tdim.length, z, ydim.length, xdim.length)) + data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + dd = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) # get vertical profiles scc = emis_util.get_scc(sector) vertical_pro = emips.vertical_alloc.read_file(z_file, scc) + if z != len(vertical_pro.get_ratios()): + print('----------Warning-----------\nThe z-dimensions in the vertical profile file and the output file do not match.\n----------Warning-----------') # read, merge and output if round(vertical_pro.get_ratios()[0], 2) != 1.0: print('Allocating: {}'.format(sector.name)) diff --git a/emips/run/run_with_script/for_WRFChem/merge.py b/emips/run/run_with_script/for_WRFChem/merge.py index 9e2cb6a..4e48669 100644 --- a/emips/run/run_with_script/for_WRFChem/merge.py +++ b/emips/run/run_with_script/for_WRFChem/merge.py @@ -74,8 +74,8 @@ def run(year, month, dir_inter, model_grid, sectors, z): print('Write variables data...') for sname, fns in dict_spec.iteritems(): print(sname) - spec_data = np.zeros((tdim.length, z, ydim.length, xdim.length)) - dd = np.zeros((tdim.length, z, ydim.length, xdim.length)) + spec_data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) + dd = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) for fn in fns: f = dataset.addfile(fn) dd = f[sname][:] diff --git a/emips/run/run_with_script/for_WRFChem/proj_2_file.py b/emips/run/run_with_script/for_WRFChem/proj_2_file.py index babbea4..e478e68 100644 --- a/emips/run/run_with_script/for_WRFChem/proj_2_file.py +++ b/emips/run/run_with_script/for_WRFChem/proj_2_file.py @@ -90,6 +90,7 @@ def run(year, month, dir_inter, model_grid, target_grid, out_species, out_specie #dd = dd.reshape(12, 1, ydim.length, xdim.length) #Set default values dd[dd==np.nan] = 0 + data[:, :, :, :] = dd else: print('{} no data!'.format(out_specie)) ''' @@ -101,7 +102,6 @@ def run(year, month, dir_inter, model_grid, target_grid, out_species, out_specie #Set default values dd[dd==np.nan] = 0 ''' - data[:, :, :, :] = dd ncfile.write(out_specie, data) ncfile.close() f_in.close() diff --git a/emips/run/run_with_script/for_WRFChem/transform.py b/emips/run/run_with_script/for_WRFChem/transform.py index 9f4db04..cc58263 100644 --- a/emips/run/run_with_script/for_WRFChem/transform.py +++ b/emips/run/run_with_script/for_WRFChem/transform.py @@ -59,7 +59,7 @@ def run(year, month, dir_inter, model_grid, out_species, out_species_aer, z): #add data to ncfile print('Process data and write to file...') for name in out_species: - data = np.zeros((tdim.length, z, ydim.length, xdim.length)) + data = np.zeros((tdim.length, zdim.length, ydim.length, xdim.length)) sname = name[2:] print(sname) if sname in f_in.varnames: From 335651c826bb08f0fd720a317ca9ff6a8ab7e388 Mon Sep 17 00:00:00 2001 From: chen-0126 <252941087@qq.com> Date: Sat, 24 Jun 2023 10:45:26 +0800 Subject: [PATCH 4/4] renew README --- README.md | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/README.md b/README.md index de2d0b7..6f73b53 100644 --- a/README.md +++ b/README.md @@ -32,7 +32,8 @@ Learn more about MeteoInfo and EMIPS in its official documentation at http://met Version ------- -EMIPS-1.1 was released (2023-06-24). Added support for CMAQ and bug fixes. +EMIPS-1.1 was released (2023-06-24). Added support for CMAQ and bug fixes. + EMIPS-1.0 was released (2023-01-26). Author