1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
#### You should not either edit or run this file ####################################
############## Below there is NOTHING to edit ##########################################
########################################################################################
case $freq in
THz) scale=33.3;;
Thz) scale=33.3;;
thz) scale=33.3;;
cm) scale=1 ;;
Cm) scale=1 ;;
CM) scale=1 ;;
meV) scale=8.0532;;
mev) scale=8.0532;;
esac
cat > Freq_plot_unit <<EOF
$scale
EOF
$PLOT_DIR/bin/k_for_bands.x < K_points
cat > matdyn.in.tmp1 <<EOF
&input
EOF
cat > matdyn.in.tmp2 <<EOF
asr='crystal',
flfrc='$FC_name.fc',
flfrq='$FC_name.freq'
&end
EOF
cat matdyn.in.tmp1 Atomic_mass matdyn.in.tmp2 ph.grid >matdyn.in
rm -f matdyn.in.tmp1 matdyn.in.tmp2 Atomic_mass
echo ' Recalculating omega(q) from C(R) ... '
$BIN_DIR/matdyn.x < matdyn.in > matdyn.out
echo ' Well done'
$PLOT_DIR/bin/bands_to_gnuplot.x <$FC_name.freq
cat > plot.GNU.begin << EOF
set term postscript enhanced color "TimesNewRoman" 18
set output '$FC_name.$freq.ps'
set nokey
set noxtics
set ylabel "Frequency, $freq"
set title "Phonon calculations for $FC_name"
set xzeroaxis lw 3
set border 15 lw 3
set encoding iso_8859_1
EOF
$PLOT_DIR/bin/E_min_max.x <$FC_name.freq
cat plot.GNU.begin plot.GNU.tmp > plot.$freq.GNU
rm -f plot.GNU.begin plot.GNU.tmp
gnuplot plot.$freq.GNU
ps2pdf $FC_name.$freq.ps
echo '###############################################################################################'
echo '###############################################################################################'
echo ' '
echo 'You have got phonon dispersion relations plotted using gnuplot (www.gnuplot.info) '
echo 'Now you can edit' plot.$freq.GNU 'in order to define E_min, E_max and Label_position '
echo 'more accurately to get publication-quality Postscript and PDF files. '
echo 'You can also redefine default parameters and add experimental data, too.'
echo ' '
echo '###############################################################################################'
echo '###############################################################################################'
|
