Skip to content
Snippets Groups Projects
Commit 7902a18d authored by Lionel GUEZ's avatar Lionel GUEZ
Browse files

Polish

parent e71d1111
No related branches found
No related tags found
No related merge requests found
......@@ -53,15 +53,17 @@ else:
# (in [llcrnrlon, llcrnrlon + 2 pi[)
lon_mask = longitude <= urcrnrlon
lat_mask = np.logical_and(latitude >= llcrnrlat,
latitude <= urcrnrlat)
lat_mask = np.logical_and(latitude >= llcrnrlat, latitude <= urcrnrlat)
longitude = longitude[lon_mask]
latitude = latitude[lat_mask]
src_crs = ccrs.PlateCarree()
# Use a conformal projection for quiver:
projection = ccrs.Stereographic(central_latitude = latitude.mean(),
central_longitude = longitude.mean())
##projection = ccrs.NorthPolarStereo()
fig = plt.figure()
ax = plt.axes(projection = projection)
......@@ -75,8 +77,7 @@ else:
quiver_return = wind_cartopy.plot(ax, longitude, latitude,
ugos[lat_mask][:, lon_mask],
vgos[lat_mask][:, lon_mask],
scale = args.scale,
scale_units = "width")
scale = args.scale, scale_units = "width")
ax.quiverkey(quiver_return, 0.9, 0.9, 1, r"1 m s$^{-1}$",
coordinates = "figure")
......
......@@ -1279,7 +1279,7 @@ Minimum cible pour le test \verb+Get_1_outerm+ en :
Indices dans le domaine global, base 1, de ce minimum : 39, 225. SSH
en ce mimimum : \np{0.2759}. Cyclone numéro 4 dans ce cadre. Contour
le plus intérieur : SSH = \np{0.293300003} m. Contour le plus
extérieur : SSH = \np{0.460744} m. Contour de vitesse maximale : SSH =
extérieur : SSH = \np{0.460504} m. Contour de vitesse maximale : SSH =
\np{0.3815} m. Cf. figures (\ref{fig:test_get_eddy}) et
(\ref{fig:test_get_snapshot}).
\begin{figure}[htbp]
......
......@@ -245,11 +245,7 @@
"&config_nml /\n"
],
"commands": [
[
"mkdir",
"-p",
"SHPC/Slice_0/Cyclones"
],
["mkdir", "-p", "SHPC/Slice_0/Cyclones"],
[
"$build_dir/Inst_eddies/Tests/test_set_max_speed",
"$src_dir/Inst_eddies/Tests/Input/Region_1/huv_region_1_2006_01_01.nc",
......
......@@ -57,7 +57,6 @@ program test_get_1_outerm
integer n_cont ! number of good contours found and stored in cont_list
TYPE(shpfileobject) hshp_cont_list
namelist /main_nml/ ind_targ_extr, innermost_level, cyclone, date, &
grid_lon_lat
......
......@@ -32,7 +32,8 @@ def loop_inst_eddies(files, bbox, d, my_slice):
if os.access("config_nml.txt", os.R_OK):
print("inst_eddies.loop_inst_eddies: Will use config_nml.txt.")
else:
sys.exit('inst_eddies.loop_inst_eddies: "config_nml.txt" not found')
sys.exit('inst_eddies.loop_inst_eddies: "config_nml.txt" not found '
'-- Aborting')
nco_instance = nco.Nco()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment