From fd2d9ba26b807b44573c04973cfc38221088a7bd Mon Sep 17 00:00:00 2001 From: David Jean du Preez Date: Mon, 7 Mar 2022 16:03:33 +0100 Subject: [PATCH] add missing variables to output --- .README.txt.swp | Bin 16384 -> 0 bytes dynamic/.palm_dynamic_output.py.swp | Bin 0 -> 40960 bytes dynamic/palm_dynamic_output.py | 11 +++++++++++ 3 files changed, 11 insertions(+) delete mode 100644 .README.txt.swp create mode 100644 dynamic/.palm_dynamic_output.py.swp diff --git a/.README.txt.swp b/.README.txt.swp deleted file mode 100644 index d5666534131a220441d401770251c97d4211c743..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 16384 zcmeHO-HT*L6)zKw*=RITQSilbR^6VY`}Tb7Y}O%(PG)9z7It>mo!x}EE|==Q)!mi7 zU%7Sf^voC&SkRYyM!|;!1QGuL!Ke=+{!UfhkDi%q$crf4 zhBN(D_tdH1`98O2u<+!CZGP5UWOzKl*we3ndT{=}m+tz)+Gki$q?wd~ps()qh-A>4 z>fXY_!t$M#vFA&M{p}|{M{ezx zEdnh9?+k(a+49oj0$p=*?vwn(AG-d|E^k-22($>a2($>a2($>a2($>a2($?NUn7vG zcd-|6g1g))w%qf(4n5!P-gn*eg+t%Fx9wAlK#M?&K#M?&K#M?&K#M?&K#M?&K#M?& zK#RctAOa#_>_N!2O$|Tp|Cjy$FF(lGPe3n&o(HXfPJ`|P{q_TleGhaK^d-y`zc6*J_Y*Idl>sBXbgG)^qY4x z_G{2>&}q;gG2Hk8==-4CpeI0YVj%Ku&&>w)$ZO|nUVRaF7^r7NhW~m>j%#TGR zIXy-ylspyrgy#vj3zAUeGBDq0W)>a#Gf6Uaki@wN`8X4Sl5ws*{-g*C$uA|bpi7yPNK=lHJ7ay7(OwQwqlNm{X+ zsDSH<$kMlDsdEtrA`5sRca_NLN@9Fo6k010^II~DWyq6JEhdj-K1l+6#&#WIPA0!P z3ZvhJO!fzv3j9#+Y6H-<4Cjxeh&8ugaeO614av<(2x{Hs=R~N>by|%hQ7-EcgL%*J zrnP<~^eqRyQImNOQ3)#u(4!xwfO%TvR#}{P=J0t^3i#2o7+5-);5MoJW}H#Gl*)&w2BIR19<&OQe*}XxtGH*2wGzW)e3^c zEP6|&?HWr`8ruojVU|RO#W+j`umh3jnHm(itXSvRdB}rI!B)Iuxjv0YZ+pA9vB4*k z)hJrk8d%|c{2H(fkuU9Rp0^J0!2s z*>s|Y6Yk`7vea@I*I@BMy<_<*E~0_VNb-8Mg#jNKVS*&05Ybs$cnjgJCW8b17>CeK zpLh$K5^(%+O!1BED%kS?YPoa$JgXS9X(h}UUe67VyhD)87!hM>mOSH6$NKE)nq4MG zGL!L8LR(1|s2CBJ>wTTWfpMrXiFLENCssLy578_N#jf3pS9$;xgtl?oqK;p`f$c)JmOBJO_ zW^z{JS}1-g;BSzJh43$izO5l_|EF=$!ru?ae&;EBhx_=azJDo zE2O=K+iQ^d%zD5+l(EQEQo(~4r3?fCG2B%_fjCz!Ezu$DNg@W|lCUslQt@drfh^E1 z0&?c|hpkjsRBR2anKL^PN2R7f#rUtfQl$#nMoi#PYn#mtL zG0WhZP+`XfWBU$Ax4=%MR(gJwtzW)$?!x)&SJ$pxxO|B*N_R7P_%VZi-=IC`v8`l^ zIFog`%E|X*Ih0z9Y@f=9VUlHVws7CsU4h+Er)Vb*N4;RIQp-F?^kSuCyohrZmbHnd zl499=Xx^xvR5+XatE;pk_os^HG*)};8kUm%rwFpTNXc;d0EH5=P=?9`p8~cj@m9gd zl@!(on&ksykqRh^sJsZI_Si*1VKi4FG6$qg6w)g->~b5LaaNM2On;FDhY}|Lk_1y|ujPDRX8;p-Ewkv3@9;#U( zj;u8+j0{EeHcOGPQiyvAzedL4qa-M(_JuXCGNxh$6@=;lih-q=#Dx5a%-wJYmkF#d z#e5*S2%ae@=y;5xzF)>K;^sjbs4+(dgymS{G{%P_nVi0KIfOo2`Z6)b52g7O@RqtJ zSAeK-0-Ljum{jTlO`(x}t%uswMv;=u&G_tGY3??bQ^*MTqjrgDE*|&jjlJF-+8CX& zd*gaKO%Uu*qRcRMl`MsVG!#Q{g;cs2poWFBCfR<&abT^_7Wf>ep?1cF8550Wc~ZY% z@=FvE#-WO($Jpkch|*B5!f;UV$IH-24K<ZDWlHse0W<1g6c%+ijSG6GHPgYw_Nci!e;4=#&z!Z>bcOan)OqsHNYt>spLS z)P{&8-wfpZMsIqPvC^+-trr!yJUCcyBH}m!H!qO=JQKtIRNYejAQ|>Wn)X#fI@adr z?&eUaE@cC08;u9ZK3}0(}Ma9O%oSv!IWIUPVv; zGUx{AVbB}s>0bx^9OQ!@2mJ-|z6RO@k?iL|?W0AYMW98XMc}`WK+pcymhEth6p0*( zjPh$Nr*%z)3Oz*ajH0V|l^Uj9sy?(Xkar_IQ|fA@=h@_s;)BLTqBiE#|3fI_3Z zlCf>Sz|CVPOdWz`Jf;e9fQ*dIHfd6(F}(r{`^^;70v?F~9Z52i^L^K-lI+b|z_$et zEf+(S@N`0UqC8Mp+0C}M*m*5REjcT;o86S-#UX~IC=s;`QGVHK%$8D(0icDIM4GBF zq~#bJUZl1-=N4@AzjDyu78wRYOAAX&^JoOtWR?JAjiSkHQ3*?osQ#Op;TZO;>J9_! zNRjUB#aycc$p@%}!OwI~GT)~{@X(;xo z*q-&SP_Hv8a%x+rB8$m7O`8w#0YBds&!EXzcLrq#a@mMK%uUf&cAlk;Q6Ej1RsZJ* zA9bV8ndz7teY#THO%Gh(P}ZsrEtJY7i@mO!nfhq_9E&~Ue@0mZr)9+gEJe{g^Dy(a zp@d}_Z6DdYY&BM ziWt0%O{3FI4S%!+$BJ}~EW}OeOK#*txJ6JM-HEG4X#1Nv<=v3F(xJvR_ zPcQHG9MSZ#06?k6dhEPD9xaK(*X^ub+*zwzON_i?n)#7%BlLUJxpw)`4sS55u@f7& zpqkvUOpZgFob$2Enw1r&5xr?`tpwq*h)4(K}^>YPwrtqx3fYCyE; zt(=j)N4orsw*-*z!h+SsP!-^<;Jof?-)&(;YusEcg)r%=%$nV5#LOB;sa|&G;snR7 z-Zx&r4P3#pih&T~X0k+omb*NS78knw%J!1K;%{7i8W`ZiA(l+c*AsfRWL8g>D?E%c c9D5wV-lRN|j%~~Cy0rur)hyVZp diff --git a/dynamic/.palm_dynamic_output.py.swp b/dynamic/.palm_dynamic_output.py.swp new file mode 100644 index 0000000000000000000000000000000000000000..71ec2f95fd6bf73e6686fd3220998f1544c77edb GIT binary patch literal 40960 zcmeHQeUvQ4Rc}5NLV`im0LH}14q;|@XFng$@D`Q;Zy!;^g9i^DAInSIJJY+{y!m3f zXLe>WK#X`GipCQ{$nj`AVq(M?H5zhENbn>fhKL6BB#I;&l87Kjay0&7&Ph)E-TLUR z{+yo2$^q<#-*!)T)vc;qRk!NiTUBpu{?O0tXXnf=7`RR`jLFxYwbpZ6$NtN+XBqWg zr)%4Fi#}ytO}jp`<<88{&+p0#X56d@&h5Ww-@b!$2QIv1|J+u0-K=feO>@27H9J#a4xV1cw{{&jG&zTmxJU%mD`Q z_b>wd9`H-Rb-*6rOyFr44?YGAfdOz0@EYLnV05?_=mBp6UIsi4L&dKHo4}dCcVYOr zAGiZp22KIK2xG#Z03F~8;27W$7*PHUxC!_Uc=~6=TXSd24j_bJ{yH;hChbrsZ8x5yfWmx6G z7kS{Ni5HIN>e=gcwwzX-^{qzFb{Q6e@Jt{vhJ4iQ7PIh`^2jKYJnPQ73S&yzOZ!X5o50 zB=c7(56abv^|L3ctTHi}kd?~G4a`OiqgB2ClvtPPcL_?4cp_GLuETUx4HiAnk#rLy`ZkEccn@hhX4 zT4JHJJoG`%&wvJhWx94PRD=Aj3Pktx#3*Qx9{o~D`d!~|@T@}Lc_*l~+ugd;f*8?y zW4Ij#^2LPkSe-vy)q2t@f46Y>u;|OwZ|MzsAlr626+B`g7fA&?xz4t2@}HBUqo$=( z7KNcW^dxiUmjkbbz;m+-B8h`#E!*3>a`D--E0$-uwpUT>$?9FFZ+E%asNm+rfgVoWems!kUYu2SAh;?`+xuoba7jjeD3OOUt;Bt=v|)5t`vzD`{nYuFfx z%0T{7%}z55fWZyrVu6+fGVOhyGilfgmjGN2;iuQ`h~OoLWsx0Z7Y%*A(HqQBYr~zhmr+BP!z*tXVsT^8{g#>KHQO#k zco+4q#zh)41;_l4%WSpW@ygZmjeY#-0d$8DD|5lzR1aQRWs4K;Y}#F?hOty1;AEW- zQNY!5!(Jz{U1wvH?k=*wiDq1GwY&6CMnQd&L1W>$W4JmPhxreCP_sLp=*{$@?ap@X zX1i}Ebu!Upz3xziM;~bL;oe+tv@H+$;`hqT9$d5Ys{c?hIW6ExQUOr+R@Z86P%X;XzVYokE|7;`QpL~Vb=?%H+yqU^G>8HP2?qS{SrddwgSCD!)PXfW6EhXei| zIpY!-IioJipJ|j*m4oj>6tng^0|iDQG%B2k6rb0(Fkh_Mu1g|ROc~n^3_Nu%4PEB7 zl^A70GM7avl)vcJwd$tGOb0YE$HZ?jbu`5bw~1+N%XUHDthbvMCc{B1jZzCuU8n^B zXjhGP%f!qRb6Y%Ho_aQ-(baFC85r)Aa0{L^8B#Us_Ih|`EyjfkDnw1A!;|WS0jA2C z?NW7A4Z~@6(09W`V@|#9pcX`=DJ=x$pyxRaSAFF1=B;tu z?r`-L7tied>3#cy=1q&;BW-+^i0pX2*v!-ndLbru8%vC@`lUBU>2-5Gw>t*~Y|d`= zS!d{Nwp+Ul(*HjIJ^WdK^nc|azYn_pmB6vUBo!GVGU1qTWa6dZU#ae%aSverSloUJw5u3ZVt4C;N?CVR{w-vQHb zX=`!~>XZq~n0a9T!aOW8k++>D-i}3!TiX(-iTC5v<7T`BdgA@~^tgF2fu49jK0UHJ zN^(r+9)F*coYQhV)38a22H%e=wnSQZLVB#ZBiq^r=x2H zi8NJH7@MkU3yD-!gAi0n|Nmv^@s9(f|I>}G2cY-A9cTkH!12I?(Er~D+zQ+Rlz{Ib z!dGPcZ%I4C)wq8bhJ=TK&jSm z6(B$U`M?D541D%q1pXAb3m{+q^}ri}KZo!B4qy#1fVabU{|nN8{|J2c?+5k*uLn*7 zo`KK)VcjyC zK+GlBv%4NEVf7- zK+4nb4`FS>G%K4j+sK0Q?~~yaj!5#M23(z>bsZCvRfA_s{wYNYJzENWrcgX%sXvxbfiBgczqEr_ zDJ-cpKrgX125Y68RFiO5>Mdla>XFzrW-U_IOL+3f_f6-k##A27oL z4*!+#1}Dh{m0TpJp;APWBcb6_=KkkUTEZ2rrA13gcz6tqO6Ug#F9}a$-R{`vL@iGn zNN1E_)K<8mPc@XZm17=tD0(y~({x3{=$)X4S*MH9r!?2xr3s^~DEgp9vr~sgd%NrL zVPDF(JVRG@8j+q84V=**qBG|BShPJ+RU=QeCaKT6`L-^9M1^Bx;)9JNX zY3N7gOc7E3`dMb38V!Z(ua^J*Ex&g zV+UBzCX$AVTon<|m!&EptpXCF93_`xLO+4@bKsILq<4TzI*&=`7mE$miZR#2%C8Nt zUE8#}m1xjl;Ea9|#lv5@Ql)DRE30bMJzTAO!;Xz_<%5g5s>t~iC6!VMY(jQd%y`v& z!cSFuFtle{L}NPj!D@@vX|<-s+7k*uqP9f}h*(oIkL3Ukl*tnmH>YJIM9bw1 zq}Af&fQZV2&UY#NS-Ghx_9|coE97cyKDAeCvqMX=edW(DV|eDdqV%jvxoHisfP`8j z7Nw#|(rT^6-9f%nBPQ@XQoG#$zYogoKGJca|EvA^e+7O2kAa^8UJZNzy8Z(2ZRq-6 z0DcpAHSkU7`ZoakfiFPM|3hF3I2m{c^!ou|0N;Rq|7qZzz|FuVz^TBOpw~YL{338E zPzAmXeg4b9UjTOlhk$9|KIrx9z#iaZ(CuB|V&JpT=Whkx0(=#E{Ck0WfUAHh;DgZV z>%bd;`%%s=K;>PbGVjENZj^fQ(G8W&Fg9WJUU+azjfD?c(mpF5edFRJFi~c-(;$%P zE2m{VBF8bsbb`xVuC!!jAXMW1OQFQQsQw)-KIR{)ihCIa=9~#8!rT&94$Xx9S+@~Q zBJKAJwV7ZyLD*&yEdaI^*e~whK&!@~uZ=9H);eAa0$|u5tH>_}2{5vfwJOdM00AwTU1_M5I>R3$L;b_h#Sr z&6foyDeE!=hP;Pvd<5eP6Y`NTAuD$}`ixA)O)vgPf=4!>-)oaRMN(&va3yKXIj~>( zSR(wP*r`c2+7Z>ALMLZ^1dYXD7DzOZSE@G5=KW0u(MqAsaAa9`)_Pb^Yfvox3Z<7Z zZ;0x$m^>h5Ic2Dnx69SZYTltKN6C*iC)hQ&EyXnQMpKvknUc&*o44%1R0igI`CxtOvx=b7uIPDfVN$eHrcF%b}wSIBmX~x z4tyTz!qEROz>O~2^S=a~4!j3?|3P3ka1wAk^!^p#G~nIP`!5I11nz~-ZvxA}15)=t zi2G^4lhE~l8@Li6U;k6k^B)D?3)~1afFB3`9eVz~z;(dOfv2GBKM1@7xEweYApia+ zp!auxtAG=LhoSp_2xtNafFB0F1>OInzy`1v_#X8CCxO2K{u;OgI0F0x@G#1MH{b%q zi*tcdm+A`YE*XCr+!B0u6+|VxnB{5OKts>&s^`kRvMLJQGujt@rv8jQCQ%EWzC}p4 z(oq!I+|clpF*AGnYj+yb2Y~;qW!a?c zMWyo-4`^p+cn+4F+L0<-!vAER)zW1C>nVn9ZLs7ch@mK`J2jCbdJ72+RkQA<)v>X6 zTN@-OX2?z>C;suAhDh=B6LN`cEyS}HN?_y^C4)96A1-!?_Q>f)Hth&UvqxHe;N1!h z!<>aNMX3Z)fSqd?lId^@Lsgl&huEu)6#kgJT27v7huJ~oRI}kNh_~}GiHtW7j}qrL z!=uSnLPWg^a;c=$17yb=HR&NX!5`~v^L0WLOFIdK__)K;Q^=eIS*%NNezL+uH&YCU z$92CLG{~z$$X1pre#48`4)G?Ly?zZR21K@N@{ugqlBtrU8WLMgaPjbVX*C2^I&IVsg-_CHEjb1U=_TE6p0A)6#p%c$Nu8 z1r(^)&~Th?Gf;u)o{@E7)QA`)^tf2G=5^ET4thZS8LchXIy^m&dX8jPldgG|=%so= zhDvoaR6^()k|*?g4HD zehzpg@JW>U2tZ{&oyxCYddH>b$Y5gjOzx}I1{F04b8*B@uxChY;^@f}b#$t6irgVQ zlA_qak&hx@SviuT*wd4bB6+R!v&fQ5&3E-A5j1oT(DN1@tc3gYGG9K6Rhiwc**EBu ze)Wff-6Mku9?oP_>kR~BMo{YwL}^D*lN0m(v~sZSxNBC!!l{Kz?56ppb2i?}n!Dz3 zSeeBQq&F>yS99)Uc~)3lDtsAW5U;B9A!1%2j_}}uOQ_gLWHuK)pjnkCR!Mqs-C40Mofwgf z$5Wh($Lq;VMDig{!J>R55mQ}0YKYK?%^wboZ>qsWgKv@8kQO*iYQ%yih{|Y~em+t- zDl671qojKINF~X#lBSIM1zPG=Ssw?Uf|=ATnUDrh(qlP^y-&^w6oQ5#r|!Jwm6sElO7wb?n+)xN4tm6<84D*pj*KCE zGdSuQei=j+o*5Wbk1hkrjd4^a5(b-5QIJr4KpYh*A)OMA{pL>SxHnjFKA#yHoX3KS zY&*1~D5kv=2ZDp|LdW{S+JQOnL^WjQ<_^$UrL<~vhH{52ExBnT$`?MNE*3 zI)6icXdETlmHz(~=)HB)f1&?hD0N_3`~OEk51@1ZpM<`DGjKX^8t^&j{2nj?JPMuv zR^S$326zZM{~f@7;NPI{KMByje>%(mKcVlx54a4t6!;!={=0z{;H%L0-w*sQa2wzN zbHMSyalq5i{T~Ky0>~b467VV50PX_%z>fm|4gLQy;3EL71DFSX40slHfZKshU_Wpc z@Opsk1CPNL@M}O#+66uiUJyTO#FHJkR6ABpBjbb-p@|2p3PT5th>xA-*ylmekc`7J z)KW!My=I``AM+B8&Ym@*F+ZhF6GI%%;ZPYMXG6N1Ze0JHmFhLa3MW_-?3x#xxj^oC zI{I3ixr=2A{9H~XiIcf~xuhye9TDSypvg9O(MQglU0pnk_j<5dzHtbhuzdAgbFm!y zl(lyqDIk)c0q*lUkTdliN=%KdQ>Y3|5E{?-2}r( zxb11%(1PfE6InPnnxRTB6vYWF9^7k;w#OarewUC5CJb0|@@a+->re0RkTSG+F_JTO zKARUh5j>Ds{T#5|X;H7|83|m{Szd7-6y?k76o-6~rSnV6aQQ5bEKj0ubjdkkO#Fzi zarJwPPNoW(2DnknP~p`*I$_xoHtbN$sD$T6u!wd#LWhqh!5b!s$VX9_YDeXehsv2pURm0( zmRx2ry~&qHs+XKz2qHMLC-h@IDB2IW`Bt~x@_!LPTGc2Eze$5}jvU52 zk+cuzh&#HwDDyNGMK9(!{iTon$s_uUo=lz&!z%`?9Mv4D!z7t|Q>tJjJ>j;?M{%4E zlaJy!9VQ>eF**zp9NDCLBRdTD4(KgO?=Wa&_~%%1q{8bVQvxe5gr1Z>K;=s(-npUz zK&68bDj&&spOb!>9L0|HBjL`C3`ab#n-vqoCI37^v{#^q7|0efzJaM0jC0V{@)$ID*)OXa2xdhtAJVH5$OIO2HpjDz!ktQ z;FHk(?*!@qo%`oO-@gGk1Na8Y%*(w4I(|Xo`Q!i^cySP;??;Xv8K_%aVUjNnV&wWJ zj*<*k<@nyD;vhzOl4fxbV`x~$ys>N2k%GlRjF7c?SFZY8^DP_^17=DWE#OR+MOcVOzAV#dbNjnFSueT`NKmJUql-S}*{ODQV0Big&nj6y$ literal 0 HcmV?d00001 diff --git a/dynamic/palm_dynamic_output.py b/dynamic/palm_dynamic_output.py index 8c08818..54fa349 100644 --- a/dynamic/palm_dynamic_output.py +++ b/dynamic/palm_dynamic_output.py @@ -83,6 +83,17 @@ def palm_dynamic_output(wrf_files, interp_files, dynamic_driver_file, times_sec, outfile.createDimension('Dmid', sum(nbin)) outfile.createDimension('composition_index', len(listspec)) outfile.createDimension('max_string_length', 25) + # variables + val_dmid = outfile.createVariable('Dmid', "f4",('Dmid',)) + val_dmid.setncattr('units', 'm') + _dmid, _bim = palm_dynamic_aerosol.define_bins(nbin, reglim) + _dmid = np.array(_dmid) + val_dmid[:] = _dmid[:] + + val_idx = outfile.createVariable('composition_index', "i4", ('composition_index',)) + for n in range(0, len(listspec)): + val_idx[n] = n+1 + # 2D vertical profile variables _val_composition_name = outfile.createVariable('composition_name',"S1", ("composition_index", "max_string_length")) _val_composition_name.setncattr('long_name', "aerosol composition name") -- GitLab