KEI先生のメモhttps://isotope.iis.u-tokyo.ac.jp/~kei/?IT%20memo/GrADS%20memo2#ea4a6ac9はうまくいかなかった。
0.MAからの切り出し．
描画高度をwgrib2で切り出す
1.指定場所の切り出し
A-A'を緯度経度で指定する
2.指定高度面についてAA'を取り出し．
ctl:
MANAL2016040218.1000.ctl
../../PICKUPh
2.mkpng
check_zy2.ge
オフライン解析装置
１．計算する
gox.pl
[65]jobs/2022/ukanri/hout_off
２．結果の確認
２．１付属ソフトで描画
[65]jobs/2022/ukanri/hout_off
mkpng.pl
２．２Gradsで描画
read_cx/simple_read(read_x)/plot_grads(rrmap)
1.extraxt data
[80]2022/tougou/rrmap/prg/read_cx [read_x, simple_read]
INP period inp.kouu >> OUT RRMAP/ccyymmdd/area_name0-ccyymmdd-hhff-G000-EL000000.dat
AREA10、ii、kkを指定。ii,kkを中央とする9ブロックを抽出。
2時間分102[sec]
[r]~/jobs/2023/xmp/sapporo/rrmap/read_x script=gox.pl
extact every 1hour @3/20'24
2.plot_grads
^^DEG.mode^^
<x>c_mkmap
INPUT: ccyymmdd,AREA10,HHFF,II,KK
^^exec.mode^^
<x>c_mkmap.x.a<y:rapper>gox.a.pl
read_cxで記載したAREA10、ii、kkをscriptに記述。
KOUU.lstはscriptに書き込み。
2時間分53 [sec]
備考：東北22で整備。
2.2. accumulate
<s>rrmap.large.sum.gs 20160301_sum test 1
[80]$home/jobs/2022/tougou/rrmap/prg/plot_grads/PNG_ONEDAY/
2.3. daily rr
[80]$home/jobs/2022/tougou/rrmap/prg/plot_grads/read_bin
精度評価
○scat_all
pdir="2022/akag"
polarとgridは入力条件が異なる。
★inpとoutを簡単にする
○[80]2022/ukanri/scat_all
地点抽出の条件設定が必要（7/25'23)
font フォントのエラー
○（おそらく）プロットルーチンによって呼び出されるフォントがインストールされていない．
○インストールが難しければ，現在インストール済みのフォントのみを用いる
◇使えるフォントと使えないフォントを整理すること．
MANAL
[one]2015\xmp\bb\MANAL\rain
file name=MANAL_2016040203.grb2
manal_CCYYMMDDHH.grb2
header 5char
spearator _
suffix .grb2
☆mesh.size from ctl.file
Y 0.0455 deg ➡ 163.6 byo ➡ 5.45 km
X 0.0633 deg ➡ 228.0 byo ➡ 5.07 km
☆tips :: get traget.time
$idx=rindex($msm,"_");
$len=length($msm);
$ccyymmddhh=substr($msm,$len-$idx-2,10);
☆plot
ame.b :2015-09-25
ランベルト図法なので以下を参照。
How to use wgrib2 command (hysk)
pdef←Grads内の定義について
GrADSのTips | 東北大学大学院 気象学・大気力学分野←ランベルト座標」を検索
pdef 721 577 lccr 30 140 490 170 60.000000 30.000000 140.000000 5000.000000 5000.000000
とする。
wgrib2:[r]~/jobs/2015/xmp/bb/wgrib2
sample script 4 wgirb2::[r]ext_ma.pl
$prg $in -order we:sn -match \":$exlev:\" -no_header -bin $out
★データは別ディレクトリへ入れる ./bin
result=./sapporo/MSM/Ztype_p
