仕事の話をしてみる

してみよう。

最近すっかりモデル屋さんですが。
1. 雲物理の方は、組み込みはどうやら完了している模様（プロットで確認したが…　全く一緒ではないが関数の出力は一緒だからいいことにする）。でも大して予報が改善しない（結果は変わるが見た目、観測とだいぶ異なる）のですが、それはどうやら、前線の後ろ側でできているCold Poolが強すぎるためってことになり、弱めるには雨粒の粒径をおっきくしようということになって、MPのIntercept Parameterを1/10、1/20と小さくしてみるも、確かにCold Poolは温まってきたがまたもや大して変わらず。ちなみに適当にやめてSnowの方も変えてみたが大して…　この結果は明日話す。

2. いろいろありすぎてちょー大変だった、3DVarですけれども、またもやバグが見つかり（というか、並列化したら勝手に観測値をカットしていて、実質入ってるデータがものすごい少なくて大して効かなかった）、巨大なドメインでもCPUが2個以上あるのに1個で同化するしかなくなっています…　これは何とかしたいが、原因がさっぱり分からず。下手に動くので優先順位が下がりそう。ただレーダーを同化しようとなるともちろん、CPUはいっぱい使いたいし、メモリ（というか、スタックサイズ）が足りなくなるなぁと思っていたら、オプションでデータ領域？に割り当てるというのを発見し、とりあえずそれをつけて、Single Threadで動かしてみている。とりあえずVrを同化していますが、さすがに3次元？だからか、iterateを35回してもまだ収束しません（ちなみに、METARだけなら10回で落ち着く）。ちなみに6時間METARを同化し続けた結果、そこそこレーダーの見た目と近づいてきています。もう少しやってから同化、というのもアリだと思うが、もう近かったらわざわざ入れなくてもいいだろうという話もあるし、とりあえずMETARだけの結果を（まだ走り終わっていないので）みてから考える。そういえば、Analysis Incrementsの意味がさっぱり（モデルのおっさん以外みんなそうでした…）でしたが、とりあえず図を出してみて分かって一安心。観測値とFirst Guessの差、ですね。やれやれ。あ、そういえば、観測値にNormalized Gaussianをかけたり、Random Errorを加えたりして同化するのもいいんですね。まだ理論を理解しておりませんができるということはわかったところです。まずは何もなしでやってみているところです。

と、書いてみたところで、大して面白くないですな。図はあるので示したらそこそこ楽しげだと思われますが…