基本概念可參考上一篇文章 https://learncropmodel.blogspot.com/2025/05/blog-post.html 上一篇文章提及碳排係數在Tier1當中的拆解方式如下,我們來看看文獻裡的參數吧 $EF_i = EF_c \cdot SF_w \cdot SF_p \cdot SF_o$ EF i : i條件的每日排放係數 EF c : 無施用有機肥料下持續灌溉水田的排放係數 SF w : 不同水分管理區域的調整因子(表 5.12) SF p : 不同種植前水分管理調整因子(表 5.13) SF o : 有機質肥料型態或施用量之調整因子(Eq. 5.3, 表5.14) Baseline emission factor 也就是EF c 的預設係數來自於表5.1 (page 5.53) 表5.1 預設水田甲烷排放係數,假設水稻栽培時間小於180天、連續淹水、未施用有機質 World Regional Emission factor (kg CH 4 ha -1 d -1 ) Error range (kg CH 4 ha -1 d -1 ) Region Emission factor (kg CH 4 ha -1 d -1 ) Error range (kg CH 4 ha -1 d -1 ) 1.19 0.80 - 1.76 Africa 1.19 0.80 - 1.76 East Asia 1.32 0.89 - 1.96 Southeast Asia 1.22 0.83 - 1.81 South Asia 0.85 0.58 - 1.26 Europe 1.56 1.06 - 2.31 ...
工作區 (workspace) STICS模式模擬大部分所需要的資料都會在工作區當中,但是植物參數則是在plant的資料夾當中。STICS的資料使用xml編輯,只需要用notepad就可以打開閱讀。 建立新工作區 打開JavaStics.exe,選擇File → New workspace 建立新的工作區,就可以看到玉米的範例maize。 總共包含了11個檔案 氣象日資料檔的附檔名為年份 (.1996)氣象站的資料為_sta.xml,初始狀態檔為_ini.xml、田間管理為_tec.xml、如果有實測資料則可建立.obs檔、所有的土壤資料都在sols.xml、模擬的資料再usms.xml、輸出資料的設定檔在var.mod當中 模擬檔案都在usms.xml當中、土壤資料都在sol.xml當中,每個working space只會有一個usms.xml 和sol.xml。 建立氣象檔 STICS模式內部有建立氣象資料的工具,可以讀取CSV檔之後直接生成模式所需要的氣象檔。 選擇 Model input → Local parameter → Climate → Formatting a climate file,即可打開氣象資料產生的工具檔 經過測試,可輸入的檔案為CSV氣象資料,但是STICS只能選擇分號或者空白分隔,無法讀取逗號分隔的檔案,目前選擇的方法是先產出CSV檔,再進入將逗號取代為分號,或許可以使用python 或R語言直接匯出符合使用的csv檔。 輸入的檔案為逐日氣象資料,我們可以將所有年度的氣象檔放在同一個CSV底下,STICS模式會自動產出每一個年度一個檔案,因此建議可以在STICS模式的資料夾當中另外產出一個專門存放氣象CSV的資料夾,再透過小工具產出程式可以讀取的氣象檔。 我們在這裡使用的範例如下 再依序的將每一行對應的內容填入視窗中,氣象資料的第一列是表頭,因此我們就設定number of headlines 為1。 依序填完後就可以按generate,程式就會自己產出氣象檔 批次化模擬 JavaSticsCmd.exe --run iceberg JavaSticsCmd.exe --run iceberg L1-M572_17J L2-M780_17J