ZHEI検出プログラム確認

R□~/jobs/2019/xmp/koden/bb_csv_k/sonde_chk/zhei

課題:時間検出で1つずれる。

5           6          21

250行目に記録されている(zhei2019_47412.txt)

処理の際、12時間間隔決め打ち←データがなくてもわからない。

対処:set_search_range2に時間チェックの個所を入れる。

 

 csv

/home/koshida/jobs/2014/wyo/data/2019_47412←PRESなし

/home/koshida/jobs/2014/wyo/get_tw/data/2019_47412←PRESあり

にある。

zhei.txtを欠測で埋めて、対処。

 

debug#ゾンデ前後2時間の平均値

□~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/com_el/sonde_chk

../kita_com0914/171011/BBH/bbh_17101120s_17101122e_00_kita.csv

DEGmodeで再計算してみる。

#26

../kita_com0914/171011/BBH/bbh_17101120s_17101122e_00_kita.csv

とまる原因はわかっていない。

↑以前は、行数が足りなかった。

確認:#26 2017 4 22 8 0 2017 4 22 10 0 OK

問題は

#31 2017 10 11 20 0 2017 10 11 22 0

../kita_com0914/171011/BBH/bbh_17101120s_17101122e_00_kita.csvの読み込みエラー

対処1:元データを再計算してみる

⇒~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/com_elに移動

bbh_17101120s_17101122e_00_kita.csvは23行:再計算後も行数は同じ。

At line 1521 of file /home/koshida/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/com_el/src/main.f90

で0読み。

⇒if(k-1<1)cycle !! 9/15'1 で対処

②sonde_time.ishi.txt.tmp

#49まで読み込みOK

51 2018 3 4 8 0 2018 3 4 10 0

上記①(kita)と異なり、

../ishi_com0915/180304/BBH/mxh_18030402s_18030412e_00_ishi.csv は

120行ある。

>  5  337,  29, -99,  -9, -99,  -9,     277,*****, -99,  -9, -99,  -9,     337,18.72, -99,  -9, -99,  -9,

アスタリスク*が記録されている。

xmin_phvが問題。←-203.677が記録されている。

⇒最小値を求める際の条件に、-99を加える。(L633)

さらに、xmin_d2zdrに大きな値が記録されることがあり、99に補正(L565)

 

  5  337,  29, -99,  -9, -99,  -9,     307,-3.50, -99,  -9, -99,  -9,     337,18.72, -99,  -9, -99, 

 

⇒修正後#77まで記録。

 

関数

関数を作成する。

Binのレンジを与える⇒Zの最大、BB_TOP,BB_BOT、BB_HEIを得る。

 real function height(int)
 !convert binnumber 2 height
 integer :: int
 logical :: ldeg=.true.

  if(int>0)then
   height=int*d_range +sheight -d_range/2.
  else
   height=-99.9
  endif
 end function height

Contains内に記述しているのでMainのパラメータは引き継ぐことができる。

★BB検出の関数★

IN:探査範囲、OUT:BB_TOP、BB_BOT、BB_HEI

 

 

WBTが正の領域の描画方法

ZHEIより上で正の領域がない場合を0

ひとつある場合を1

二つある場合を2

とする。

コンター図の上に重ね描きする。

★gradsファイルはひとつ

要素をひとつ追加する。

 

 

データリスト

5/23'19DL

保存DIR

C:\home\koshida\jobs\2019\xmp\koden\data\F20190523

⇒F:\船舶レーダーデータ\data\F20190523

・時間高度2019→height3km(画像)とCSVデータ←280M

 地上観測

 F:\船舶レーダーデータ\data\F20190523\地上降水強度2015-2019\2015

 

 

##

5/31'19 append mistake