Excel to Web Prototype
Developed by 方煒 臺灣大學生物機電工程學系 名譽教授｜Web implementation assisted by ChatGPT 5.4 and Codex
本頁聚焦於與 Pn simulation 直接相關的兩張工作表：Pn_Simulator 與 Pn_24 h_Simulation。
Worksheet View
Light Factor
第一張工作表整理光反應因子 GL。它以 Ki 與 PAR 為輸入，用來描述光照對最大光合速率的限制程度。
GL = 1 / (1 + Ki / PAR)Ki : irradiance rate constant, W/m2PAR : photosynthetically active radiation, W/m2P = Pmax x GC x GL x GTLCO2 Factor
第二張工作表整理 CO2 因子 GC。它先把 CO2 ppm 轉換為葉綠體中的 CO2 濃度，再用 Kc 與 CO2ch 求得 GC。
CO2ch = CO2ppm x 44.01 / 22.4 x 273 / (273 + T)GC = 1 / (1 + KC / CO2ch)KC : CO2 concentration rate constant, mg/m3CO2ch : chloroplast CO2 concentration, mg/m3Temperature Factor
第三張工作表整理溫度因子 GTL。它以葉溫 TL、最適葉溫 Tm 與溫度反應常數 a 估算溫度對光合速率的限制程度。
GTL = [2(TL + a)2(Tm + a)2 - (TL + a)4] / (Tm + a)4GTL,bound = max(GTL, 0)TL : leaf temperature, °CTm : optimum leaf temperature for Pn, °Ca : temperature response constant, °CPhotorespiration
第四張工作表整理光呼吸 Rp。它依植物型態、溫度與 CO2 濃度，估算 CO2 條件對光呼吸的抑制程度。
CO2mg = CO2ppm x 44.01 / 22.4 x 273 / (273 + T)Rp = IF(Ptype = C4, 0, IF(CO2mg ≥ 2352, 0, (-Rpmax / Rhocstop) x (CO2mg - Rhocstop)))Rpmax = 0.5, Rhocstop = 2352 mg/m3Dark Respiration
第五張工作表整理暗呼吸 Rd。它使用 Q10、Rd20 與溫度 T，估算溫度升高時暗呼吸的指數增加。
Rd = Rd20 x Q10((T - 20) / 10)Q10 : respiration temperature coefficientRd20 : dark respiration at 20°C, mg/m2/sMain Model
這一頁是單點模擬器。給定葉溫、相對濕度、PAR、風速與 CO2，可估算光合速率、光呼吸、暗呼吸、淨光合速率與蒸散相關阻力。
GTL = [2(TL + a)2(Tm + a)2 - (TL + a)4] / (Tm + a)4Tm = 25, a = 5GL = 1 / (1 + Ki / PAR)GC = CO2ch / (KC + CO2ch)Pm = Pmax x GTL x GL x GCA = RcB = CO2mg + KC + Rc x PmC = CO2mg x PmPs = (B - sqrt(B² - 4AC)) / (2A)Rc = Rac + Rlc + Rlc,incCO2mg = CO2 x 1.83025Pn = Ps - Rp - RdTr = VDD / Rv × 1000Profile Model
第七張工作表把 24 小時的溫度、相對濕度、風速、PAR 與 CO2 資料逐時帶入，計算逐時蒸散、光合速率、光呼吸、暗呼吸與淨光合速率。
Worksheet Switch
Workbook Reference
The first sheet organizes the intermediate variables behind the transpiration and photosynthesis model so the later two sheets are easier to interpret.
GTL = f(TL)GI = 1 / (1 + Ki / PAR)Rv = Rav + Rlv + RvcTr = VDD / Rv x 1000Variable Map
| Symbol | Meaning | Unit |
|---|---|---|
Tl | Leaf temperature | °C |
PAR | Photosynthetically active radiation | W/m2 |
Rav | Aerodynamic resistance | s/m |
Rlv | Leaf vapor resistance | s/m |
VDD | Vapor density difference | g/m3 |
Tr | Transpiration rate | mg/m2/s |
Sheet 1 reference is ready. Adjust the variables above to inspect the resistance structure.
Derived Values
| Symbol | Meaning | Unit | Value |
|---|
Reference Visual
Places Rav, Rlv, Rvc, and Rv in one view so the transpiration resistance structure is easy to inspect.
Inputs
GL worksheet is ready. Adjust Ki and PAR to inspect the light factor.
Results
Reference Curve
Compares the workbook Ki = 100 and Ki = 200 curves, and marks the current Ki/PAR condition.
Inputs
GC worksheet is ready. Adjust temperature, CO2, and Kc to inspect the CO2 factor.
Results
Reference Curve
Shows the GC response curve for the current temperature and Kc, and marks the current CO2 point.
Inputs
GTL worksheet is ready. Adjust TL, Tm, and a to inspect the temperature factor.
Results
Reference Curve
Compares the workbook a = 5 and a = 10 curves, and marks the current TL/Tm/a condition.
Inputs
Rp worksheet is ready. Adjust plant type, temperature, and CO2 to inspect photorespiration.
Results
Reference Curve
Shows the Rp response to CO2 for the current temperature, and compares C3 with C4.
Inputs
Rd worksheet is ready. Adjust Q10, Rd20, and temperature to inspect dark respiration.
Results
Reference Curve
Shows the exponential Rd response to temperature and marks the current temperature point.