-
Notifications
You must be signed in to change notification settings - Fork 2
/
run_plot_VALSO.bash
executable file
·93 lines (78 loc) · 4.44 KB
/
run_plot_VALSO.bash
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
#!/bin/bash
if [ $# -eq 0 ] ; then echo 'need a [KEYWORD] (will be inserted inside the output name), [FREQ] (1y or 1m) and a list of id [RUNIDS RUNID ...] (definition of line style need to be done in RUNID.db)'; exit; fi
KEY=${1}
FREQ=${2}
RUNIDS=${@:3}
. ~/.bashrc
. PARAM/param_arch.bash
load_python
# ACC
# Drake
echo 'plot ACC time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *ACC*${FREQ}*.nc -var vtrp -sf -1 -title "ACC transport (Sv)" -dir ${WRKPATH} -o "${KEY}_ACC" -obs OBS/ACC_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# GYRE
# ROSS GYRE
echo 'plot Ross Gyre time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *RG*${FREQ}*psi.nc -var max_sobarstf -title "Ross Gyre (Sv)" -dir ${WRKPATH} -o ${KEY}_RG -sf 0.000001 -obs OBS/RG_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# WED GYRE
echo 'plot Weddell Gyre time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *WG*${FREQ}*psi.nc -var max_sobarstf -title "Weddell Gyre (Sv)" -dir ${WRKPATH} -o ${KEY}_WG -sf 0.000001 -obs OBS/WG_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# HSSW
# mean S WROSS
echo 'plot mean bot S (WROSS) time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *WROSS*so*${FREQ}*T.nc -var '(mean_sosbs|mean_vosaline)' -title "Mean bot. sal. WROSS (PSU)" -dir ${WRKPATH} -o ${KEY}_WROSS_mean_bot_so -obs OBS/WROSS_botS_mean_obs.txt
# mean S WWED
if [[ $? -ne 0 ]]; then exit 42; fi
echo 'plot mean bot S (WWED) time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *WWED*so*${FREQ}*T.nc -var '(mean_sosbs|mean_vosaline)' -title "Mean bot. sal. WWED (PSU)" -dir ${WRKPATH} -o ${KEY}_WWED_mean_bot_so -obs OBS/WWED_botS_mean_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# CDW
# mean T AMU
echo 'plot mean bot T (AMU) time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *AMU_*thetao*${FREQ}*T.nc -var '(mean_sosbt|mean_votemper)' -title "Mean bot. temp. AMU (C)" -dir ${WRKPATH} -o ${KEY}_AMU_mean_bot_thetao -obs OBS/AMU_botT_mean_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# mean T EROSS
echo 'plot mean bot T (EROSS) time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *EROSS*thetao*${FREQ}*T.nc -var '(mean_sosbt|mean_votemper)' -title "Mean bot. temp. EROSS (C)" -dir ${WRKPATH} -o ${KEY}_EROSS_mean_bot_thetao -obs OBS/EROSS_botT_mean_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# mean T EROSS
echo 'plot mean bot T (EWED) time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *EWED*thetao*${FREQ}*T.nc -var '(mean_sosbt|mean_votemper)' -title "Mean bot. temp. EWED (C)" -dir ${WRKPATH} -o ${KEY}_EWED_mean_bot_thetao -obs OBS/EWED_botT_mean_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# MLD
# max mld in WEDDELL GYRE
echo 'plot max mld in Weddell Gyre time series'
python SCRIPT/plot_time_series.py -noshow -runid $RUNIDS -f *WMXL*1m*m09_*T.nc -var '(max_sokaraml|max_somxzint1|max_somxl010)' -title "Max Kara mld WG (m)" -dir ${WRKPATH} -o ${KEY}_WG_max_karamld -obs OBS/WG_karamld_max_obs.txt
if [[ $? -ne 0 ]]; then exit 42; fi
# crop figure (rm legend)
convert ${KEY}_WG.png -crop 1240x1040+0+0 tmp01.png
convert ${KEY}_RG.png -crop 1240x1040+0+0 tmp02.png
convert ${KEY}_ACC.png -crop 1240x1040+0+0 tmp03.png
convert ${KEY}_WWED_mean_bot_so.png -crop 1240x1040+0+0 tmp04.png
convert ${KEY}_WROSS_mean_bot_so.png -crop 1240x1040+0+0 tmp05.png
convert ${KEY}_EWED_mean_bot_thetao.png -crop 1240x1040+0+0 tmp06.png
convert ${KEY}_EROSS_mean_bot_thetao.png -crop 1240x1040+0+0 tmp07.png
convert ${KEY}_AMU_mean_bot_thetao.png -crop 1240x1040+0+0 tmp08.png
convert ${KEY}_WG_max_karamld.png -crop 1240x1040+0+0 tmp09.png
# trim figure (remove white area)
convert FIGURES/box_VALSO.png -trim -bordercolor White -border 40 tmp10.png
convert legend.png -trim -bordercolor White -border 20 tmp11.png
convert runidname.png -trim -bordercolor White -border 20 tmp12.png
# compose the image
convert \( tmp01.png tmp02.png tmp03.png +append \) \
\( tmp04.png tmp05.png tmp10.png +append \) \
\( tmp06.png tmp07.png tmp08.png +append \) \
\( tmp09.png +append \) \
tmp11.png tmp12.png -append -trim -bordercolor White -border 40 $KEY.png
# save figure
mv ${KEY}_*.png FIGURES/.
mv ${KEY}_*.txt FIGURES/.
mv tmp10.png FIGURES/${KEY}_legend.png
mv tmp11.png FIGURES/${KEY}_runidname.png
# clean
rm tmp??.png
#display
#display -resize 30% $KEY.png