#!/bin/csh -f
########################################################################
# Output turbospectrum/moog/multi format compatible 
# A control plot (interpol_check.ps) is displayed at the end.
# Extrapolation is not advised, even if allowed by this program.
# Requires an "cubic" set of 8 MARCS binary format models,
# in other words 
# !!!!!   MODELS MUST DIFFER 2 BY 2 BY ONLY ONE PARAMETER !!!!!!  
# !!!!!!! ORDER OF THE INPUT MODELS MATTERS !!!!!!!
# here is the order of the files
# model1: Tefflow logglow zlow
# model2: Tefflow logglow zup
# model3: Tefflow loggup zlow
# model4: Tefflow loggup zup
# model5: Teffup logglow zlow
# model6: Teffup logglow zup
# model7: Teffup loggup zlow
# model8: Teffup loggup zup
########################################################################
set Tinf  = 6250 
set Tsup  = 6500
set ginf  =   +4.0
set gsup  =   +4.5
set zinf  =   -2.50
set zsup  =   -2.00

set ainf  =   +0.40
set asup  =   +0.40
set mass  =    0.0 
set vtinf =   01
set vtsup =   01

set geom  = p
set comp  = st

if ($geom == p)  set path_geom = pp
if ($geom == s)  set path_geom = sph
if ($comp == ap) set path_comp = alpha_poor
if ($comp == st) set path_comp = stand_comp

set model_path = ../marcs/
set path_Tinf  = ${path_comp}/${path_geom}/${Tinf}K
set path_Tsup  = ${path_comp}/${path_geom}/${Tsup}K

echo $path_Tinf
echo $path_Tsup

#MARCS binary format (.true.) or MARCS ASCII web format (.false.)?
set marcs_binary = '.false.'

set model1 = ${geom}${Tinf}_g${ginf}_m${mass}_t${vtinf}_${comp}_z${zinf}_a${ainf}_c+0.00_n+0.00_o${ainf}_r+0.00_s+0.00.mod
set model2 = ${geom}${Tinf}_g${ginf}_m${mass}_t${vtsup}_${comp}_z${zsup}_a${asup}_c+0.00_n+0.00_o${asup}_r+0.00_s+0.00.mod
set model3 = ${geom}${Tinf}_g${gsup}_m${mass}_t${vtinf}_${comp}_z${zinf}_a${ainf}_c+0.00_n+0.00_o${ainf}_r+0.00_s+0.00.mod
set model4 = ${geom}${Tinf}_g${gsup}_m${mass}_t${vtsup}_${comp}_z${zsup}_a${asup}_c+0.00_n+0.00_o${asup}_r+0.00_s+0.00.mod
set model5 = ${geom}${Tsup}_g${ginf}_m${mass}_t${vtinf}_${comp}_z${zinf}_a${ainf}_c+0.00_n+0.00_o${ainf}_r+0.00_s+0.00.mod
set model6 = ${geom}${Tsup}_g${ginf}_m${mass}_t${vtsup}_${comp}_z${zsup}_a${asup}_c+0.00_n+0.00_o${asup}_r+0.00_s+0.00.mod
set model7 = ${geom}${Tsup}_g${gsup}_m${mass}_t${vtinf}_${comp}_z${zinf}_a${ainf}_c+0.00_n+0.00_o${ainf}_r+0.00_s+0.00.mod
set model8 = ${geom}${Tsup}_g${gsup}_m${mass}_t${vtsup}_${comp}_z${zsup}_a${asup}_c+0.00_n+0.00_o${asup}_r+0.00_s+0.00.mod

#enter here the values requested for the interpolated model 
foreach Tref    ( 6275    )
foreach loggref (   +4.11 )
foreach zref    (   -2.20 )
foreach name    ( marcs_hd84937 )
#format output possibilities: TURBO, ATLAS, MOOG, MULTI, COMP
set format = MULTI
#set model_out = s${Tref}_g${loggref}_m${mass}_t${vtinf}_${comp}_z${zref}_a${ainf}_c+0.00_n+0.00_o${ainf}_r+0.00_s+0.00${format}
set model_out = ${name}_marcs_${Tref}_${loggref}_${zref}.${format}
#set model_out = ${format} 

#### the test option is set to .true. if you want to plot comparison model (model_test)
set test = '.false.'
set model_test = stand_comp/pp/5750K/p5750_g+4.5_m0.0_t01_st_z-0.25_a+0.10_c+0.00_n+0.00_o+0.10_r+0.00_s+0.00.mod

# interpolation program (for further details see interpol_modeles.f)
echo 'Teff = '${Tref}', log g = '${loggref}', [Fe/H] = '${zref}

./intermarcs << EOF

'${model_path}${path_Tinf}/${model1}'
'${model_path}${path_Tinf}/${model2}'
'${model_path}${path_Tinf}/${model3}'
'${model_path}${path_Tinf}/${model4}'
'${model_path}${path_Tsup}/${model5}'
'${model_path}${path_Tsup}/${model6}'
'${model_path}${path_Tsup}/${model7}'
'${model_path}${path_Tsup}/${model8}'
'${model_out}'
${Tref}
${Tref}
${loggref}
${zref}
${comp}
${test}
${marcs_binary}
'${model_path}${model_test}'
EOF