-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrun_tests.m
101 lines (77 loc) · 2.44 KB
/
run_tests.m
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
% Newton Nguyen
% 05/14/2018
% A script to run tests of the inversions located in DriverScript
% Exporting results of test to test_results.mat for python plotting because pythong plots are being used in the paper
% Reference Case from AJT 2017 paper
DriverScript
case0.concentrations = out_anal;
case0.ems = anal_soln;
case0.nh_ch4_ems = anal_soln{1}(:,1);
case0.sh_ch4_ems = anal_soln{1}(:,2);
case0.nh_co_ems = anal_soln{1}(:,13);
case0.sh_co_ems = anal_soln{1}(:,14);
case0.nh_oh = anal_soln{1}(:,11);
case0.sh_oh = anal_soln{1}(:,12);
save('case0_test_tau0.mat')
%%%1. Run the basecase
%%% Alex's PNAS Paper
% ignoreCO = true
% interactiveOH = false
% fixedOH = true
DriverScript_case1
% Give variables names
case1.concentrations = out_anal;
case1.ems = anal_soln;
case1.nh_ch4_ems = anal_soln{1}(:,1);
case1.sh_ch4_ems = anal_soln{1}(:,2);
case1.nh_co_ems = anal_soln{1}(:,13);
case1.sh_co_ems = anal_soln{1}(:,14);
case1.nh_oh_ems = anal_soln{1}(:,11);
case1.sh_oh_ems = anal_soln{1}(:,12);
save('case1_test_tau0.mat')
DriverScript_case2
% Give variables names
case2.concentrations = out_anal;
case2.ems = anal_soln;
case2.nh_ch4_ems = anal_soln{1}(:,1);
case2.sh_ch4_ems = anal_soln{1}(:,2);
case2.nh_co_ems = anal_soln{1}(:,13);
case2.sh_co_ems = anal_soln{1}(:,14);
case2.nh_oh_ems = anal_soln{1}(:,11);
case2.sh_oh_ems = anal_soln{1}(:,12);
save('case2_test_tau0.mat')
%%% Case 3:
% ignoreCO = false
% interactiveOH = true
% fixedOH = true
DriverScript_case3
% Give variables names
case3.concentrations = out_anal;
case3.ems = anal_soln;
case3.nh_ch4_ems = anal_soln{1}(:,1);
case3.sh_ch4_ems = anal_soln{1}(:,2);
case3.nh_co_ems = anal_soln{1}(:,13);
case3.sh_co_ems = anal_soln{1}(:,14);
case3.nh_oh_ems = anal_soln{1}(:,11);
case3.sh_oh_ems = anal_soln{1}(:,12);
save('case3_test_tau0.mat')
DriverScript_case4
case4.concentrations = out_anal;
case4.ems = anal_soln;
case4.nh_ch4_ems = anal_soln{1}(:,1);
case4.sh_ch4_ems = anal_soln{1}(:,2);
case4.nh_co_ems = anal_soln{1}(:,13);
case4.sh_co_ems = anal_soln{1}(:,14);
case4.nh_oh_ems = anal_soln{1}(:,11);
case4.sh_oh_ems = anal_soln{1}(:,12);
save('case4_test_tau0.mat')
DriverScript_case5
case5.concentrations = out_anal;
case5.ems = anal_soln;
case5.nh_ch4_ems = anal_soln{1}(:,1);
case5.sh_ch4_ems = anal_soln{1}(:,2);
case5.nh_co_ems = anal_soln{1}(:,13);
case5.sh_co_ems = anal_soln{1}(:,14);
case5.nh_oh_ems = anal_soln{1}(:,11);
case5.sh_oh_ems = anal_soln{1}(:,12);
save('case5_test_tau0.mat')