BBH/で
はTOP&BOT
mxh_1710300855_12.kita.csv
はMAX(MIN)とVAL
昨日、発見したように、移動平均をかけすぎると、ピーク位置がずれる*。
一方で、移動平均しないとをBB位置が不定になる場合がある**。
*EL8:移動平均後で抽出していたBB位置と観測結果を動径方向のみ移動平均した生値に近いGrads分布図で不一致が見られた。
**EL6:移動平均時には検出できていたBBが未検出となった。
⇒mxh_201710300900_kita.1009
3386 1
Heightが高い。ZDRが定義されていない。
RHIではEL6でも定義できてほしい。
↑
求めることができない理由を確認する。
☆RMLはEL1~EL11まで検出できた。
★EL6を検証。
kmarkc(3)= 195←3387m、ZDR=1dB
↑
検出範囲を0~3000mとする。下限は8BINとしているので、上限を3,000m/dhで求める。
673.mへ改善。
MXHとBBHをMXH12に保存する。
⇒PPTに張り付けた表の根拠を整理する。
kmarkc(3)の検出を確認する。
system "./main.run 17103008 58 15 08 kita";では37
C:\home\koshida\jobs\2019\xmp\mkgrads2010\bb_csv\prg3\all_sec\ELV12\fort.505.xlsx
peak val 0.450 thresh: 0.040 kmarkc(3)= 29 hei= 768. L1252
移動平均の回数を抑えることで、見た目と検出BB高度が一致する。
RunMean回数修正後
修正前
★BB計算手順
1.Grads.dat作成
□vprofile確認:~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/all_sec/ELV12/vprof
-rw-rw-rw- 1 koshida hydrolab 2563200 4月 11 2019 ../../../../all_sec/kita/20171030/kita1710300855_01.dat
↑
Grads.datの変更はしていない
⇒grads make soruce :: modules/m_elvs.f90
2.EL12それぞれのBB
~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/all_sec/ELV12
※el_check.pl:単独1仰角実施- main.run
▶(result):mlh171030_0832_02.kita.csv:bbinfo
▶(result):bbh171030_0832_02.kita.csv :top~bot
▶(result):bin171030_0832_02.kita.csv :top,bot,pek bin
▶(result):mxh171030_0832_02.kita.csv:pek bin & value
※exec.pl:12仰角実施:EX =go.a.pl
⇒bin_1710300855_12.kita.csv 等
▶mlh*_EL.RADN.csv in kita1025/171030
▶mxh*_EL.RADN.csv & bbh*_EL.RADN.csv in kita1025/171030/BBH
▶bin*_EL.RADN.csv in kita1025/171030/BIN
▶mxh_*_RADN.csv & bbh_*_RADN.csv in current
3.仰角ごとデータを5分ごとにまとめる。
◇cat_height.a.pl exec in test.pl --- make subdir ABN/
※テスト版:test.pl
※実行版:merge12el.x.pl INPKOUU
▶(result)bbh_201710300835_kita.csv -top~bot height
▶(result)mxh_201710300835_kita.csv -max height & value
▶(result)ABN/bin_201710300835_kita.csv -top&bot max bin used for Grads
4.calc & compile 12elevation BBH info in sonde time
※ext_bb/go.x.pl
▶(result)ABN/ahe_STARTTIME000_FINISHTIME00_kita.csv
5.RHI図化
□~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/all_sec/ELV12/mkrhi/cgraph
※plot_tes.pl
★sub_set_bbinfo ABNデータを読んでBB高度(bbinfoh)を返す。
◀inp.kita.1 :: period - written in script.
◀~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/all_sec/ELV12/kitaMMDD :: bbinfo
◀gdat:~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/all_sec/ELV12/mkrhi/kita/20171030
▶modules/m_elvs.f90 出力Dirを設定
★8:30~8:34のRHIは bin_201710300835_kita.csvを参照:: Time proceed 5min.
⇒図は開始時。参照データは終了時。
・bbh_201710300855_kita.csv 3タイプのBB。3つまで。12仰角。
1.TOP&BOT
./171030/ABN/bin_201710300900_kita.csv
2.BBinfo値
./mxh_201710300900_kita.csv
3.BIN番号
●MXHのZML-EL1⇒-99,1.43 、EL2⇒551,1.17
from
kita0928/171030/BBH/mxh_1710300858_02.kita.csv
●BINのZML-EL1 103、EL2⇒74
from
kita0928/171030/BIN/bin_1710300858_02.kita.csv
☆元の値はそのまま引用されている。
◆ src
MXH
write(25,'(i4,",",f4.1,",")',advance='no') int( bb_hei_rml), xmax_dbz
write(25,'(4(i4 ,","))',advance='no')int(bb_hei_rml2),ibb_val2,int(bb_hei_rml3),ibb_val3
!pml
write(25,'(i8,",",f5.2,",")',advance='no') int(bb_hei_pml),xmin_phv !***9/15'19***
write(25,'(2(i4,",",f5.2,","))',advance='no')int(bb_hei_pml2),bb_val_pml2,int(bb_hei_pml3),bb_val_pml3
!zml
write(25,'(i8,",",f5.2,",")',advance='no') int(bb_hei_zml),xmax_zdr
write(25,'(4(i4,","))',advance='yes')int(bb_hei_zml2), ibb_val2,int(bb_hei_zml3),ibb_val3
BIN
write(26,'(3(i4,","))',advance='no') ibb_top_rml,ibb_bot_rml,kmarkc(1)
write(26,'(3(i4,","))',advance='no') ibb_top_rml2,ibb_bot_rml2,kmarkc(7)
write(26,'(3(i4,","))',advance='no') ibb_top_rml3,ibb_bot_rml3,kmarkc(10)
!pml
write(26,'(3(i4,","))',advance='no') ibb_top_pml,ibb_bot_pml,kmarkc(2)
write(26,'(3(i4,","))',advance='no') ibb_top_pml2,ibb_bot_pml2,kmarkc(8)
write(26,'(3(i4,","))',advance='no') ibb_top_pml3,ibb_bot_pml3,kmarkc(11)
!zml
write(26,'(3(i4,","))',advance='no') ibb_top_zml,ibb_bot_zml,kmarkc(3)
write(26,'(3(i4,","))',advance='no') ibb_top_zml2,ibb_bot_zml2,kmarkc(9)
write(26,'(3(i4,","))',advance='no') ibb_top_zml3,ibb_bot_zml3,kmarkc(12)
☆BINとMXHの不一致が問題なのは、評価が不定となるから
⇒bb_hei_zmlとkmarkc(3)の関係を調べる
L259:bb_hei_rml=height(kmarkc(1))←RML(参考)
L1530:bb_hei_zml=height(kmark(3))←kmarkc(3)の誤植!?
518 kmark(3)=k !! insert Phv
prg修正後の結果:
rainbow% diff mxh_201710300900_kita.csv mxh_201710300900_kita.csv.3
で、mxh_201710300900_kita.csvにはEL11,EL12を除く仰角に、値が記録された。ZMLのRHI図と整合!
3要素を横並びにする(RHIをつなげる)
評価(BBの○)のある版、ない版を整理。
RHI版で走る。