j1dZddlmZddlZddlmZdZgdZgdZdZ d Z d!d Z d"d Z d#d Z d$d Zd%dZd&dZd'dZd(dZd)dZd*dZd+dZd,dZd-dZd.dZd/dZd0dZd1dZd2dZd1dZd3dZd4d5dZd6dZ d4d7d Z!y)8a}Python port of packages/engine (the CMAH guidance engine). Must produce output identical to the JS engine for every vector in packages/engine/vectors/cases.json (see tests/test_engine_vectors.py). The server uses it to recompute runs on finalize (never trusting the browser), to backfill suggestions for imported forecasts, and to replay proposed config versions against history. ) annotationsN)Anyz1.0.0)north northeasteast southeastsouth southwestwest northwest)uppermiddlelowerc2tj|dzS)z5JavaScript Math.round (half rounds toward +infinity).?mathfloorvs '/var/www/avy-guidance/api/app/engine.pyjs_roundrs ::a#g cD|ytj|dzdzdz S)Ng@rrrs rround4rs&y ::a#gm $s **rc||dk(st|try tjt | S#t t f$rYywxYw)NF) isinstanceboolrisnanfloat TypeError ValueErrorrs ris_numr&"sKyAGz!T2::eAh''' z "s9A  A c2t|r t|SdSN)r&r#rs r_numr)+say58*d*rc.t|t||Sr()maxmin)rlohis rclampr//s r3r1: rct|jjj}t |dk7s|dt vs |dt vry|d|dfS)Nrrr)strstriprsplitlenASPECTSBANDS)locpartss rparse_locationr96s[ HNN  " " $ * * ,E 5zQ%('1U1XU5J 8U1X rctt}|xsgD]$}t|}|s|j|d&|S)Nr)setr9add) locationsoutr7ps rbands_ofr@=s? %CB 3  GGAaDM Jrcnt|tr |dks|dkDryt|dzdt|dzS)Nr )r intr5r6)cells rcell_to_locationrGFs@ dC D1Hr dQh % "2!3 44rct|}|sytj|ddztj|dzS)NrrCr)r9r6indexr5)r7r?s rlocation_to_cellrJLs<sA  ;;qt q 7==1#6 66rc |d|d |d tt|dt dz }tt|dtdz }tj|}t t dz |dz}||z }tj|}t tdz |dz}||z }d fd } | ||d|z z| |||zz} | ||d|z z| |||zz} | d|z z| |zzS)N sensitivity distributionlikelihoodMatrixrrc~tj|dxsij|dxsdS)Nkeyr)r#get)disidistmsenss rrFzinterp_likelihood..cellcs<aeeDHUO,277RHMANNr)rRrErSrEreturnr#)r/r#r4rrr,)dist_possens_posconfigd0d1fds0s1fsrFbottomtoprTrUrVs @@@rinterp_likelihoodrcVs - D . !D!"AU8_aTQ7HU8_aTQ7H H B SY]BF #B BB H B SY]BF #B BBO"b\QV $tB|b'8 8F r2,!b& !DRL2$5 5C QV sRx ''rc z|jdxsi}|jdxsi}td|jd|jd|jd|jdfDsyt|d|d|}t|d|d|}t||t ||dS)NrLrMc32K|]}t|ywr(r&.0xs r z#likelihood_range..ns[Qvay[r,r+r,r+)rQallrcr,r+)problemrZsdr-r.s rlikelihood_rangerqks M"(bA N#)rA [155<uquuU|QUUSX\"Z[ [ 1U8QuXv 6B 1U8QuXv 6Br2;s2r{ 33rc|sgStttj|ddd}tttj|ddd}t t ||dzS)Nr,rr+)rEr/rrceillistrange)rr-r.s rlikelihood_optionsrxus`  U4::ah'A . /B U499QuX&1 - .B b"q&! ""rc  t||}|sy|jd}t|r7tt t t |dd}|t|vr|Stt t |dddS)NlikelihoodChoicerrsr+)rqrQr&rEr/rr#rx)rnrZrwchoicecs rlikelihood_forr}}sv&)A  [[+ ,F f~ huV}-q!4 5 "1% %H uXah'A. //rc t|dd}ttdt|t j t |dz S)N dangerTable1rg& .>)r4rEr+r,rrtr#)sizerZcolss r size_classrsA vm$S) *D s1c$ %+*< =>? @@rcx|djt|}|syt|t||dz S)Nrrr)rQr1rEr) likelihoodrrZrows r danger_forrs=   # #C O 4C  s:dF+a/0 11rc Tddddd}tD]}d}|D]w}t||}|jdxsi}|rt|jdsA|t |jdvr^t |t ||d|}y|||<t |d||d<|S)Nr)r rroverallrr+r=r)r6r}rQr&r@r+r)problemsrZr>bandbestr?likers rdanger_suggestionsrsQ1 =C 3 DA!!V,D55=&BDvdhhuo68AEE+$677tZd5k6BCD DD S^T2I 3 Jrc~tt|jdt|jdzS)Nr=)r!r@rQ)abs r _shares_bandrs. {+,xk8J/KK LLrcxt||t||}}|r|sy|d|ddzkxr|d|ddzkS)NFr,r+r)rq)rrrZrarbs r_likelihood_overlaprsO a (*:1f*EB R e95 A % D"U)r%y1}*DDrc t|jdxsi|jdxsi}}td|jd|jd|jd|jdfDsyt|dt|dkxrt|dt|dkS)Nrc32K|]}t|ywr(rfrgs rrjz _size_overlap..s_Qvay_rkr,r+F)rQrmr#)rrsasbs r _size_overlaprs UU6] b!%%-"52B _266%="&&-PRPVPVW\P]"^_ _ E uRY/ / XE"U)4DbQViHX4XXrc Ni}|D]%}|jd||t|d<'g}|jdxsgD]}|jt|d|jt|d}}|r|sBt||t|||t ||d}|jdduxrt |j } |j|dd|dt|dt|d|| |jd d |S) Ntype pairingRulesrr)locationrrrequireJustificationF-source)rPrroverlapstrongr)rQrErrrrmvaluesappend) rrZby_typer?flagsrulepapbrrs r pairing_flagsrs1!G ( 55= $&'GC& N #( E >*0b Sc^,gkk#d3i..IB $R,-b"f=!"b)  01>X3w~~GWCX s)Ad3i[1c^c^" ((8,    ( Lrcg|jdxsg}t|jdxsg}d5d6fd }t|D]g\}}|jdt|dnd}|jdxsi} t | jdt | jd} } t |jdxsijd} t |jd xsijd} | | | | z d kDr |d d d d||dvr| | dkr| | t kr |dddd||dk(r*| | dkr |dd dd|| | t kDr |dddd||dk(r |dddd||jdrXgr\|dd d d|j|Dchc]"}|jdt|d$}}d!|vrd"|vrd#|vr |d$dd%dd&|vr |r |d'dd(d|rtfd)|D}tDcgc] }t|j|xsd*"}}t||r0d"k\r+tfd+|Drdk(rd,nd-}|d.d |d/d0t}|D]}|t|jdz}!tD];}t|j|xsd*dk\s%|s(||vs-|d1dd2|d3d4=Scc}wcc}w)7Nr conditionsc4j|||||dy)N)idleveltextr problemIndex)r)id_rrr problem_indexr>s rr<zguidance_checks..adds #tv_lmnrrrr+r,rLrMr size_rangewarnzXKeep the range of sizes to two or less and describe outliers in the problem description.uFG §3.4)rrg?small_loose_or_storminfozTable 1 lists Dry Loose and Storm Slab when there is potential for ~D1.5 or larger avalanches, or a widespread pattern of small avalanches.z FG Table 1rsdps_sizezODeep Persistent Slab is for slabs capable of high-consequence avalanches (D3+).dps_sensitivityzDeep Persistent Slab is reserved for low likelihood / high consequence scenarios and is normally stubborn to triggers. Temporary increases are acceptable if likelihood is expected to return to low. glide_durationzyGlide timing is highly uncertain: consider keeping the problem for at least 48 hours after the onset of glide avalanches.r= no_locationzSMark the aspects and elevations where the problem is most likely to be encountered. new_snow_24hrfirst_storm_on_pwlzThe first storm on a persistent weak layer should be a Storm Slab unless you expect the terrain-management considerations of a Persistent Slab. no_hazardstructure_not_problemz|Assign a problem only when there is an associated avalanche hazard; poor structure alone belongs in the Forecast Discussion.c3K|]W}t|jdxsijdduxr"t|dddk\xr t|Yyw)rr+Nr)r)rQr#r})rhr?rZs rrjz"guidance_checks..sl  !%%-%2**51 2$ > *ai&'1, *q&) * sAA rc3(K|] }|k ywr()rhlvlrbs rrjz"guidance_checks.. s#@#C#I#@sEXTREMEHIGHhigh_all_bandsz danger with D3+ avalanches should be assigned to all elevation bands unless there is a specific, explainable reason these avalanches will not affect the other elevations.uFG §2.1band_without_problemzThe uu band is rated Considerable or higher but no problem is located there — check the locations or explain in the text.uFG §3.2r() rr1rr1rr1rr1r int | NonerWNone) rQr; enumeraterEr) WIDESPREADSTUBBORNanyr6r+r@)inputsrZratingsrrr<ir?trs_maxs_minsens_mindist_maxtypesbigrlevelsnametouchedr>rbs ` @@rguidance_checksrsCzz*%+HVZZ -34Jo(#M1eeFm7C& NTuuV}"DHHUO,d488E?.Cu}-388?@~.4"99%@A  !2uu}q7H  f'ACMOP Q ;5,AUZbeoZo &1~@LNO P 6 UQYJ(y|HJKL#8(;%v0wyEGHI 6  &+fhtvw xk"b  v'|IKL M'M*&. Kv1JS6^ KE K#U q~ &+|~J Kj X #V.lnx y     5::q#gkk!n)*::&k 3!8#@#@ @ #q9fD  &TF3^+_ak lE 4A xk 23 3G 4 AA7;;q>&Q'1,aw>N*Fd1#>s5tv@A A J1 L;s:LL)%Lc t|jdxsgd|jdxsd}g}|D]b}t||}|jt |d|rt |dt |ddndt |t||ddt||}t||t||ti|d|i||d S) Nr maxProblemsrrr,r+rl)rlikelihoodRangelikelihoodOptionsr) engineVersionr suggestions pairingFlagschecks) rurQrqrrErrxr}rENGINE_VERSIONrr)rrZr per_problemr?rwrs r compute_runrsFJJz*0b12RFJJ}4M4RQRSHK    Q 'AfIYZ6!E(+;F1U8DT#U`d%7%:,Q7     %Xv6K'"%h7!"BV"BZ"BFKX  rcP|xsi}g}t|jdxsgD]\}}|t|dkr|d|ni}|jds|j|ddd|jdxsi}|jddvr|j|ddd|jd rgr|j|d d d|d D]G}|d s |j|dxsdj r0|jd|dddI|S)Nrrz,Set the sensitivity and distribution ranges.)rcoderrr+)NrzSet the expected size range.r=rz#Mark at least one aspect/elevation.rrrPr justificationzrs#  `$    + +5 7(*4# 0A 2(MEY@5v.r