forked from wangyf/AH-RJMCMC
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathMake_all_figures
executable file
·125 lines (99 loc) · 3.31 KB
/
Make_all_figures
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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
#!/bin/bash
# Script to run all models and make all figures
if [ "$#" -ne 1 ]; then
echo "Run: Make_all_figures <Fortran/Python>"
exit 0
fi
if [ "$1" == "Fortran" ]; then
echo "Using Fortran..."
Run_command="AH"
make all
else
echo "Using Python..."
Run_command="python AH_RJMCMC.py"
fi
echo "Running RJMCMC...."
$Run_command input_Synt_Paris700 > Screen_output_Synt_Paris700
$Run_command input_Paris700 > Screen_output_Paris700
$Run_command input_Synt_Hawaii > Screen_output_Synt_Hawaii
$Run_command input_Hawaii > Screen_output_Hawaii
$Run_command input_Lubeck_Paris700 > Screen_output_Lubeck_Paris700
$Run_command input_Paris700_prior_sampling > Screen_output_Paris700_prior_sampling
$Run_command input_Hawaii_prior_sampling > Screen_output_Hawaii_prior_sampling
$Run_command input_Lubeck_Paris700_prior_sampling > Screen_output_Lubeck_Paris700_prior_sampling
$Run_command input_Synt_Hawaii_prior_sampling > Screen_output_Synt_Hawaii_prior_sampling
$Run_command input_Synt_Paris700_prior_sampling > Screen_output_Synt_Paris700_prior_sampling
$Run_command input_Paris700_no_age_errors > Screen_output_Paris700_no_age_errors
$Run_command input_Paris700_no_age_errors_twice_F_error > Screen_output_Paris700_no_age_errors_twice_F_error
$Run_command input_Paris700_no_age_errors_min_5muT > Screen_output_Paris700_no_age_errors_min_5muT
$Run_command input_Tema_etal_2017 > Screen_output_Tema_etal_2017
# Remove outlying data from Paris700 and rerun the model..
cd Outputs_Paris700
python ../remove_outlying_data.py
cp data_within_credible_bounds.dat ../archeomagnetic_data/Paris700_within_credible_interval.txt
cd ..
$Run_command input_Paris700_pass_2 > Screen_output_Paris700_pass_2
# only Fortran version..
if [ "$1" == "Fortran" ]; then
IL input_Synt_Paris700_IL > Screen_output_Synt_Paris_IL
cd Outputs_Synt_Paris700_IL
python ../make_plots.py
cd ..
fi
#
echo "Making all the figures..."
# Now make all the plots...
# .......
cd Outputs_Synt_Paris700
python ../make_plots.py joint_highlight_points 119 154
python ../make_joint_distribution_plots.py 0154
python ../make_joint_distribution_plots.py 0119
# only with Fortran
if [ "$1" == "Fortran" ]; then
python ../make_AH_IL_comparison.py ../Outputs_Synt_Paris700_IL
fi
#
cd ..
cd Outputs_Paris700
python ../make_plots.py
python ../period_analysis_Paris700.py
cd ..
cd Outputs_Hawaii
python ../make_plots.py
python ../period_analysis_Hawaii.py
cd ..
cd Outputs_Synt_Hawaii
python ../make_plots.py
cd ..
cd Outputs_Lubeck_Paris700
python ../make_plots.py
python ../make_Stratified_marginal_age_plot.py
python ../make_joint_distribution_plots.py 13
python ../make_joint_distribution_plots.py 22
cd ..
cd Outputs_Paris700_pass_2
python ../make_plots.py
cd ..
cd Outputs_Paris700_no_age_errors
python ../make_plots.py
cd ..
cd Outputs_Paris700_no_age_errors_min_5muT
python ../make_plots.py
cd ..
cd Outputs_Paris700_no_age_errors_twice_F_error
python ../make_plots.py
cd ..
python make_age_sampling_comparison.py
cd Outputs_Tema_etal_2017_dataset
python ../make_Posterior_with_shifted_ages.py
cd ..
$Run_command input_Herve_etal_2017 > Screen_output_Herve_etal_2017
cd Comparisons/Herve_2017_comparison
python make_Herve_comparison.py
cd ..
cd Comparisons/LeGoff_Paris700_comparison
python make_Paris700_comparison_plot.py
cd ..
cd Comparisons/Tema_etal_2017_comparison/
python make_Tema_comparison.py
cd ..