From bf022aaad0d5b3402b653e0a7bd2808d425cdb16 Mon Sep 17 00:00:00 2001 From: David Jean du Preez Date: Mon, 9 May 2022 12:05:34 +0200 Subject: [PATCH] fix aerosol concen. variables --- configurations/mbees_les_test_10.conf | 2 +- dynamic/.palm_dynamic_output.py.swp | Bin 0 -> 53248 bytes .../palm_dynamic_aerosol.cpython-38.pyc | Bin 6015 -> 6864 bytes .../palm_dynamic_output.cpython-38.pyc | Bin 11989 -> 12138 bytes .../__pycache__/palm_wrf_utils.cpython-38.pyc | Bin 17463 -> 17463 bytes dynamic/palm_dynamic_aerosol.py | 164 +++++++++++------- dynamic/palm_dynamic_output.py | 36 ++-- 7 files changed, 124 insertions(+), 78 deletions(-) create mode 100644 dynamic/.palm_dynamic_output.py.swp diff --git a/configurations/mbees_les_test_10.conf b/configurations/mbees_les_test_10.conf index e1387f4..c95765d 100644 --- a/configurations/mbees_les_test_10.conf +++ b/configurations/mbees_les_test_10.conf @@ -52,7 +52,7 @@ nbin = [1,7] reglim = [3.9e-8, 5.0e-8, 2.5e-6] nf2a = 1.0 # Either SO4 or OC must be included in listspec -listspec = ['DU', 'SO4'] +listspec = ['DU', 'SO4', 'OC', 'NH','NO', 'BC'] #listspec = ['DU','SO4', 'OC','BC','SS','NH','NO'] # radiation diff --git a/dynamic/.palm_dynamic_output.py.swp b/dynamic/.palm_dynamic_output.py.swp new file mode 100644 index 0000000000000000000000000000000000000000..6d364c98281a8aaffcfbba3945cb82b22282ffc7 GIT binary patch literal 53248 zcmeI53zTF9^GBQTXPWDw*z*%i_+E2;RqVM5PA@43QVygXNI8&lAmu>H zfs_L&2T~5C97s8ka^NB1K%tn)+)T=|q8PLOZ}s-CqKJPA{1&(xJQF+veC4MznR~%? zU==Kb%fJQT8DJMU0et>(naoGQt-uGf;HxO%Zw6O@Gr?294^Zr{0#5*+!^rUG;0?eB zF9iF*)4?wARPY511s?;i1~u?2U<#Z79s%w`QGXA(8e9bCz)9d<6#U!4?}KG<^UyaVLHJeUJtLixWL+ypw{GH@w)GI#*K zz6&(LuL0rn&w?*wAGi;^7rffo3Ow^E7dP{jQq5cR_?xSg{g&S-txR~6-kBL2XWT2- zRx0gc$-{5LTPygbdfWF3rDol)SA1`^S+9EHR>8|nd?xH`xKXT@E2Wv0W~tCB?JrkL zHNRZ1<+A&$FNqgk%Ca@li>#d2-U681U+9v7PB!g8hLXD3vWgMb8mAzb;yiS3O~+ zUT@M7>9f9$EpJ2{W7wJVbfcx&qf<_IT2f(iM~~{aMP^6{Ld4rHH*$Q%=c?qZSN!sJ zDX(rPbZ_mIPF5F2l96vFXF)j#*3!lAR}e!SmeoqFef!VfJF~ygLY%d7*7d8{EO$!H zyvVI7Z)}53M)s?P*5cT~^A5`!W1cRR5n200wu#W%uN}#}{ffW>;w^z>jQM0@rMN1r z#of%4G(+JjJfHLmkoWW3%e%tzEoC>k$sUaaZKS* zW)MYFZswO}XQl!Xpd=But$IU8ANTNJY>-_V`f8=Uxm#KrSH)gW8oD1=uWB6H7cm%k zNT$%L)_p{IQ?UslThhIPO49AA&nq++Gw$VO2v-ys6NRmA_=I;@zsgd$hw$aK? zWna|cPaDV;*_c%VziPjP`8Qg?j+LHKhrE#>zz_k&(-~axrmNW$Eq;;woCcd zO1;n`U+S8heg>bJS=%Wk8bu+VPEra$gDNx&wKYj5BI2^97ZkT9)srnvWi$HDX7nA~ zaW^N6@?=xLL)Q(Zpkms#^kxq(ZEwgo`6AU|zAGw(ZMIE-$_fEn9SKkq4B8?Ntqb1v z1_ZZh-gAnJJXf+-r=kc`$A-LoVQK+g`NkwbC~9(274|a@J#W1v73vk%Dk9aZbb}2~ zI@DUJSt@cT+4c6)_BwiYNCWgTLSfviw`eq;GsdWJa;C&OIgOD~IP=zI*jm z4+@^daHuFEMkZfTJJhVNlzd<0G84;Evr(@wniZ9Z^}SZz%4ScCdozKQ{KY@bLNTv1 z(}3obC-=IIj*R>(xTJW(Kl2CAp7GiKGPp+Et}q%deCqqkJ#4?DtF( z2?-0B!5C^KX86q!Ot@8-R#j7o`8JD0P=g#4G>7b6!l*T+nC(k-Ol30VYJ`9ktBozM;L|N_kkpO+UiUNE zY&O%P!}>C6Q{ghPkF)^jPm5}xo5D=y@Ol{qhLP7yh_o)S(JV9hQSHp|!?YBcH^`8! z-dJYlQ1t3GCHJ-HFRC1j@0{)+l*Tjs^DbU#mK!ZH%XCVbo@w87hf6Oye|q2Z&N~n` zZ*g{ywecAtvg7&krgu%#3q@^P@RZlDgk;7~+uin?yO+zg-KAQ`YizaF>$N=@(f@Bl zPkldntmyxC{rbJ=`Zt1K2j4}Xe+M`RJQAD)?n9sd1b7`Nfg<=-unYVvdi`I3E5R>; z`_b(`1KQvcumDa4??a!z7IeUc;H&8G{{U_PzYAUgM8{tM|A`KND`~8Z=+mNw!f6}rJs}oDF+_z9FXDGXpJ)~mAYTb=^2R|-Ym(a zDzwj3eA!>Lf=tXaz0DswFgMGzDi+$P5;_tsHDykqro%_4r{>}b^mO>>^wbg`TUvr{TE%H z)ChQVY^@-Xrfmu%Q?+d&k*aMFBU804B9W?X79&$N?IIwW)HsGC8gB-N4g^HA@n?-{ zNK>~F8@GQ$F~S4AW;l~mA;IKim7tY+vsqecXwo+)dnz?MH zE@u6?P^$5qmx2tab$}y*9b1z1H}d5g*1=U?v_@^2)8~0Oh80MZa;CG+!)iNWZf;4LyP) zxC8y&*zRIE3$1tKY85MFdATONJ|mf6d5h*KhL0o4o~4asG25Y)D)h>jk}TS-qX<)( zViXg*i7A~*L=uuuz~b2SNI&%uk3uh)J!NdTy{hVVW=9(`Yg%eynONYz*pD~F$U4M( zYWR!18JFT4OS6h>XhB&rmOQbMNpeFZqsGIbl99xQp{T~fp(w%q9f;TYpivIsubno9l9~b@K`tN@Pz5g2UO3(uH;AC(s zy8oYomxGsq95@kt8NL5La0_@A_&R$39pDD=8c+onfJtyKX}t-&5lFgslm17*hrvIg z>u-bSg7d+R==z7j+2B<0-{||dfDeH`0Z#`vqw~KCWWl%5`R@R?gFgl@1s8%n;3?o9 z>;SidcY|v{6}$lK1vg*^xf;mY{r#-bZ-Nc*P1fiC4!jAR1@;5Wf4S0Int7IG+))=r zLuLRXx&z}ohPPr_M6oXx(<*j!EU|?ZC7dxINtaykmOV@~{#Yn3d!-N?*WC=suq<7; zUFtnThh&|$qE{fp>3|~wvXiWr;BGMKW%H8tyky3Z5k~6xr_JmTV+v;jC>QpXHvvQU6U1- zrkb4&okpzZ#52^VKg1^r^(6ZRCCTnx1 zi73G>aywrL8iDF-wbl88clM$;H)EEe6Prks=xjiOhOP)K#t1~}M!!6$y-2KK2E_!@ zI?KJqVZ8cPhPnqa+CcJ{)o)M@SkaxYf~^GAR3y<@1th$3j9ls(eG*7N2rjv!#W8V7 zr;REkf+*yK4VFpk5S{1%EHs@ce3NhEF|t$P)9yzqt`>_!BpK_o;dEk@-rH_$F#WNL zn+SW!1m2paKP%l;xT4yB*tbKj#IoJ0udElExp>r2g2bN$3^Gm? z&B9i0N$gi@>Uwml*xG88sFdu*^R7iSh%r-W$Y7MqUbye#BaX>5DZ;`25e`0YPhSGh zJ9y!M2=fP9$%Z(ZYc6V+GI^k zr`o*+P1$%NI};-FGH>>pcDm(K^L-g$%gQ%4oAAs!%i5e38Gm)_&YQ0pc_bedAbA!K zAZe^BAZ*|27ENFo$HpF_1mBBttwfBPpL08JNJI)K3A?SqOoqH{)oEciVjot>T*6K- ztc!KTj+IbDQMT^48cr_T)c~M8e`c8Ro-ZaaJ*Zfv|Gy5s_%OP&=>NOS>fl?^^&bh| zgKobGHo!gT^)Cliuph`;|2NR*Mc2O*l)#tJ=Whmw!D--rboy_DkAZiC7Pt(Y4Nd~G z=KqgC*8E=&I^c3}C%XRI!Eb_}1LC)T9e5nLo4kGiNS;4W-ai6PlKiKilmjUTQVtw% z9MHqx5PzB1#4VBTTM3UtoqgvFl@CFZG3#s0NZ@9upBm@J)4b|4xDPh3?VN2?bV&Ip zGq!&Eis(kEk<_V1Rg0;Sw4d!YX0=rIpwMoq7S^JM%j4Hb8Zp&q^JK<0CR$QzEM-Gr z;1Z%tK%(J{ZY{;^J3x@722;Jr78}>P4eWuzXuv9Br1{EA`DNX*dn@v(!4xhP9t4+u z{0ydcM~O9vO(^~UF!N7qbbp2WIKj`xfAp7(m173+fzaN|dzK1@)3HE|VftR7jKO3A1u0^ju z0A$_%HLTA+A7sHhSf4L|vq1*j&KmtSU=chGe2TK%41OD24lV;yKEVSdoe{eGi3LY!1=bhiag&k2(B}s`PG5(g)pkw$0kCyLE#vl1P3^B&Sp(}}6bQ@$+3<(=1tt;19VRX!DNwj1o6-Ik3Fd0W+sv``r zUZ0McGEM~w62De&iG4HJWEn26#$V`@b?s6)A6nkx@3K5383X zy&uc$h)PQ~k8cOu35RxZs*`1~E%~GTHIh=bVZ*fr8~dRz=S-Vri)Ic_R@b|Al=7O^ zL`$*@wzA~WyD{0a+v`Kmg9lx9>53cKhk21ZL$+^I-%{qz!@bO^$Zhe#vaK=c>Xf+M zr^FApuHW74l@{|8B;tlG@9`+&w&m=QVa7A3^*`KY+=%skJcu}5e($;7u;o2|1)-{J@&A7TI`NCpjYa=I*XYN0p!feCI3LLV z|4#$)1?Yf9@PFw3p9R-}XM-6a=li`4ycEd!e&59oa2I$pxEeeQ`~vtQHi16?F9MU` ziQxO#1-=Jv0&fI|!2e<&_yl+tkbVDqfSlv^7H}9$ffK=(uoc`3-VB}(&H%SyGdKip z$6l}ro(t~5R`9n#><0_rvEUSNld&IUxjq0N-VTHx7lRQ#ZV&JEy;SM-M8aLaXU`na zGkA&})Bd<^cJC;WK)ug~?J8^}>6ZnyVUB-l+L$|q5J zXQ;ndbi5eKW(mDTyT!>#K5oy{TJKDyz9HIR!rmH}ICLl!f2~C6$w3(|RBW_`^?AR( zSFN59f-n1et-dF`<;bP=$ZdV4)^Q)#SNx9aA?{|Vza(byf+Nt{rNmsiIW%WUiD$aa zvNcotbEQl$K_1oJMp6!?*5qYWaaG&&+s)O&3P(+u69DWzXM_2CFVS0+D;#cR_^THiFoM%ejWL^%dBWL<`o7ePFQe!L)*@Jd8Dd!-m1#fIBPJbfm{U~gz zt;lht!S+W9az9-HZH;1d)-rA!`5uqgmU!*K@j4RE+HIYR2(7aD&RPQ36VvNYr=3Kn zy+&`5JSbNJIKlBwIn9%v2VVt;T6NcpeiL%>nik4TbvZ**RdcqPVCFUFgI$Otu^7|T zc}WXib~5sAcbengTB(jZY`wA0VRzzK9jP4=KWvV7U0G#EvD-8zV%B2A$An@Sk(7^~ zL)n}h*l2^Z&(76{=<2=72^)CW#k)UpR;9xzgtb_$BgdJA<8j|zLH%WE*ez@miXWyD z+a;XX{Oy3NU3P=JS>vxVQ!N#0Igyf!T1G~yiL#kNg6dWqUwYe7RpSrnU#YeXD}R_E zjOk4fyvySdIF|F1_WDZWxYBJK5z;gMdZAH5l*ihKnwv>(MMXBFW>^oFHxW`QeaQli z%FQL&n(NM5B!MI|24>|e~s8s5-xUzLf9$glyS8UvAo^u1dJV)`om3n z3v`}@WBU(WP@7We_W z{k7nQU?2E-a6h{JC&7n61&Gi8$I zj#+T2$@gxG%F%8hMSbp*9NaLNqS|0QkfJ{GX*h~v3>bX0gv3&X9+<~Y*Rdi>f;iH< zBSwpGoBvW)ox>}KLvzAF!XPSKiSxFjYb#1?UG86t<>CyR8pRM-&j4n9LfmCwuv z++vP$h)9PGVFL~kag!P56cJT1urhIkbsD&wBEosZL$G3$4l^9ZQ98_U6i4YW!%-Ze z!wA8Fd~6TwFztRFCShhsC3kh0c0z}VcTANI^xQ_P&)>Fvb8xF@^V1caJs#be3`cQ< z{*pwntBTbrRn*JmqB^F8(=VFMe5Wg#84U(U(>$E+PtCNqLDY2H^1K1m%uG;rcF@`U z*TuV&Wi@HUwSF{A-8;`?p zTn$4iX$BA{#jj}%Ln=v@jXNH|fy%dF);SnsrCEAf8lV++^^dLq+so{=YnI$+)zEEK z?gTu!mD&{K8eP$3y4c5=%)U2D&SFh%K$EvAbmvSBq&uIt&a{d^TcSayz7~ z|GxyC{WCz;|55kV=L>wf1zZp0{68NY2A{(Qa08I@|0cn!u>rgaEPxY%(D^11KY-r> zawfoD@Hf~2t_QooUt~o1u>l+erwy_N zbXA_1qK*e|SnUwg8SDXMjETpMtJpxNC#S$kwsYGe#E=>Fn-Du@5Ta(wNo(cuzSuN- zpZ#W2NH#gk)1K@7BU@L68ZY%(bCMXrkG%C<(~?+xGU{*akCO0etfccK>_EvTZ$Rl? z{DA~}?J{xMqmP(3#${>U;o0M{RSGFi&m_A}+?LqMSsy?nw#9_Jimd6qe?^}@Ea!gr z*n1L5kC};GK8YXu8K1^P_J`!M+u5m!sbPH>hbZ~TIk9Ucp*uIiYm_X7tn9z=xFZbX zL?`We9BiiSjwycwmn=uF=1BG0-ALlsdrbLD65S*5&-4Pj zX-NMfs_@Nvm5-FIz(;JTH7}7N@);?^j#rmtasRxU-s2&mmcHk*fQ_r_Y#Lbh*o%O zv*%9U-1C3!`Nc9Z)n}d$bMkL*$~$$S583I{6&w+Pz*O^bIquw?|JdDq)%n%oN$2J; zjEoe+eg>V|EAbRF)nUE%jGwt_6;{Z|WBXCIl$<`QS5E46yH4@|52mW=U(pL2jq7Ew}EgvLS#ateKcaSK5Jz??tW?_kw3=8=*9v~A-sYDc^>2g7 zkAc*y`nWl^yKRM}qZUW3Tl$yK?9U9aL9?oU;lOe73T<3Awy*U6e`YQGovf3K{y%Lh z@YCr09|X?=r-K_=>t6uhXN~_(@J?_#_z%|heQ*Ta$=d!0z#i~qa2;#=Rq*eu?cW8K z!Fk|x@GaK&?~yfr&;m~fr-A!e>whnJEw}`n1;)SwtozIF{JjR80&c@l@g^W=1KvsA zRsK!h(@V;Mlmo{T2dq3zx<>Xk9v8=I?t`xuui`Ori@ zq8+{WDU-vU#N!k=X(iWYSS;Lu%Z)hcs!gXsNa}=R85LbhcGq#tuN|dMICWWdN}X`B z$?ia@6OQRhN-kP&wmwoZnT?@OfOWtmL=L1JR#ZpEzofBXXY;3gMp7pne#*MslI7M1 zUE!${4*p1^`oqNLBKe;TC!DbNc;Zyb2Kk79}*A`m*Z1in@rGb8K&FGs)sGa&ju zuJG#fJsW23(@V;MlmjUTQVygXNI8&lAmu>Hfs_L&2T~5C97s9vBjA7-8+Yml>%v`B Gng0jsr!=ep literal 0 HcmV?d00001 diff --git a/dynamic/__pycache__/palm_dynamic_aerosol.cpython-38.pyc b/dynamic/__pycache__/palm_dynamic_aerosol.cpython-38.pyc index c9f5a0074b33719da08c1f656b4b6652362d7051..ee15406aaa8acad3e5e8511334bdd9353953911b 100644 GIT binary patch delta 3289 zcma)8O^h5z74EA3nVz5Soxh#gU+)GRCu3t{JBe9{vypK~Fc8N9g5!yu$?R0mdS+*O z*3~_m^-`@tStFD~q>UN@2NbPO;3IMn;()|~Bar|h4kZKbG%`lv4I832F3NGUlE?T!6mqZAHhth zwTaLenba6%Iy2tYBPwW{G7~hQjdV4#Kr)n+0LfI64U!e=+lH{(W@rhrNt*N?i4tLg z5!7uC==0hZy?Hvag#jv?X-xY*9b(e}^Q15zsZW7hnNUjlnZYgDx%7 zuoXl^WDHH$Ry6)6mR=$J-@ZXKeE5@>^l89(H@d>Mf6Zwffo%_tp`Dq*BiCEG&W#D zn*G{aK9lB}dL%%uu@;v~sIM(cbCox1>oR`DSyr-*vLXldmk0Cwkp95E6mt{@kB?daP@@ei;H_QN0Sb&pMd$!y`)n2&RW{W} z*bK0342`y-ZR~BtM`Zu0GCCd}+GwQ7?QUI!c9GrZ;R9eQwu`@v+QlIcc=Q1S2XNwJ z|IY?0=#T={9@e69R)|Wd3k#SorM()JRhEk)Ey|Na>SH}jhza0$Q4|!*&0Tq*6}4zm zOp1&c7v&}_h4GXqiK(&{m4MVUAUm*pd_;2V_78;vBp<(Wu+A^_{00ZwYCM5pBP0>h z2pI$i!3DTaq`u%%9SY;fE*z!X@3$y$P z!W@81Gq2h5z1mXC4|dOqqu!dwJM~p*1w9^mOb+zG3#HlenYYg8(f%=n&mbH__$LNpy;-=buA36~CXGw=g68^YQO;X9gJcstrM@1*=}8YWeEP zFQ7W1)cx*(aLK7rAIAoMZJ5X}_xe7o^E-FJF2MY0@ctQn3L*&?piQ2DxZ31AL~s%A z;nQ|t;6j@}M^cn~55AZGm=2x;YyNG7N0H`LoKVcUh|+fut|2H+yoSbb<nME{;W7gLqw%jJ zd;{Sn1gtH-gwOz}CLdOCgELgB^1lkwJbxL-uOKKod=n*o7tHvb(h>5r_~)e;Ey|mq z|HY3JltOWA)46%oDRw@EUnMzwGnIGsdRQ9qKj_ia-bv$vzA(PuIc|fI`eXb{{gh7m zrFh1C`}DruVNN>gv!XVlQQn6iUaac1`k9Bv z+6y=yOX}FH`XKSH8U`U(DT*tJIa9?hzLt2d_`e5H>HJUp>%{DUD!JlF+6z?ElcwMA zuHI2GkcHKHr(0up{Cc<5s6pFkww5^_S$$e4j&cu&7|vm(sI#bY+1pb|>QlzP(Cx8) e$6G`vK}iD^+fd_hQ9EltX`i=`+9&KO+xjnXhtenj delta 2415 zcma)7&2Jl35Z|}nU2pucjvaqw<2VjYzg%ePM_SSkT0zqfkZ22SDK$yAu6LbmVz1q| zYd4LaWr+wOfhwh1B)A}w14un_$eA$q^GKs<6_?b~FK=hun*T{=Iu z&VzjYIXs0vkGC3Lz;4z&anYc4+Y@QYnlap`=`qC| zAAlFA8WzR6cmdOL=m;*veAad-q-!|tH)=lSMK|1PwB(p|_eb^I*pBTJn~_4I8_Fl$ z;i23vLU)71vIwyGi~@n1o`=6(3ApF@R7Ygu9IuNirxTTHOFkz774)efP_j?SK=nC+ z6@m55iH^j`5~=X3q#bgCIzibIKglig)6d#rhLW((>D7G?@;6$A5_zy-K@j{dBN=Lvmk352r-od{FrqL_dl$#P4i z-LI2&-10N1dtcmWLlWECI&R(36&7=! z1hSnSLHR*~CWN;8|9FDLg>I!G`2^Ov%VDLbS$lvHHUQp|EUCxm`DNZ71}d>d4Rp51 zR*jFqwC;FY1^*aqXKmS&H7!EnNu{~J&(rfJzKZQZ2vJvR5tvvE9zby}`M<=;9 zt_`taHd5hL&gXM1pW)g$2>Kbzu`Fm0ZzS>C{ku$DNv(1Olgn0km}BY{O5ls~2=#(j zF5Eajd2##+@uZ6Tw=k3}#!bs~w31V++l^U+8d|gHj)+59a~6!&4}u}ZZ`3@E>?g$XuCMgNXL8(PJz8n!1_P20htNBhw{juwdrv{IfmYB)QQ0oJm@(WZM;5Dr=>pum~$fv8g?<^rG1EUqPTX z-K0gID79vaVuMz+g`&GCO}P1RK1Ro2_oD1OE7vV|F`OIkPV3&tRyM}pHCo4bE`f!8 z4}jy7Fn&Y$B=`jBCjv>3G#Lge0+dYX&9!*sEB=i=U{3cU z>_d1J;V{Axgrf+@06d{;R2+&EA1)LN{`TeDgy^e zq7&%qB*L2rrx5(3yWn1pjY=1hzvMP!hY!C6WKnLFx))MWviysQzJubRFK*Vhm&W}y zmJfdht2Tzq@$U0jD#xFIVw{D4_bmghKuL~5k?73Yk^VpV&T&wt?;~JG(;bz#g|Y&| zYX~18+(rl>6cIi|m_nFFxC2m7daX$D>hb%KP9n9{$K9%QPNNX}cw24!vRnXJAx-b% zQV9VMnfetfBem77^fIzn5UwIzL%>F&I)VXE?AvUqGOeHxt4rSl7^Uwbc^%#3M+mF^0t!Az9(9g*mXq2ioO@C;WYGzr3?yZptR6Y=wE1|z&4z)P3Oi2V+>~4pHlj3p-|fG_U<2^ z%HlNO!wlvudxofq$sTm+i;2;g_~e6$4>Qp|nfc(K850v@Of=EoIn;t^l;nQ%xxfGC z{CfHP()X9ECQYUa3IA?>eRI|sdsyWpWFsFU8}g($Sx?Y_Cd_p z3nA!%UW{^47e*awF;aSTmfsHYus1P0+L_NCIYvMJ}>Yo>+0 za2WeMy=H`=yGqVOJ)+7muR`6_0}@ED$Xl9_2XRWwII<(dOT8lFsOZ}a#}?&vMLi5Q z>fN22ZqvA*pnnEd?9|;L{WD0jZrIk!wHmEXYuQwCKGE(XmD$8HePQTSBe+DRs8h777x+Tj^dodFe+;E%hu!Yk8v5H?W#7cnhuCnNb3f* z8tgQPo7fR?;qQINcJ%Em9jb@+?C8rm3}=dGa{P5B8YphkNPD&?aAW;-RHzjtoto8T=oI_8ThPON#p(gx*Y>UcznfjHsXQ6Vf z3?{(|`|s*rHHpt?YAc3MBR-NTT*?%CP-e5A=TuxvAMM*d0X$FpB5mxOG8-E*D)#?~ zp;cgKPs;Y$r~ju`#MrS?nhP$_x<$jfQ9XhaIpHE)Do*Ew%MdB?3`9#j3$YRd#7mrj zInFR&tawjUOu|AbVi8hXJczB9AYF1~;7W*`{CAdqVt5z-cV^#XuG;tIWe@wwXJNnA28opS)VXAxG5b54XmFD^^SOrmvV)QZ z9dz!Kay5e7NM=DEBnz!!<1JP;))XaP_OK~Je5|JVK+uoACMsiQ0I69}7m`g-H&TnB zAd;P0X%lFi$nN&kXi-xBDDz$BRK{2v54(scR#+ji}wK}m!Jbk9fA%b zbqYGfMp~XIgQLq`?4Z4wU9?{zA@-@gmHlkr%)jAyM9K>KSSpi^Guq~DJAv0zFeKB7 zWGb$n!x!~5%Cs?=PHGX{f;!4RX=^$#!ab+B9O814%NWYE1%pzNR5p^xs4?82$cb?G z43}u$=6piPI4f@-8Ja*#L6=>M7Rqm^i3o3*Ef|+0srdDHDyvTNgee|V9R-La(^_1; zG8Ii-i)U*IF7aD-uRWKyxSDljKkIf68D@Dxj9qu%AaVAiyK&fD{P~=B=*5=YsE)C7K40~Do`0JQze;tIt@|2}oZ&dmh4>x@x0>Rc<5k^SahB(rQTFi zVGCU{6S9+-?PY__hdubR^v%U(i9Re@_U5wLx6P+%mL-dC?%z4bTARh#3Ngtj5@l^e+*=!Q%ulwfmd}02pdN)F*Vi>K9ld|%pU|u(aS2d}PD`-PRf{Ikj z3WDwyL38q)4?9_$a}W5-+yGwiKwv|sxv*=!34#zJUTH}7K`-=?xa9Ih=#B>HhcFD# zDi8M1st9c5XIL~Z_O7?V;08hh`9kWhQ^kcA>^+NMJB*NDq>i-OsM^%_4f&*~`Lv#b zd7%Mz&=cP%7-96fr1`O*`|V_Y*oXZffbhQfOq1dW?U(KttJ)F30ls4w@4ybb7sWM6 z*$H+WAThZq&mxFM=D~q|tpX0b)N!0T?kV>W#6jw~7aY6?#u3<8;UK)pxkn2Yksb zj&TGgxgDFHye2Cjf?W$|EGklP;{UN~`j!akXyMoRnc24!>vAA~=F=D~LQWxO7KghWu@4^G~R?3v~qp z90fPLc3t;@nP2Bw6=CzUYtJw1ES=yi88g})#~77oIEV*d^aXQxFzO_l(e6I}>Ka@R z4swlpy~2mqxM+=yG}M^$mEPb4wk_({49YFqRVTa!C)pG@pMujBj>4G=$6&s~aexXZ zAjufcmL31$F6SUssaSx;XFNiWO2bkmk%4T5&%=cZE08O*YeS&c!i?Zki$B$LpkhPq z9T63ZKbZ4Uq#&)!m3;5i755lGDs$Eesij0Fm&#^xPBrTsznqCJrQ*(UCAFMTidjnu zAw8L}4jD>vJe|!Yir1_|NdG1M_sZQS-!y+MS`fKx+l!3FZ*3k?SCf;-Z>>IbtJu_b zL+n%q=w6|JDHtv2z}F^djhuQYSvdL0-7Qy9fQ+_BQILGrzRee+xh+_v^}UqZIrUMp zbLywm!6{6T{jP6-uAO`xq2%JUl~O0CK}v2;L*#c_JxteKe7%j5htunnyqvaE>gF^; z9&|jGc1&H_O1^Wn%g*Ik`XYJgxPXSqtkWg4WrAFEeqMi~ZgOm9Cc10v;MmmeJ;fH+ zeIzUU>2S!Oo&L@V3MG9ilT4)($}C+dQ&dhEQ<;<+U5+WrAu{UT(tMC*CYfw!a)8M+ z$-3KW&oGQJnJ*IeV}y>7eAfYVlziXSp+CmrbL5Y%W_>N9pW_iRc{I_8imRS>9deVK zzMc9ws{!)JcWFFHv80ckKcC1zsV=7~^y919__>%;s*fd%c!bH@OfEB-VRD7ZRVJ9pQ6_wl8l!iZ z= 1.0): - _var_b = 0 + _var_b[:,:,:,:] = 0 elif (nf2a < 1.0): - _var_b = 1-_var_a + _var_b[:,:,:,:] = 1-_var_a + # aerosol concen.# _var_con = outfile.createVariable('ls_forcing_'+ side + '_aerosol', "f4", ("time", _dim1, _dim2, "Dmid"), fill_value=fillvalue_float) _var_con.setncattr('lod', 1) _var_con.setncattr('long_name',"boundary condition of aerosol concentration") - _var_con = palm_dynamic_aerosol.aerosolConWrfchemBoundary(dimensions, _dim1, _dim2, interp_files, side, nbin, reglim, wrfchem_bin_limits) - + _var_con[:,:,:,:] = palm_dynamic_aerosol.aerosolConWrfchemBoundary(dimensions, _dim1, _dim2, interp_files, side, nbin, reglim, wrfchem_bin_limits) + #--------------------------------------------------------------------------- # create dynamical & chemical variables in output def add_interpDim(dynam_chem_variables): - print('Adding dimensions to Dynamic Driver') + print('Adding dynamical variables') # surface pressure _val_surface_forcing_surface_pressure = outfile.createVariable('surface_forcing_surface_pressure', "f4",("time")) # geostrophic wind @@ -259,13 +267,15 @@ def palm_dynamic_output(wrf_files, interp_files, dynamic_driver_file, times_sec, _val_ls_forcing.setncattr('units', 'kg/m3') else: _val_ls_forcing.setncattr('units', 'ppm') - outfile.close() + #outfile.close() # function - create dynamical & chemical variables add_interpDim(dynam_chem_variables) + + outfile.close() # read interpolated files and write values for dynamical & chemical variables def add_interpValues(dynam_chem_variables): - print('Adding initializing variable values to Dynamic Driver') + print('Adding initializing variables to Dynamic Driver') infile = netCDF4.Dataset(interp_files[0], "r", format="NETCDF4") outfile = netCDF4.Dataset(dynamic_driver_file, "r+", format="NETCDF4") # initialization variables @@ -296,7 +306,7 @@ def palm_dynamic_output(wrf_files, interp_files, dynamic_driver_file, times_sec, # time dependent variables - dynamical & chemical variables if not nested_domain: - print('\nAdding time dependent variables value to Dynamic Driver') + print('\nAdding time dependent variables to Dynamic Driver') outfile = netCDF4.Dataset(dynamic_driver_file, "r+", format="NETCDF4") for ts in range(0, len(interp_files)): # geostrophic wind -- GitLab