-
Notifications
You must be signed in to change notification settings - Fork 2
/
40_analysis_figs.yml
104 lines (80 loc) · 3.6 KB
/
40_analysis_figs.yml
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
target_default: 40_analysis_figs
packages:
- EGRET
- EGRETci
include: 34_analyze_data_groups_monthly.yml
sources:
- 30_analyze_data/code/plot_conc_flux.R
- 40_analysis_figs/src/plot_monthly_change.R
- 40_analysis_figs/src/plot_annuals.R
- 40_analysis_figs/src/plot_change_conditions.R
targets:
40_analysis_figs:
depends:
- figures/monthly_fig_conc.png
- figures/series/tp_contour_plot.png
- figures/series/tp_conditions_panel_plot.png
- figures/series/ss_contour_plot.png
- figures/series/ss_conditions_panel_plot.png
- figures/series/ortho_p_contour_plot.png
- figures/series/ortho_p_conditions_panel_plot.png
all_targets:
command: list_targets(remake_file = '34_analyze_data_groups_monthly.yml')
monthly_targets:
command: grep(I('tp_groups_'), all_targets, value = TRUE)
figures/monthly_fig_conc.png:
command: create_monthly_fig(fig_dat = monthly_targets, fig_name = target_name, fig_height = 4, fig_width = 6, type = I('conc'))
figures/monthly_fig_percent.png:
command: create_monthly_fig(fig_dat = monthly_targets, fig_name = target_name, fig_height = 4, fig_width = 6, type = I('perc'))
figures/monthly_fig_flux.png:
command: create_monthly_fig(fig_dat = monthly_targets, fig_name = target_name, fig_height = 4, fig_width = 6, type = I('flux'))
figures/annual_tp_flux.png:
command: plot_flux_annuals(target_name, tp_wy_out, tp_wy_out_boot, unit = 8)
figures/annual_orthop_flux.png:
command: plot_flux_annuals(target_name, ortho_p_filt_wy_out, ortho_p_filt_wy_out_boot, unit = 8)
figures/annual_ss_flux.png:
command: plot_flux_annuals(target_name, ss_wy_out, ss_wy_out_boot, unit = 8)
# under which conditions figures
# this combines the contour plot + slices of that plot into multiple panels
# for each constituent
# TP
tp_breaks:
command: seq(0,0.8,0.1)
figures/series/tp_contour_plot.png:
command: plot_contours(tp_wy_out, target_name, conc_breaks = tp_breaks)
qval:
command: calc_probs(tp_wy_out, probs = I(c(0.5, 0.9, 0.99)))
figures/series/tp_conditions_panel_plot.png:
command: plot_conditions(
file_name = target_name,
eList = tp_wy_out,
year1 = I(1992), year2 = I(2005), year3 = I(2018),
date1 = I('03-15'), date2 = I('06-15'),
qlegend = I(c(0.3, 1.5)), tlegend = I(c(1991, 0.7)),
qvals = qval, conc_range1 = I(c(0.02, 2)), conc_range2 = I(c(0, 1.4)))
# SS
ss_breaks:
command: seq(0,400,50)
figures/series/ss_contour_plot.png:
command: plot_contours(ss_wy_out, target_name, ss_breaks)
figures/series/ss_conditions_panel_plot.png:
command: plot_conditions(
file_name = target_name,
eList = ss_wy_out,
year1 = I(1992), year2 = I(2005), year3 = I(2018),
date1 = I('03-15'), date2 = I('06-15'),
qlegend = I(c(0.25, 1000)), tlegend = I(c(2005, 700)),
qvals = qval, conc_range1 = I(c(20, 2000)), conc_range2 = I(c(0,800)))
# ortho_p
ortho_p_breaks:
command: seq(0,0.6,0.05)
figures/series/ortho_p_contour_plot.png:
command: plot_contours(ortho_p_filt_wy_out, target_name, ortho_p_breaks, yearStart = I(1993), yearEnd = I(2018))
figures/series/ortho_p_conditions_panel_plot.png:
command: plot_conditions(
file_name = target_name,
eList = ortho_p_filt_wy_out,
year1 = I(1993), year2 = I(2005), year3 = I(2018),
date1 = I('03-15'), date2 = I('06-15'),
qlegend = I(c(0.21, 0.9)), tlegend = I(c(2006.5, .68)),
qvals = qval, conc_range1 = I(c(0.01, 1)), conc_range2 = I(c(0,0.7)))