#!/bin/bash


# ----------------
# Input Parameters
# ----------------


# [EXE] FileName Freq ExcitationType WaveformType Tdelay Twidth FinalTime FluxType PolyOrder MatrixFree ScalaStudy 

# 1. Filename: Simulation Filename
FileName="PATCH" 

# 2. Freq: Central Modulation Frequency (MHz)
Freq=2350

# 3. Planewave Waveform: (0)Time Harmonic (1)Gauss (2)Neumann (3) Modulated Gauss
waveform=1

# 4. Port Waveform (0)TH (1)Gauss (2) Neumann (3) Linear Chirp with hanning window (4) Linear Chirp 
port_waveform=4

# 5. Tdelay (sec): Delay of excitation pulse
# Add 5e-10
Tdelay=2e-10 

# 6. Tau (sec): Parameter for pulse width
Tau=7e-9

# 7. FinalTime (sec): Termination Time (Set to -1 to activate energy convergence termination)
FinalTime=1e-8

# 8. Gamma
Gamma=0.025

# 9. VTU Output Flag (0) True and (1) False
vtu=0

# 10. Flag to Write Currents / Fields on Surface (0) Off (1) On
WriteSurfBCFlag=1

# 11. Poly Order: (1) First order (2) Second order
PolyOrder=2

# 12. Free Space Comparison: (0) True and (1) False
FreeSpaceComparison=1

# 13. Sampling Rate: sampling frequency / nyquist frequency
#     Usually set to 12.5
SamplingRate=12.5

# 14. Padé Acceleration: Number of points used to calculate the Pade Aprox
# Padé Points: Number of points used to calculate the Padé convergence (starting from the top of the list)
#     Set to 0 if don’t want to use Padé Approximants
Pade=0

# 15. Time to perform Pade in ns (set to -1 if want to do convergence study)
# Padé Simulated Time: Instead of doing convergence study, once the Padé simulated time is achieved
#     Padé will be performed at the end of that Padé period. (ns)
#     Set to -1 if you want to perform a convergence study.
PadeTime=-1

# 16. Energy Convergence Factor (%): factor in % required to the energy to decay finish the simulation
# 17. Energy Convergence Factor (%): factor in % required to the energy to decay to finish the simulation (Won’t work unless FinalTime = -1)
EnergyConv=20

# 17. Energy Convergence Number of Points: number of Points used to the determine energy convergence (set to 0 for all)
#     Set to 0 to use all the probes
EnergyConvPoints=5



echo "================================================================================"
echo "EXECUTABLE"
EXE="/media/digitalstorm/edf2449b-a566-40a2-8d0f-393bbe8ea1c0/maxwell_td_hybrid_mesher/Maxwell_PML_Regular/MaxwellTD_GPU_Port/build/maxwelltd_CUDA_LTS_DOUBLEpre_FLOATpro-"
echo $EXE
echo $FileName
echo "================================================================================"


# Remove old log files
rm -f *.log

# Check FreeSpaceComparison (needs to be defined or passed in)
if [ "$FreeSpaceComparison" -neq  "1" ]; then
    mkdir -p Analytic
    cd Analytic || exit
    rm -f *
    cd ..
fi

# Create VTU output directory
mkdir -p VTU_LTS
cd VTU_LTS || exit
rm -f *
cd ..

# Create PROBES output directory
mkdir -p PROBES
cd PROBES || exit
rm -f *
cd ..

# Run DGTD simulation with tee logging
"$EXE" "$FileName" "$Freq" "$waveform" "$port_waveform" "$Tdelay" "$Tau" "$FinalTime" \
       "$Gamma" "$vtu" "$WriteSurfBCFlag" "$PolyOrder" "$FreeSpaceComparison" \
       "$SamplingRate" "$Pade" "$PadeTime" "$EnergyConv" "$EnergyConvPoints" | tee -a computeDGTD.log

# Move current files to CURRENT, and analytical outputs to Analytic
mkdir -p ./CURRENT
mv -f *.curJ ./CURRENT/
mv -f *.curM ./CURRENT/
mv -f *.csv ./Analytic/

# Post-process RCS
python PORT_POSTPROCESSING_FINAL.py