Drought model workflows

[rburghol]

• MLLR:
  • (NEW) Meta-model workflow:
    • Run all sbatch run_model drought 2026 all auto mllr analyze
    • Run one sbatch run_model drought 2026 02035000 auto mllr analyze
  • Run stat calculation for eventual map display
    • Script: https://github.com/HARPgroup/vahydro/blob/master/drupal/dh_drought/src/r/call_all_gages_mllr_v03.R
    • Ex (all gages): sudo -u www-data Rscript /opt/model/vahydro/drupal/dh_drought/src/r/call_all_gages_mllr_v03.R
  • Run chart creation (MLLR and context) for a single gage
    • Script: analyze_winter_one.R
    • Ex: https://deq1.bse.vt.edu:81/usgs/mllr/usgs_02035000_mllr_bar_winterflows_2025_.png
  • TODO:
    • Add context charts to show distribution of:
      • flows during the recharge period
      • Precip during the recharge period
• Baseflow Analysis
  • Extract baseflow events from time series
    • Flag baseflow
  • Project baseflow into future given zero (or minimal) recharge.

MLLR Risk map:

Sample gage: 02035000

New gage winter flow distribution

Create the Meta Model

META_MODEL_ROOT=/opt/model/meta_model
export META_MODEL_ROOT
/opt/model/meta_model/scripts/make_model_shell drought "mllr" "download analyze"

outflow

Run for all gages

filename="/var/www/html/tmp/mllr_drought_probs_2025.tsv"
while IFS= read -r line; do
    #echo "Text read from file: $line"
    IFS=$'\t' read hydrocode tsvalue varkey tstime hcode <<< "$line"
    gageno=${hydrocode//"usgs_"/}
    echo "Handling $gageno"
    Rscript /opt/model/meta_model/scripts/usgs/analyze_winter_one.R $gageno
done < $filename