2012-05-01から1ヶ月間の記事一覧

SW

☆A New Hope 1.Wind shake the stars.○ 2.Point of Origin. 3.Black Knight, White Princess and Pawns 9.Rogues, Rebels and robots 10.The Luke Skywalker Initiative 11.The Jedi Nexus 12.The Case for Relbellion 13.The force and counter force☆Empir…

MLと強度の関係

H:\jobs\2010mpx\mkgrads\malti_elv WP2Dm.xls WR2Dm.xls 他にもある。Zと高度差の関係はエクセル不明。

malti−上端と下端を調べる。

上端、下端がそろっていると言うことは意味がある。ドップラ速度を見ることができると良い。(つぎの課題)

malti

dhと層厚の関係を調べる。 2010/10/10で1:45には値が記録されている。 幅の記録は正しいが、高度がRで20、Pで0となっている。原因究明。←直接高度を記録していない。Rと思ったのは高度差、Pと思ったのは幅の差。

malti

高度の計算をしているので、相関を出力する。欲張らない。ρhvの強度は次の課題。 ◆gradsのLineは計算を間違えているとしても、出力を間違えない。必要な値は、PMLとRMLがともに観測された場合。 ◇ih_2ddr2>0.and.ih_2ddp2>0を出力。

malti Zh_max高度と0度高度の確認

9時21時のみ確認。 Zmaxと0度高度は抽出済み。Zmax高度を出力。

malti-確認

0228:比較的きれいにBBが確認できる。 たとえば201102281100 は比較的強い値が出ているが、RMLの範囲が811.5 511.5 となっている。←要調査。 ★228を描画する予定なのに、BBが4000くらいに出ている。読み込みの値がおかしい。 133-135を対象としているはず。…

malti-画像

まとめる。たぶん整理していない。

malti-bb_csv

とまっているので、進める。 計算はできたはずなので文章を修正。⇒雪が降っているときはPMLが定義できていないはず。 02/28 03/07 03/21 この3日は除く。除くことが適正かどうかは、TH図を見て判断する。

malti TH図の整理からcsvへ

計算結果を整理。grdデータについて2011はどのように表記するか? NG=TH図で時間軸が変更になっていない。←修正。 bb_csv 2011エラーあり。←pathnameにCCYYを挟むときは、元の文字列の長さを変えない。

malti PMLの検出

閾値がおかしい可能性あり。 kmarkは番号を1要素増えたことにあわせる。←番号をGRDファイルの変数入力並びと同じにする。 RMLの検出を合わせる。★ ☆リスト:H:\jobs\2010mpx\mkgrads\elv8\mlh_zdr_d100.xls H:\jobs\2010mpx\mkgrads\malti_elv\mlh_zdr_d100_…

malti bb_csv

出力項目を4つにしたので、プログラムを修正する。 cal2d/malti/RunMean -- shin??_12s.dat 2011と2010を分けていない。 @9/22'14 YYMMDDを記録するように修正。 (3:~/jobs/2010/mpx/bb/cal2d/malti/RunMean) MMDDs.grdも同じところに記録しておく。elv…

malti次の計算

N:の計算と M:の結果の整理 他のドライブの計算のために! 「F:\jobs\2010mpx\mkgrads\malti_elv」 ⇒「Z:\Pub\user\koshida\MPX\mkgrads\malti_elv」 80c[C:\home\jobs\2010\mpx\grads\malti_elv] ■:d100「c_go.x」の実行elv12_scanave ⇒/f/jobs/2010mpx/mk…

malti計算を実施

M:から ■注意点:TMPの下のキャッシュを消す。 ↑debug=.false.でキャッシュを消す。debug=.true.であれば、再計算のために解凍したデータを残す。@11/29'13 ■複数日動くか?

bb_csv malti

1030.grdからデータを切り出したときに、grads描画用に、top&bottomを記載している。そこでの値が、RMLでひとつ、PMLで2つずれている。原因を調査。。また、PMLの高度は最小ρhvで定義するのが良い。とすればMMDD.grdデータの作成を修正する必要がある。 ■計…

d100実施後の1030.grdを確認

なぜ、エラーが出るか? ■F:\jobs\2010mpx\mkgrads\malti_elv\mk100\d100data div_write.gs確認。 >ひとつにはデータがt=1〜288であるのに、set t 1 289としているから。 run F:\jobs\2010mpx\mkgrads\malti_elv\mk100\d100data\div_write.gs 1030 でローカ…

cal2d malti

5:20@2010/10/30のρhvが0となっている。elv8の計算では0.25.きちんと整理しよう。 >直前の計算ではd100で実施したので、elv#12もd100で計算する。←保留。 ◇ F:\jobs\2010mpx\mkgrads\cal2d\malti(確認のみ) >> check /home/kisho/koshida/jobs/2010/mpx/b…

elv#12の確認ML malti

elv8と比較。当たり前だがsurfaceは同じ。PMLやRMLの高度がどうなっているか? ρhvの値がおかしい。Zはそれほど違わないが、ρhvの値が異なる。

cal2d on F

Fから2次微分を計算。スクリプトがUNIX用なので書き換えが必要。 #libraryの場所が違う。 #gradsのディレクトリが違う C:\Program Files\PCGrADS\win32e\Grads.exe に変更 ⇒perlのなかでは、\が理解されていない。空白なしで記述すべき。 どうしてもProgra…

MLの検出

記憶では、移動平均を実施。スキップするか? k2010/9/1の日記はsyakaのもの。←PMLで検索 XMPは?elv8:http://d.hatena.ne.jp/koshida_tomoki/20110211 記憶があいまいだが、100m移動平均を実施したのはあくまでTH図あるいは鉛直プロファイルを描くためにgrad…

MLmalti 1日の統計値

見た目も重要であるが、検出される統計量を比較した方がわかりやすい。 目標:館野との0度高度比較。5分値の比較。 そのまま比較しても違いがわからないから計算結果を示すのが良い。

malti ML

BB内のZ強度と地上強度は明らかに関係がある。さて、層厚、層の高度はどうであるか?

ML malti

この図自体はよく見るもの。違いを示すためには、意味を良く考えること。 思い出したこと;ゾンデが存在する9時21時に限定する。と言う考えもある。 図を縮小することも大切。 見た目がよいと言うことは大切。2010年10月30日21Z(6時JST) 論文では使わないか…