7. Menggabungkan ke NetCDF + plot gaya basin-atlas (vflood-netcdf)

Untuk output “satu file untuk seluruh deret waktu” bukan ratusan .out/.tif — dan untuk plotting dengan cara yang sama seperti basin-atlas memplot NetCDF-nya:

vflood-rri RRI_Input.txt                              # jalankan simulasi
vflood-netcdf RRI_Input.txt --variables hs,hr,qr --plot-last --plot-max --geotiff-max   # gabung + preview

Menulis <out folder>/simulation.nc — dimensi (time, lat, lon) dengan satu variabel per output yang diminta, metadata CF-ish (units, standard_name bila berlaku).

Baca lewat netCDF4.Dataset (pilihan library sama dengan writer vflood-netcdf; bukan xarray). Koordinat lat/lon adalah cell-center dengan lat menurun (row 0 adalah baris paling utara). Ini tidak dibalik paksa menjadi lat naik — cek urutan lat[:].

Jika --plot-last ditambahkan, akan dibuat satu PNG per variabel (timestep terakhir) menggunakan atlas_plot.py.

--plot-max merender frame puncak, bukan timestep terakhir — timestep terakhir belum tentu yang paling menarik (ekor resesi run bisa jauh melampaui puncak sebenarnya). Untuk tiap variabel, cari timestep di mana nilai maksimum domain-wide-nya terjadi (postprocess.domain_max_series + hydrograph_stats), lalu cetak step dan tanggal aslinya:

vflood-netcdf: wrote simulation_qr_max.png (qr_max=405.4195 at step 120, t=432000s, date=2023-09-22 00:00)

Butuh start_date.txt untuk tanggalnya — fallback ke step/detik saja kalau tidak ada.

--geotiff-max menulis frame puncak yang sama sebagai GeoTIFF georeferensi asli — simulation_<var>_max.tif, langsung bisa dibuka di QGIS/ArcGIS. Independen dari --plot-max — bisa pakai salah satu atau keduanya. Tombol "Bundle NetCDF" di VFloodControl selalu memakai ketiganya (--plot-last --plot-max --geotiff-max).

Opsinya:

--variables hs,hr,hg,qr,gampt_ff   variabel yang digabung (default: hs,hr,qr)
--out <path>                       default: <out folder>/simulation.nc
--plot-last                        juga render preview PNG per variabel (timestep terakhir)
--plot-max                         juga render + laporkan frame puncak per variabel
--geotiff-max                      juga tulis frame puncak sebagai GeoTIFF georeferensi asli