學習記錄十:VIC模型—全局參數(shù)文件的制作及陸面模型的運行
1.建立文件夾,將①"官網(wǎng)的小例子\params\vic\global_param.STEHE"(文件后綴改為“txt”)+②"官網(wǎng)的小例子\params\vic\veglib\LDAS"+③制備好的"土壤參數(shù)文件"+④"植被參數(shù)文件"(文件后綴改為“txt”)拷貝至param文件夾中;將制備好的"氣象驅動數(shù)據(jù)"拷貝至forcing文件夾中,result為輸出文件存儲位置。

2.修改全局參數(shù)文件:打開“global_param.txt”,修改輸入數(shù)據(jù)路徑、起止時間、輸出數(shù)據(jù)路徑等。(注意:文件路徑中的"\"均需改為"/")




3.使用Cygwin軟件,運行VIC模型

①將VIC模型源代碼"VIC_4.1.2.g",拷貝至“vic”文件夾中。

②修改"VIC_4.1.2.g"文件夾中的"Makefile"文件。

③打開Cygwin,進行編譯,輸入make后回車運行即可。

④將Cygwin運行出的結果,拷貝至Cygwin安裝軟件目錄bin里


⑤打開Cygwin,繼續(xù)編譯。

Error1:?layer 0 mineral bulk density (0.000000) must be > 0

原因:土壤參數(shù)文件錯誤,但代碼中6對應的bulk density本來就是0,而vic模型要求此指標值必須大于0,重新制備土壤數(shù)據(jù)也依然如此


解決方法:將6對應的bulk density值改為了0.1(治標不治本,因為這個錯誤停滯了一周,實在不想繼續(xù)耗時間了,所以直接粗暴的改了數(shù)據(jù),不知道這樣改會不會對后面造成較大影響,但模型確實能運行了,后期如果結果有問題再糾正吧)
(備注:代碼中的參數(shù)是教程作者自己算出來的,說可以直接用,但至于為什么會與vic模型有如此出入,我也還沒搞懂,目前還沒找到真正原因和解決方法,后期有時間自己算一遍參數(shù),好像還挺復雜的......)
Error2:HHY_global_param.txt文件中輸出文件路徑錯誤,多了引號" ",刪除即可


終于,運行成功了!?。。。。≈鹁W(wǎng)格輸出,可能有些慢,耐心等待即可!??!




太難啦,真的太難啦,每一步走的都很不容易,一個小問題,對我來說就是大問題,我需要很久去理解去解決,今天從15點到24點不吃不喝,就耗在上面兩個小問題,最后還是請了外援才得以解決,其中一個問題還沒有徹底解決。代碼順利運行不報錯的時候,我真的要哭出來了!?。IC模型真的好折磨人?。?!有沒有大佬帶帶我啊,求求了!?。?/p>