Changeset 3458 for trunk

Timestamp:

Oct 11, 2024, 11:32:49 AM (6 weeks ago)

Author:

jbclement

Message:

PCM:
Addition of a python script to visualize the layerings over time and on a reference grid for elevation which is specified by the user.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• deftank/visu_evol_layering.py (added)
• deftank/visu_layering.py (modified) (3 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -442 +442 @@
 == 11/10/2024 == JBC
 Modification of dimension detection for the variables written in "diagpem.nc": in particular for 'nb_str_max' which can evolve and match the value of other dimensions.
+
+== 11/10/2024 == JBC
+Addition of a python script to visualize the layerings over time and on a reference grid for elevation which is specified by the user.

trunk/LMDZ.COMMON/libf/evolution/deftank/visu_layering.py

@@ -11 +11 @@
 ##############################
 ### Parameters to fill in
-filename = 'startpem9.nc' # File name
+filename = 'startpem.nc' #'starts/restartpem226.nc' # File name
 igrid = 1  # Grid point
 islope = 1 # Slope number
@@ -17 +17 @@
 ##############################
 
+
+####################################################################################
 ### Open the NetCDF file
 if os.path.isfile(filename):
@@ -83 +85 @@
 plt.savefig('layering_simpleprofiles.png')
 
-# Stackplot for a layering
+# Content profiles for a layering
+plt.figure()
+plt.step(contents[0,:],height,where = 'post',color = 'r',label = labels[0])
+#plt.plot(contents[0,:],height,'o--',color = 'r',alpha = 0.3)
+plt.step(contents[1,:],height,where = 'post',color = 'b',label = labels[1])
+#plt.plot(contents[1,:],height,'o--',color = 'b',alpha = 0.3)
+plt.step(contents[2,:],height,where = 'post',color = 'y',label = labels[2])
+#plt.plot(contents[2,:],height,'o--',color = 'y',alpha = 0.3)
+plt.step(contents[3,:],height,where = 'post',color = 'g',label = labels[3])
+#plt.plot(contents[3,:],height,'o--',color = 'g',alpha = 0.3)
+plt.grid(axis='x', color='0.95')
+plt.grid(axis='y', color='0.95')
+plt.xlim(0,1)
+plt.xlabel('Volume fraction [m3/m3]')
+plt.ylabel('Elevation [m]')
+plt.title('Content profiles for the layering')
+plt.savefig('layering_simpleprofiles.png')
+
+# Stack content profiles for a layering
 fig, ax = plt.subplots()
 ax.fill_betweenx(height,0,contents[0,:],label = labels[0],color = 'r',step = 'pre')