2017-01-01から1年間の記事一覧

作業リスト

1年ではインパクトが弱い。複数年評価するためには、データのDLが必要。 →リストアップ 2004年データはZ圧縮されているので展開が必要。 →全データ必要となるわけではないので、リストアップする。 ○ スペクトルの表示 →精度評価に加えて、スペクトルを評価…

作業リスト

・現象論 ・スペクトル強度の確認 ・

精度評価

実施できるところを続ける−アメダスで実施

過去データの描画MRR

心構えのみ ○TH図を見る理由 →どういう現象が起きているかを記述する ○精度評価 →結果を示す

元に戻せるかやってみる

現状は2度変更(オブジェクトから) ・make clean-再コンパイル ・オブジェクトを戻して−再コンパイル

コンパイルがうまく行かない時はオブジェクトを見直す

make時に参照moduleが見えない理由として、記録済みのオブジェクトを参照していることがある(考えられる)

降水日は1mm/日で。

あまり小さい雨量はブライトバンドが表現しきれていないと考える。

2015〜2017を整理

2013のprgを利用。 - XRAIN missing

mkhyeto

2014~2017の降水日を整理する。

mkhyetoの仕様

ccyymmddを入力し、5分平均の降水量を抽出する。 ★2010,2011を調べた時は、日降水量を計算し1mm以上の降水日を抽出した。 /home/koshida/jobs/2013/prg/get_aws/mkhyeto

ブライドバンド判定に、榊原・武田を利用する

・移動成分と停滞成分を分ける

修正点

・精度評価をきちんと実施すること ・計画が大切。 ・「ブライトバンド」の定義を広げることができるように。 ↑ ρHVで観測する融解層は「ダークバンド」と呼ぶべき現象であるが、”融解層付近に出現する特異な信号”という意味で「ブラントバンド」と呼んでも…

文章の修正

できるところをしっかり inp.tmpが作成されない。

締め切りに向けて

・C合成の評価は追加する。 工程:13日提出(夜) 12日計算 13日文章修正 file:///H:\jobs\2017\doc\171215_koshida_D.doc OGO(H:\jobs\2017\develop)追加

鉛直エコーのタイプわけ(10年)

タイプわけのためのスクリプトが必要。 12/15原稿では削除。 ・同様の観測例を示す。 P29

融解層内Z頭打ちの理由

・早い粒子落下によって融解の効果が遅れる←×:Zには関連しない ・初期粒子があまり大きくない←融解中の最大粒子が大きくない←落下速度が大きいはず ↑上空のドップラ速度が大きければ、傍証にはなる。 ・Zは数の効果と大きさの効果があるが大きさの効果が小…

エクセルのグラフからZR関係を作成する方法

横軸:dBZ 縦軸:Rg 縦軸を対数で表示すると Y=C*exp(D*X) の近似式が得られる。 (C,D:係数、Y:Rg、X:dBZ) 両辺の常用対数をとると ln[Rg]=ln[C]+D*dBZ ln[Rg]=ln[C]+D*10*log[Z]

合わない理由を整理

5dBZの出現の理由。 →⊿Zが大きいほど凝集成長による粒子形成が顕著。 MRRでは、雲粒付着成長と凝集成長の区別がつきにくい。 それほど明確ではない。大粒子がとらえにくい(過小評価する)。

netcdfをgradsで表示

arare C:\home\koshida\jobs\2017\develop\ai\mkinputsdfopenでなく、cntlを作成。 ◆◆ D論:おごの結果を入れる。 10分雨量で整理しているが、1分での比較観測も可能。 アメダスでの効果も評価可能。 ↑ 合成雨量でBBが目立たない理由として2仰角合成が考えら…

ローカルでTH図の描画

mk6.plの実行 確認したい点は、「地表のZが大きいとき融解層の⊿Zが小さい」こと。 20:00に比較的大きな値が観測されている。 ■時間高度断面図 ◎スクリプト ga-> cbar_line -x 3. -y 6. -c 1 2 -m 0 -l 1 -t "20:07" "19:52" ◇ ◇ ◇ 弱エコー時 強エコー時 Zの…

20171122 2000

降雨の強い時はBB幅が大きい。

LIVE(Komaba)の修正

0400のグラフが表示されない。 PS以下のLIVE(過去の確認)を削除。 /export/raid2/koshida/MRR/data/prg/read_mrr2/LIVE 0420に実施。 live04.png ::空白。 o_err:: 'Can't open description file' 原因:sedで200mを100mに変える設定をコメントアウトしたと…

他との比較

比較が大切

Today'sのY軸を修正

sedが問題かもしれない。

作業のリストアップ

10年分計算する ・DTCの計算漏れ確認 ・bb_csv以降←直近 ・比較 ・Zが大きいときのスペクトル #必ずしも雪の話になっていない

解析期間(これまでMRR@東京で観測してきたデータから抽出)

R:~/jobs/2017/xmp/tokyo/mrr/check_mode ①50m内挿 ②BB抽出 ③rg取り出し ④比較compbb

OFFLINE見直し

けりをつける。 11/15'17 delta_h extarct←ihei(1)を利用。sapporoでは実施済みのはず。 Today'sのリカバリ。c_mk24を実施。

ensemble

test arare C:\home\koshida\jobs\2017\yokoyama\grads

OFFLINEMRR

◇repair オフライン操作の際 /home/koshida/jobs/2003/prg/read_mrr2/fig9.png ↑ oneday.pngとしてリンク が更新されることに注意。 ■2/23'19 20150122 repair 20150115 ■ ■20150609 7/17'22 □R:~/jobs/2013/prg/read_mrr2/OFFLINE X:go.pl #sub_script c_m…

作業内容

・平均時間について 5分→10分 ①5分で計算した値を平均する ②プロファイルを10分平均し、特徴量を計算する。 ※まず、縦軸を修正する。 ⇒OFFLINE作業:gradsデータはあるはず。 _↑オンラインではonline.dataにデータを保存しているので日付データはない。 _b…