plopm -i results/MODEL -v 'depth * 0.001' -s ,,1:217 -agg min -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -cbl 'Depth [km]' -cbn 5 -cbf .2f -c managua_r -hide 0,0,0,1 -fn figure5b &
plopm -i 'results/MODEL results/MODEL results/MODEL results/FILLED_MODEL results/FILLED_MODEL results/FILLED_MODEL' -v 'porv' -s ',,1:30 ,,57:111 ,,117:217 ,,1:30 ,,57:111 ,,117:217' -r 0 -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -sg 2,3 -fs 24,24 -cbp 0.1,0.95,0.8,0.02 -cbf .1f -fz 24 -cbn 5 -st 0 -cbl 'Pore volume [m$^3$]' -clog 1 -t 'Cook, porv = 1.06e10 [m$^3$]  Johansen, porv = 3.76e10 [m$^3$]  Statfjord, porv = 8.17e10 [m$^3$]  Cook, porv = 1.21e10 [m$^3$] (+14\%)  Johansen, porv = 5.75e10 [m$^3$] (+53\%)  Statfjord, porv = 1.43e11 [m$^3$] (+75\%)' -rdl 1 -c cet_diverging_rainbow_bgymr_45_85_c67 -fn figure6 &
plopm -i results/FILLED_MODEL -v pvtnum -s ,,: -r 0 -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -u opm -hide 0,0,1,1 -c '203;203;203' -fn figure7 &
plopm -i 'results/FILLED_MODEL results/STANDARD_COARSENED_MODEL results/DUAL_VERTICAL_TF_COARSENED_MODEL results/DUAL_NO_VERTICAL_TF_COARSENED_MODEL' -v 'pressure - 0pressure' -s ',,1:217 ,,1:5 ,,1:5 ,,1:5' -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -sg 1,4 -fs 32,12 -cbp 0.1,0.95,0.8,0.02 -fz 24 -cbn 5 -st 0 -cbl 'Pore-volume-weighted-average pressure increase ($\bar{p}^m_{(i,j)}$) [bar], end of simulation (t=525 years)' -rdl 1 -dg 0,0,1,1 -r 30 -fn year_525 -t 'Filled model  Standard coarsened model  Dual (vertical TF) coarsened model   Dual (no vertical TF) coarsened model' -fn figure8 &
plopm -i 'results/FILLED_MODEL' -v 'pressure - 0pressure' -s ',,1:217 ,,1:5 ' -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -fs 12,12 -fz 24 -cbn 5 -st 0 -cbl 'Pore-volume-weighted-average pressure difference [bar] (coarsened)' -r 30 -di results/STANDARD_COARSENED_MODEL -hide 0,0,1,0 -t '$E^{standard}_{global}$=1.83 [bar]' -cl '[-18,18]' -fn figure9a &
plopm -i 'results/FILLED_MODEL' -v 'pressure - 0pressure' -s ',,1:217 ,,1:5' -dg 0,1 -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -fs 12,12 -fz 24 -cbn 5 -st 0 -cbl 'Pore-volume-weighted-average pressure difference [bar] (coarsened)' -r 30 -di results/DUAL_VERTICAL_TF_COARSENED_MODEL -hide 1,0,1,0 -t '$E^{dual\;(vertical\;TF)}_{global}$=1.52 [bar]' -cl '[-18,18]' -fn figure9b &
plopm -i 'results/FILLED_MODEL' -v 'pressure - 0pressure' -s ',,1:217 ,,1:5' -dg 0,1 -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -fs 12,12 -fz 24 -cbn 5 -st 0 -cbl 'Pore-volume-weighted-average pressure difference [bar] (coarsened)' -r 30 -di results/DUAL_NO_VERTICAL_TF_COARSENED_MODEL -hide 1,0,1,0 -t '$E^{dual\;(no\;vertical\;TF)}_{global}$=1.39 [bar]'  -cl '[-18,18]' -fn figure9c &
plopm -i 'results/STANDARD_COARSENED_MODEL results/STANDARD_COARSENED_MODEL results/STANDARD_COARSENED_MODEL' -v 'pressure - 0pressure' -s ',,1:5' -tr '[-495000,-6.605e6]' -x '[0,95000]' -xu km -xnt 2 -y '[0,160000]' -yu km -yf .0f -ynt 2 -xf .0f -fz 24 -cbn 5 -st 0 -cbl 'Pore-volume-weighted-average pressure difference ($E^{m}_{diff,(i,j)}$) [bar], end of simulation (t=525 years)' -r 30 -cl '[-18,18]' -sg 1,3 -fs 24,12 -cbp 0.1,0.95,0.8,0.02 -c RdYlGn -fn to_extract_colorbar_for_figure9 &
wait