2019-01-01から1年間の記事一覧

BB高度

BB高度を探査範囲の、上下端に一致しないように設定。 ※下の画像は実体ではなくてリンクと思われる。

Sonde時間30分ー平均の結果を参照している

平均時間を参照しているのであれば、個別の時間を評価する。

bbh mxhの抽出

最初にエクセルベースで、必要データを取り出す。 ⇒ピーク高度は整理できる。上端下端はそろっていないと利用しない。

vprofile作成手順

1.データDL 2.合成マップ(仰角分布) 3.地点ごとのプロファイル作成

v_profile check -- one by one

30分平均 □~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/com_el/sonde_chk/mkgrads ⇒5分ごと、シンプルに X:$prg $its_ite_type $radn $mtype $zhei $bbtop $bbpek $bbbot $wbt $sonde_time $nko $son_cou $DEG ↑exec.x.plで実行。 1.設計 ctlファイルは、…

cron list check

yosoku_c /export/raid2/koshida/MRR/data/reame/TENKIZU/c_cp TENKIZU :: cron DL 5/22'22 re.start save about 5month--- script error saved month 12/13'23

ZHEI検出プログラム確認(point)/aero

R□~/jobs/2019/xmp/koden/bb_csv_k/sonde_chk/zhei X=main.exe S=exec.py @1/16'23➡~/jobs/2022/akagi/zhei 「radn指定が可能なようにする」 課題:時間検出で1つずれる。 5 6 21 250行目に記録されている(zhei2019_47412.txt) ↑ 処理の際、12時間間隔決め…

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 とまる原因はわかってい…

bbh計算(XMP)

①合成 注:合成の前にelevation data convertが必要 □~/jobs/2019/xmp/mkgrads2010/all_sec X:c_go.a INP.KOUU □~/jobs/2019/xmp/mkgrads2010/all_sec/comishi ②BBH計算(毎5分) □~/jobs/2019/xmp/mkgrads2010/bb_csv/prg3/com_el @10/26'19 ZDR condtion…

関数

関数を作成する。 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…

WBTが正の領域の描画方法

ZHEIより上で正の領域がない場合を0 ひとつある場合を1 二つある場合を2 とする。 コンター図の上に重ね描きする。 ★gradsファイルはひとつ 要素をひとつ追加する。

0℃高度再挑戦

3時間毎のMSM初期値を利用。 AERO版は2015年2月2日を参照。

データリスト

5/23'19DL 保存DIR C:\home\koshida\jobs\2019\xmp\koden\data\F20190523 ⇒F:\船舶レーダーデータ\data\F20190523 ・時間高度2019→height3km(画像)とCSVデータ←280M ・ time.info added H:\koden\New\20160403\03 @5/24'21 地上観測 F:\船舶レーダーデ…

TH図の作成(キタ)

▽:計算結果.or.Result ーーーーー 0.データの取得 DL:~/jobs/2018/python/dias DATA:/export/raid2/koshida/DIAS 1.mkgrads (PPI画像データ) □~/jobs/2019/xmp/mkgrads2010/all_sec ▽./RADN/CCYYMMDD Y:c_go.a 注:対象レーダはgo_a.x内で選択 2.セクタ…

perl技 9時と21時を取り出す

時間でループしてit(4)が9時、21時であれば,printする @11/21’22 JSTをUTCに変換して,1,0で求める. ★前後のSonde時刻を求めるとき itが12Zより前(it[3]/12<1⇒int(it[3]/12)=0)であれば,0と1が前後のSonde時刻 itが12Zより後(it[3]/12>=1⇒in…

DIASファイルの展開

/export/raid2/koshida/DIAS/ex.pl DLしたデータの展開。 ひまわり DIASシステムへのログイン から申請が必要。 (ウ)データセットの利用を予定する期間 いつまで利用するかを申請すること。

作業イメージ

gradsdatを作成する 2次差分を計算する BB高度の抽出

mkgrads作業

1時間ごとに、PMH0を取り出す。 □~/jobs/2018/python/dias/ex_test tarの展開テスト ★過去の作業記録 A~/jobs/2015/xmp/mkgrads/multi_elv Rへコピー 4/1'19 3からコピー R~/jobs/2019/xmp/mkgrads2010/all_sec makeはライブラリのディレクトリを変更して実…

単体レーダの累積雨量

一冬-雨量取り出し⇒目標1日(設定;3時間) 累加⇒1日(設定;2時間) ⇒累加の際に、ZRとKDPを分ける↑取り出しの際にQFも記録する。 作図⇒1日(設定;2時間)

記者発表資料【国土交通省】

関東地整 「XRAIN(エックスレイン)」の配信エリアを拡大します! http://www.ktr.mlit.go.jp/kisha/river_00000346.html 氏家、八斗、 http://www.ktr.mlit.go.jp/ktr_content/content/000083183.pdf 東北地整 平成29年 8月10日 XバンドMPレーダ盛岡局…

0℃高度の時間変化

館野の高層観測による 0℃高度時間変化(館野)

ZDRによるBB判定

下端:融け切ったところと考え、ZDRが極大から一定値へ移行した高度の一定値側⇒粒子が球形になっている上端高度。 上端:0より大となった高度。 ↑ まとめるとZDR>=0の領域をBBと判定する。 ZDRの鉛直分布

アメダスデータのwget

マクロでメッセージを表示させない

BB判定方法

ビーム中心高度とBBTOP、BBBOTの差で評価 code judge bb_top-cnter bb_bot-center 1 bb_in 正 負 0 snow 負 負 2 rain 正 正 -1 NoBB N.A. N.A. 現時点では、TOP、BOTを目視で判断

pathの設定 IMProtoo

http://www.bnote.net/linux/path.shtmletc/cshrc ←リンク切れ に記述。 参考: http://www.bnote.net/linux/path.shtml 管理者権限のない中での http://keisanbutsuriya.hateblo.jp/entry/2015/01/26/224127 管理者

repair

201511は実施済み 201522を修正

0℃高度

取り出し準備。 データはとったら削除。

0℃高度調査(Grid)

F90プログラムはある。[MSM pall.version src.check@1/16'23](A/c/home/koshida/jobs/2018/xmp/rgpolar/search0/get0kod/c_go) (仕様)ido,kdo,cccyymmddhhffを与えることでZHEIを得る。 perlでどのように使うか、を考え、関数likeにした。 ↑時間がかかり…

作業整理

・極座標の表示 ・地形とのマッチング(ただの趣味) ・昨日(1/31'19)の状況について ・定数変更後の雨量換算 ・GZ展開方法 ・csv作成 ・0度高度抽出 ・リンクについて調査 st 10:56--意外とリストアップできない

grads 複数データの描画

期間をかえるにあたって、読み込みファイルを変える ☆ファイル名は同じでよい ↑*.ctlは変える必要がない。 描画して、PNGに保存する際に、出力ファイル名を指定する。 条件:計算期間を与える(its,ite,または、ccyymmdd_ccyymmdd) ⇒its,iteをもらうとccyymm…