Changeset 3824 for trunk/LMDZ.MARS/util

Timestamp:

Jul 2, 2025, 2:58:14 PM (9 days ago)

Author:

jbclement

Message:

Mars PCM:
Small corrections to the "display_netcdf.py" script in the "util" folder.
JBC

Location:

trunk/LMDZ.MARS/util

Files:

• display_netcdf.py (modified) (13 diffs)
• display_netcdf.yml (modified) (2 diffs)

trunk/LMDZ.MARS/util/display_netcdf.py

@@ -210 +210 @@
             linewidth=0.5
         )
-        #gl.top_labels = False
-        #gl.right_labels = False
 
         # Plot data in PlateCarree projection
@@ -229 +227 @@
         label = varname + (f" ({units})" if units else "")
         cbar.set_label(label)
-        ax.set_title(f"{varname} — {pole.capitalize()} Pole", pad=50)
+        ax.set_title(f"{varname} — {pole.capitalize()} polar region", pad=20, y=1.05, fontsize=12, fontweight='bold')
 
         figs.append(fig)
@@ -362 +360 @@
 
     # Create plotter
-    plotter = Plotter(title="3D topography view", bg="bb", axes=0)
+    plotter = Plotter(title="3D globe view", bg="bb", axes=0)
 
     # Configure camera
@@ -409 +407 @@
         plt.xlabel(time_var or "Time Index")
         plt.ylabel(varname + (f" ({var.units})" if hasattr(var, "units") else ""))
-        plt.title(f"{varname} vs {time_var or 'Index'}")
+        plt.title(f"{varname} vs {time_var or 'Index'}", fontweight='bold')
         if output_path:
             plt.savefig(output_path, bbox_inches="tight")
@@ -442 +440 @@
         cbar = plt.colorbar()
         cbar.set_label(varname + (f" ({var.units})" if hasattr(var, "units") else ""))
-        plt.title(f"{varname} averaged over latitude")
+        plt.title(f"{varname} averaged over latitude", fontweight='bold')
         if output_path:
             plt.savefig(output_path, bbox_inches="tight")
@@ -564 +562 @@
                     cbar = fig.colorbar(cf, ax=ax, pad=0.02)
                     cbar.set_label(varname + (f" ({var.units})" if hasattr(var, "units") else ""))
-                    ax.set_title(f"{varname} (interpolated map over physical_points)")
+                    ax.set_title(f"{varname} (interpolated map over physical_points)", fontweight='bold')
                     ax.set_xlabel(f"Longitude ({getattr(dataset.variables[lonv], 'units', 'deg')})")
                     ax.set_ylabel(f"Latitude ({getattr(dataset.variables[latv], 'units', 'deg')})")
@@ -600 +598 @@
                 plt.xlabel(varname + (f" ({var.units})" if hasattr(var, "units") else ""))
                 plt.ylabel(coord + (f" ({dataset.variables[coord].units})" if hasattr(dataset.variables[coord], "units") else ""))
-                plt.title(f"{varname} vs {coord}")
+                plt.title(f"{varname} vs {coord}", fontweight='bold')
                 if output_path:
                     plt.savefig(output_path, bbox_inches="tight")
@@ -612 +610 @@
         plt.xlabel("Index")
         plt.ylabel(varname + (f" ({var.units})" if hasattr(var, "units") else ""))
-        plt.title(f"{varname} (1D)")
+        plt.title(f"{varname} (1D)", fontweight='bold')
         if output_path:
             plt.savefig(output_path, bbox_inches="tight")
@@ -620 +618 @@
         return
 
-    # if dslice.ndim == 2:
+    if dslice.ndim == 2:
         lat_idx2 = find_dim_index(dims, LAT_DIMS)
         lon_idx2 = find_dim_index(dims, LON_DIMS)
@@ -630 +628 @@
             lats = dataset.variables[latv][:]
             lons = dataset.variables[lonv][:]
+
+            # Correct latitudes order
+            if lats[0] > lats[-1]:
+                lats = lats[::-1]
+                dslice = np.flipud(dslice)
 
             # Handle masked arrays
@@ -675 +678 @@
             cbar.set_label(varname + (f" ({dataset.variables[varname].units})"
                                       if hasattr(dataset.variables[varname], "units") else ""))
-            ax.set_title(f"{varname} (lat × lon)")
+            ax.set_title(f"{varname} (lat × lon)", fontweight='bold')
             ax.set_xlabel(f"Longitude ({getattr(dataset.variables[lonv], 'units', 'deg')})")
             ax.set_ylabel(f"Latitude ({getattr(dataset.variables[latv], 'units', 'deg')})")
@@ -702 +705 @@
         plt.xlabel("Dim 2 index")
         plt.ylabel("Dim 1 index")
-        plt.title(f"{varname} (2D)")
+        plt.title(f"{varname} (2D)", fontweight='bold')
         if output_path:
             plt.savefig(output_path, bbox_inches="tight")
@@ -891 +894 @@
         plt.xlabel(time_var or "Time Index")
         plt.ylabel(varname + (f" ({var_obj.units})" if hasattr(var_obj, "units") else ""))
-        plt.title(f"{varname} vs {time_var or 'Index'}")
+        plt.title(f"{varname} vs {time_var or 'Index'}", fontweight='bold')
         if output_path:
             plt.savefig(output_path, bbox_inches="tight")

trunk/LMDZ.MARS/util/display_netcdf.yml

@@ -163 +163 @@
   - mpg123=1.32.9=hc50e24c_0
   - msgpack-python=1.1.1=py312h68727a3_0
-  - multidict=6.6.2=py312h178313f_0
+  - multidict=6.6.3=py312h178313f_0
   - munkres=1.1.4=pyhd8ed1ab_1
   - ncurses=6.5=h2d0b736_3
   - netcdf4=1.7.2=nompi_py312h3805cb1_102
   - nlohmann_json=3.12.0=h3f2d84a_0
-  - numpy=2.3.0=py312h6cf2f7f_0
+  - numpy=2.3.1=py312h6cf2f7f_0
   - ocl-icd=2.3.3=hb9d3cd8_0
   - opencl-headers=2025.06.13=h5888daf_0
@@ -174 +174 @@
   - openjpeg=2.5.3=h5fbd93e_0
   - openldap=2.6.10=he970967_0
-  - openssl=3.5.0=h7b32b05_1
+  - openssl=3.5.1=h7b32b05_0
   - packaging=25.0=pyh29332c3_1
   - pandas=2.3.0=py312hf9745cd_0