#&Lc@sddlZddlZddlZddlTdZdZdZdZdZeZ eZ eZ dZ d Z d Zd Zd Zd ZdZdZdZdZdZdZejZejZdZdZeZdZdZ dZ!dZ"dS(iN(t*g?cCsdS(Ng((twavetT((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pyt NoContinuumVscCs|dkrV|dkrVtjdd|d|dd|dd |d }n\|d kr|d krtjd d|d|dd|dd|d }nd}d|d|S(Ng9@g0@g~j!goKS8?gf,Nig>igh@g@gQȀgBfj?g# H?gGQ>gbN+=ggr@g333333?(tmathtexp(RRtkappa((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytCO2ForeignContinuumZs88cCs|dkrV|dkrVtjdd|d|dd|dd |d }n\|d kr|d krtjd d|d|dd|dd|d }nd}dd|d|S(Ng9@g0@g~j!goKS8?gf,Nig>igh@g@gQȀgBfj?g# H?gGQ>gbN+=gg?gr@g333333?(RR(RRR((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytCO2SelfContinuumhs88cCs|dkrV|dkrVtjdd|d||d|dd|d }nr|d kr|d kr|d }tjd d|d|dd|dd|d d|d}nd}|d|dS(Ng@@g@g/$U(@gzM?g{ `?gor>igc=igh@gp@g@gʡg._MT%agVZҩ>igYPgj靫+=g+\ga4=iggr@g@(RR(RRtkappaCtx((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytH2OSelfContinuumvs8 DcCs-dtjtjtjd|tj|S(NgY@(RtpitphystctB(twavenumR((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytPlancksc Cs[t|}g}|j}g|D]}d|kr%|^q%}g|D]}d|krJ|^qJ}t|}xt|D]}||}||} tj|jddtj|jddg} t} tj|jdd| _tj|jdd| _ d| _ t j ||| _ || | _|j| qW|S(NtdHtlogKappat.iitExponentialSum(t readTablet listVariablestlentrangetstringtatoftsplittDummytnu1tnu2tbandTypetNumericRRRtappend( tfileNameRtbandListtheaderstht binHeaderstlogKappaHeaderstnHeaderstith1twavest bandParams((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytloadExpSumTables$  %%   8   cCs|j|jd}t|d|t|d|}tj| }|t|}tj|j |d}|t |j | tj| |j | S(Ng@gr@g.@( RRtForeignContinuumt SelfContinuumRRtSelfRatR!t searchsortedRtsumR(tpathtpathqR-tnut continOpacitytcontintpathTottimax((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pyt TransEsumss $cCst|}|j|j}|j|jd}tjtj}tj|tj} tj|tj} xYt |D]K} d} d} d}xt | |dD]}t |||dt |||}||||d|||}t |ddkr+t |||t |d}n$t j t d||dt}t|ddkrt|||t|d}n$t j t d||dt}| |d|||t||t||d||tt} |||||t||t||d||tt}t| ||}| |||7} |dkrPqqWt ||t ||dt| ||}| t ||| ||| | } d} d} d}xt | d dD]}t |||dt |||}||||d|||}t |ddkrt |||t |d}n$t j t d||dt}t|ddkrt|||t|d}n$t j t d||dt}| |d|||t||t||d||tt} |||||t||t||d||tt}t| ||}| |||7} |dkrPqqWt ||d  t| ||}| t ||| ||| | g?g-C6?iiii(RRRtGHGtMolecularWeightt BackgroundGasR!tzerostFloatRRR/RRtTstartTrefR0tpreftabstgt cosThetaBartTrans(tpRtTgtqR-tntDeltaRtmassrattIptImR*tSumR4R5tjtdBtmoleContwTFtwTStTrans1tBddTermtfluxtheattdelPlustdelMinustARtC((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pyt radCompBandsx ("#$#$HD 0'("#$#$HD $'  0 cCst|}|j|j}|j|jd}tjtj}tdtjt d} | |t |t } t j |t j } | | dg?iii(RRRR<R=R>tLstellarR tsigmatTstellarRR!R?R@RR/RRRARBR0RCRDRERFRG(RHRRIRJR-RKRLRRMt OrbitFacttswFluxRXR4R5RQRSRTRURVRYRZR[R\RR]((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytstellarRadCompBand@sB  "#$#$HD 0 c Cst|}tj|tj}tj|tj}xRtD]J}|jtkr=t|||||\}} ||7}|| 7}q=q=W||fS(N(RR!R?R@tbandDataRtLWCutoffR^( RHRRIRJRKRXRYtbandtfluxBandtheatBand((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pyt radCompLWs   c Cst|}tj|tj}tj|tj}x@tD]8}t|||||\}} ||7}|| 7}q=W||fS(N(RR!R?R@ReRd( RHRRIRJRKRXRYRgRhRi((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytradCompStellars   cCsd}|dd|}d|tjdd|tj}t||d}|||dd}d|tj|}|tjdtj}t||||\} } tdj } tdj } t t } | d | g|| | d g|}| d|S( NgY@g?gư>i(g@g@iigMbP?gj@( R<R=R>tsetpLinR!tonesR@tradCompReRRtrombergR(tpsAirRIt GHGmolarcontptoptpstMbarRHRRJRXRYRRtmtflux1((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pyt OLRDryAirs$   (cCs@gt|D] }|||||d^q }tj|S(Ni(RR!tarray(RsRrRKR*RH((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pyRls3cCsJ||d|d}gt|D]}|||^q#}tj|S(Ng?i(RR!Rx(RsRrRKtratR*RH((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytsetpLogs'cCsd|d}||}gt|D]$}|d||||d^q%}gtd|dD] }d|d|||^qc}||}tj|S(Niig?ig?(RR!Rx(RsRrRKtn1tn2R*RHtp1((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytsetpExtraGroundRess  7: g#@g@g@p@g?g.AgX@g@g@g@(#RRR tClimateUtilitiesRFRRRR tH2OForeignContinuumR/R0RR.R;R^RdRjRkRwRlRzR~REtairR>tCO2R<RCRBRGR1RfRaR_(((s/Users/rtp1/Havsornen/PlanetaryClimateBook/ClimateWorkBook/CoursewareSandbox/PlanetaryClimateCourseware/ChapterScripts/Chapter4Scripts/miniClimtFancy.pytKs@            p A