観測がある高度に対して、求めたい高度が粗度長を基準として、何倍になっているかで求める。
求める風速＝観測の風速＊（LN(求める高度/粗度長）／LN(観測のある高度/粗度長）
sample reference
Randall D, Houck MH, Wright JR. 1990. Drought
management of existing water supply system.
Journal of Water Resources Planning and
Management 116(1): 1-20. DOI: 10.1061/(ASCE)
0733-9496(1990)116:1(1).
「苗字 名前略,（コンマ）」が1セット
苗字 名前略,コンマ, （人数分続く）．（ピリオド）年号．（論文名）．（雑誌名）：Page，（コンマ）:DOI．（ピリオド）
rrgrr.rrとrgrr.rgのcombine perl技
2つのファイルを読み込んで結合する
~/jobs/2019/xmp/rgpolar/script/combine.pl
rrは別CPUで計算するのでrgと組み合わせて保存。
★手法整理
１．・uryo.csvをもとにデータ抽出する
※uryo.csvはjcod,kcod,kname,i_dfb,k_dfb,jn1,jname,refのみで構成
rain_a.csvを参照。
◆rgpolarをエクセルに取り込む１２／２２’２３
old.2016.version(~/work_kos/jobs/2022/ukanri/hout/rgpolar)
kouu,stn,ccyymmddhhff,jcod,kcod,rg,rr,rr_on,iflag1,kosu,f_calc,status1,status2,dist,shade1,shade2,bh1,bh2,jname,kname
一括処理
◆plotppi,ppiplot, plot_ppi,(plppi -secave>)
・平均条件を変える
・データの抽出
・RT変換（TMPにできる）
・描画（RHI）
を一括処理する
●mkrhi--EXWIND
1.extract target(2016/4/3 3:00~7:00)-elements :: ~/jobs/2021/xmp/all_sec/EXWIND （PV00,1要素抽出）
☒X=test.pl
[r]~/jobs/2023/xmp/sapporo/mkgrads (modified for PHM0,etc for 1elem)
☒X=exe.pl 6/1'23
[odir]./ishi/PHM0/20160403 (RADN/KIND/CCYYMMDD)
2.rtconv
[r]~/jobs/2023/xmp/sapporo/plot_ppi/rtconv1
☒X=exec.x.pl
[odir]./ishi_g/20160403/PHM0
[e.g.]g_ishi1604030415_01.CSRV.dat
3.plot_ppi(skip OK)--just check
[r]~/jobs/2023/xmp/sapporo/plot_ppi
☒X=exec.v.pl [PV00] ／exec.z.pl [PHN0]
[odir]PNG/radn/CCYYMMDD/KIND/EL##
＞＞~/jobs/2021/xmp/mkrhi
"set clevs -30 -20 -15 -10 -5 -2 2 5 10 15 20 30"
"set ccols 14 4 11 5 13 3 0 10 7 12 8 2 6"
cost 115 [sec]for 4hour.12elv.plot
★2,3SKIP可for.just.plot_RHI
★3-2mkzdr
~/jobs/2023/xmp/sapporo/plot_ppi/mkzdr[ppi向けzdr]
☒exe.x.pl
ODIR=~/jobs/2023/xmp/sapporo/plot_ppi/rtconv1/ishi_g/20160403/PZDR
★3-3mkcsr
~/jobs/2023/xmp/sapporo/plot_ppi/mkcsr
SDIR=<rtconv>
☒exe.x.pl
ODIR=~/jobs/2023/xmp/sapporo/plot_ppi/rtconv1/ishi_g/20160403/CSR[H/V]
4.make RHI data set(combine ppi data to Selected sector)
[r]~/jobs/2023/xmp/sapporo/mkrhi
☒X=exe.pl
SDIR=../mkgrads/ishi/PHM0/20160403/
ODIR=./RADN_w/exCCYYMMDD
☆skip ZDR calc
★4-2mkzdr
[r]./cgraph/mkzdr
☒X=exe.x.pl
USE: inp.RADN
SDIR=../ishi/20160403/P(H/V)(M/N)0
ODIR=../ishi/20160403/(P)ZDR(N)
==OLD==
☒X=go.v.pl (vector)
OUT=RADN/CCYYMMDD
5.plot RHI ~/jobs/2021/xmp/mkrhi/cgraph
[r]~/jobs/2023/xmp/sapporo/mkrhi
☒X=exe.pl
period_$inp="$refd/inp.$radn";
ODIR=./ishi/20160403/PHM0/sec149
==OLD==
☒X=plot_rhi.pl
