EL8のZDRをチェック

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回数修正後

http://hydro.iis.u-tokyo.ac.jp/~koshida/obs_mrr/sapporo/XMP/kita/rhi_sonde/20171030/PVM0/kita20171030_0855.zdr.png

修正前

http://hydro.iis.u-tokyo.ac.jp/~koshida/obs_mrr/sapporo/XMP/kita/rhi_sonde/20171030/PVM0.esc/kita20171030_0855.zdr.png

 

★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.csvbbinfo

▶(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.

⇒図は開始時。参照データは終了時。