#! /bin/bash
# take UNSAT-H input file process to .bin, run UNSAT-H, 
# and call DOUT_THETA to get results in form for comparison
# to WEPS run

# USE: unsat_run  [file name preface (not including (.inp) part)]

cmdfile="unsat.cmd"

# create DIN301.EXE command input file
# input file name prefix
echo "${1}" > $cmdfile
# specify list file output (full=3)
echo "3" >> $cmdfile
# specify output to list file
echo "1" >> $cmdfile

# run DIN301.exe to create .BIN file
./din301.exe < $cmdfile > NULL

# view list file created
#less ${1}.lis

# create UNS301.EXE command input file
# input file name prefix
echo "${1}" > $cmdfile
# specify simulation progress interval (0=none; 1=yearly; 2=daily)
echo "1" >> $cmdfile

# run UNS301.exe for simulation
./uns301.exe < $cmdfile

# remove command file
rm unsat.cmd

# run DOUT_HLAYER script to create single file with THETA values
# for all years and layers
dout_hlayer ${1}