-
Notifications
You must be signed in to change notification settings - Fork 0
/
horizontal_section_plot.py
141 lines (118 loc) · 5.82 KB
/
horizontal_section_plot.py
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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
"""
Plot horizontal section (lon/lat plot) of specified variable from specified file
Fill variables in "PARAM" section:
- file_msk : [string format], /path/to/file.nc
- file_in : [string format], /path/to/file.nc
- file_msk_2 : [string format], /path/to/file.nc
- file_in_2 : [string format], /path/to/file.nc
Note: If you want to make a difference between 2 files, just fill in
variables file_msk_2 and file_in_2, otherwise leave them
empty strings (e.g file_msk_2 = '' and file_in_2 = '')
- varname : [string format], variable to visualize
- month : [list format], month(s) at which variable will be plotted
Example:
month = [3] plot variable at march
month = [7, 9] plot variable averaged from july to september
month = [] plot variable averaged over all time steps
- depth_in : [list format], depth indice(s) at which variable will be plotted
Example:
depth_ind = [0] plot variable at surface
depth_ind = [5, 13] plot variable integrated from depth indices 5 to 12
(Remember in python when slicing, indice defining the end is excluded)
depth_ind = [] plot variable integrated over all depth indices
- cmap_color : [string format], colormap color
- cmap_range : [list format], defines bounds and steps of the different
intervals in the custom colormap
[[val_min, val_1, step_1], [...], [val_n, val_max, step_n]]
Example:
cmap_range = [[0, 10, 1],[10,31,5]]
cmap_range = [[0, 10, 1]]
- contour_lines : [list format], defines contour lines to display for the plot
[contour_min, contour_max, step]
- projection : [string format], projection to use of the plot
- lon_lim : [float format], longitude limits for plot
- lat_lim : [float format], latitude limits for plot
"""
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from tools.process_grid import disp_grid_coastline
from tools.process_oce import (get_coords_h,
get_unit,
get_proj,
get_var_h,
compute_land_mask,
param4plot,
plot_var_h,
plot_cbar,
plot_cbar_legend_h)
# PARAMS
#-----------------------------------------------------------------------#
# Files to load
file_msk = ''
file_in = ''
file_msk_2 = '' # default: ''
file_in_2 = '' # default: ''
# Variable name
varname = ''
# time and depth selection
month = [7,9] # ex: [1] => jan / [3,6] => ave mar to Jun / [] => ave over all time steps
depth_ind = [0] # ex: [0] => ind 0 / [6,20] => ave ind 6 to 19 / [] => ave over all depth indices
# Color colormap
cmap_color = 'Spectral_r' # default: 'viridis', 'Spectral_r', 'RdBu_r', 'BrBG_r'
# Custom colormap and contour
cmap_range = [] # default: [] - format: [[v_min, v_1, stp_1], [...], [v_n, v_max, stp_n]]
contour_lines = [] # default: [] - format: [contour_min, contour_max, step]
# Type of projection
projection = 'polar_S' # default: 'robinson', 'polar_N', 'polar_S'
# Zoom
lon_lim = [-180,180] # default: [-180,180]
lat_lim = [-90,0] # default: [-90,90]
#-----------------------------------------------------------------------#
# Extract time indice(s) from month variable
#-----------------------------------------------------------------------#
# Switch to indices for time instead of month
time_ind = [i-1 for i in month]
# In python when slicing, last indices is excluded
# So we add this lines to take in account last month specified by user
if len(time_ind) == 2:
time_ind[-1] = time_ind[-1] + 1
#-----------------------------------------------------------------------#
# Compute variables
#-----------------------------------------------------------------------#
coords = get_coords_h(file_in)
var = get_var_h(file_in,
file_msk,
varname,
time_ind,
depth_ind)
coords_msk, land_mask = compute_land_mask(file_msk)
unit = get_unit(file_in, varname)
if file_in_2:
coords_2 = get_coords_h(file_in_2)
var_2 = get_var_h(file_in_2,
file_msk_2,
varname,
time_ind,
depth_ind)
coords_msk_2, land_mask_2 = compute_land_mask(file_msk_2)
var2plot = var - var_2
else:
var2plot = var
#-----------------------------------------------------------------------#
# PLOT PARAMS
#-----------------------------------------------------------------------#
projDataOut = get_proj(projection)
cmap_param = param4plot(var2plot, cmap_color, cmap_range)
projDataIn = ccrs.PlateCarree()
plot_limits = [lon_lim, lat_lim]
#-----------------------------------------------------------------------#
# PLOT
#-----------------------------------------------------------------------#
fig, ax = plt.subplots(figsize=(10, 10), subplot_kw={'projection': projDataOut})
mapp = plot_var_h(ax, projDataIn, projDataOut, plot_limits, coords, var2plot, contour_lines, cmap_param)
cbar = plot_cbar(ax, mapp, cmap_param)
plot_cbar_legend_h(varname, unit, cbar, time_ind, depth_ind)
disp_grid_coastline(ax, coords_msk[0], coords_msk[1], projDataIn, land_mask[0, ::], 1.25)
if file_in_2:
disp_grid_coastline(ax, coords_msk_2[0], coords_msk_2[1], projDataIn, land_mask_2[0, ::], 1.25, line_color=(0.5,0.2,0.2))
#-----------------------------------------------------------------------#