set -exv

# exp will be a combination of cas, met and target

cas=MNHRAD     # tune Ts, target MNH SCM with no radiation
cas=MNHRADSURF # tune surf, target best MNHRAD

met=th_qv_nebmax
met=th_qv_nebmax_lower
met=th_qv_nebmax_lower_3h

target=MNHREF
target=MNH
target=MNHRADSCM

cat >| metrics_th_qv_nebmax_lower_3h_MNH << eod
ARMCU_MNHRAD_zav-050-200-theta_7_9,305.352022638194,0.01
ARMCU_MNHRAD_zav-050-200-qv_7_9,0.0161136349946564,1e-8
ARMCU_MNHRAD_nebmax_7_9,0.2609765,0.0002
ARMCU_MNHRAD_zav-050-200-theta_4_4,302.734175180088,0.01
ARMCU_MNHRAD_zav-050-200-qv_4_4,0.015894891734047,1e-8
ARMCU_MNHRAD_nebmax_4_4,0,0.0002
ARMCU_MNHRAD_zav-050-200-theta_11_11,306.596396593343,0.01
ARMCU_MNHRAD_zav-050-200-qv_11_11,0.0154874571279434,1e-8
ARMCU_MNHRAD_nebmax_11_11,0.1883256,0.0002
eod

cat >| metrics_th_qv_nebmax_lower_3h_MNHRADSCM.txt << eod
ARMCU_MNHRADSURF_zav-050-200-theta_7_9,305.383451642004,0.01
ARMCU_MNHRADSURF_zav-050-200-qv_7_9,0.0161136054283566,1e-8
ARMCU_MNHRADSURF_nebmax_7_9,0.2562185,0.0002
ARMCU_MNHRADSURF_zav-050-200-theta_4_4,302.596063589043,0.01
ARMCU_MNHRADSURF_zav-050-200-qv_4_4,0.0159248006066983,1e-8
ARMCU_MNHRADSURF_nebmax_4_4,0,0.0002
ARMCU_MNHRADSURF_zav-050-200-theta_11_11,306.679790494546,0.01
ARMCU_MNHRADSURF_zav-050-200-qv_11_11,0.0155678739202058,1e-8
ARMCU_MNHRADSURF_nebmax_11_11,0.1763872,0.0002
eod

case $met in 
  th_qv_nebmax_lower_3h)
    metopt="-metrics_file metrics_th_qv_nebmax_lower_3h_$target.txt"
  ;;
esac

case $cas in
  MNHRAD)
    param=param_RAD
    samples=30
    ;;
  MNHRADSURF)
    param=param_MOSAI
    samples=60
    ;;
esac

exp=${cas}_${met}_target$target

bash bench.sh -wdir $exp -param $param -opt_model "-nproc 6" -waves "1 2" -model MESONH $metopt -sample_size_next_design $samples | tee -a log/log_run_$exp.out

exit

# Est ce que je peux avoir un script qui squizz les premières étapes et pour un même jeu de simulations,
# fait une étape de history matching avec différents jeux de métriques ?
# tentative mais compliqué, vaut mieux repartir de 0...

# exp = MNHRADSURF ; ref = LES ? ou SCM-MNHRAD ?
# dans le premier, la surface (les flux) peuvent compenser des erreurs de nuages du SCM-MNH
# => sachant que j'ai aussi changé le forçage advection pour faire du rayonnement interactif...

# dans le second , revient à retrouver les flux de surface et Ts du cas forcé
# => peut être mieux comme ça on isole juste la surface ?

# est ce que je vise les nuages ou les flux de surface ?
# est ce que j'autorise à varier les paramètres de l'atmosphère ou uniquement la surface ?

# en fait je devrais d'abord caler MNHRAD sur REF : avoir les bonnes
# températures et nuages quand j'active le rayonnement, en fonction du forçage en
# XTS que j'applique


# j'ai essayé de tuner les TS RAD de MNHRAD pour avoir la bonne theta, qv et rneb_max vs REF.
# ça plante dans htune_Emulating je crois 

# BON ça me soule,
# je crois que j'ai un pb de pas de temps qui correspondent pas en tout cas
# j'arrive pas à être sure que c'est en phase LES0.nc et mes configs MESONH
# EXS_lower_3hEG