littlepadawan
commented
1 month ago GitHub does not support uploading the script files, but here is the content of the interpolation script for p5865_g+4.0_z-0.26_a+0.10_mg+0.51_ca-0.12:
#!/bin/csh -f
set model_path = /home/ida/Documents/thesis/model_atmospheres/1D
set marcs_binary = '.false.'
#enter here the values requested for the interpolated model
foreach Tref ( 5865 )
foreach loggref ( 4.01 )
foreach zref ( -0.257 )
set modele_out = /home/ida/Documents/thesis/wrapper/TASS/output/2024-06-04_0839/temp/p5865_g+4.0_z-0.26_a+0.10_mg+0.51_ca-0.12.interpol
set modele_out2 = /home/ida/Documents/thesis/wrapper/TASS/output/2024-06-04_0839/temp/p5865_g+4.0_z-0.26_a+0.10_mg+0.51_ca-0.12.alt
#plane-parallel models
set model1 = p5750_g+4.0_m0.0_t01_st_z-0.50_a+0.20_c+0.00_n+0.00_o+0.20_r+0.00_s+0.00.mod
set model2 = p5750_g+4.0_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
set model3 = p5750_g+4.5_m0.0_t01_st_z-0.50_a+0.20_c+0.00_n+0.00_o+0.20_r+0.00_s+0.00.mod
set model4 = 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
set model5 = p6000_g+4.0_m0.0_t01_st_z-0.50_a+0.20_c+0.00_n+0.00_o+0.20_r+0.00_s+0.00.mod
set model6 = p6000_g+4.0_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
set model7 = p6000_g+4.5_m0.0_t01_st_z-0.50_a+0.20_c+0.00_n+0.00_o+0.20_r+0.00_s+0.00.mod
set model8 = p6000_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
#### the test option is set to .true. if you want to plot comparison model (model_test)
set test = '.false.'
set model_test = 'Testwebformat/p5750_g+4.5_m0.0_t01_ap_z-0.25_a+0.00_c+0.00_n+0.00_o+0.00_r+0.00_s+0.00.mod'
# interpolation program (for further details see interpol_modeles.f)
./interpol_modeles <<EOF
'${model_path}/${model1}'
'${model_path}/${model2}'
'${model_path}/${model3}'
'${model_path}/${model4}'
'${model_path}/${model5}'
'${model_path}/${model6}'
'${model_path}/${model7}'
'${model_path}/${model8}'
'${modele_out}'
'${modele_out2}'
${Tref}
${loggref}
${zref}
${test}
${marcs_binary}
'${model_test}'
EOF
end
end
end
There is a bug resulting in spectra sometimes not being generated. At the moment, it is unclear why.
Logs from running the interpolator and Turbospectrum are saved in the
tempdir and reveals two different problems../interpol_modeles: Command not found.. This means that there is no model atmosphere forbabsmaandbsynto read, and no spectrum is generated.babsmascript fails with the error message . This means that no opacity file is generated, and thebsynscript lacks necessary input.Attachements:
Spectra generated: p5011_g+4.5_z+0.46_a+0.00_mg+0.55_ca-0.07 interpolate_p5011_g+4.5_z+0.46_a+0.00_mg+0.55_ca-0.07.log p5011_g+4.5_z+0.46_a+0.00_mg+0.55_ca-0.07_babsma.log p5011_g+4.5_z+0.46_a+0.00_mg+0.55_ca-0.07_bsyn.log
Babsma error: interpolate_p5232_g+4.8_z-0.42_a+0.17_mg+0.71_ca+0.25.log p5232_g+4.8_z-0.42_a+0.17_mg+0.71_ca+0.25_babsma.log p5232_g+4.8_z-0.42_a+0.17_mg+0.71_ca+0.25_bsyn.log
Interpolation error: interpolate_p5865_g+4.0_z-0.26_a+0.10_mg+0.51_ca-0.12.log