Tool 07

ET models / 蒸散模式

作者：方煒，臺灣大學生物機電工程學系名譽教授 | 網頁實作協作：ChatGPT 5.4 與 Codex

依據上傳的 MATLAB 程式 etmodel.m 整理而成，保留原始程式中的空氣狀態推導，並比較 8 種蒸散模式的結果。

原始檔 / Source files etmodel.m
Workbook sheets 2
ET models 8
Current mode Worksheet 1

Worksheet mode

Worksheet 1 / 蝛箸除?€???典?

蝚砌?撘萄極雿”?渡???MATLAB 隞??16 ???賂?? 9 ?撓?亥? 7 ?撠€潘??其??舀?敺 8 ???芋撘€?/p>

霈隤芣?

  • Tdb嚗嗾?澈摨佗?簞C
  • RH嚗撠?摨佗?%
  • U嚗?撅支???2 m ?◢??m/s
  • Tl嚗?皞恬?簞C
  • LAI嚗??Ｙ??
  • CAI嚗?撅日蝛??唳?Ｙ?瘥?/li>
  • Rn嚗?撅文?嗡?頛餃???嚗/m簡
  • G嚗?撅文?梢€?嚗/m簡
  • rw嚗?銵函征瘞?????s/m
  • rs嚗除摮??s/m
  • 弇嚗偌?掃???梧?J/kg
  • ?嚗征瘞??摨佗?kg/m糧
  • es嚗db 銝?憌賢??豢掃憯?Pa
  • e嚗征瘞?瘙賢?嚗a
  • 帤嚗ˊ?瘙賢??脩???嚗a/簞C
  • 帠嚗??虜?賂?Pa/簞C
  • 帠c嚗耨甇???虜?賂?Pa/簞C
  • ratio嚗?/ 帠
  • dT嚗?es - e) / 帠 ?€撠?銋澈撌殷?簞C
  • Twb嚗??澈摨佗?簞C嚗BD嚗??澈??簞C

Core equations

公式摘要

  • rs = 140.9 + 703 × exp(-0.013 × Rn)；rs：氣孔阻力
  • λ = 2502535.259 - 2385.76 × Tdb；λ：汽化潛熱
  • ρ = 100000 / 287 / (Tdb + 273.16)；ρ：空氣密度
  • pws = 6.11 × exp(17.4 × Tdb / (239 + Tdb))；pws：飽和蒸汽壓，單位為 mbar
  • es = pws × 100；es：Tdb 下的飽和蒸汽壓
  • e = RH / 100 × es；e：實際蒸汽壓
  • δ = 41.45 × exp(0.06088 × Tdb)；δ：飽和蒸汽壓曲線斜率
  • γ = cp / λ × Patm / 0.6219；γ：濕球常數
  • γc = γ × (1 + rs / rw) + δ；γc：修正濕球常數
  • ratio = δ / γ；ratio：斜率比
  • dT = (es - e) / γ；dT：(es - e) / γ 的溫度當量

此處 cp = 1006 J/kg/°C，Patm = 101325 Pa。頁面亦用數值法求解濕球溫度 Twb，並以 WBD = Tdb - Twb 表示濕球溫差。

Worksheet mode

Worksheet 2 / 璅∪?瘥?

蝚砌?撘萄極雿”?? MATLAB 蝔?銝剔? 8 ???芋撘??銝€?湔撘?頛詨 ET rate ??latent heat flux??/p>

璅∪?皜

  • Penman (Rn, U, VPD), 1948
  • Penman-Monteith (Rn, G, VPD), 1973
  • Penman-Monteith (Tl, Tdb, RH), 1973
  • Stanghellini (Rn, G, Tl, LAI, VPD), 1987
  • Stanghellini? rule of thumb, 1987
  • Fynn (Rn, G, LAI, CAI, VPD), 1993
  • P-M-Fang (Rn, G, WBD), 2003
  • P-M-Fang (Tl, WBD), 2003

Core equations

璅∪??砍?

  • Penman: E = [帤(Rn-G) + 帠繚0.001157繚(1+0.53U)(es-e)] / [(帤+帠)弇]
  • P-M: E = [帤(Rn-G) + ?cp(es-e)/rw] / [帠(1+rs/rw)+帤] / 弇
  • P-M by Tl: E = ?cp(estl-e) / [帠(rs+rw)弇]
  • Fynn: E = [帤繚CAI(Rn-G) + 2LAI?cp(es-e)/rw] / [帠(1+rs/rw)+帤] / 弇
  • Stanghellini: ?誑 radiative resistance 靽格迤 Rn嚗?隞? P-M ??瘙?E
  • Rule of thumb: E = (0.35Is + 30) / 弇
  • Fang (Rn,G,WBD): ?? NTU ?????? LHF ??E
  • Fang (Tl,WBD): 隞?Eff = (Tdb-Tl)/WBD 隞?? LHF 撘?/li>