現状のアルゴリズムでまとめる

昨日、発見したように、移動平均をかけすぎると、ピーク位置がずれる*。

一方で、移動平均しないとを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に張り付けた表の根拠を整理する。

 

 

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.

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

BBinfo元データ#RHI

●GradsDat

../../../../all_sec/kita/20171030/kita1710300858_08.dat

●binfile_o

bin_YYMMDDHHFF_EL.RADN.csv

UNIT=26

●main.run

 

For EL8  1  enum:           3 kmarkc          37   975.989929  in fort.803

For EL9  1  enum:           3 kmarkc          31   975.820435   in fort.803

 

 

BBinfoファイルの確認

bbh_201710300855_kita.csv 3タイプのBB。3つまで。12仰角。

1.TOP&BOT

bbh_201710300900_kita.csv

./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_zmlkmarkc(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図と整合!

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