1&&\"string\"==typeof v&&!d.checkClone&&je.test(v))return e.each((function(i){var o=e.eq(i);y&&(t[0]=v.call(this,i,o.html())),Me(o,t,n,r)}));if(p&&(a=(i=be(t,e[0].ownerDocument,!1,e,r)).firstChild,1===i.childNodes.length&&(i=a),a||r)){for(u=(s=w.map(ye(i,\"script\"),He)).length;f0&&me(a,!u&&ye(e,\"script\")),s},cleanData:function(e){for(var t,n,r,i=w.event.special,o=0;void 0!==(n=e[o]);o++)if(V(n)){if(t=n[Y.expando]){if(t.events)for(r in t.events)i[r]?w.event.remove(n,r):w.removeEvent(n,r,t.handle);n[Y.expando]=void 0}n[Q.expando]&&(n[Q.expando]=void 0)}}}),w.fn.extend({detach:function(e){return Ie(this,e,!0)},remove:function(e){return Ie(this,e)},text:function(e){return $(this,(function(e){return void 0===e?w.text(this):this.empty().each((function(){1!==this.nodeType&&11!==this.nodeType&&9!==this.nodeType||(this.textContent=e)}))}),null,e,arguments.length)},append:function(){return Me(this,arguments,(function(e){1!==this.nodeType&&11!==this.nodeType&&9!==this.nodeType||Le(this,e).appendChild(e)}))},prepend:function(){return Me(this,arguments,(function(e){if(1===this.nodeType||11===this.nodeType||9===this.nodeType){var t=Le(this,e);t.insertBefore(e,t.firstChild)}}))},before:function(){return Me(this,arguments,(function(e){this.parentNode&&this.parentNode.insertBefore(e,this)}))},after:function(){return Me(this,arguments,(function(e){this.parentNode&&this.parentNode.insertBefore(e,this.nextSibling)}))},empty:function(){for(var e,t=0;null!=(e=this[t]);t++)1===e.nodeType&&(w.cleanData(ye(e,!1)),e.textContent=\"\");return this},clone:function(e,t){return e=null!=e&&e,t=null==t?e:t,this.map((function(){return w.clone(this,e,t)}))},html:function(e){return $(this,(function(e){var t=this[0]||{},n=0,r=this.length;if(void 0===e&&1===t.nodeType)return t.innerHTML;if(\"string\"==typeof e&&!De.test(e)&&!ve[(he.exec(e)||[\"\",\"\"])[1].toLowerCase()]){e=w.htmlPrefilter(e);try{for(;n3,re.removeChild(t)),s}}))}();var Ue=[\"Webkit\",\"Moz\",\"ms\"],Xe=v.createElement(\"div\").style,Ve={};function Ge(e){var t=w.cssProps[e]||Ve[e];return t||(e in Xe?e:Ve[e]=function(e){for(var t=e[0].toUpperCase()+e.slice(1),n=Ue.length;n--;)if((e=Ue[n]+t)in Xe)return e}(e)||e)}var Ye=/^(none|table(?!-c[ea]).+)/,Qe=/^--/,Je={position:\"absolute\",visibility:\"hidden\",display:\"block\"},Ke={letterSpacing:\"0\",fontWeight:\"400\"};function Ze(e,t,n){var r=te.exec(t);return r?Math.max(0,r[2]-(n||0))+(r[3]||\"px\"):t}function et(e,t,n,r,i,o){var a=\"width\"===t?1:0,s=0,u=0;if(n===(r?\"border\":\"content\"))return 0;for(;a<4;a+=2)\"margin\"===n&&(u+=w.css(e,n+ne[a],!0,i)),r?(\"content\"===n&&(u-=w.css(e,\"padding\"+ne[a],!0,i)),\"margin\"!==n&&(u-=w.css(e,\"border\"+ne[a]+\"Width\",!0,i))):(u+=w.css(e,\"padding\"+ne[a],!0,i),\"padding\"!==n?u+=w.css(e,\"border\"+ne[a]+\"Width\",!0,i):s+=w.css(e,\"border\"+ne[a]+\"Width\",!0,i));return!r&&o>=0&&(u+=Math.max(0,Math.ceil(e[\"offset\"+t[0].toUpperCase()+t.slice(1)]-o-u-s-.5))||0),u}function tt(e,t,n){var r=Fe(e),i=(!d.boxSizingReliable()||n)&&\"border-box\"===w.css(e,\"boxSizing\",!1,r),o=i,a=_e(e,t,r),s=\"offset\"+t[0].toUpperCase()+t.slice(1);if(We.test(a)){if(!n)return a;a=\"auto\"}return(!d.boxSizingReliable()&&i||!d.reliableTrDimensions()&&A(e,\"tr\")||\"auto\"===a||!parseFloat(a)&&\"inline\"===w.css(e,\"display\",!1,r))&&e.getClientRects().length&&(i=\"border-box\"===w.css(e,\"boxSizing\",!1,r),(o=s in e)&&(a=e[s])),(a=parseFloat(a)||0)+et(e,t,n||(i?\"border\":\"content\"),o,r,a)+\"px\"}function nt(e,t,n,r,i){return new nt.prototype.init(e,t,n,r,i)}w.extend({cssHooks:{opacity:{get:function(e,t){if(t){var n=_e(e,\"opacity\");return\"\"===n?\"1\":n}}}},cssNumber:{animationIterationCount:!0,columnCount:!0,fillOpacity:!0,flexGrow:!0,flexShrink:!0,fontWeight:!0,gridArea:!0,gridColumn:!0,gridColumnEnd:!0,gridColumnStart:!0,gridRow:!0,gridRowEnd:!0,gridRowStart:!0,lineHeight:!0,opacity:!0,order:!0,orphans:!0,widows:!0,zIndex:!0,zoom:!0},cssProps:{},style:function(e,t,n,r){if(e&&3!==e.nodeType&&8!==e.nodeType&&e.style){var i,o,a,s=X(t),u=Qe.test(t),l=e.style;if(u||(t=Ge(s)),a=w.cssHooks[t]||w.cssHooks[s],void 0===n)return a&&\"get\"in a&&void 0!==(i=a.get(e,!1,r))?i:l[t];\"string\"===(o=typeof n)&&(i=te.exec(n))&&i[1]&&(n=se(e,t,i),o=\"number\"),null!=n&&n==n&&(\"number\"!==o||u||(n+=i&&i[3]||(w.cssNumber[s]?\"\":\"px\")),d.clearCloneStyle||\"\"!==n||0!==t.indexOf(\"background\")||(l[t]=\"inherit\"),a&&\"set\"in a&&void 0===(n=a.set(e,n,r))||(u?l.setProperty(t,n):l[t]=n))}},css:function(e,t,n,r){var i,o,a,s=X(t);return Qe.test(t)||(t=Ge(s)),(a=w.cssHooks[t]||w.cssHooks[s])&&\"get\"in a&&(i=a.get(e,!0,n)),void 0===i&&(i=_e(e,t,r)),\"normal\"===i&&t in Ke&&(i=Ke[t]),\"\"===n||n?(o=parseFloat(i),!0===n||isFinite(o)?o||0:i):i}}),w.each([\"height\",\"width\"],(function(e,t){w.cssHooks[t]={get:function(e,n,r){if(n)return!Ye.test(w.css(e,\"display\"))||e.getClientRects().length&&e.getBoundingClientRect().width?tt(e,t,r):Be(e,Je,(function(){return tt(e,t,r)}))},set:function(e,n,r){var i,o=Fe(e),a=!d.scrollboxSize()&&\"absolute\"===o.position,s=(a||r)&&\"border-box\"===w.css(e,\"boxSizing\",!1,o),u=r?et(e,t,r,s,o):0;return s&&a&&(u-=Math.ceil(e[\"offset\"+t[0].toUpperCase()+t.slice(1)]-parseFloat(o[t])-et(e,t,\"border\",!1,o)-.5)),u&&(i=te.exec(n))&&\"px\"!==(i[3]||\"px\")&&(e.style[t]=n,n=w.css(e,t)),Ze(0,n,u)}}})),w.cssHooks.marginLeft=ze(d.reliableMarginLeft,(function(e,t){if(t)return(parseFloat(_e(e,\"marginLeft\"))||e.getBoundingClientRect().left-Be(e,{marginLeft:0},(function(){return e.getBoundingClientRect().left})))+\"px\"})),w.each({margin:\"\",padding:\"\",border:\"Width\"},(function(e,t){w.cssHooks[e+t]={expand:function(n){for(var r=0,i={},o=\"string\"==typeof n?n.split(\" \"):[n];r<4;r++)i[e+ne[r]+t]=o[r]||o[r-2]||o[0];return i}},\"margin\"!==e&&(w.cssHooks[e+t].set=Ze)})),w.fn.extend({css:function(e,t){return $(this,(function(e,t,n){var r,i,o={},a=0;if(Array.isArray(t)){for(r=Fe(e),i=t.length;a1)}}),w.Tween=nt,nt.prototype={constructor:nt,init:function(e,t,n,r,i,o){this.elem=e,this.prop=n,this.easing=i||w.easing._default,this.options=t,this.start=this.now=this.cur(),this.end=r,this.unit=o||(w.cssNumber[n]?\"\":\"px\")},cur:function(){var e=nt.propHooks[this.prop];return e&&e.get?e.get(this):nt.propHooks._default.get(this)},run:function(e){var t,n=nt.propHooks[this.prop];return this.options.duration?this.pos=t=w.easing[this.easing](e,this.options.duration*e,0,1,this.options.duration):this.pos=t=e,this.now=(this.end-this.start)*t+this.start,this.options.step&&this.options.step.call(this.elem,this.now,this),n&&n.set?n.set(this):nt.propHooks._default.set(this),this}},nt.prototype.init.prototype=nt.prototype,nt.propHooks={_default:{get:function(e){var t;return 1!==e.elem.nodeType||null!=e.elem[e.prop]&&null==e.elem.style[e.prop]?e.elem[e.prop]:(t=w.css(e.elem,e.prop,\"\"))&&\"auto\"!==t?t:0},set:function(e){w.fx.step[e.prop]?w.fx.step[e.prop](e):1!==e.elem.nodeType||!w.cssHooks[e.prop]&&null==e.elem.style[Ge(e.prop)]?e.elem[e.prop]=e.now:w.style(e.elem,e.prop,e.now+e.unit)}}},nt.propHooks.scrollTop=nt.propHooks.scrollLeft={set:function(e){e.elem.nodeType&&e.elem.parentNode&&(e.elem[e.prop]=e.now)}},w.easing={linear:function(e){return e},swing:function(e){return.5-Math.cos(e*Math.PI)/2},_default:\"swing\"},w.fx=nt.prototype.init,w.fx.step={};var rt,it,ot=/^(?:toggle|show|hide)$/,at=/queueHooks$/;function st(){it&&(!1===v.hidden&&e.requestAnimationFrame?e.requestAnimationFrame(st):e.setTimeout(st,w.fx.interval),w.fx.tick())}function ut(){return e.setTimeout((function(){rt=void 0})),rt=Date.now()}function lt(e,t){var n,r=0,i={height:e};for(t=t?1:0;r<4;r+=2-t)i[\"margin\"+(n=ne[r])]=i[\"padding\"+n]=e;return t&&(i.opacity=i.width=e),i}function ct(e,t,n){for(var r,i=(ft.tweeners[t]||[]).concat(ft.tweeners[\"*\"]),o=0,a=i.length;o1)},removeAttr:function(e){return this.each((function(){w.removeAttr(this,e)}))}}),w.extend({attr:function(e,t,n){var r,i,o=e.nodeType;if(3!==o&&8!==o&&2!==o)return void 0===e.getAttribute?w.prop(e,t,n):(1===o&&w.isXMLDoc(e)||(i=w.attrHooks[t.toLowerCase()]||(w.expr.match.bool.test(t)?pt:void 0)),void 0!==n?null===n?void w.removeAttr(e,t):i&&\"set\"in i&&void 0!==(r=i.set(e,n,t))?r:(e.setAttribute(t,n+\"\"),n):i&&\"get\"in i&&null!==(r=i.get(e,t))?r:null==(r=w.find.attr(e,t))?void 0:r)},attrHooks:{type:{set:function(e,t){if(!d.radioValue&&\"radio\"===t&&A(e,\"input\")){var n=e.value;return e.setAttribute(\"type\",t),n&&(e.value=n),t}}}},removeAttr:function(e,t){var n,r=0,i=t&&t.match(P);if(i&&1===e.nodeType)for(;n=i[r++];)e.removeAttribute(n)}}),pt={set:function(e,t,n){return!1===t?w.removeAttr(e,n):e.setAttribute(n,n),n}},w.each(w.expr.match.bool.source.match(/\\w+/g),(function(e,t){var n=dt[t]||w.find.attr;dt[t]=function(e,t,r){var i,o,a=t.toLowerCase();return r||(o=dt[a],dt[a]=i,i=null!=n(e,t,r)?a:null,dt[a]=o),i}}));var ht=/^(?:input|select|textarea|button)$/i,gt=/^(?:a|area)$/i;function vt(e){return(e.match(P)||[]).join(\" \")}function yt(e){return e.getAttribute&&e.getAttribute(\"class\")||\"\"}function mt(e){return Array.isArray(e)?e:\"string\"==typeof e&&e.match(P)||[]}w.fn.extend({prop:function(e,t){return $(this,w.prop,e,t,arguments.length>1)},removeProp:function(e){return this.each((function(){delete this[w.propFix[e]||e]}))}}),w.extend({prop:function(e,t,n){var r,i,o=e.nodeType;if(3!==o&&8!==o&&2!==o)return 1===o&&w.isXMLDoc(e)||(t=w.propFix[t]||t,i=w.propHooks[t]),void 0!==n?i&&\"set\"in i&&void 0!==(r=i.set(e,n,t))?r:e[t]=n:i&&\"get\"in i&&null!==(r=i.get(e,t))?r:e[t]},propHooks:{tabIndex:{get:function(e){var t=w.find.attr(e,\"tabindex\");return t?parseInt(t,10):ht.test(e.nodeName)||gt.test(e.nodeName)&&e.href?0:-1}}},propFix:{for:\"htmlFor\",class:\"className\"}}),d.optSelected||(w.propHooks.selected={get:function(e){var t=e.parentNode;return t&&t.parentNode&&t.parentNode.selectedIndex,null},set:function(e){var t=e.parentNode;t&&(t.selectedIndex,t.parentNode&&t.parentNode.selectedIndex)}}),w.each([\"tabIndex\",\"readOnly\",\"maxLength\",\"cellSpacing\",\"cellPadding\",\"rowSpan\",\"colSpan\",\"useMap\",\"frameBorder\",\"contentEditable\"],(function(){w.propFix[this.toLowerCase()]=this})),w.fn.extend({addClass:function(e){var t,n,r,i,o,a,s,u=0;if(h(e))return this.each((function(t){w(this).addClass(e.call(this,t,yt(this)))}));if((t=mt(e)).length)for(;n=this[u++];)if(i=yt(n),r=1===n.nodeType&&\" \"+vt(i)+\" \"){for(a=0;o=t[a++];)r.indexOf(\" \"+o+\" \")<0&&(r+=o+\" \");i!==(s=vt(r))&&n.setAttribute(\"class\",s)}return this},removeClass:function(e){var t,n,r,i,o,a,s,u=0;if(h(e))return this.each((function(t){w(this).removeClass(e.call(this,t,yt(this)))}));if(!arguments.length)return this.attr(\"class\",\"\");if((t=mt(e)).length)for(;n=this[u++];)if(i=yt(n),r=1===n.nodeType&&\" \"+vt(i)+\" \"){for(a=0;o=t[a++];)for(;r.indexOf(\" \"+o+\" \")>-1;)r=r.replace(\" \"+o+\" \",\" \");i!==(s=vt(r))&&n.setAttribute(\"class\",s)}return this},toggleClass:function(e,t){var n=typeof e,r=\"string\"===n||Array.isArray(e);return\"boolean\"==typeof t&&r?t?this.addClass(e):this.removeClass(e):h(e)?this.each((function(n){w(this).toggleClass(e.call(this,n,yt(this),t),t)})):this.each((function(){var t,i,o,a;if(r)for(i=0,o=w(this),a=mt(e);t=a[i++];)o.hasClass(t)?o.removeClass(t):o.addClass(t);else void 0!==e&&\"boolean\"!==n||((t=yt(this))&&Y.set(this,\"__className__\",t),this.setAttribute&&this.setAttribute(\"class\",t||!1===e?\"\":Y.get(this,\"__className__\")||\"\"))}))},hasClass:function(e){var t,n,r=0;for(t=\" \"+e+\" \";n=this[r++];)if(1===n.nodeType&&(\" \"+vt(yt(n))+\" \").indexOf(t)>-1)return!0;return!1}});var xt=/\\r/g;w.fn.extend({val:function(e){var t,n,r,i=this[0];return arguments.length?(r=h(e),this.each((function(n){var i;1===this.nodeType&&(null==(i=r?e.call(this,n,w(this).val()):e)?i=\"\":\"number\"==typeof i?i+=\"\":Array.isArray(i)&&(i=w.map(i,(function(e){return null==e?\"\":e+\"\"}))),(t=w.valHooks[this.type]||w.valHooks[this.nodeName.toLowerCase()])&&\"set\"in t&&void 0!==t.set(this,i,\"value\")||(this.value=i))}))):i?(t=w.valHooks[i.type]||w.valHooks[i.nodeName.toLowerCase()])&&\"get\"in t&&void 0!==(n=t.get(i,\"value\"))?n:\"string\"==typeof(n=i.value)?n.replace(xt,\"\"):null==n?\"\":n:void 0}}),w.extend({valHooks:{option:{get:function(e){var t=w.find.attr(e,\"value\");return null!=t?t:vt(w.text(e))}},select:{get:function(e){var t,n,r,i=e.options,o=e.selectedIndex,a=\"select-one\"===e.type,s=a?null:[],u=a?o+1:i.length;for(r=o<0?u:a?o:0;r-1)&&(n=!0);return n||(e.selectedIndex=-1),o}}}}),w.each([\"radio\",\"checkbox\"],(function(){w.valHooks[this]={set:function(e,t){if(Array.isArray(t))return e.checked=w.inArray(w(e).val(),t)>-1}},d.checkOn||(w.valHooks[this].get=function(e){return null===e.getAttribute(\"value\")?\"on\":e.value})})),d.focusin=\"onfocusin\"in e;var bt=/^(?:focusinfocus|focusoutblur)$/,wt=function(e){e.stopPropagation()};w.extend(w.event,{trigger:function(t,n,r,i){var o,a,s,u,l,f,p,d,y=[r||v],m=c.call(t,\"type\")?t.type:t,x=c.call(t,\"namespace\")?t.namespace.split(\".\"):[];if(a=d=s=r=r||v,3!==r.nodeType&&8!==r.nodeType&&!bt.test(m+w.event.triggered)&&(m.indexOf(\".\")>-1&&(x=m.split(\".\"),m=x.shift(),x.sort()),l=m.indexOf(\":\")<0&&\"on\"+m,(t=t[w.expando]?t:new w.Event(m,\"object\"==typeof t&&t)).isTrigger=i?2:3,t.namespace=x.join(\".\"),t.rnamespace=t.namespace?new RegExp(\"(^|\\\\.)\"+x.join(\"\\\\.(?:.*\\\\.|)\")+\"(\\\\.|$)\"):null,t.result=void 0,t.target||(t.target=r),n=null==n?[t]:w.makeArray(n,[t]),p=w.event.special[m]||{},i||!p.trigger||!1!==p.trigger.apply(r,n))){if(!i&&!p.noBubble&&!g(r)){for(u=p.delegateType||m,bt.test(u+m)||(a=a.parentNode);a;a=a.parentNode)y.push(a),s=a;s===(r.ownerDocument||v)&&y.push(s.defaultView||s.parentWindow||e)}for(o=0;(a=y[o++])&&!t.isPropagationStopped();)d=a,t.type=o>1?u:p.bindType||m,(f=(Y.get(a,\"events\")||Object.create(null))[t.type]&&Y.get(a,\"handle\"))&&f.apply(a,n),(f=l&&a[l])&&f.apply&&V(a)&&(t.result=f.apply(a,n),!1===t.result&&t.preventDefault());return t.type=m,i||t.isDefaultPrevented()||p._default&&!1!==p._default.apply(y.pop(),n)||!V(r)||l&&h(r[m])&&!g(r)&&((s=r[l])&&(r[l]=null),w.event.triggered=m,t.isPropagationStopped()&&d.addEventListener(m,wt),r[m](),t.isPropagationStopped()&&d.removeEventListener(m,wt),w.event.triggered=void 0,s&&(r[l]=s)),t.result}},simulate:function(e,t,n){var r=w.extend(new w.Event,n,{type:e,isSimulated:!0});w.event.trigger(r,null,t)}}),w.fn.extend({trigger:function(e,t){return this.each((function(){w.event.trigger(e,t,this)}))},triggerHandler:function(e,t){var n=this[0];if(n)return w.event.trigger(e,t,n,!0)}}),d.focusin||w.each({focus:\"focusin\",blur:\"focusout\"},(function(e,t){var n=function(e){w.event.simulate(t,e.target,w.event.fix(e))};w.event.special[t]={setup:function(){var r=this.ownerDocument||this.document||this,i=Y.access(r,t);i||r.addEventListener(e,n,!0),Y.access(r,t,(i||0)+1)},teardown:function(){var r=this.ownerDocument||this.document||this,i=Y.access(r,t)-1;i?Y.access(r,t,i):(r.removeEventListener(e,n,!0),Y.remove(r,t))}}}));var Tt=e.location,Ct={guid:Date.now()},Et=/\\?/;w.parseXML=function(t){var n;if(!t||\"string\"!=typeof t)return null;try{n=(new e.DOMParser).parseFromString(t,\"text/xml\")}catch(e){n=void 0}return n&&!n.getElementsByTagName(\"parsererror\").length||w.error(\"Invalid XML: \"+t),n};var St=/\\[\\]$/,kt=/\\r?\\n/g,At=/^(?:submit|button|image|reset|file)$/i,Nt=/^(?:input|select|textarea|keygen)/i;function Dt(e,t,n,r){var i;if(Array.isArray(t))w.each(t,(function(t,i){n||St.test(e)?r(e,i):Dt(e+\"[\"+(\"object\"==typeof i&&null!=i?t:\"\")+\"]\",i,n,r)}));else if(n||\"object\"!==x(t))r(e,t);else for(i in t)Dt(e+\"[\"+i+\"]\",t[i],n,r)}w.param=function(e,t){var n,r=[],i=function(e,t){var n=h(t)?t():t;r[r.length]=encodeURIComponent(e)+\"=\"+encodeURIComponent(null==n?\"\":n)};if(null==e)return\"\";if(Array.isArray(e)||e.jquery&&!w.isPlainObject(e))w.each(e,(function(){i(this.name,this.value)}));else for(n in e)Dt(n,e[n],t,i);return r.join(\"&\")},w.fn.extend({serialize:function(){return w.param(this.serializeArray())},serializeArray:function(){return this.map((function(){var e=w.prop(this,\"elements\");return e?w.makeArray(e):this})).filter((function(){var e=this.type;return this.name&&!w(this).is(\":disabled\")&&Nt.test(this.nodeName)&&!At.test(e)&&(this.checked||!de.test(e))})).map((function(e,t){var n=w(this).val();return null==n?null:Array.isArray(n)?w.map(n,(function(e){return{name:t.name,value:e.replace(kt,\"\\r\\n\")}})):{name:t.name,value:n.replace(kt,\"\\r\\n\")}})).get()}});var jt=/%20/g,qt=/#.*$/,Lt=/([?&])_=[^&]*/,Ht=/^(.*?):[ \\t]*([^\\r\\n]*)$/gm,Ot=/^(?:GET|HEAD)$/,Pt=/^\\/\\//,Rt={},Mt={},It=\"*/\".concat(\"*\"),Wt=v.createElement(\"a\");function Ft(e){return function(t,n){\"string\"!=typeof t&&(n=t,t=\"*\");var r,i=0,o=t.toLowerCase().match(P)||[];if(h(n))for(;r=o[i++];)\"+\"===r[0]?(r=r.slice(1)||\"*\",(e[r]=e[r]||[]).unshift(n)):(e[r]=e[r]||[]).push(n)}}function Bt(e,t,n,r){var i={},o=e===Mt;function a(s){var u;return i[s]=!0,w.each(e[s]||[],(function(e,s){var l=s(t,n,r);return\"string\"!=typeof l||o||i[l]?o?!(u=l):void 0:(t.dataTypes.unshift(l),a(l),!1)})),u}return a(t.dataTypes[0])||!i[\"*\"]&&a(\"*\")}function $t(e,t){var n,r,i=w.ajaxSettings.flatOptions||{};for(n in t)void 0!==t[n]&&((i[n]?e:r||(r={}))[n]=t[n]);return r&&w.extend(!0,e,r),e}Wt.href=Tt.href,w.extend({active:0,lastModified:{},etag:{},ajaxSettings:{url:Tt.href,type:\"GET\",isLocal:/^(?:about|app|app-storage|.+-extension|file|res|widget):$/.test(Tt.protocol),global:!0,processData:!0,async:!0,contentType:\"application/x-www-form-urlencoded; charset=UTF-8\",accepts:{\"*\":It,text:\"text/plain\",html:\"text/html\",xml:\"application/xml, text/xml\",json:\"application/json, text/javascript\"},contents:{xml:/\\bxml\\b/,html:/\\bhtml/,json:/\\bjson\\b/},responseFields:{xml:\"responseXML\",text:\"responseText\",json:\"responseJSON\"},converters:{\"* text\":String,\"text html\":!0,\"text json\":JSON.parse,\"text xml\":w.parseXML},flatOptions:{url:!0,context:!0}},ajaxSetup:function(e,t){return t?$t($t(e,w.ajaxSettings),t):$t(w.ajaxSettings,e)},ajaxPrefilter:Ft(Rt),ajaxTransport:Ft(Mt),ajax:function(t,n){\"object\"==typeof t&&(n=t,t=void 0),n=n||{};var r,i,o,a,s,u,l,c,f,p,d=w.ajaxSetup({},n),h=d.context||d,g=d.context&&(h.nodeType||h.jquery)?w(h):w.event,y=w.Deferred(),m=w.Callbacks(\"once memory\"),x=d.statusCode||{},b={},T={},C=\"canceled\",E={readyState:0,getResponseHeader:function(e){var t;if(l){if(!a)for(a={};t=Ht.exec(o);)a[t[1].toLowerCase()+\" \"]=(a[t[1].toLowerCase()+\" \"]||[]).concat(t[2]);t=a[e.toLowerCase()+\" \"]}return null==t?null:t.join(\", \")},getAllResponseHeaders:function(){return l?o:null},setRequestHeader:function(e,t){return null==l&&(e=T[e.toLowerCase()]=T[e.toLowerCase()]||e,b[e]=t),this},overrideMimeType:function(e){return null==l&&(d.mimeType=e),this},statusCode:function(e){var t;if(e)if(l)E.always(e[E.status]);else for(t in e)x[t]=[x[t],e[t]];return this},abort:function(e){var t=e||C;return r&&r.abort(t),S(0,t),this}};if(y.promise(E),d.url=((t||d.url||Tt.href)+\"\").replace(Pt,Tt.protocol+\"//\"),d.type=n.method||n.type||d.method||d.type,d.dataTypes=(d.dataType||\"*\").toLowerCase().match(P)||[\"\"],null==d.crossDomain){u=v.createElement(\"a\");try{u.href=d.url,u.href=u.href,d.crossDomain=Wt.protocol+\"//\"+Wt.host!=u.protocol+\"//\"+u.host}catch(e){d.crossDomain=!0}}if(d.data&&d.processData&&\"string\"!=typeof d.data&&(d.data=w.param(d.data,d.traditional)),Bt(Rt,d,n,E),l)return E;for(f in(c=w.event&&d.global)&&0==w.active++&&w.event.trigger(\"ajaxStart\"),d.type=d.type.toUpperCase(),d.hasContent=!Ot.test(d.type),i=d.url.replace(qt,\"\"),d.hasContent?d.data&&d.processData&&0===(d.contentType||\"\").indexOf(\"application/x-www-form-urlencoded\")&&(d.data=d.data.replace(jt,\"+\")):(p=d.url.slice(i.length),d.data&&(d.processData||\"string\"==typeof d.data)&&(i+=(Et.test(i)?\"&\":\"?\")+d.data,delete d.data),!1===d.cache&&(i=i.replace(Lt,\"$1\"),p=(Et.test(i)?\"&\":\"?\")+\"_=\"+Ct.guid+++p),d.url=i+p),d.ifModified&&(w.lastModified[i]&&E.setRequestHeader(\"If-Modified-Since\",w.lastModified[i]),w.etag[i]&&E.setRequestHeader(\"If-None-Match\",w.etag[i])),(d.data&&d.hasContent&&!1!==d.contentType||n.contentType)&&E.setRequestHeader(\"Content-Type\",d.contentType),E.setRequestHeader(\"Accept\",d.dataTypes[0]&&d.accepts[d.dataTypes[0]]?d.accepts[d.dataTypes[0]]+(\"*\"!==d.dataTypes[0]?\", \"+It+\"; q=0.01\":\"\"):d.accepts[\"*\"]),d.headers)E.setRequestHeader(f,d.headers[f]);if(d.beforeSend&&(!1===d.beforeSend.call(h,E,d)||l))return E.abort();if(C=\"abort\",m.add(d.complete),E.done(d.success),E.fail(d.error),r=Bt(Mt,d,n,E)){if(E.readyState=1,c&&g.trigger(\"ajaxSend\",[E,d]),l)return E;d.async&&d.timeout>0&&(s=e.setTimeout((function(){E.abort(\"timeout\")}),d.timeout));try{l=!1,r.send(b,S)}catch(e){if(l)throw e;S(-1,e)}}else S(-1,\"No Transport\");function S(t,n,a,u){var f,p,v,b,T,C=n;l||(l=!0,s&&e.clearTimeout(s),r=void 0,o=u||\"\",E.readyState=t>0?4:0,f=t>=200&&t<300||304===t,a&&(b=function(e,t,n){for(var r,i,o,a,s=e.contents,u=e.dataTypes;\"*\"===u[0];)u.shift(),void 0===r&&(r=e.mimeType||t.getResponseHeader(\"Content-Type\"));if(r)for(i in s)if(s[i]&&s[i].test(r)){u.unshift(i);break}if(u[0]in n)o=u[0];else{for(i in n){if(!u[0]||e.converters[i+\" \"+u[0]]){o=i;break}a||(a=i)}o=o||a}if(o)return o!==u[0]&&u.unshift(o),n[o]}(d,E,a)),!f&&w.inArray(\"script\",d.dataTypes)>-1&&(d.converters[\"text script\"]=function(){}),b=function(e,t,n,r){var i,o,a,s,u,l={},c=e.dataTypes.slice();if(c[1])for(a in e.converters)l[a.toLowerCase()]=e.converters[a];for(o=c.shift();o;)if(e.responseFields[o]&&(n[e.responseFields[o]]=t),!u&&r&&e.dataFilter&&(t=e.dataFilter(t,e.dataType)),u=o,o=c.shift())if(\"*\"===o)o=u;else if(\"*\"!==u&&u!==o){if(!(a=l[u+\" \"+o]||l[\"* \"+o]))for(i in l)if((s=i.split(\" \"))[1]===o&&(a=l[u+\" \"+s[0]]||l[\"* \"+s[0]])){!0===a?a=l[i]:!0!==l[i]&&(o=s[0],c.unshift(s[1]));break}if(!0!==a)if(a&&e.throws)t=a(t);else try{t=a(t)}catch(e){return{state:\"parsererror\",error:a?e:\"No conversion from \"+u+\" to \"+o}}}return{state:\"success\",data:t}}(d,b,E,f),f?(d.ifModified&&((T=E.getResponseHeader(\"Last-Modified\"))&&(w.lastModified[i]=T),(T=E.getResponseHeader(\"etag\"))&&(w.etag[i]=T)),204===t||\"HEAD\"===d.type?C=\"nocontent\":304===t?C=\"notmodified\":(C=b.state,p=b.data,f=!(v=b.error))):(v=C,!t&&C||(C=\"error\",t<0&&(t=0))),E.status=t,E.statusText=(n||C)+\"\",f?y.resolveWith(h,[p,C,E]):y.rejectWith(h,[E,C,v]),E.statusCode(x),x=void 0,c&&g.trigger(f?\"ajaxSuccess\":\"ajaxError\",[E,d,f?p:v]),m.fireWith(h,[E,C]),c&&(g.trigger(\"ajaxComplete\",[E,d]),--w.active||w.event.trigger(\"ajaxStop\")))}return E},getJSON:function(e,t,n){return w.get(e,t,n,\"json\")},getScript:function(e,t){return w.get(e,void 0,t,\"script\")}}),w.each([\"get\",\"post\"],(function(e,t){w[t]=function(e,n,r,i){return h(n)&&(i=i||r,r=n,n=void 0),w.ajax(w.extend({url:e,type:t,dataType:i,data:n,success:r},w.isPlainObject(e)&&e))}})),w.ajaxPrefilter((function(e){var t;for(t in e.headers)\"content-type\"===t.toLowerCase()&&(e.contentType=e.headers[t]||\"\")})),w._evalUrl=function(e,t,n){return w.ajax({url:e,type:\"GET\",dataType:\"script\",cache:!0,async:!1,global:!1,converters:{\"text script\":function(){}},dataFilter:function(e){w.globalEval(e,t,n)}})},w.fn.extend({wrapAll:function(e){var t;return this[0]&&(h(e)&&(e=e.call(this[0])),t=w(e,this[0].ownerDocument).eq(0).clone(!0),this[0].parentNode&&t.insertBefore(this[0]),t.map((function(){for(var e=this;e.firstElementChild;)e=e.firstElementChild;return e})).append(this)),this},wrapInner:function(e){return h(e)?this.each((function(t){w(this).wrapInner(e.call(this,t))})):this.each((function(){var t=w(this),n=t.contents();n.length?n.wrapAll(e):t.append(e)}))},wrap:function(e){var t=h(e);return this.each((function(n){w(this).wrapAll(t?e.call(this,n):e)}))},unwrap:function(e){return this.parent(e).not(\"body\").each((function(){w(this).replaceWith(this.childNodes)})),this}}),w.expr.pseudos.hidden=function(e){return!w.expr.pseudos.visible(e)},w.expr.pseudos.visible=function(e){return!!(e.offsetWidth||e.offsetHeight||e.getClientRects().length)},w.ajaxSettings.xhr=function(){try{return new e.XMLHttpRequest}catch(e){}};var _t={0:200,1223:204},zt=w.ajaxSettings.xhr();d.cors=!!zt&&\"withCredentials\"in zt,d.ajax=zt=!!zt,w.ajaxTransport((function(t){var n,r;if(d.cors||zt&&!t.crossDomain)return{send:function(i,o){var a,s=t.xhr();if(s.open(t.type,t.url,t.async,t.username,t.password),t.xhrFields)for(a in t.xhrFields)s[a]=t.xhrFields[a];for(a in t.mimeType&&s.overrideMimeType&&s.overrideMimeType(t.mimeType),t.crossDomain||i[\"X-Requested-With\"]||(i[\"X-Requested-With\"]=\"XMLHttpRequest\"),i)s.setRequestHeader(a,i[a]);n=function(e){return function(){n&&(n=r=s.onload=s.onerror=s.onabort=s.ontimeout=s.onreadystatechange=null,\"abort\"===e?s.abort():\"error\"===e?\"number\"!=typeof s.status?o(0,\"error\"):o(s.status,s.statusText):o(_t[s.status]||s.status,s.statusText,\"text\"!==(s.responseType||\"text\")||\"string\"!=typeof s.responseText?{binary:s.response}:{text:s.responseText},s.getAllResponseHeaders()))}},s.onload=n(),r=s.onerror=s.ontimeout=n(\"error\"),void 0!==s.onabort?s.onabort=r:s.onreadystatechange=function(){4===s.readyState&&e.setTimeout((function(){n&&r()}))},n=n(\"abort\");try{s.send(t.hasContent&&t.data||null)}catch(e){if(n)throw e}},abort:function(){n&&n()}}})),w.ajaxPrefilter((function(e){e.crossDomain&&(e.contents.script=!1)})),w.ajaxSetup({accepts:{script:\"text/javascript, application/javascript, application/ecmascript, application/x-ecmascript\"},contents:{script:/\\b(?:java|ecma)script\\b/},converters:{\"text script\":function(e){return w.globalEval(e),e}}}),w.ajaxPrefilter(\"script\",(function(e){void 0===e.cache&&(e.cache=!1),e.crossDomain&&(e.type=\"GET\")})),w.ajaxTransport(\"script\",(function(e){var t,n;if(e.crossDomain||e.scriptAttrs)return{send:function(r,i){t=w(\""
+ ],
+ "text/plain": [
+ ":DynamicMap []\n",
+ " :Image [x,y] (z)"
+ ]
+ },
+ "execution_count": 13,
+ "metadata": {
+ "application/vnd.holoviews_exec.v0+json": {
+ "id": "1002"
+ }
+ },
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "# plot composite polygon boundary on a latlon map \n",
+ "def callback(x_range, y_range):\n",
+ " cvs = ds.Canvas(plot_width=650, plot_height=400, x_range=x_range, y_range=y_range)\n",
+ " agg = cvs.polygons(geo_poly, geometry='geometry',)\n",
+ " return hv.Image(agg).opts(cmap=\"viridis\").opts(width=650, height=400)\n",
+ "\n",
+ "polymap = hv.DynamicMap(\n",
+ " callback, streams=[hv.streams.RangeXY()]\n",
+ ")\n",
+ "polymap"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 14,
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stderr",
+ "output_type": "stream",
+ "text": [
+ "/Users/ms/anaconda3/envs/vortexse_analysis/lib/python3.8/site-packages/holoviews/operation/datashader.py:413: SettingWithCopyWarning: \n",
+ "A value is trying to be set on a copy of a slice from a DataFrame.\n",
+ "Try using .loc[row_indexer,col_indexer] = value instead\n",
+ "\n",
+ "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
+ " df[category] = df[category].astype('category')\n"
+ ]
+ },
+ {
+ "data": {},
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "application/vnd.holoviews_exec.v0+json": "",
+ "text/html": [
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "
\n",
+ "
\n",
+ ""
+ ],
+ "text/plain": [
+ ":DynamicMap []\n",
+ " :RGB [lon,lat] (R,G,B,A)"
+ ]
+ },
+ "execution_count": 14,
+ "metadata": {
+ "application/vnd.holoviews_exec.v0+json": {
+ "id": "1163"
+ }
+ },
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "# plot a dynamic map of NLDN flashes grouped as per IC and CG category\n",
+ "\n",
+ "tmp = pd.concat([final_cg_nldn,final_ic_nldn])\n",
+ "\n",
+ "# Lasso boundaries + storm specific NLDN IC flashes\n",
+ "\n",
+ "from holoviews.streams import RangeXY\n",
+ "import holoviews.operation.datashader as hd\n",
+ "\n",
+ "points = hv.Points(tmp[['lat','lon','type']],kdims=['lon','lat'])\n",
+ "pts = hd.datashade(points, aggregator=ds.count_cat('type'),width=400, height=400)\n",
+ "\n",
+ "dynamic = hv.util.Dynamic(hd.aggregate(points, width=20, height=12, streams=[RangeXY]), \n",
+ " operation=hv.QuadMesh) \\\n",
+ " .opts(tools=['hover'], alpha=0, hover_alpha=0.2)\n",
+ "\n",
+ "# polymap * (pts * dynamic).relabel(\"Dynamic square hover\")\n",
+ "\n",
+ "hd.dynspread(pts, threshold=0.10).opts(height=500,width=850)\n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "In the plot above, you will notice that IC flashes are shaded red while CG flashes are shaded blue.\\\n",
+ "Also, the saturation of individual pixel corresponds to the count of flashes within each spatial bin"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 16,
+ "metadata": {},
+ "outputs": [
+ {
+ "name": "stderr",
+ "output_type": "stream",
+ "text": [
+ "/Users/ms/anaconda3/envs/vortexse_analysis/lib/python3.8/site-packages/holoviews/operation/datashader.py:413: SettingWithCopyWarning: \n",
+ "A value is trying to be set on a copy of a slice from a DataFrame.\n",
+ "Try using .loc[row_indexer,col_indexer] = value instead\n",
+ "\n",
+ "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
+ " df[category] = df[category].astype('category')\n"
+ ]
+ },
+ {
+ "data": {},
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "data": {
+ "application/vnd.holoviews_exec.v0+json": "",
+ "text/html": [
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "\n",
+ "
\n",
+ "
\n",
+ ""
+ ],
+ "text/plain": [
+ ":DynamicMap []\n",
+ " :Overlay\n",
+ " .Image.I :Image [x,y] (z)\n",
+ " .RGB.I :RGB [lon,lat] (R,G,B,A)"
+ ]
+ },
+ "execution_count": 16,
+ "metadata": {
+ "application/vnd.holoviews_exec.v0+json": {
+ "id": "1329"
+ }
+ },
+ "output_type": "execute_result"
+ }
+ ],
+ "source": [
+ "# Next, we overlay the NLDN flashes (CG+IC) on top of the composite polygon boundary \n",
+ "# This is more of a sanicty check to confirm if our NLDN subset conforms to lasso polygons\n",
+ "\n",
+ "from holoviews import dim, opts\n",
+ "\n",
+ "\n",
+ "def callback(x_range, y_range):\n",
+ " cvs = ds.Canvas(plot_width=650, plot_height=400, x_range=x_range, y_range=y_range)\n",
+ " agg = cvs.polygons(geo_poly, geometry='geometry',)\n",
+ " return hv.Image(agg).opts(cmap=\"viridis\").opts(width=650, height=400)\n",
+ "\n",
+ "polymap = hv.DynamicMap(\n",
+ " callback, streams=[hv.streams.RangeXY()]\n",
+ ")\n",
+ "\n",
+ "\n",
+ "points = hv.Points(tmp[['lat','lon','type']],kdims=['lon','lat'])\n",
+ "points.opts(color='k', marker='+',size=50)\n",
+ "\n",
+ "pts = hd.datashade(points, aggregator=ds.count_cat('type'),width=400, height=400)\n",
+ "\n",
+ "dynamic = hv.util.Dynamic(hd.aggregate(points, width=20, height=12, streams=[RangeXY]), \n",
+ " operation=hv.QuadMesh) \\\n",
+ " .opts(tools=['hover'], alpha=0, hover_alpha=0.2)\n",
+ "\n",
+ "# polymap * (pts*dynamic)\n",
+ "polymap*hd.dynspread(pts, threshold=1).opts(height=500,width=850)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Combine LMA and NLDN flash data to plot 1 minute flash count in the storm of interest"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 17,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "# set the time column in NLDN dataframe as index\n",
+ "final_cg_nldn.index = final_cg_nldn['time']\n",
+ "# drop the time column since it's already added as index\n",
+ "final_cg_nldn = final_cg_nldn.drop(['time'],axis=1)\n",
+ "\n",
+ "# bifurcate NLDN CG data into positive and negative CG dataframes\n",
+ "# we use a standard value of (+/-)15 kA current to remove noisy data \n",
+ "pos_cg_nldn = final_cg_nldn[(final_cg_nldn['current'] > 15)].index.value_counts().sort_index()\n",
+ "neg_cg_nldn = final_cg_nldn[(final_cg_nldn['current'] < -15)].index.value_counts().sort_index()\n",
+ "\n",
+ "# Also, count the total number of CG flashes\n",
+ "total = final_cg_nldn[(final_cg_nldn['current'] > 15) | (final_cg_nldn['current'] < -15)].index.value_counts().sort_index()\n",
+ "\n",
+ "# Resample total CG, +CG, and -CG dataframes to bin flashes\n",
+ "# in 1 minute time interval for the entire duration\n",
+ "total1 = total.resample(\"Min\").agg({'sum'})\n",
+ "neg1 = neg_cg_nldn.resample(\"Min\").agg({'sum'})\n",
+ "pos1 = pos_cg_nldn.resample(\"Min\").agg({'sum'})\n",
+ "\n",
+ "# Load data from lasso analysis generated flash stats csv file using lmatools\n",
+ "lasso = '/path/to/flashsort_lmatools/figures-length/testofflashstats/flash_stats.csv'\n",
+ "\n",
+ "# read the first and third columns in csv file containing \n",
+ "# info about time and number of flashes per 1 minute interval\n",
+ "lasso_data = pd.read_csv(lasso, usecols=[0, 2])\n",
+ "lasso_data = lasso_data[0:-1]\n",
+ "lasso_data.columns = ['time', 'number']\n",
+ "\n",
+ "# convert time values to datetime objects and use as dataframe index \n",
+ "lasso_data.index = pd.to_datetime(lasso_data['time'])\n",
+ "lasso_data = lasso_data.drop(columns=['time'])"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 18,
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ ""
+ ]
+ },
+ "execution_count": 18,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAABF0AAAHGCAYAAACmfvR1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABoSklEQVR4nO3deZhcVZn48e+bfSGhQgiboQg0iyBqVBZxgXQEdQwgIgryQ3YQR5yRcQUcFkXEBVFRVGSPo7JHcRwQMBhEAVFAkTUNoUFQ2YqEhCWdPr8/qjp2uqs6XdW1dfX38zz1VPe5557z1qnlVr997j2RUkKSJEmSJEnVNarRAUiSJEmSJLUiky6SJEmSJEk1YNJFkiRJkiSpBky6SJIkSZIk1YBJF0mSJEmSpBow6SJJkiRJklQDYxodQC2MGjUqTZw4sdFhSJIkSZKkPlasWJFSSiNiEkhLJl0mTpzI8uXLGx2GJEmSJEnqIyJebHQM9TIiMkuSJEmSJEn1ZtJFkiRJkiSpBky6SJIkSZIk1YBJF0mSJEmSpBow6SJJkiRJklQDJl0kSZIkSZJqwKSLJEmSJElSDYxpdACSJEmSpLV7/vnnefrpp3nllVcaHYqGudGjRzNlyhTWW289xo8f3+hwWppJF0mSJElqci+99BL/+Mc/mDlzJhMnTiQiGh2ShqmUEitXrmTp0qV0dnaSzWZNvNSQpxdJkiRJUpN76qmnmDFjBpMmTTLhoiGJCMaNG8f666/PtGnTePbZZxsdUksz6SJJkiRJTe6ll15inXXWaXQYajFTp05l2bJljQ6jqIhYJyJOiIi/RMSyiHg6In4XEYfGMMo8enqRJEmSJDW5rq4uxozxzzdV19ixY1m1alWjw+gnIkYB/we8BbgYOBuYBHwIuBDYFvhswwIsg+9aSZIkSRoGhtE/9zVMNPFramfgbcA3U0rH9RRGxDnA/cBHMOkiSZIkSZJUtqmF+yd6F6aUXomIp4Fhc+Vfr+kiSZIkSdIAIoJDDz200WGMJLcDOeAzEfGBiMhGxDYR8WXgTcApjQyuHM50kdQSzv/1bVzTMaXotr3alnHE3J3rHJGkVtHR0UFnZ2fZ+2WzWdra2moQkSS1pnJOdXnkkUeYNWvWgHWWLFnCRRddxD777MPs2bOHFlyZurq6uOSSS/jJT37CXXfdxfPPP8+UKVN43etex7777ssRRxzBpEmT1tino6ODs88+mxtuuIHOzk5eeukl1l9/fWbPns1ee+3FwQcfzOTJk+v6OBolpfRcROwNnAdc1mvTMuD9KaUFDQmsAiZdJLWEazqm8OCLM9l64uNrlD/44kyu6XicI+Y2KDBJw15nZye5XI5MJjPofXK5HIBJF0kqw/z589f4/eabb+bcc8/l6KOP5u1vf/sa22bMmLHW9pYsWcKpp57KrFmz6pp0eeqpp9h777259dZb2XnnnfnEJz7BxhtvTC6XY9GiRRx33HHcfPPNXHbZv3IJF198MR/5yEcYM2YMH/zgBznmmGOYNGkSf//731m0aBHHHnssCxYs4Lrrrqvb46ixMRFxR6/fz00pndunzgvAPcDPgd8B6wEfA34cEe9NKV1fn1CHxqSLpJax9cTHWXDUdmuU7fPDexsUjaRWkslkaG9vH3T9hQsX1jAaSWpNBx100Bq/d3V1ce6557LLLrv029asUkrst99+3HrrrXz729/m4x//+Brb/+u//osHH3yQK664YnXZjTfeyOGHH87222/PL3/5S171qletsc8JJ5zAww8/zKWXXlqXx1AnXSmlHUptjIjXkk+0HJdS+n6v8p+QT8T8MCLaUkrNt/RSH17TRZIkSZI0bCxfvpzjjz+etrY2xo8fz0YbbcTBBx/Mo48+urrORRddtDpZfthhhxERRARz5swBoLu7my996UvsuuuubLTRRowbN45sNstHP/pRnnnmmYpj+8UvfsGiRYvYf//9+yVcemy99daccMIJq3//zGc+A8Bll13WL+HSY4sttuD444+vOK5h6DhgAnB578KU0grgf4HNgFn1D6t8znSRJEmSJA0LXV1dvOtd7+KWW25hv/3245Of/CQPPfQQ3/ve9/jVr37FHXfcwcyZM9l111054YQTOP3009c4PWnDDTcE4JVXXuFrX/sa73//+3nve9/L5MmT+cMf/sD555/Pb3/7W/74xz8ybty4suPrmcFy9NFHD6r+kiVL+NOf/sSuu+7KNttsU3Z/Lawn+zS6yLYxfe6b2rAIUpIkSZK0pod+vHmjQxjQVgc+UvU2L7zwQm655RY+/elP89WvfnV1+e67786ee+7J8ccfz/z589liiy3YY489OP3004uenjR+/HiefPJJJk6cuLrsmGOO4S1veQtHHnkkCxYs4IMf/GDZ8d1zzz0Ag76GzF/+8peS9VesWMGKFSvWKJs+fXpZFxwexu4F3gkcCqx+oiMiA7wXeA7oaERg5fL0IkmSJEnSsHD11VczatSofqfazJs3j9mzZ/Ozn/2M7u7utbYTEasTLqtWrSKXy/H0008zd25+9YXbbrutoviWLl0KwNSpU4dc/6STTmLGjBlr3IZy6tMw803gWeCMiJgfEcdExAnAncDGwOdTSl2NDHCwTLpIkiRJkoaFRx55hE022YRp06b12/aa17yGZcuW8fTTTw+qrcsuu4ydd96ZiRMnMm3aNGbMmMEWW2wBwHPPPVdRfD3Jk2XLlpVVvyf50ttHPvIRrr/+eq6//nre+c53VhTPcJVSehTYCZgPtANnA58DHiO/ZPQ5DQyvLJ5eJEmSRryOjg46OzuLbit3uWhJUu2klKrSzlVXXcX+++/PTjvtxLe+9S023XRTJkyYwKpVq3j3u989qNkyxWy//fb86U9/4s4771w9a2Zt9QHuuuuuftu22morttpqKwB+9KMfVRTPcJZS6gAOaXQcQ2XSRZIkjXidnZ0lkyuZTIZsNlv/oCRpLWpxzZRm19bWxrXXXlv0M/vee+9l6tSprL/++gADXvtk/vz5TJgwgYULFzJp0qTV5ffff/+Q4nv/+9/PJZdcwnnnnTeopMvmm2/OG9/4Rn7729/ywAMPeDHdFuTpRZIkSeSTK+3t7UVvbW1tjQ5PkgTss88+dHd3c8YZZ6xR/n//93/ceeed7L333owalf8zd5111gHg2Wef7dfO6NGjiYg1ZrSklDjttNOGFN9ee+3Frrvuyk9+8hPOOaf4GTCLFy/my1/+8urfv/KVrwDwwQ9+kCeeeKLoPtWa4aP6c6aLJEmSJGlYOPTQQ7n44ov5yle+wpIlS9h1111ZvHgx55xzDhtuuCGnn3766rrbbbcdU6ZM4ZxzzmHSpElkMhk22GAD5s6dy3777ceVV17J3LlzOfjgg1m5ciULFizot1pQuSKCK664gr322ouPfexjzJ8/n7333puNNtqIXC7Hb3/7W37+85/z/ve/f/U+u+++O+effz4f+chH2HrrrfnABz7Am970JiZNmsQ//vEPFi1axK9+9Ss23nhjJkyYMKT4VH8mXSRJkiRJw8LYsWO57rrrOO2007j00ku56qqryGQyfOADH+C0005j0003XV134sSJ/PSnP+Xzn/88n/jEJ3j55ZfZbbfdmDt3LgcccADLli3jrLPO4lOf+hTTpk1jr7324owzzmD69OlDinHGjBncfPPNXHLJJfzkJz/hzDPP5Pnnn2fq1Km8/vWv51vf+haHHXbYGvsceuihvP3tb+fb3/42N9xwA5dffjmvvPIK66+/Pq9//ev57ne/y4c//GEmT548pNhUf9GK05QmT56cli9f3ugwJNXRPj+8F4AFR203qHJJ6m3hwoUAtLe3N2V7knTfffex7bbbNjoMtaBGvLYiYkVKaURkkLymiyRJkiRJUg2YdJEkSZIkSaoBky6SJEmSJEk1YNJFkiRJkiSpBky6SJIkSZIk1YBJF0mSJEmSpBpoSNIlIiZFxCMRkSLiO0W2bxMRCyLiuYhYHhE3R8TcRsQqSZIkSZJUiUbNdPkCsH6xDRHRBvwO2AX4KvBpYB3guojYvW4RSpIkSZIkDUHdky4R8UbgE8DJJap8GcgA70opfTmldA7wduAJ4LsREfWIU5IkSZIkaSjqmnSJiNHAD4FrgauKbJ8M7A3clFK6q6c8pfQCcB6wNbBjXYKVJEmSJEkagnrPdDkOeDVwbIntrwPGA78vsu3Wwr1JF0mSJEmS1PTqlnSJiM2BU4EvpJSWlKi2SeH+b0W29ZS9qsqhSZIkSZIkVd2YOvb1PeAR4BsD1JlUuH+5yLaX+tRZQ0QcDRwNMG7cuApDlCRJkiRJqo66zHSJiIOAdwLHpJRWDlB1ReF+fJFtE/rUWUNK6dyU0g4ppR3GjKlnLkmSJEmSVC933HEHhx56KFtssQUTJ05k8uTJbL/99nziE5/g/vvv71f/pZde4pxzzmHu3LnMmDGDsWPHkslk2HHHHfnsZz9bdB+pWmqenYiI8eRnt/wS+HtEbFnY1HOa0LqFsqfJr1DUe1tvPWXFTj2SJEmSJLW4L3zhC5xyyimsv/76HHjggWy77baklPjrX//KpZdeyne+8x2ee+45pkyZAsDDDz/MnnvuyX333cduu+3Gcccdx8Ybb8wLL7zAXXfdxQUXXMDXv/51Ojs7edWrvJKFqq8eU0ImAjOAeYVbXwcVbp8Gvk/+1KJditR7c+H+jhrEKEmSJElqYhdccAEnn3wy7e3tXH311ay77rprbP/qV7/KqaeeSkoJgBdffJF58+bR0dHBVVddxfve975+bb700kucddZZRERdHoNGnnokXZYDHyhSPgM4h/zy0ecDf04pvRAR1wD7RsTrU0p3A0TEOsCRwEPA7XWIWZIkSZLUJF555RVOPPFE1llnHS699NJ+CReAiRMncsYZZ6z+/bzzzuP+++/nhBNOKJpwAZgwYQLHH398zeKWap50KVzD5Yq+5RExq/BjR0qp9/bjgXcAv4qIs4ClwFHkTy+al3rSlpIkSZKkEeGWW27h73//Ox/+8IeZMWPGoPa54or8n5lHHnlkLUOTBtR0V5xNKS2OiLcCZwCfA8YBfwLenVK6oaHBSZIkSVKT2OxbjY5gYI/+Z/XauueeewCYPXt2WftMnTqVzTfffI3yVatW8dxzz61RNnnyZCZOnDjkOKW+GpZ0SSktAYqeOJdSug94b10DkiRJkiQ1paVLlwIwderUsvbZaKON+pXfd999vPa1r12j7Gtf+xqf+tSnhhakhpWI2AbYjfw1ZTchfwmUCcAzwFPAfcBvgFtSSssr7afpZrpIkiRJkkamVatW8dRTT61RNnHixNXJlmXLlg26ralTp65O1vS2+eabc/311wNw9913m2wZQSJiQ+Bw8pcw2az3piLV30v+7JuuwrVnf5BSur7cPk26SJIkSZKawmOPPdbvdKBDDjmEQw45BIA777xz0G1tv/32LFq0iEceeWSNNidPnszuu+8OwJgx/kk8EkTExsAXgQ+Tz4P0JFmeAO4EngaeBV4E1ivcNgdeT/6SJ/sC74uIxcDnU0qXD7ZvX2GSJEmSNAxV85opzWKjjTZaPQulxyabbMKWW27JRhttxIIFC3jmmWeYPn36Wtvab7/9WLRoEeeddx5f+tKXahWymlxEnAJ8EpgMdAG/BH4MLEop/W0t+44D3kB+1suHgK2An0bEccBHUkp/WVv/o4YUvSRJkiRJVTJhwgR23333NW7bbbcd48aN40tf+hLLli1j//33L3qa0UsvvcQJJ5yw+pSiI488kle/+tV87Wtf4+qrry7an4vjjggnkZ/BcjywcUppr5TST9aWcAFIKb2SUrotpXRCSmlzYA7wK+DNQPF1yPtwposkSZIkqekdfvjhPPbYY5x66qlsueWWHHjggWy33XZ0d3dz3333cfnll/PPf/6T448/HshfC+Z///d/2XPPPdl3332ZM2cO73znO9loo41YunQp999/P5deeimjR49m0003bfCjUw0dD3xnKBfD7ZFSWgQsioidyF94d61MukiSJEmShoWTTz6ZefPmcfbZZ7NgwQK+973vMWrUKNra2th///356Ec/ypQpU1bX32KLLfjjH//IBRdcwBVXXMGZZ57J888/z+TJk9lyyy058sgjOeKII9hmm20a+KhUSymlr9SgzdsHW9ekiyRJkiRp2Nhhhx24+OKLB11/4sSJfOxjH+NjH/tYDaNSNRWuw3LyAFW6Ukpj6xTOkJh0kSRJI0ZHRwednZ39ynO5HJlMpv4BSZKkYq4CFhcpfx3waeCa+oZTOZMukiRpxOjs7CyaYMlkMmSz2cYEJUmS1pBS+jPw577lEfGDwo/nV6OfiJgEHAm8C9gMmJhSauu1fV1gXj6k9JNK+jDpIkmSRpRMJkN7e3ujw5AkSWUoJEgOAP4GXFuF9mYDPwNmAlEo7ruc1VLgRODVEfFwSum2cvtxyWhJkiRJktTsPghMBS5MKa0aSkMRMR34X2BT4E/Ap8gnWNaQ8muKn08+KbNPJX2ZdJEkSZIkSfU0JiLu6HU7ehD7HEF+JsoFVej/E8DGwI3AzimlbwAvlqj7v4X7iqbJenqRJEmSJEmqp66U0g6DrRwR2wBvA25MKT1Shf73Ip/A+UxKqXugiimlByJiJbBlJR0500WSJEmSJDWzIwr351WpvTbgFeCuQdZfSv7UprKZdJEkSZIkSU0pIsYABwPPAldXqdlR5Gfb9L1wbrH+RwFTgOWVdiRJkiRJktSM9gI2BOanlF6uUptPApMiYoNB1H0zMA5YUklHJl0kSZIkSVKz6jm16PwqtrmocH/IQJUKs1y+SP76LzdW0pFJF0mSJEmS1HQiYhPg3cDtKaW/VLHpbxfuT4yIOSX63hRYQH7VolXAOZV0ZNJFkiRJkiQ1o0OB0VTvAroApJTuAk4jf3HcGyLi1sLPRMR3I2IR0AHMK+xyfErp4Ur6csloSZIkSZLUdFJKpwOn16jtkyPiOfLJl516bToGiMLPLwKfSymdXWk/znSRJEmSJLWMFStW8M1vfpO3v/3trLfeeowdO5YNN9yQ97znPVx00UV0dXX12+fuu+/mqKOOYuutt2by5MlMmDCBzTbbjH333Zf58+ezcuXKBjwS1VpK6ZvAZsAngCuA3wG3Az8DPgNsMZSECzjTRZIkSZLUIhYvXsy8efN48MEH2X333Tn++ONZf/31+ec//8kNN9zAYYcdxr333stXv/rV1fucfvrp/Pd//zfTpk1j//33Z/vtt2fcuHE8/vjj3HDDDRx88MH89re/5Qc/+EEDH5lqJaX0DPlrvHx7bXUrYdJFkiRJktTUbrrpJtrb21m4cCFz5swpWufFF19kzz335OGHH+bKK69k3333XWP7Zz/7Wf7whz/whz/8YXXZRRddxIknnsg73vEOrrzyStZdd9019jn55JO56667uOWWW6r+mNQ4EZEFVqWU/jbI+psAY1JKneX2ZdJFkiRJkjTsnXfeeTzwwAN89rOf7Zdw6bHjjjuy4447AvDKK69wwgknMGXKFC677LJ+CZces2fPZvbs2bUKW42xBHgSeNUg698CbEoFORSv6SJJkiRJGvauuOIKAI4++uhB1b/lllt48skned/73sd6661Xy9DUnGLtVYZUH3CmiyRJkiQNT1HR34D1k1Jdu7vnnnuYMmUKW2yxxaDrA0VnsSxbtoyXX3559e8RwfTp06sSp4alCUD/KzAPgkkXSZIkSVLT6O7u5tlnn12j7Pnnn199//TTT6+xbb311mPUqFEsXbqUDTfccND9LF26FICpU6f223bYYYdx5ZVXrv598uTJvPDCC4NuW62jcD2XGcA/K9nfpIskSZIkqWl0dnay+eabF922zz779Ct75JFHmDVrFlOnTmXZsmWD7qcn2dKTfOnt5JNP5phjjgHgk5/8JB0dHYNuV80nInYF5vQpXiciThpoNyADvKfw822V9G3SRZIkSZLUNDbaaCOuv/76NcruvvtuPvWpT/H1r3+d17/+9f3qA2y//fYsWrSIhx9+eFCnGG2//fYA3HXXXf22vfa1r+W1r30tANOmTavkYai5tAMnA73PeZtcKFubAF4CvlxJxyZdJEmSJGk4qvM1U+plwoQJ7L777muUjRmT/9P1TW96U8klo9///vezaNEizjvvPE4//fS19vPWt76VjTbaiKuvvppvfOMbXrOltS0BftPr992AlcDvB9inG1gK3ANcnFJaXEnHrl4kSZIkSRr2jjzySLbZZhu+/vWv87Of/axonT/+8Y+cc845AIwbN47TTz+dZcuWsf/++6++bkxfqUWTWyNJSunilFJ7z61Q/GzvsiK3d6SU3pdS+u9KEy7gTBdJkiRJUguYNGkSv/jFL5g3bx777LMP73znO9ljjz2YPn06Tz31FAsXLuS6667jM5/5zOp9DjvsMJ544gn++7//m7a2Nvbff3+23357xo4dy5NPPsmvfvUrfvvb364+FUkt4zDgxXp0ZNJFkiRJktQSttxyS+68805+8IMfcOWVV/KlL32JF154gfXWW48ddtiBiy++mAMPPHCNfU488UTmzZvHd77zHX71q19x4YUX0t3dzQYbbMCb3vQm5s+fz/7779+gR6RaSCldXK++TLpIkiRJkpranDlzBn2az6RJkzjuuOM47rjjBt3+7NmzOe+88yoNTyrJpIskSZIkSRoxIuLXFeyWUkrvKHcnky6SJEmSJGkkmTPIej3Tq4I1l5seNJMukiRJkiRpJDl1LdvXBXYGdgGeAb4HrKqkI5MukiRJkiRpxEgprS3pAkBEzAWuArZLKe1XSV+jKtlJkiRJkiSplaWUfg38J/C+iDiykjZMukiSJEmSJBV3KflTiypKunh6kSRJUg3kcjkWLlxY1j7ZbJa2trYaRSRJksqVUnopIpYD21ayv0kXSZKkKstms2Xvk8vlAEy6SJLURCLiVeQvrPtCJfubdJEkSaqytra2spMn5c6KkSRJtRURE4FzCr/+pZI2TLpIkiRJkqQRIyJOWkuVCcCmwLuA6UACvltJXyZdJEmSJEnSSHIK+UTK2gTQDXwppfTjSjoy6SJJkiRJkkaSRQycdOkCngPuBi5LKT1UaUcmXSRJkiRJGkBEcMghh3DRRRc1OhRVQUppTr36GlWvjiRJkiRJKiUiBn1bsmTJWttbsmQJp5xyCnfddVfNY+/dZ0Rw7LHHDqpeRLDnnnsWrbNy5UpmzJhBRDBr1qySbZ1zzjlEBOuuuy4rVqwYSviqAWe6SJIkSZIabv78+Wv8fvPNN3Puuedy9NFH8/a3v32NbTNmzFhre0uWLOHUU09l1qxZzJ49u5qhVs2ECRO49tprefLJJ9l4443X2Pbzn/+cp59+mgkTJgzYxgUXXEBbWxsdHR1cfvnlHHLIIbUMWWUy6SJJkiRJariDDjpojd+7uro499xz2WWXXfptaxV77bUXV199NfPnz+czn/nMGtsuuOACXve617Fq1SpeeOGFovvffffd/PGPf+SSSy7hrLPO4oILLjDp0mQ8vUiSJEmSNGwsX76c448/nra2NsaPH89GG23EwQcfzKOPPrq6zkUXXUR7ezsAhx122OpTeebMmQNAd3c3X/rSl9h1113ZaKONGDduHNlslo9+9KM888wzdXssG2ywAe95z3u48MIL1yh/8sknue666zjssMMG3P/8889nnXXWYd999+XQQw9l0aJFPPRQxdd8bToRsV5EfD0iFkfESxHxVEQsjIi3r33vtbYdEXFoRFwXEU9GxMsRsWqAW1cl/Zh0kSRJkiQNC11dXbzrXe/ijDPO4I1vfCNnnXUWH/rQh7jsssvYeeedefzxxwHYddddOeGEEwA4+uijmT9/PvPnz+fEE08E4JVXXuFrX/saW221FZ/+9Kf59re/zR577MH555/PnDlzeOWVV+r2mA4//HDuv/9+fv/7368uu/jiixk9evSAM3xefvll/ud//of99tuPyZMnc+CBBzJ27Nh+CZzhKiI2A/4IHAJcAfw7cDqwBHjVENseD1wPnA/sAWwIjCW/RPRAt7J5epEkSZIkDUOXXXZZo0MY0Ac/+MGqt3nhhRdyyy238OlPf5qvfvWrq8t333139txzT44//njmz5/PFltswR577MHpp59e9PSk8ePH8+STTzJx4sTVZccccwxvectbOPLII1mwYEFN4i9m3rx5bLjhhlx44YXssssuQH6mzl577cX6669fcr8FCxbw7LPPrj6daP3112fevHlcfPHFfPGLX2T06NF1ib+GfkQ+Z/G6lNKTVW77s8Dcws9XAT8DniC/VHRVOdNFkiRJkjQsXH311YwaNYrjjz9+jfJ58+Yxe/Zsfvazn9Hd3b3WdiJidcJl1apV5HI5nn76aebOzf8dftttt1U/+BLGjBnDhz/8YS699FJWrFjBLbfcwgMPPMDhhx8+4H7nn38+s2bNYrfddltdduihh/LEE09w7bXX1jrsmoqIXYG3AV9NKT0ZEWMjYlIVuzgASMCpKaX9UkrzU0o3ppR+M9Ctko5MukiSJEmShoVHHnmETTbZhGnTpvXb9prXvIZly5bx9NNPD6qtnlOSJk6cyLRp05gxYwZbbLEFAM8991xV416bww8/nKVLl3LVVVdxwQUXsMkmm/Cud72rZP1HH32UG2+8kXe+8510dHSwePFiFi9ezDbbbMPUqVM5//zz6xh9TbyncN8ZEdcALwLLI+LBiKjGVZU3J590ObMKbQ3I04skSZIkScNCSqkq7Vx11VXsv//+7LTTTnzrW99i0003ZcKECaxatYp3v/vdg5otU03bbrstO++8M9/97ne55557OPbYYwc8PejCCy+ku7ubc889l3PPPbff9l/84hf885//ZIMNNqhl2LW0TeH+h8BD5K/rMh74L2B+RIxNKQ3l4jXLgNEppeLLQlWRSRdJkiRJGobqdc2RZtLW1sa1115LLpcjk8msse3ee+9l6tSpq6+DElH6uqfz589nwoQJLFy4kEmT/nXWyv3331+TuAfj8MMP5yMf+QjAgKsWpZS46KKLmD179uoLA/f297//nY9//OPMnz+fT37ykzWLd4jGRMQdvX4/N6XUO3s0pXC/DGhPKb0CEBFXAw8Dp0fExSmlSrNjfwDeFRHrpZSerbCNQTHpIkmSJEkaFvbZZx9++ctfcsYZZ3DGGWesLv+///s/7rzzTg466CBGjcpfRWOdddYB4Nln+/9NPXr0aCJijRktKSVOO+20Gj+C0g444ACeeOIJ1ltvPbbeeuuS9W644QYeffRR/uM//oP99tuvaJ1vfOMbXHDBBc2cdOlKKe0wwPYXC/c/6Um4AKSUnouInwMHk58Nc1+F/X8DeDdwHPDfFbYxKCZdJEmSJEnDwqGHHsrFF1/MV77yFZYsWcKuu+7K4sWLOeecc9hwww05/fTTV9fdbrvtmDJlCueccw6TJk0ik8mwwQYbMHfuXPbbbz+uvPJK5s6dy8EHH8zKlStZsGABK1asqEqcd9xxR9EEzpgxY/jc5z5XdJ+pU6dyyimnrLXtnuu17LvvviXr7Lvvvpx55pnceuutvPnNbx5c0M3l8cL934ts61nJqP+FfQYppXRjRHwW+HJEvAKcmVKqzpPfh0kXSZIkSdKwMHbsWK677jpOO+00Lr30Uq666ioymQwf+MAHOO2009h0001X1504cSI//elP+fznP88nPvEJXn75ZXbbbTfmzp3LAQccwLJlyzjrrLP41Kc+xbRp09hrr70444wzmD59+pDjvO2224qugDR+/PiSSZfBePbZZ1mwYAFvfOMbmTVrVsl673//+znzzDO54IILhmvS5XbgGGBmkW09Zf+stPGI+HXhx2XAKcDxEfHXwu+lpJTSO8ruq1oXImomkydPTsuXL290GJLqaJ8f3gvAgqO2G1S5pJFp4cKFALS3tzc4kv6aOTZJjXffffex7bbbNjoMtaBGvLYiYkVKafIA26cBjwJLgVf3XPA2IjYmf2HdJ1JKpc/BWnv/lVwLJqWUSl/duARnukiSJEmSpKZRuHbLp4AfALdGxAXAOOCjhftjh9jFqUPcf9BMukiSJEmSpKaSUjo3Ip4GPgN8EegGfg8cmFK6ZYhtm3SRJEmSJEkjV0rpKuCqRscxFKMaHYAkSZIkSVIrMukiSZIkSZJUA55eJEmSJEmSWlJEPFz4cXFK6Z19ysqRUkpt5e5k0kWSJEmSJLWqWYX7l4qUlSNV0rlJF0mSJEkaBlJKRESjw1ALSamiPMJwc1jh/vkiZTVXl6RLRGwDnAS8EdgEGAt0Ar8EvpZSerJI/a8Au5Ffg/tPwMkppV/XI15JkiRJaiZjxoyhq6uLsWPHNjoUtZCVK1cyevToRodRUymliwdTViv1mukyE9gYuBp4HOgCXgscDRwQEbNTSv8EiIg24HeFOl8ln406CrguIv4tpXRDnWKWJEmSpKYwYcIEXnjhBaZNm9boUNRCli5dypQpUxodRkurS9IlpXQjcGPf8ohYBFwGHEo+wQLwZSADvCmldFeh3iXAX4HvRsSr0wiZAyVJkiRJADNmzKCzs5Px48czceJETzNSxVJKrFy5kqVLl/Lcc8+RzWYbHVJLa/Q1XR4t3E8DiIjJwN7ATT0JF4CU0gsRcR7wBWBH4PY6xylJkiRJDTNhwgQ23HBD/v73v/Pyyy83OhwNc6NHj2bKlClks1nGjx/f6HAaLiImkp/8MeD5eymlznLbrmvSJSImAOsAE4DtyF+3BfLXdgF4HTAe+H2R3W8t3Jt0kSRJkjTirLvuuqy77rqNDkNqCRGxDvAZ4ABgMEtBJyrIodR7psuRwNm9fl8CHJRSurnw+yaF+78V2ben7FW1CU2SJEmSJLW6iNgAWARsBQz2XL2Kzumrd9JlAXA/+dkubyB/KtGMXtsnFe6LzZd7qU+dNUTE0eQvzMu4ceOqEKokSZIkSWpBXwK2BlYAZwLXAf8gv6BPVdU16ZJSepz86kUACyLiSuAPETExpfRl8g8Y8qcY9TWhcL+iyDZSSucC5wJMnjzZC+1KkiRJkqRi9iR/utChKaUratnRqFo2vjYppT8DdwL/Xih6onBf7BSinrJipx5JkiRJkiQNxrrAK8DVte6ooUmXgonAeoWf/0L+1KJditR7c+H+jnoEJUmSJEmSWtJjQFdKaVWtO6pL0iUiNipR3g5sT2FlopTSC8A1wJyIeH2veuuQvwjvQ7hykSRJkiRJqtwCYFJE7Fjrjup1TZfvRcTGwK+BR8lfn+VN5JdmWgZ8slfd44F3AL+KiLOApcBR5E8vmpdS8notkiRJkiSpUl8FPgB8PyLekVLK1aqjeiVdfgIcAnyY/GpFiXzy5QfA11JKnT0VU0qLI+KtwBnA54BxwJ+Ad6eUbqhTvJIkSZIkqQWllJ6JiN2BHwP3RsQPyF/KZNla9ltUbl91SbqklC4DLiuj/n3Ae2sXkSRJaiYdHR10dnYW3ZbNZmlra6tKW7lcjkwmU0mIkiSptXQBS4CdgJMGUT9RQQ6lGS6kK0mSRrjOzk5yuVy/8lwuVzKBUm5bAJlMhmw2W0GEkiSpVUTELOB35E8xAohB3CrKn9Tr9CJJkqQBZTIZ2tvb1yhbuHBh1dqSJEkq+AKwCfA08FngOuAftVjNyKSLJEmSJEkaSd5B/nShD6WUbqxlR55eJEmSJEmSRpIM8CL5FZZryqSLJEmSJEkaSR4FIqWUat2RSRdJkiRJkjSSXAZMiIi5te7IpIskSZIkSRpJvgLcA/wwIjavZUdeSFeSJEmSJI0kHwDOB04B/hIRVwK3A8sG2imldEm5HZl0kSRJkiRJI8lF5FcvAgjgoMJtIAkw6SJJkiRJkjSATv6VdKkpky6SJEmSJGnESCnNqldfXkhXkiRJkiSpBpzpIkmSJEmSmkpElDr9Z3lKaZ26BjMEJl0kSZIkSVIzuhk4t0/ZykYEUimTLpIkSZIkqRk9nFL6UaODGAqv6SJJkiRJkppSRIyLiGFzOlFfznSRJElSxZ5f/GOWLfl50W1TZu3NulseWOeIJEktZD/gIGB0RDwFXAp8PqX0fGPDGjyTLpIkSarYsiU/5+Xn7mX8tO3WKH/5uXsBTLpIkooZExF39Pr93JRS32u33A5cDiwGpgLvAY4FdouIt6SUXqhPqENj0kWSJElDMn7adszc/adrlD1+wwENikaSNAx0pZR2GKhCSmnnPkWXRMSfgS8B/1m4b3pe00WSJEmSJA0HXwNeAeY1OpDBMukiSZIkSZKaXkppJfAEsP5Q2omIX0fEjRHRVp3ISvP0IkmSJEmS1PQiYgIwE7h1iE29DViZUuoYelQDc6aLJEmSJElqGhExvcSmL5KfPHLNELv4B/nTlGrOmS6SJEmSJKmZfD4i3gwsBDqBdcivXtQO3AacPcT2FwEHRMRWKaWHhtjWgJzpIkmSJEmSmslNwFLgEOCbwKnAesCJwJyU0otDbP/rQBdwZkTEENsakDNdJEmSJElS00gp/Qz4WQ3bvzMiPgRcBNwSEV8Dfgf8M6WUqtmXSRdJkiRJkjRiRMSqXr/uDFzRa1up3VJKqewcikkXSZIkSZI0ktT0lKLeTLpIkiRJkqSRpL1eHZl0kSRJkiRJI0ZK6Tf16svViyRJkiRJkmrApIskSZIkSRqxIm/9iMhWu22TLpIkSZIkacSJiDdGxFXA88A/gIf7bJ8WET+IiO9HxLhK+jDpIkmSJEmSRpSI+DDwe2AfYB3yKxqtsapRSuk5YHPgKGCPSvox6SJJkiRJkkaMiNgW+CEwFvg2sAPwdInql5BPxry3kr5cvUiSJEmSJI0k/wWMA76bUvoEQESsKlH314X7XSrpyJkukiRJkiRpJJkLJOAra6uYUnoCWAFUdJFdky6SJEmSJGkk2QRYnlJ6fJD1XwQmVtKRSRdJkiRJkjSSvAyMi4hYW8WImAhkyK9wVDaTLpIkSZIkaSRZQv4iulsNou57gNHAvZV0ZNJFkiRJkiSNJNeSX5HoPweqFBHTga+Sv/7L/1bSkUkXSZIkSZI0kpwFvAAcExEnR8SU3hsjYmJEHAjcAWwOPAN8v5KOTLpIkiRJkqQRI6X0D+BAYCVwEvAUMB0gIv4KPAvMBzYjf/2XD6WUllbSl0kXSZIkSZI0oqSUfgHsCvwRGAeMIX/K0bbA+MLPdwK7ppRurLSfMUMPVZIkSZIkaXhJKd0O7BQRrwPeRn4p6dHA34FbUkp3DLUPky6SJEmSJGnESin9GfhzLdr29CJJkiRJkqQacKaLJEmSJEkakSJiHLAHsAOwQaH4n+RXLro+pfTKUNo36SJJkiRJkkaciDgWOBlYr0SVZyPiCymlsyvtw6SLJElSk8jlcixcuLDotmw2S1tbW50jkiSpNUXEecBh5FcpAngc+Fvh51cBM8kvI/3NiHhDSunwSvrxmi6SJElNIJvNkslkim7L5XJ0dnbWNyBJklpURBwIHE4+4fIjYOuUUjaltEvhlgW2Ai4p1DmksE/ZnOkiSZLUBNra2krOZCk1+0WSJFXko0ACvpNS+s9iFVJKHcChEfE88HHg34Efl9uRM10kSZIkSdJI8jrySZcvDKLuFwp1X1tJRyZdJEmSJEnSSJKAXErpmbVWzNfJFfYpm0kXSZIkSZI0kjwIrBsR66ytYqHOVOCBSjoy6SJJkiRJkkaSC4DR5K/VsjbHFuqeX0lHJl0kSZIkSVJTi4hJEfFIRKSI+M5Q2kopfR+4FPhiRJxcbMZLREyOiJOALwI/TSmdW0lfrl6kltXR0VFyec1sNltyhQhJUnPJ5XJlrd6Ty+VKLr0sSZKGrS8A61ejoYi4AHgRWAacBHw6Iu4A/kb+2i0zgR2AicDzwEuFffpKKaUjBurLpItaVmdnZ9Ev3rlcDsCkiyQNA9lstux9MplMRftJkqTmFBFvBD4BfAY4swpNHko+uRKF3ycBu5aomwEO6VWXXvsmwKSLRq5MJkN7e/saZeX8t1SS1FhtbW0mySVJGsEiYjTwQ+Ba4Cqqk3S5hApXIyqXSRdJkiRJktSsjgNeDby/Wg2mlA6tVltrY9JFkiRpBDn/17dxTceUktv3alvGEXN3rmNEagW+riTVQkRsDpwKfCGltCQiZjU4pLK5epEkSdIIck3HFB58cWbRbQ++OHPAP5ylUnxdSSrTmIi4o9ft6BL1vgc8AnyjjrFVlTNdJEmSRpitJz7OgqO261e+zw/vbUA0ahW+riSVoSultMNAFSLiIOCdwK4ppZX1Cav6nOkiSZIkSZKaRkSMJz+75ZfA3yNiy4jYEtisUGXdQllmEG39OSL2rWJsm0TEdyLis4Op70wXDWsdHR10dnYW3VZsuWhJanUDfS5CfgnmvqsBVbKPJElSDU0EZgDzCre+DircPg18fS1tbQtcHhF/Ab4PXJZSerbcgCJiV+Bg4P8B44DPD2Y/ky4a1jo7O0smVzKZDNlstv5BSVIDDfS5mMvlAPolUCrZR5IkqYaWAx8oUj4DOIf88tHnA38eRFuvBb5GPnnzXeCbEXE9cDNwO3B3Sum53jtExDhgC2AnYGdgT2AmEMBT5C/u+4PBPBCTLhr2MpkM7e3tjQ5DkppGqc/FhQsXVnUfSZKkWihcw+WKvuW9Vi/qSCn1216irfuBvSKiHfgssAf5BMx7erW7CsgBLwPTyM+0WaNr4G/AucA3U0rLBvtY6nJNl4jYOiK+EBG3RsRTEbEsIu6KiBMjYnKR+ttExIKIeC4ilkfEzRExtx6xSpIkSZKk1pJSWphSejfQBnwZuBPoJp9QGQOsD7wKmFQoC+CfwOXA+4DNUkpfLCfhAvWb6XI48DHg58D/ACuBduA04IMR8eaU0osAEdEG/A7oAr4KPA8cBVwXEf+WUrqhTjFLkiRJkqQmkVJaQj4ZMtQ2TgROjIgpwA7AJuRPXZoAPEP+FKJ7U0oPDqUvqF/S5Qrgyyml53uVfT8iHiL/YI8AvlMo/zKQAd6UUroLICIuAf4KfDciXp1SSnWKW5IkSZIktaDCrJWanktdl9OLUkp39Em49Li0cL89QOFUo72Bm3oSLoX9XwDOA7YGdqxttJIkSZIkSUNXl6TLAGYW7v9RuH8dMB74fZG6txbuTbpIkiRJkqSm17DViyJiNHAS+Wu3/LhQvEnh/m9Fdukpe1WNQ1MT6ujooLOzs195qSVOJUlqNblcruhqUtls1iW9JUlqUo2c6fJN4M3ASSmlBwplkwr3Lxep/1KfOmuIiKMj4o6IuKOrq6uqgarxOjs7yeVy/cozmQzZbLb+AUmSVEfZbLboPxlyuVzRf0pIkqTm0JCZLhHxReBY4NyU0pd7bVpRuB9fZLcJfeqsIaV0Lvk1s5k8ebIX2m1BmUyG9vb2RochSVLdtbW1FZ3NUmzmiyRJah51n+kSEacAnwcuBI7ps/mJwn2xU4h6yoqdeiRJkiRJktRU6pp0iYiTgZOBS4Ajiyz9/BfypxbtUmT3Nxfu76hdhJIkSZIkSdVRt6RLRJwEnALMBw5LKXX3rVNYGvoaYE5EvL7XvusARwIPAbfXJWBJkiRJkqQhqMs1XSLiY8CpQCdwA3BgRPSu8o+U0vWFn48H3gH8KiLOApYCR5E/vWhekdkxkiRJkiRJTadeF9LdsXCfBS4usv03wPUAKaXFEfFW4Azgc8A44E/Au1NKN9QhVkmSJEmSNIJFxDRgVUpp6VDaqUvSJaV0KHBoGfXvA95bq3gkSZIkSdLIFBGbALsD/0wpXdtn22vITxZ5Q+H33wFHpJQerKSvuq9eJEmSJEmS1ECHk19ReU7vwoiYCPySfMIlCre3AjdExNRKOjLpIkmSJEmSRpLdC/eX9ik/BNgUeJb8tWUPAh4nf43Zj1XSkUkXSZIkSZI0kswq3N/fp3xfIAEnpJTOTyn9mHzyJYC9K+nIpIskSZIkSRpJ1geWppRe7CmIiFHAW8gnXa7oVfd6oBvYppKOTLpIkiRJkqSRZDQwvk/Za4FJwF9TSs/1FKaUuoHngMmVdGTSRZIkSZIkjSRPAuMjYvNeZe8q3P+uSP11yF/npWwmXSRJkiRJ0kjy+8L9yRExKiJmAB8lf2rRdb0rFhIz48knaspm0kWSJEmSJI0k3yrcfxjIAY8BmwGPAL/oU3ePwv2fKuloTCU7SZIkldLR0UFnZ2dZ++RyOTKZTG0CkiRJ6iWldHtEHA58G5hSKL4fOCCl1NWn+sGF+4WV9OVMF0mSVFWdnZ3kcrmy9slkMmSz2doEJEmS1EdK6WJgI2Bn8isTbZ9S+nPvOhExDjgXOAz430r6caaLJEmqukwmQ3t7e6PDkCRJKqmwZPQfBtj+CnDJUPpwposkSZIkSVINONNFkiRJkiSNOBExCTiS/HLRmwETU0ptvbavC8wDUkrpJ5X0YdJFkiRJkiSNKBExG/gZMBOIQnHqU20pcCLw6oh4OKV0W7n9eHqRJEmSJEkaMSJiOvkL425KfinoT5FPsKwhpZSA88knZfappC+TLpIkSZIkaST5BLAxcCOwc0rpG8CLJer2rFpU0QoBnl4kSdIIksvlWLhwYb+yTCbTmIA0ZMWe0x7ZbJa2trai2yRJGsH2In8q0WdSSt0DVUwpPRARK4EtK+nIpIskSSNENpstWp7JZEpuU3Mb6HnL5XIAJl0kSeqvDXgFuGuQ9ZcC61bSkUkXSZJGiLa2Nv8AbzEDPaelZr9IkiRGAV2Fa7YMKCJGAVOA5ZV2JEmSJEmS1BQiYpuI+J+IuC8ino+IFRFxf0R8IyI2rkIXTwKTImKDQdR9MzAOWFJJRyZdJEmSJElSM5lJ/kK3VwPHk7/w7fXA0cAfB5ksGciiwv0hA1UqzHL5Ivnrv9xYSUeeXiRJkiRJkppGSulGiiQ5ImIRcBlwKPDVIXTxbfIJlxMj4g8ppZuK9LUp8F3yqxZ1AedU0pEzXSRJkiRJ0nDwaOF+2lAaSSndBZwGTAVuiIhbCz8TEd8tJHc6gHmFXY5PKT1cSV/OdNGwcMUjqwBor2hl9NZw/q9v45qOKUW37dW2jCPm7lzniCRJgp89+w6uf/6tjL9izfKXnzuJPda9hY81JixJagnPL/4xy5b8vOT28dO2Y8abTqpjRPUVEROAdYAJwHbAVwqbfjnUtlNKJ0fEc+STLzv12nQMEIWfXwQ+l1I6u9J+nOmiYeHRlzfg0ZeHetre8HZNxxQefHFmv/IHX5xZMhkjSVKtXf/8W3nopc36lT/00mZc//xbGxCRJLWOZUt+zsvP3dvoMGphTETc0et2dIl6RwJPAY8B1wEZ4KCU0s3VCCKl9E1gM/LXjLkC+B1wO/Az4DPAFkNJuIAzXaRhZeuJj7PgqO3WKNvnhy35ISxJGka2mvAol+7X9/j0aInakqRyjJ+2HTN3/2mjw6i2rpTSDoOotwC4n/xslzcAewMzqhlISukZ8td4+XY12+1h0kWSJEmSJDWdlNLjwOOFXxdExJXAHyJiYkrpyw0MbdA8vUiSJEmSJDW9lNKfgTuBf290LIPlTBdJkiRJkjRcTATWq0ZDETEF2BN4XaHNsQNUTymlI8rtw6SLJEmSJElqGhGxUUrp70XK24HtgZuq0MehwLfIXy9mdXGRqqlQngCTLpIkSZIkaVj7XkRsDPwaeJT8ktFvAg4AlgGfHErjEfEu4HzyyZSXgN8DTwBdQ2m3GJMukiRJkiSpmfwEOAT4MPnVihL55MsPgK+llDqH2P5nyCdcfg+8N6X09BDbK8mkiyRJkiRJahoppcuAy2rYxZvIJ3IOrWXCBUy6aIg6Ojro7CyeZMxms7S1tdU5IkmSJEmSBjQGeCGl9FCtO3LJaA1JZ2cnuVyuX3kulyuZjJEkSZIkqYE6gPERMbrWHTnTRUOWyWRob29fo2zhwoUNikaSJEmSpAH9CDgD+DfgF7XsyJkukiRJkiRpJPkm8AfgnIjYqpYdOdNFkiRJkiS1pIg4uMSm+cAXgLsj4grgNvLLUZeUUrqk3P5NukiSJEmSpFZ1EfmVikoJ4P8VbgNJgEkXSZIkSZKkgk4GTrrUlEkXSZIkSZLUklJKsxrZv0kXNY2Ojo6Sy0yntB4RUbW+crlcyRWWstksbW1tVetLkiRJkjQyuXqRmkZnZye5XK7otoggRlUn6ZLNZslkMkW35XK5kokfSZIkSZLK4UwXNZVMJkN7e3u/8m8tvrdqfbS1tZWcyVJq9oskSZIkaWSIiD2BPYBVwC9TSjdU2pYzXSRJkiRJ0ogREftGxMMR8f0i274B/Aw4FvhP4LqI+FqlfZl0kSRJkiRJI8newGbAzb0LI+KNwCfILyP9GNBR+Pm/ImJOJR2ZdJEkSZIkSSPJjoX7G/uUH164vxrYIqW0NfBd8omXoyrpyKSLJEmSJEkaSWYAXSmlv/cpfyeQgK+klLoLZacX7neppCOTLpIkSZIkaSTJAC/0LoiI6cCWQC6ldHtPeUrpSWA5sHElHbl6kWoml8uVtRpQLpcruZSzJI1EHR0dAy5jn81mS67GVmsDxebnuSRJanIvAOtGxNiU0spC2dsK978vUn9lkbJBcaaLaiKbzZb9hTuTyZDNZmsTkCQNQ52dneRyuaLbcrncgAmZWhsoNj/PJUlSk7uf/HVa3tOrbH/ypxb1vbjuJGBdoO+pSIPiTBfVRFtbW8P++ypJrSSTydDe3t6vvJyZhLVSKjZJkqQmdxXwZuC8iHg1+VOH9ge6gcv71N2RfILmkUo6MukiSZIkSZJGku8ABwGvI3+h3CiUn51SerhP3X3Jz4BZVElHJl0kSZIkSdKIkVJ6KSLeBnyC/KpEOeAXKaWf9K4XEeOA3YBO4FeV9GXSRZIkSZIkjSgppReA09ZS5xVg9lD68UK6kiRJkiRJNWDSRZIkSZIkqQY8vUiSJEmSJI1IETETeAswE5jMvy6q209K6Qvltm/SRZIkDSiXy/VbojqXy5HJZBoTkCRJ0hBFxPrA94F9GCDR0lOd/ApGJl0kSVL1ZLPZouWZTKbkNkmSpGYWEZOBm4BtgVeAu4GdCj/fDmwEbFmo/izwl0r7MukiSZJKamtro62trdFhSJIkVdPHgO2A+4F3pJSejIhu4NmU0q4AEbEZ8BVgP+DalNJXKunIC+lKkiRJkqSR5H3kTxc6PqX0ZLEKKaVHU0oHAJcBp0fEOyrpyKSLJEmSJEkaSV5duL+2T/nYInVPJH9Nl49X0pFJF0mSJEmS1DQiYuuI+EJE3BoRT0XEsoi4KyJOLFyPZagmALmU0su9yl4C1ulbMaX0CPA8+Wu+lM2kiyRJkiRJaiaHA8cBHeRXDPo08ABwGvC7iJg4xPb/AUyNiN45kaeAcYUlpFeLiNHkl5KeXklHJl0kSZIkSVIzuQKYmVL6fymls1NK308p7Q98CXgdcMQQ23+UfD5kk15ldxXu39en7t7kFyH6ZyUduXqR1Ecul2PhwoX9yrPZrCt4SFqrjo4OOjs7i24r9TlSap9cLkcmk6l2iFLTGuj9MxCP0ZLUWlJKd5TYdCn5a6xsP8QubgTeBswFLunV9l7AlyNiAvkkzOuBz5O/6O7/VdKRM12kXrLZbNE/cHK5XEVfAiWNPJ2dneRyuX7lA32OlNonk8mQzWarHKHUvEq9FwbiMVqSRpSeU3/+McR2riF/cdz/16vsJ8BNwCTgDPIX2f0KMLXQ3ymVdORMF6mXtra2ov8pKzbzpVWd/+vbuKZjStFte7Ut44i5O9c5Imn4yWQytLe3r1G2ts+RYvuMJCtOO5LRl/+86LZVH9ibSZ8/r84R1c5Aj3XF3M1YPq/4jI0ps/Zm3S0PXKPs+cU/ZtmS4m3d8bdd+M2Kdr61+N41yh98cSZbT3y8gsiLe+eiy9j19l/C/0xao/y0J1ewaKf3wFGnlNVeue+FkXSMrpeBvguA3wc0sJH0eV4Pk/+3g0m/fhSmzSleYfZs+OY36xdQAxWurXIS0AX8eChtpZTupM8klJRSioh55Ge27A9sSv4CutcCn08pPVFJX3WZ6RIRx0fE5RHxcESkiFiylvrbRMSCiHguIpZHxM0RMbcesUoj3TUdU3jwxZn9yh98ceaAX8AkaShGX/5zxix+pl/5mMXPlPzyPlyVeqxjO3L5L9ZFvPzcvUWTK8uW/JyXn7u3yB7wmxXtPLxyi37lW098nL3alpUZdWm73v5LNn/s/n7lmz92fz4Zo2Gn1HcB8PuA1m4kfZ7Xw6RfP8rYjlyjw6iFMRFxR6/b0YPY55vAm4GTUkoP1CKolNKLKaUTU0pbppTGp5Q2SCkdnFKqeEplvWa6nA48C/wJyAxUMSLagN+Rz159lXxm6Sjguoj4t5TSDbUNVdLWEx9nwVHbrVG2zw+Lf6mXpGrp2nI64+9e8xp1L79+gwZFU1vFHitz5jAemLn7T/vVf/yGA0q2NX7adkX3eWXJc8wcu5IFR25XZK/qemTTV7P9TTetWbZNRStrqkkU+y4Afh/Q4Iykz/N6WNmWYXyfz9gW0JVS2mGwlSPii8CxwLkppS/XLqzqq1fSpS2l9DBARNxDkbWve/ky+cTMm1JKdxX2uQT4K/DdiHh1SinVNlxJkiRJktRoEXEK+VN+LgSOaWw05atL0qUn4bI2ETGZ/HJMN/UkXAr7vxAR55Ffn3tH4PZaxClJkiRJkppDRJwMnEx+haEjK5mAERFVW5WgktOMmu1Cuq8DxgO/L7Lt1sK9SRdJkiRJklpYRJxEfsWg+cBhKaXuCpt6pEohJSrIoTRb0mWTwv3fimzrKXtVnWJRQUdHR8mlGHO5XNEllltRLpcre4WEbDZbdDUkSbUx0OcVFH9PVnufen4uFvtcGkmfy6rMQK/frq7pxKgoa5/1C/907PtaXN+zwSVJFYqIjwGnAp3ADcCBEWscn/6RUrp+sM1VK6xKdmq2pEvPWoMvF9n2Up86ayhc7fhogHHjxlU/shGss7Oz5Jf4TCZDNlu12VpNq5LHmMvlAEy6SHU00OdVqfdktfep1+diqT5GyueyKjfQ6zelBEX+jzjQPqUkElG177mSpBFmx8J9Fri4yPbfAINNumxelYgq1GxJlxWF+/FFtk3oU2cNKaVzgXMBJk+e7L9WqiyTydDe3t7oMBqmra2t7ORJubNiJFVHqc+rgd6T1dynXir5XJJ6lHr9nnn/n8ve557Cfx77bvuLCRdJUoVSSocCh1aprUer0U6lRjWy8yKeKNwXO4Wop6zYqUeSJEmSJElriIjnIuLePmUHR8QH6tF/s810+Qv5U4t2KbLtzYX7O+oXjiRJkiRJGsbWBV7sU3YR8CRwea07b6qZLimlF4BrgDkR8fqe8ohYBzgSeAhXLpIkSZIkSYOzkn9drqS3upwHW5eZLhHxYWCzwq8zgHER8fnC74+mlOb3qn488A7gVxFxFrAUOIr86UXzKlmXW5IkSZIkjUh/B2ZGxFYppYfq3Xm9Ti86AtitT9kXC/e/Ib/uNgAppcUR8VbgDOBzwDjgT8C7U0o31CFWSZIkSZLUGn4NHAL8PiJ+DbxQKF83Ii4oo52UUjqi3M7rknRJKc0ps/59wHtrE42K6ejooLOzs+i2cpeIlKRqGOhzKZvNunKPNAwkUlmr+fmdQ5JUAycBu5M/e2a/XuUTGdwKSYn8qUiJ/ISSsjTbhXTVIJ2dnSW/6GQyGbLZbP2DkjSilfpcyuVyACZdpGEgyjxd3u8ckqRqSyk9FhGvAT4EvBaYRD7Z8iJwWa37N+mi1TKZDO3t7Y0OQ5JWK/a5VM5/zSU1TkQ+4eJ3C0lSo6WUlgI/6Pk9Ig4Fnk8pHVbrvk26SJIkSZKkkWQR8HQ9OjLpIkmSJEmSRoxyrzs7FKPq1ZEkSZIkSdJIYtJFkiRJkiSpBky6SJIkSZIk1YBJF0mSJEmSpBrwQrqSpGEnl8sVXTo6l8uRyWRq3lct+pGaRUrJ17wkSVVi0kWSNKxks9mS2zKZzIDbq9VXtfuRmkWMCujuX+5rXpKkyph0kSQNK21tbbS1tbVcX1IzGD1qNIyC9vb2RociSVJL8JoukiRJkiRJNWDSRZIkSZIkqQY8vUiSJEmSJLW0iJgKbAR0pZQe7rPtLGDdAXY/PaW0uJJ+TbpIkiRJkqRWNx/YEzgT+EyfbQcAGwywbwbYt5JOTbpIkiRJkqSWFRFbAXsBTwMnDVB1UZGyNwHvjYjNU0qPlNu3SZdhrKOjg87OzqLbstmsK240sUY/dwP139U1Pb9kqGpuoOcBir8WKtlH/5LL5Vi4cGG/skwm07D+6x2DRpaurq6GvuYbrdHH22bgcUOSANivcP/DlNJLpSqllPot3xcRXwROBD4EnF5ux15Idxjr7Owkl8v1K8/lcgMeXNV4jX7uSvUPkFIidaeax6CBn4dSr4VK9lFeNpst+odmJpMhm802rP96xqCRZcL4CYwZ0///ayPp9dbo420z8LghSQC8DUjA/1aw72W92iibM12GuUwmQ3v7msm4Yv9FVfNp9HNXrH+AM+//c91iUOnnYaDXQiX7CNra2hr639xG96+RZ+LEiUycOLHo58VI0ujjbTPwuCFJbF+4/0MF+94LdAGvqaRjky6SJEmSJKmVTQeWppRWltj+J2C9YhtSSqsiYmmp7Wtj0kWSJEmSJLWyMeRnqxSVUpq3lv3HUmH+xGu6SJIkSZKkVvYMMCUiJpS7Y0RMBKYAz1bSsUkXSZIkSZLUynqWet6lgn3fUrh/uJKOTbpIkiRJkqRW9hsggKMq2Pco8isf/aaSjr2mywjT0dFRdGnAXC5XcilTSUPT6Pddqf57ZLPZhq2qM1Bsfi5JqqVcLldy9Z5Gfi4OV6u6V5G6U9Ex7eqaToyKBkQlSavNBz4HfDAifppS+vlgdoqI9wIfIJ90+VElHTvTZYTp7Owkl8v1K89kMmSz2foHJI0AjX7fleof8n90DJSQqbWBYvNzSVKtZLPZkkndRn8uDlepO5FSKr4tJVJ38W2SVEpEHB8Rl0fEwxGRImJJpW2llO4HfkI+B3JZRHwyIsYP0Pf4iPgUcGmh6MeFNsrmTJcRKJPJ0N7e3ugwpBGl0e+7Uv2X+i9vPTV6bCSNPG1tbSVnsjTD5+JwFRFFP8/PvP/PDYhGUgs4nfzFa/8EZKrQ3jHAG4FXA18FPhcRNwB38a+L5K4HzAZ2L/wcwL3ARyvt1KSLJEmSJElqNm0ppYcBIuIeYJ2hNJZSeiEi3gZcBswFpgMfLNz66jkn8kZg/5TS8kr79fQiSZIkSZLUVHoSLlVu89mU0u7A3sC1wAryCZbetxWFbXullPZIKVW0VHQPZ7po2OjuWs7jNxywRtm65/yJcVPbGH/R/9W07/N/fRvXdEwpa5+urukAfGvxvUW37bLuEmp9QsWK045k9OX9rxG1U1cX/5z7FmjgKR2lYuux6gN7M+nz59Uxon9Z2/O9V9syjpi7cx0jqp3fPbmM3z8/q+jrtB6Pc6D+B9JKzwEM/H4o970wXN9bX3xxU9Yb+wKvqnNMQ1VqvJ/umsbP3vghfrVr/3+elXqs/1gO69x7F0teM6ffPqlrBU/ttg0zd1+z/KHLxjLjNw/w/Jj++3xqJdzZfiDsd3QZj6g+BnotlHp/l9qnGY63qr5KXiPN3E8pzfyZ3YqG4/H26a5pPNe1Ll+8ovj27WbAybsNuZsRJ6X0C+AXETEG2Iz8rBeAZ4BHU0pd1erLmS4aFkZPmM6oMZP7lY954Am484817/+ajik8+OLMqrX38MpZ/P75WVVrr5TRl/+cMYuf6Vc+YcnzbPDr39W8/4GUig1gzOJnBjyI1dpAz/eDL84sOwHXzH7//CweXjmrX3m9Hmep/gfSas8BlH4/VPJeGK7vrRfTJJ7rnl50WzMrNd6ZR55kh1tvLLpPqcf6szceyL2vml10n+xjDzLjNw/0K5/xmwfIPvZg0X22+9tdvPdPPx4g+sYp9VoY6P1dybG4XsdbVV8lr5Fm7qeUZv7MbkXD8Xj7XNe6rOgueb3X4WxMRNzR69aQ/xCklLpSSh0ppdsLt45qJlzAmS4aJsZM3JAxEzdk5u4/XaP85TEb1C2GrSc+zoKjtht0/Z4L8RW7oNye36vfBeW6tpzO+Lv/uUbZ8tesV7f+B1IsNoCXX1+/57WUUs/3Pj8sb0bGcLDF2CUsOOp1a5TV83EW638grfgcQPH3Q6XvheH43vrrN0cDo+sfUBUUG+/Fr5nDZCjrsd64x9HcuMfRXLpf/z7u2Wankv13bro12//1pv4b5uRjaFbFXgtre38X26dZjreqvkpeI83cTynN/Jndiobj8XbSqJeLHhuGua6U0g717LAwq2USQEpp6SD3mVr4cXlKaVUl/TrTRZIkSZIktbqfAs8BF5WxzwUV7LMGZ7pIkvpZ1b2K1J1KLp2azWb7Lbfa0dFBZ2dnyTaL7SNJ1ZDL5Yp+XqW0HhFRZI/K2hvIcP2MK/ZYu7qmE6PKH7eBpFT6mFJMz7V6isVWrHxtmvX5WdXdTeru5ndFHs9OXV3EKP9HLlVDRLwG2Bd4Hji8jF2PAt4BfCgivpBSeqjcvn0XS5L6Sd2JlFLRbblcrmhypbOzk1wuV9Y+kjRU2WyWTCZTdFtElJ08GKi9UobrZ1ypx5pSInUXPwZUIkZFRcmvamnm5yd1d5c83uafh+46RyS1rP9XuD8npZQb7E4ppeeAs8nnTg6qpGNnukiSioqIotdIGOi/i5lMpux9JGko2traSs5gKHdltLW1V8pw/Ywr9VjPvL+618IZPWo0jCp+3Z1Sep67vvuUKh9Isz8/pY63yxuYqJJa0NuBBFxZwb5XAZ8H5lTSsUkXSZIkSZLUVCLiw+SXcwaYAYyLiM8Xfn80pTS/jOa2BrqBOysI5c+FfV9dwb4mXSRJkiRJUtM5AtitT9kXC/e/AcpJumSAXCp1Pt8AUkrdEZED1i13XzDpIkmSJEmSmkxKaU4Vm1sBTBnC/usAL1ayoxfSlSRJkiRJreyfwNiIKHsZs8I+4wptlM2kiyRJkiRJamW3Fu73rWDf9xfub6ukY08vanIdHR0ll7jL5XIllzTM5XJFr9Q+0D6qvlLPQ0rrlb104qruVaTu1K+9rq7pZS+HuTYpldfPQK/TbDZb9ioQq7q7Sd3d/K7I2FXS3nBV7PXje1iSVA/FvgtA/vsA9F8RqJLvNpJUR78APgx8OiJ+lFJ6cjA7RcQmwKfIr3z0i0o6dqZLk+vs7CSXyxXdlslkyGaz/cqz2WzJP8pK7aPqG+h5iIiyEyWpO1Hsuk8pJVJ32deDKilGRdEvTQP1U+p1msvlSiZjBpK6u4s+1krbG45KvX58D0uSaq3Ud4EB96ngu40k1dGVwEPAdOC6wZxmFBFbAtcC6wOLgcsr6diZLsNAJpOhvb190PXb2tpGzEyAZjbQ8/CtxfdW1GZE9HstnHn/nytqq5TRo0bDKMrup9jrtNh/yAar2GMdSnvDje9jSVKjlPouAP/6DtN3W6XfbSSpHgorEB0CLAReA/w5In4ELCC/jPSzharrAW8A3gccCEwCXgYOrWTlIzDpIkmSJEmSWlxK6daI+CD5paanAkcWbqUE8ALw4ZTS7yvt19OLJEmSJElSy0spXQPsQP5UoUQ+sVLslgp13pRS+tlQ+nSmiyRJkiRJGhFSSouB/SNiA6Cd/OlG08knW54G/gosTClVtER0XyZdJEmSJEnSiFJIqlxa6348vUiSJEmSJKkGnOkitYCUUtGVfXbq6gLgd3227ZhS2UtBjiSruleRuouPaVfX9KZeEjOXy/WLO5fLlVy+HIq/flJab8DXSCX9FNtnbf1I0nBR7DMOIJvNFl2NrqOjg87OzpLtldqvlFLtre2zeaQo9fx0dU0H+q+QWKp8bcp93qqt2q8rSUNn0kUa5mJUQHeZ+0QQo5zoVkrqTpRaES6lVPZ410s2my1anslkSm4r9frJv0aKJ0Mq6adk/wP0I0nDRanPuFwuB1D0j9zOzs6SCZGB9iulVHsDfTaPFPV6/JU8b9VW7deVpKEz6SINc6NHjYZR0N7e3m/by2Pyb/G+23rKVVpEFB3TM+//cwOiGZy2trayv0iVev18a/G9Ve2n1D4D9SNJw0Wpz7i1zZLIZDJFjzXlzq5YW3sj3UDHrZ7jUKnjYDnjWenzVm3Vfl1JGhr/1S1JkiRJklQDJl0kSZIkSZJqwKSLJEmSJElSDZh0kSRJkiRJqgGTLpIkSZIkSTXgEiZNoqOjg87Ozn7lpZZ8kxohpVT0yvfN8Dpd1b2K1N0/vq6u6S5JLElqmFwuV9Gxs9R+2WzWJX9HkFLfvaB6r4VSf4dUux9VJqUElF59KpPJ8IY3vKGeIalMznRpEp2dneRyuX7lmUyGbDZb/4CkPmJUEFE8edEMr9PUnVYflNYoT4nU3b9ckqRay2azJRMrAx07S+2Xy+UG/ONYrSVGjSr53auar4VSf4dUux9VLuF32eHMmS5NJJPJ0N7e3ugwpKJGjxoNo2jq12hE9IvvzPv/3KBoJEkjXVtbW0UzBErtV+o/3WpNo0eNglGjin73qvZrodTfIb7mmkPQ/zuuhg9nukiSJEmSJNWASRdJkiRJkqQaMOkiSZIkSZJUAyZdJEmSJEmSasCkiyRJkiRJUg2YdJEkSZIkSaoBl4yWmkwulyu6PF9K6xERDYioNlZ1d5O6u/ldkce6Y0olH2s1x6ejo4POzs6qtCVJUj2UOg7mcjkymUz9A9IaSj0/ANlstqIlxKupWHyVvHZKfYcC2KmrC6Dfd7ydurqIUcX/5z9Qe9Uct1LfPweKrdy2Km1PrcukS5PwICnIH1RKiQhiVOskAlJ3Nymlotvyj7X/gara49PZ2Vn0i0arjbUkqTUMdBzMZDIDblftDTT+uVwOoKFJl1LxVfLaKfUdaiApJejuLqu9ao9bqe+fA8VWbluVtqfWZdKlSbzhDW9odAhqAm1tbSUPKt9afG+do6m9iKC9vb1f+ctjin801WJ8MplMvxhacawlScPfQMdBNd5Az0+p2S/1VO3XT7HvUPCv73F9ty1fyyziYu3VYtyKff9cW2zltDWU9tSamnLOU0SMiojjIuL+iHgpIh6LiDMjYnKjY5MkSZIkSbXVKnmBpky6AGcB3wDuBT4OXA78B3BNRDRrzJIkSZIkqTpaIi/QdKcXRcRryA/oVSml9/cqfwT4NnAA8OMGhSdJkiRJkmqolfICzZgd+hAQwDf7lP8QWAEcVO+AJEmSJElS3bRMXqAZky47At3A7b0LU0ovAXcVtkuSJEmSpNbUMnmBKLXMVaNExF+ADVJKGxbZdhnwAWB8SumVUm1Mnjw5LV++vIZRqt72vwLufQq2m7Fm+X+fuBOzHruPJZtuW9P+V3SPZ9Kol9l+40lVae+eJ1esbnOoMQzUVs/YfPFLa3xWDThulfRTSkqJRCLofwX3zR+/n0dmvprTTv9Dv22VPK+l4v7LE8t5MU1gYrzUP75CbNHnCvMDPd8DtVdKqX6aQTWf72r2P5B6xVZPlbxXy22r0vaqaaDne/kf7uLeV83m66fftEZ5o2Nem1Lj/akT5rDd3+5i8o6z++1T6rH2HOcu3a9/P/dssxPZxx6kc9Ot1yjvKdv+gf7PN3PmwF13wez+MZRSyfurVAylYobKPnv8vKifgca61JhWsk8lMVT7Oa1XP6W+E1X6fahUe9X+zlFJP6U+Fz9/wo6rH2+/fkq0N9B3yUqUGu+BYiu3rUrbK2Wzxx/i0Zlb8doH+/cznEXEipRSyQviViMv0CyacabLJKDUp9tLveqsISKOjog7IuKOrq6umgWnxnjvNv0TLgC/fesH6/JFfNKol5k+sXqvq+kTu8o+iJeKYaC2lmy6Lb996wf7lQ80bpX0M5BSB8lHZr6aRbsU+ctiLfGVUiruaeNeKZkgKRXbQM/3QO2VUq0vCrVQ7ee7Wv0PpF6x1VMl79Vy26q0vWoa6Pl+YbvZ3Nl+YL/yRse8NqXG+872A3lhu9lF9yn1WLebkT/eFZOb946iyYvOTbcmN+8dxXc68MCyEi5Q2furVAylYobKPnv8vKifgca61JhWsk8lMVT7Oa3ncbDY94GhfB8q1l4tvnOU20+pz8VFu+xXMgkxUHvVfEylxnug2Mptq9L2Snl05lb8Y/e3VaWtJjOm5+/3wu3oPtsrygs0I2e6SJIkSZKkunGmS2M9AawfEeOLbHsV8PRwGFhJkiRJklSRlskLNGPS5Q/k49qpd2FETABmA3c0ICZJkiRJklQfLZMXaMaky6VAAj7Rp/wo8uds/U+9A5IkSZIkSXXTMnmBprumC0BEnA0cC1wN/BLYFvgP4BZgbkqpe6D9vaaLJEmSJEnNaW3XdCnUGVJeoFk0a9JlNPmM1tHALOBp8pmuk1JKL6xtf5MukiRJkiQ1p0EmXYaUF2gWTZl0GSqTLpIkSZIkNafBJF1aRTNe00WSJEmSJGnYM+kiSZIkSZJUAyZdJEmSJEmSasCkiyRJkiRJUg2YdJEkSZIkSaoBky6SJEmSJEk1YNJFkiRJkiSpBky6SJIkSZIk1YBJF0mSJEmSpBow6SJJkiRJklQDkVJqdAxVFxHdwIuNjqPBxgBdjQ6iCTgOjgE4BuAY9HAcHANwDMAxAMegh+PgGIBjAI5Bj3qNw8SU0oiYBDKm0QHUyJ9SSjs0OohGiog7RvoYgOMAjgE4BuAY9HAcHANwDMAxAMegh+PgGIBjAI5BD8eh+kZEZkmSJEmSJKneTLpIkiRJkiTVQKsmXc5tdABNwDHIcxwcA3AMwDHo4Tg4BuAYgGMAjkEPx8ExAMcAHIMejkOVteSFdCVJkiRJkhqtVWe6SJIkSZIkNZRJF0mSJEmSpBoYVNIlIkZFxHERcX9EvBQRj0XEmRExuUjdbSJiQUQ8FxHLI+LmiJhbpN7WEfGFiLg1Ip6KiGURcVdEnFis3XLarpUajcM2EfE/EXFfRDwfESsK7X8jIjYuEUfDxqEWY1Bkv0kR8UhEpIj4Tok6LTcGhcdb7PZCifpNPwYRsVNEfDsibomIFwqP59ChtturfiuOwfERcXlEPFyou2QtcQyLz8VyxiGG2fGhRmPQkseGct8PffZtiWNDBZ8Jw+bYUOh/reMQeQdFxE8jYnHh9d0ZET+PiJ0rbbdP/aZ+LVQ4BsPm+FCLMYgWPDZUMAYtd2yo5L3Qp4+mPjYU+q/VZ8KwOj6MaCmltd6AbwEJuAo4CvgGsBL4NTCqV7024BngH8DxwL8Ddxbq7t6nzTOAZcD/AB8HjgEuLfRzNzCxT/1Bt12rW43G4R2F/U8v1DsaOBt4AXgC2KCZxqEWY1Ckj68XXhsJ+E6R7S05BoU2FwEH9bntP4zH4BRgFfBX4JbCPocOtd0WH4NUeFzXA88CSwaoO5w+Fwc9Dgyz40ONxqBVjw1lvR/69NEqx4ayxoBhdGwY7DgAEwp17gROA44APg88DnQDB1U6vs0wDjUcg2FzfKjFGNCCx4YKxqDljg2VvBf69NHUx4ZajgPD7Pgwkm+DeZG8pvBEX9mn/OOFJ/rAXmWXkf8iMbtX2TrAo8ADFC7cWyjfAVi3SH+nFdo9tk/5oNuuyUDVaBwG6O8DhXY/0yzjUI8xAN4IdAH/RekPz5Ycg8L+Fw0yjuEyBhsCkws/78fAf2QOut1WHYNCnS16/XwPA3+pHk6fi+W8FobN8aGWr4US/Q33Y0NFY0BrHRvK/UwYFseGcsYBGAPsVmT/DYGnyf9RMKrcdpthHGo1BoVtw+L4UMPXQcsdGyp5HZTob9geG4YyBjT5saHW48AwOj6M9NtgnqCeD7K39ymfACwHfln4fTLwEnBjkTb+u9DGToPo77WFut/vVVaVtoc0UPUfh50Kdb/cLONQ6zEARgN/BH4BzCr24dnKY9DzwQmMA9YZIIZhMQZF9lvbH9qDbrdVx6BI/ZJfqhs9BvUch177Nd3xoQFjMGyPDZWOAS10bKhkDBgmx4ahjEOfulcW2tioknYbPQ61GoMidZr2+FCvMehVt2WODUMYg5Y5Ngx2DBgGx4ZajwPD6Pgw0m+DuabLjuSzc7f3LkwpvQTcVdgO8DpgPPD7Im3c2quttZlZuP9Hr7JqtT0UNR2HiJgQEetHxMyIeCfwg8KmX/aq1uhxqPVr4Tjg1cCxA8TQ6mOwH7ACWBYR/4yIsyNi3T51hssY1LLdVh2DcjR6DHrar+c4NOPxoaZj0GLHhkq10rGhUsPh2NDT/lDHYSbwCpCrsN1Gj0OtxqAcI20MWvXYUHIMRtCxYaDXwXA4NvS0X8txGC7HhxFtMEmXTYCnU0ovF9n2N2D9iBhXqNdTVqwewKsG6igiRgMnkZ8m9uM+MQyp7Sqo9TgcCTwFPAZcB2TIn7t3c58YKmm7Wmo2BhGxOXAq8IWU0pK1xFBW21VWy9fB7eTP9d8POIT8eZ7HAjdHxDp9Yii37Woa7BjUst1WHYNyY+jpr1gM0Dyfi0PWxMeHWo9BKx0bytaCx4ZKDJdjQ08MFY9DRLyH/H/sLy38QVJJu40eh1qNQbkx9PRXLAZokTFo1WPDIMag5Y8NA43BMDo29MRQq9fCcDo+jGhjBlFnElDsRQL5aUo9dSYVfi5Wt3e9gXwTeDNwQkrpgT4xDLXtoar1OCwA7id/bt0bgL2BGUViqKTtaqnlGHwPeIT8haXWFkO5bVdTzcYgpdT3yuSXRMSfgS8B/1m4771fs4/BKzVst1XHoNwYKBFHs30uVmMcvklzHh9qPQYLaJ1jQyVj0GrHhrLHYBgdG3rar2gcImIrYD75PwA+OYR2Gz0OtRqDcmOgRBytNgbfpMWODYMcgwW08LFhEGMwXI4NPe3XZByG2fFhRBvMTJcV5KcjFTOhV50VhZ+L1e1dr6iI+CL5zNy5KaUvF4mh4rarpKbjkFJ6PKV0Q0ppQUrpZPLZyq9ExPF9Yii77SqqyRhExEHAO4FjUkorBxHDoNuugbq8H3r5GvkP4Xl9YqhG25Ua7BjUst1WHYNyY6BEHM32uTgkTX58qOkYtNixoSwtemyolmY8NvS0X/Y4FP5rfSP56wr8W0rpqSG02+hxqNUYlBsDJeJomTFoxWPDYMeglY8NaxuDYXZs6Gm/np8JzXp8GNEGk3R5gvy0p2JP0qvIT5d6pVCvp6xYPSg+pYmIOIX8klgXkl/+rVgMFbVdRTUfh95SSn8mv4zXv/eJYchtD0HVx6DQ1jfIn4P694jYMiK2BDYr1Fu3UJbpFcOg2q6Rer8OVvb02SeGIbc9BIMdg1q226pjUG4MPf0ViwGa53OxYsPg+FDX18IwPzYMWgsfG6qiSY8NPTGUNQ4RMQtYSP4/9nuklP4yxHYbPQ61GoNyY+jpr1gMMMzHoBWPDUN5HbTKsWFtYzAMjw09MdTztdCsx4cRbTBJlz8U6u3UuzAiJgCzgTsKRX8hP2VplyJtvLlwf0ffDRFxMnAycAlwZEr5yyj3UVHbVVbTcShhIrBer98bPQ61GIOJ5KdDzgMe6nW7qbD9oMLvR1bQdi3U9XVQaHcma14cbriMQS3bbdUxKEejxwBqPA7D5PjQiNfCcD02lKNVjw1V0aTHBihzHCJiM/J/WKxL/g+LO6vQbqPHoVZjUI6WHoNWPDZU6XUwrI8NgxyD4XZsgDq/Fpr4+DCyrW15I/LLsA20tvhBvcouJ7/+9+t7lfWs//0gfdb/Jn/hq0T+Q3PANejLbbvat1qNAyWWgQPaC23c2Ke8YeNQizEAxpK/+FPf20cLbf5f4fetW3UMCuXTS/T3tUK7nxmOr4M+29e2ZHRZ7bbiGBSpX3JJ0EaPQa3HgWFyfKjh+6Eljw2DHQNa+NhQ5utg2Bwbyh0H8v+ZfoT8ahw7VnN8h8troZwxKNJP0x4fajkGtOCxocz3QkseGwY7BgyzY0ONXwvD6vgw0m89f/gOKCLOJn/O5NXkp3NtC/wHcAswN6XUXai3JfmrKK8EzgKWAkcVXmzzUkrX9WrzY8B3gE7y64N39+n2Hyml63vVH3TbtVKjcbga2Jj81aYfJX9e3ZuAA8ifWzcnpXRXr/oNHYdajEGJfmaR/9D5bkrp2D7bWm4MIuIs8pnmheTfE+sA7yF/EL0NaE8pvdir/nAZg82ADxd2ew351/VV5KfAAsxPKT1abruFuq06Bh/mX9NkPw6MA84s/P5oSml+r7rD6XNx0OMw3I4PNRqDVj02lPV+KNLPLIb/saGc18GwOjYUYljrOETEFOBuYHPgbPospVpwfUpp9X9qW+34UOEYDJvjQy3GoBWPDRWMQcsdGyp5LxTpZxZNemwoxFCL18KwOz6MaIPJzACjyV8x+QHyU5P+Rv58unWK1N0W+Bn5DN0K4LfA7kXqXUQ+C1fqdlOlbdfqVqNx+CDwv+SXfHsJeJH81cjPBrIl4mjYONRiDEr0M6vwOvjOSBgD4L3kl/z7W+F1sBy4CzgBmDBcxwCYs5b3+ZxKx7aFx+CmAere1ExjUKtxYJgdH2o0Bi15bCj3/VCkn1kM82NDma+DYXVsGOw49HoeR+zxocIxuGmAuje1+hjQgseGCsag5Y4NlbwXivTT00bTHRtq+FoYdseHkXwb1EwXSZIkSZIklWcwF9KVJEmSJElSmUy6SJIkSZIk1YBJF0mSJEmSpBow6SJJkiRJklQDJl0kSZIkSZJqwKSLJEmSJElSDZh0kSRJkiRJqgGTLpIkqS4iYk5EpIg4pdGxSJIk1YNJF0mSVDURMauQWLmo0bFIkiQ12phGByBJkkaM24FtgacbHYgkSVI9mHSRJEl1kVJaAdzf6DgkSZLqxdOLJElSVRSu1fJI4ddDCqcZ9dwOLXVNl4i4qVA+NiJOioiOiHgpIu6PiKN61TsmIv4SES9GxOMRcWpEFP0uExE7R8QVEfH3iHglIh6LiB9ExCa1evySJEl9OdNFkiRVy01ABvhP4G5gQa9tdxW2DeSnwM7AL4GVwH7AuRGxEngdcAjwC+BGYG/gJGAF8JXejUTEYcAPgZeBnwOPAVsBRwJ7RcSbU0qdlTxASZKkckRKqdExSJKkFhERs8jPdrk4pXRon21zgIXAqSmlU3qV3wTsBtwB7JFSyhXKtyB/OtJyIAe8LaX0t8K2DLAYSMDGKaWuQvnWwD1AJ7BbT/3CtrnA9cDPU0rvq96jliRJKs7TiyRJUrP4XE/CBSCl9DDwW/IzZL7YO4FSqHcNsD7wql5tfBQYC/xn7/qFfX5NfubLXhExpTYPQZIk6V88vUiSJDWLO4qUPVG4/2ORbT1JlZnAo4Wfdync7xYROxbZZwNgNLB1iTYlSZKqxqSLJElqCiml54sUdxXuB9o2tlfZ9ML9p9fS3TplhCZJklQRky6SJKmV9CRn1k0pLW1oJJIkacTzmi6SJKmaVhXuRzeo/1sL929vUP+SJEmrmXSRJEnV9Bz5FYWyDer/O+SXmz6rsJLRGiJiXESYkJEkSXXh6UWSJKlqUkovRMRtwNsj4n+AB8nPfvl5nfq/PyIOBy4A/hoR1xZiGEs+EfR24Cng1fWIR5IkjWwmXSRJUrV9GDgLeDfwISCAx4El9eg8pfSjiLgb+CTQDrwTWE5+JaQrgEvrEYckSVKklBodgyRJkiRJUsvxmi6SJEmSJEk1YNJFkiRJkiSpBky6SJIkSZIk1YBJF0mSJEmSpBow6SJJkiRJklQDJl0kSZIkSZJqwKSLJEmSJElSDZh0kSRJkiRJqgGTLpIkSZIkSTVg0kWSJEmSJKkG/j+iXFzTaF5cMwAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ ""
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# overlay total flash rate (LMA), total CG flash rate (NLDN), \n",
+ "# +CG flash rate (NLDN), and -CG flash rate (NLDN) \n",
+ "\n",
+ "from matplotlib.lines import Line2D\n",
+ "\n",
+ "fig = plt.figure(figsize=(18,8))\n",
+ "ax1 = fig.add_subplot(111) # axis 1 for CG flash rates\n",
+ "\n",
+ "with pd.plotting.plot_params.use('x_compat', True):\n",
+ " lasso_data.plot(drawstyle='steps', color='darkgray', ax=ax1)\n",
+ "\n",
+ "plt.xticks(rotation= 0)\n",
+ "ax1.legend([],frameon=False)\n",
+ "\n",
+ "ax2 = ax1.twinx()\n",
+ "\n",
+ "with pd.plotting.plot_params.use('x_compat', True):\n",
+ " total1.plot(drawstyle='steps', color='goldenrod',\n",
+ " ax=ax2, rot=0,label='total CG')\n",
+ " neg1.plot(drawstyle='steps', color='dodgerblue',\n",
+ " ax=ax2, rot=0,label='-CG')\n",
+ " pos1.plot(drawstyle='steps', color='red', ax=ax2, rot=0,label='+CG')\n",
+ "\n",
+ "# adjust and set the location of ticks and ticklabels\n",
+ "ax2.xaxis.set_major_locator(\n",
+ " mdates.MinuteLocator(byminute=[0, 10, 20, 30, 40, 50]))\n",
+ "ax2.xaxis.set_major_formatter(mdates.DateFormatter('%H%M'))\n",
+ "ax2.set_xlabel('Time (UTC)', fontsize=25)\n",
+ "ax2.set_ylabel('CG flashes (per minute)',\n",
+ " rotation=90, labelpad=10, fontsize=25)\n",
+ "\n",
+ "legend_elements = [Line2D([0], [0], color='goldenrod', lw=4, label='total CG'),\n",
+ " Line2D([0], [0], color='dodgerblue', lw=4, label='-CG'),\n",
+ " Line2D([0], [0], color='red', lw=4, label='+CG'),\n",
+ " Line2D([0], [0], color='darkgray', lw=4, label='total LMA')\n",
+ " ]\n",
+ "\n",
+ "ax2.legend((legend_elements[0], legend_elements[1],legend_elements[2],legend_elements[3]),\n",
+ " ('Total CG', '-CG', '+CG','Total LMA'),\n",
+ " bbox_to_anchor=(0.48, 0.5, 0.5, 0.5), prop={'size': 18})"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": []
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "vortexse",
+ "language": "python",
+ "name": "vortexse_analysis"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 3
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython3",
+ "version": "3.8.8"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 4
+}