Changeset 3673 for trunk/LMDZ.COMMON/libf/evolution/deftank

Timestamp:

Mar 7, 2025, 10:01:32 AM (4 months ago)

Author:

jbclement

Message:

PEM:

• Correction for the update of average pressure at the end.
• Addition of pore-filled ice information for 'stratum' in the layering algorithm.
• Deactivating the soil temperature shifting because of bad values given by interpolation.

JBC

Location:

trunk/LMDZ.COMMON/libf/evolution/deftank

Files:

• visu_evol_layering.py (modified) (4 diffs)
• visu_layering.py (modified) (3 diffs)

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

@@ -67 +67 @@
     stratif_h2oice = np.zeros((ngrid,nfile,nslope,max_nb_str))
     stratif_dust = np.zeros((ngrid,nfile,nslope,max_nb_str))
-    stratif_air = np.zeros((ngrid,nfile,nslope,max_nb_str))
+    stratif_pore = np.zeros((ngrid,nfile,nslope,max_nb_str))
     for var_name in datasets[0].variables:
         if 'top_elevation' in var_name:
@@ -97 +97 @@
                             stratif_dust[i,j,k,:datasets[j].variables[var_name].shape[1]] = datasets[j].variables[var_name][0,:,i]
             print(f"Processed variable: {var_name}")
-        elif 'air_volfrac' in var_name:
+        elif 'pore_volfrac' in var_name:
             for i in range(0,ngrid):
                 for j in range(0,nfile):
                     for k in range(0,nslope):
                         if f'slope{k + 1:02d}' in var_name:
-                            stratif_air[i,j,k,:datasets[j].variables[var_name].shape[1]] = datasets[j].variables[var_name][0,:,i]
+                            stratif_pore[i,j,k,:datasets[j].variables[var_name].shape[1]] = datasets[j].variables[var_name][0,:,i]
             print(f"Processed variable: {var_name}")
 
@@ -109 +109 @@
         ds.close()
 
-    stratif_data = [stratif_heights,stratif_co2ice,stratif_h2oice,stratif_dust,stratif_air]
+    stratif_data = [stratif_heights,stratif_co2ice,stratif_h2oice,stratif_dust,stratif_pore]
     return stratif_data, max_top_elevation, longitude, latitude
 
@@ -158 +158 @@
 
 ### Figures plotting
-subtitle = ['CO2 ice','H2O ice','Dust','Air']
+subtitle = ['CO2 ice','H2O ice','Dust','Pore']
 cmap = plt.get_cmap('viridis').copy()
 cmap.set_under('white')

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

@@ -43 +43 @@
 stratif_h2oice_volfrac = []
 stratif_dust_volfrac = []
-stratif_air_volfrac = []
+stratif_pore_volfrac = []
 for i in range(1,nslope + 1):
     stratif_thickness.append(nc_file.variables['stratif_slope' + str(i).zfill(2) + '_thickness'][:])
@@ -50 +50 @@
     stratif_h2oice_volfrac.append(nc_file.variables['stratif_slope' + str(i).zfill(2) + '_h2oice_volfrac'][:])
     stratif_dust_volfrac.append(nc_file.variables['stratif_slope' + str(i).zfill(2) + '_dust_volfrac'][:])
-    stratif_air_volfrac.append(nc_file.variables['stratif_slope' + str(i).zfill(2) + '_air_volfrac'][:])
+    stratif_pore_volfrac.append(nc_file.variables['stratif_slope' + str(i).zfill(2) + '_pore_volfrac'][:])
 
 ### Display the data
 igrid = igrid - 1
 islope = islope - 1
-labels = 'CO2 ice', 'H2O ice', 'Dust', 'Air'
+labels = 'CO2 ice', 'H2O ice', 'Dust', 'Pore'
 contents = np.zeros([4,len(stratif_top_elevation[islope][0,:,igrid]) + 1])
 height = np.zeros(len(stratif_top_elevation[islope][0,:,igrid]) + 1)
@@ -64 +64 @@
     contents[1,1 + i] = stratif_h2oice_volfrac[islope][0,i,igrid]
     contents[2,1 + i] = stratif_dust_volfrac[islope][0,i,igrid]
-    contents[3,1 + i] = stratif_air_volfrac[islope][0,i,igrid]
+    contents[3,1 + i] = stratif_pore_volfrac[islope][0,i,igrid]
 contents[:,0] = contents[:,1]