#!/bin/sh
for i in 400
do
    bcl=$(printf "%04d" $i)
   
    # eg 1 high spont and wavelets, low wave
    spont_rate=$(printf "%.2f" 100.0)  # spontaneous sparks 
    wavelet_spont_rate=$(printf "%.2f" 5.0) # spont->wavelet
    wave_rate=$(printf "%.2f" 1.0) # wavelet -> wave
    CaSR50=$(printf "%.5f" 0.0) # shift of CaSR50 for waves and spont etc
    kCaSR_spark=$(printf "%.5f" 1.00) # steepness of sparks

    # also increase Jup, ICaL and IK here as an ISO like model to promote SR loading
    for model in hAM_CRN hAM_GB hVM_GB
    do
        ./model_single_cell_spark Model ${model} BCL ${bcl} Reference SCRE_${model} spark_alpha_spont_scale ${spont_rate} spark_alpha_wavelet_spont_scale ${wavelet_spont_rate} spark_spont_k_CaSR_scale ${kCaSR_spark} spark_wave_CaSR50_shift ${CaSR50} spark_spont_CaSR50_shift ${CaSR50} spark_pwave_scale ${wave_rate} Results_Reference ${bcl}_spnt_${spont_rate}_wvlt_${wavelet_spont_rate}_wave_${wave_rate}_k_${kCaSR_spark}_CaSR50_${CaSR50} Beats 200 Total_time 85000 Jup_scale 1.75 ICaL_scale 1.25 IKr_scale 1.25 IKs_scale 1.25 IKur_scale 1.5
    done
    
    # eg 2 high waves
    spont_rate=$(printf "%.2f" 10.0)  # spontaneous sparks 
    wavelet_spont_rate=$(printf "%.2f" 1.0) # spont->wavelet
    wave_rate=$(printf "%.2f" 50.0) # wavelet -> wave
    CaSR50=$(printf "%.5f" 0.00) # shift of CaSR50 for wave
    kCaSR_spark=$(printf "%.5f" 1.00) # steepness of sparks
    # can also explore thresholds (set higher/lower to drastically change spont)

    for model in hAM_CRN hAM_GB hVM_GB
    do
        ./model_single_cell_spark Model ${model} BCL ${bcl} Reference SCRE_${model} spark_alpha_spont_scale ${spont_rate} spark_alpha_wavelet_spont_scale ${wavelet_spont_rate} spark_spont_k_CaSR_scale ${kCaSR_spark} spark_wave_CaSR50_shift ${CaSR50} spark_spont_CaSR50_shift ${CaSR50} spark_pwave_scale ${wave_rate} Results_Reference ${bcl}_spnt_${spont_rate}_wvlt_${wavelet_spont_rate}_wave_${wave_rate}_k_${kCaSR_spark}_CaSR50_${CaSR50} Beats 200 Total_time 85000
    done
done