2010-11-01から1ヶ月間の記事一覧

hikaku

hikakusuru

TODO:軸を修正する

2010/10/10のみ描画可能。 軸を変えることが必要。kouu DTCを実施した降雨で決まっている。 20101010_j.dat 20101028_j.dat 20101030_j.dat 20101101_j.datTH図作成16降雨 F:\jobs\doc\20101118\BKlst1.xls

get_mpx1の後。

おそらく読み込み成功。 TODO 地点データを書き出すので、ran,secを計算し 配列に組み込む。

kdpgkdprを作るために

DTCをRainで実施できるようにする必要がある。 計算は途中のはず。

kdpgkdpr

何をやっているか忘れないようにする。Level1データを読み込む。 rgkdpでget_mpx1を使えるようにしている。 /f/jobs/2010mpx/seido/rgkdp [地上とレーダのKdpの比較] m_mpxdataの修正。 ファイル名を作るところでSTOP.iran,isecを与えれば、データは作成でき…

平均の高度増分

あらかじめ計算しておく。

単一仰角でTH図を描く

3km〜4kmをIIS上空で観測している仰角を選択する。 新横 仰角 8.1 ° レンジビン(高度差)213.7m

仰角の選定

厳密には高度によって異なるが、平均を与える。

バイト変換(rgkdbr)

0が14になる。(どこで?) UNIX⇒PCの特徴。 2バイトの整数に注意。endianの問題。 違いを確認。

仰角分布のファイルを調査

G:にはなかった。 Z:にはなかった。 71c:mydocになし、Homeになし。 81:C:mydocになし、Homeになし。 F:doc\20101013にあり。

20度でやってみる

結構上まで見えるのではないか? 90度にこだわる必要はない。 新横でやればよい。関東でもやっておく。 ToDo ・200m間隔データの作成−100mにしておく。−距離重みで内挿。 ・11.8°、9.7°データで鉛直分布を比較する。←こっちのほうが楽。 ・Level1データ(kdp)…

データ確認−TH図

上空に弱いエコーがあるところでρhvが大きくなる(相対的な値) そのようなところでは、ρhvが下方ほど大きい。

C:\home\jobs\2010\syaka\grads\verpを参考に

gradsのスクリプトを書く

強度別にρhv評価

やってみる。 比較。 25dBZアップ。ρhvが合わない。 TH図の確認。

プログラム修正

g95ではsubrouitineで再度、integer定義してもエラーが出ない。 作りたい画像。作業。 THZ図 散布図-rr1rr2←全降雨必要。 #現在は1日のみ。リストを与える?1日95K. ⇒rr1rr2に降雨番号、時間番号を加える。降雨リスト 2010 07 29 0 5 NG 2010 08 09 0 5 201…

gswinの変換がうまくいかない。

エラーは特にないが、pngファイルができていない。

値の読み込み 精度評価

エラーが出る。 rgrr(j)でJが記入されていない。1日5分の配列が取れない。7/1を基点に半日で考える。

コンパイルはできた XMP-Kdpの比較 精度評価

参照プログラムのコンパイルはできた。実行は、エラー無し。取り出す値が大きい? ρhv・Z比較用のscatはコンパイルまで。読み込み、比較、の動作確認は未。 確認必要。 XMP,RD80の雨←依頼。ダウンロード。DTC. 対象日: 20101004 KdpR Go−4時のみ 20101010 ○…

XMP kdp計算&TODO

・DTC ・Kdp_g--ease..←考え甘い!「Floating exception」 > How to Make 5min.ave. ⇒Kdp_gはR内で計算するつもりであった。PCでもよい。 ・Kdp_r--pickup ##before the calc Kpd_r, you must pickup Kpd_r from Level 1 data, or calc from Level 0. インパ…

XMP Kdp比較

RD80からKdpが計算できるようにしておく。 main文を忘れた。。 ~/jobs/2010/mpx/calc_kdp

Zとρhvの比較

2つのレーダであわせることを考える。 まずZをあわせることが必要?11月までHDDに入れる。Mへの移動メモ ・SHINYOKO 20101021 ファイルなし? ・1031、1101PAST DL未 mkgradsの実施11/7'10->11/10'10 ・平均が必要

XMP 2レーダの比較

まず比較してみる。 合わなければ、合わせる方法を考える。 7/29はエラーでとまる?理由が不明。ファイルがなければ、欠測とするはず。rrrrの構成 CCYYMMDDFF elv1 rr1.1 rr1.2 rr1.3 rr2.1 rr2.2 rr2.3 CCYYMMDDFF elv2 rr1.1 rr1.2 rr1.3 rr2.1 rr2.2 rr2.…