Index: trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh	(revision 3777)
+++ trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh	(revision 3790)
@@ -333,4 +333,7 @@
         cleanfiles starts/restartpem .nc $(($iPEM - 1))
         rm -f startpem.nc
+        if [ -f "starts/startpem.nc" ]; then
+            cp starts/startpem.nc .
+        fi
         if [ $irelaunch -eq $(($nPCM_ini - 1)) ]; then
             cp diags/data2reshape${irelaunch}.nc data2reshape_Y1.nc
Index: trunk/LMDZ.COMMON/libf/evolution/deftank/visu_evol_layering.py
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/deftank/visu_evol_layering.py	(revision 3777)
+++ trunk/LMDZ.COMMON/libf/evolution/deftank/visu_evol_layering.py	(revision 3790)
@@ -106,5 +106,5 @@
         'dust':      'h_dust',
         'pore':      'h_pore',
-        'pore_ice':  'icepore_volfrac'
+        'pore_ice':  'poreice_volfrac'
     }
     var_info = {prop: set() for prop in property_markers}
Index: trunk/LMDZ.COMMON/libf/evolution/deftank/visu_layering.py
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/deftank/visu_layering.py	(revision 3777)
+++ trunk/LMDZ.COMMON/libf/evolution/deftank/visu_layering.py	(revision 3790)
@@ -55,6 +55,6 @@
         'h_h2o': f"stratif_slope{idx_str}_h_h2oice",
         'h_dust': f"stratif_slope{idx_str}_h_dust",
-        'h_air': f"stratif_slope{idx_str}_h_pore",
-        'volfrac': f"stratif_slope{idx_str}_icepore_volfrac",
+        'h_pore': f"stratif_slope{idx_str}_h_pore",
+        'poreice_volfrac': f"stratif_slope{idx_str}_poreice_volfrac",
     }
 
@@ -82,6 +82,6 @@
     h_h2o = data['h_h2o']
     h_dust = data['h_dust']
-    h_air = data['h_air']
-    volfrac = data['volfrac']
+    h_pore = data['h_pore']
+    poreice_volfrac = data['poreice_volfrac']
 
     layers = top.shape[1]
@@ -94,5 +94,5 @@
         + h_h2o[0, 0, grid_index]
         + h_dust[0, 0, grid_index]
-        + h_air[0, 0, grid_index]
+        + h_pore[0, 0, grid_index]
     )
     raw_height[0] = top[0, 0, grid_index] - total_thickness0
@@ -104,5 +104,5 @@
             + h_h2o[0, i, grid_index]
             + h_dust[0, i, grid_index]
-            + h_air[0, i, grid_index]
+            + h_pore[0, i, grid_index]
         )
 
@@ -135,7 +135,7 @@
         h2o = h_h2o[0, layer_idx, grid_index] / thickness
         dust = h_dust[0, layer_idx, grid_index] / thickness
-        air = h_air[0, layer_idx, grid_index] * (1.0 - volfrac[0, layer_idx, grid_index]) / thickness
-        pore = h_air[0, layer_idx, grid_index] * volfrac[0, layer_idx, grid_index] / thickness
-        contents[:, idx + 1] = [co2, h2o, dust, air, pore]
+        air = h_pore[0, layer_idx, grid_index] * (1.0 - poreice_volfrac[0, layer_idx, grid_index]) / thickness
+        poreice = h_pore[0, layer_idx, grid_index] * poreice_volfrac[0, layer_idx, grid_index] / thickness
+        contents[:, idx + 1] = [co2, h2o, dust, air, poreice]
 
     if num_filt > 0:
@@ -276,10 +276,10 @@
         h_h2o = slope_data['h_h2o']
         h_dust = slope_data['h_dust']
-        h_air = slope_data['h_air']
+        h_pore = slope_data['h_pore']
         total_thickness0 = (
             h_co2[0, 0, igrid_input]
             + h_h2o[0, 0, igrid_input]
             + h_dust[0, 0, igrid_input]
-            + h_air[0, 0, igrid_input]
+            + h_pore[0, 0, igrid_input]
         )
         raw_height[0] = top[0, 0, igrid_input] - total_thickness0