Inverse Solution その2

昨日の雨量のエラーは、マクロのミスでした。変換したデータをチェックすると、重複した時間のデータが1か所あり、削除することですんなり動きました。

修正したら計算は早い！
θr、α、n、Ks の4つを同定する設定で、答えを出すまで2分弱。さすが1次元です。ちなみにPCは2年前の Dual Core でそれほど良いスペックではありません。3次元はストレスになりますが、1次元はストレス解消になりますね。
今回は事前に実施していた Instantaneous profile method によって、Van Genuchten モデルではKsが室内透水試験結果の10倍くらいないと合わないことが分かっていました。現場データを使用した今回の逆解析結果も似たような傾向を示しています。
ある程度の妥当な範囲を計算前に設定しておかないと、答えが出ないことも再確認しました。

とりあえず、問題なく計算が回ることはわかりました。明日からはモデルを層分けし、観測データを増やし、より詳細な条件にて計算してみましょう。