AMATERASS日射量プロダクトを用いたデータ解析実習

mkdir -p ~/amaterass_jica
cd ~/amaterass_jica

java -version
javac -version
bzip2 --version
date --version
gnuplot --version

wget1 --version
wget --version

cd ~/amaterass_jica
wget -c https://downloads.amaterass.science/tmap_v3.0cc_20260815d.tar.gz
tar xzf tmap_v3.0cc_20260815d.tar.gz
cd tmap
javac *.java

wget -O download_day1.sh https://amaterass.science/ja/jica-training/download_day1.sh
chmod +x download_day1.sh
./download_day1.sh

~/amaterass_jica/data/
├── reference/
│   ├── standard_4000x4000.iodc.42.5.lat.fld.3km.bin
│   └── standard_4000x4000.iodc.42.5.lng.fld.3km.bin
└── 20260501/
    └── 202605010912.dwn.sw.flx.sfc.fld.3km.bin

wget -O analysis_iodc_day1.tar.gz https://amaterass.science/ja/jica-training/analysis_iodc_day1.tar.gz
tar xzf analysis_iodc_day1.tar.gz
javac analysis.iodc.day1/*.java

java -cp analysis.iodc.day1 IODCCrop \
  data/20260501/202605010912.dwn.sw.flx.sfc.fld.3km.bin \
  data/reference/standard_4000x4000.iodc.42.5.lat.fld.3km.bin \
  data/reference/standard_4000x4000.iodc.42.5.lng.fld.3km.bin \
  51.0 73.0 34.02 47.02 \
  202605010912.central_asia.bin

202605010912.central_asia.bin と、Tmap用の格子情報 202605010912.central_asia.bin.grid.txt が生成されます。

切り出し処理で作成した 202605010912.central_asia.bin.grid.txt には、Tmapで必要な格子サイズと緯度・経度範囲が保存されています。まずその内容を確認します。

cat 202605010912.central_asia.bin.grid.txt

GRID=202605010912.central_asia.bin.grid.txt

WIDTH=$(sed -n 's/^width=//p' "$GRID")
HEIGHT=$(sed -n 's/^height=//p' "$GRID")
NORTH=$(sed -n 's/^north=//p' "$GRID")
LAT_SPAN=$(sed -n 's/^lat_span=//p' "$GRID")
LAT_DIV=$(sed -n 's/^lat_div=//p' "$GRID")
WEST=$(sed -n 's/^west=//p' "$GRID")
LON_SPAN=$(sed -n 's/^lon_span=//p' "$GRID")
LON_DIV=$(sed -n 's/^lon_div=//p' "$GRID")

cd tmap
java tmap ../202605010912.central_asia.bin 1400 0 8 "W/m2" auto png   "$NORTH" "$LAT_SPAN" "$LAT_DIV"   "$WEST" "$LON_SPAN" "$LON_DIV"   "$WIDTH" "$HEIGHT"
cd ..

wget -O download_day2.sh https://amaterass.science/ja/jica-training/download_day2.sh
chmod +x download_day2.sh
./download_day2.sh

wget -O analysis_iodc_day2.tar.gz https://amaterass.science/ja/jica-training/analysis_iodc_day2.tar.gz
tar xzf analysis_iodc_day2.tar.gz
javac analysis.iodc.day2/*.java

java -cp analysis.iodc.day2 IODCAverage hourly data/20260501 2026050109 \
  2026050109.hourlymean.dwn.sw.flx.sfc.fld.3km.bin

java -cp analysis.iodc.day2 IODCAverage daily data/20260501 20260501 \
  20260501.dailymean.dwn.sw.flx.sfc.fld.3km.bin

java -cp analysis.iodc.day2 IODCAverage monthly data/daily 202605 \
  202605.monthlymean.dwn.sw.flx.sfc.fld.3km.bin

wget -O download_day3.sh https://amaterass.science/ja/jica-training/download_day3.sh
chmod +x download_day3.sh
./download_day3.sh

wget -O analysis_iodc_day3.tar.gz https://amaterass.science/ja/jica-training/analysis_iodc_day3.tar.gz
tar xzf analysis_iodc_day3.tar.gz
javac analysis.iodc.day3/*.java

java -cp analysis.iodc.day3 IODCPointSeries data/20260501 \
  data/reference/standard_4000x4000.iodc.42.5.lat.fld.3km.bin \
  data/reference/standard_4000x4000.iodc.42.5.lng.fld.3km.bin \
  20260501 45.02 56.48 point_A_20260501.dat

java -cp analysis.iodc.day3 IODCPointSeries data/20260501 \
  data/reference/standard_4000x4000.iodc.42.5.lat.fld.3km.bin \
  data/reference/standard_4000x4000.iodc.42.5.lng.fld.3km.bin \
  20260501 41.02 58.48 point_B_20260501.dat

出力
├── point_A_20260501.dat
└── point_B_20260501.dat

gnuplotでは、出力形式、軸、時刻書式、入力ファイルを順に指定します。ここではコマンドをそのまま入力して2地点を同じグラフに描きます。

gnuplot

set terminal pngcairo size 1200,700 enhanced font ",14"
set output "point_timeseries_20260501.png"
set title "Surface downward shortwave flux: 1 May 2026"
set xlabel "UTC"
set ylabel "W/m^2"
set xdata time
set timefmt "%Y-%m-%dT%H:%M:%S"
set format x "%H:%M"
set grid
set key top left

plot "point_A_20260501.dat" using 1:2 with linespoints title "45.02N, 56.48E",      "point_B_20260501.dat" using 1:2 with linespoints title "41.02N, 58.48E"

exit

point_timeseries_20260501.png が生成されます。