Index: trunk/LMDZ.PLUTO/libf/phypluto/phyredem.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/phyredem.F90	(revision 3732)
+++ trunk/LMDZ.PLUTO/libf/phypluto/phyredem.F90	(revision 3735)
@@ -14,5 +14,5 @@
                        emisice, emissiv,             &
                        iceradius, dtemisice, phisfi
-  use iostart, only : open_restartphy, close_restartphy, & 
+  use iostart, only : open_restartphy, close_restartphy, &
                       put_var, put_field, length
   use mod_grid_phy_lmdz, only : klon_glo
@@ -41,7 +41,43 @@
   real,intent(in) :: pzgam(ngrid)
   real,intent(in) :: pzthe(ngrid)
- 
+
+  character(ndscrpt), dimension(ldscrpt), parameter :: dscrpt_tab_cntrl = (/ &
+      "(1)  Number of atmospheric columns in physics     ", &
+      "(2)  Number of atmospheric layers                 ", &
+      "(3)  Final day                                    ", &
+      "(4)  Final time of day                            ", &
+      "(5)  Planet radius (m)                            ", &
+      "(6)  Rotation rate (rad.s-1)                      ", &
+      "(7)  Gravity (m.s-2)                              ", &
+      "(8)  Molar mass of the atmosphere (g.mol-1)       ", &
+      "(9)  = r/Cp           (=kappa in the dynamics)    ", &
+      "(10) Length of a solar day (s)                    ", &
+      "(11) Physics time step (s)                        ", &
+      "(12) -                                            ", &
+      "(13) -                                            ", &
+      "(14) Length of year (in solar days)               ", &
+      "(15) Minimum star-planet distance (AU)            ", &
+      "(16) Maximum star-planet distance (AU)            ", &
+      "(17) Date of periastro (sols since N. spring)     ", &
+      "(18) Obliquity of the planet (deg)                ", &
+      "(19) Default surface roughness (m)                ", &
+      "(20) -                                            ", &
+      "(21) -                                            ", &
+      "(22) -                                            ", &
+      "(23) -                                            ", &
+      "(24) Emissivity of northern cap ~0.95             ", &
+      "(25) Emissivity of southern cap ~0.95             ", &
+      "(26) Emissivity of martian soil ~.95              ", &
+      "(27) -                                            ", &
+      "(28) -                                            ", &
+      "(29) -                                            ", &
+      "(30) -                                            ", &
+      "(31) Mean scat radius of CO2 snow (north)         ", &
+      "(32) Mean scat radius of CO2 snow (south)         ", &
+      "(33) Time scale for snow metamorphism (north)     ", &
+      "(34) Time scale for snow metamorphism (south)     ", &
+      "(35) Soil volumetric heat capacity                "/)
   real :: tab_cntrl(length) ! nb "length=100" defined in iostart module
-  
+
   ! Create physics start file
   call open_restartphy(filename)
@@ -55,5 +91,5 @@
   tab_cntrl(1) = float(klon_glo)  ! number of nodes on physics grid
   tab_cntrl(2) = float(nlay) ! number of atmospheric layers
-  tab_cntrl(3) = day_ini + int(time)         ! final day 
+  tab_cntrl(3) = day_ini + int(time)         ! final day
   tab_cntrl(4) = time -int(time)            ! final time of day
 
@@ -73,5 +109,5 @@
   tab_cntrl(14) = year_day  ! length of year (sols) ~668.6
   tab_cntrl(15) = periastr  ! min. star-planet distance (AU)
-  tab_cntrl(16) = apoastr   ! max. star-planet distance (AU) 
+  tab_cntrl(16) = apoastr   ! max. star-planet distance (AU)
   tab_cntrl(17) = peri_day  ! date of periastron (sols since N. spring)
   tab_cntrl(18) = obliquit  ! Obliquity of the planet (deg) ~23.98
@@ -91,5 +127,5 @@
   tab_cntrl(34) = dtemisice(2) ! time scale for snow metamorphism (south)
 
-  tab_cntrl(28) = 0. 
+  tab_cntrl(28) = 0.
   tab_cntrl(29) = 0.
   tab_cntrl(30) = 0.
@@ -98,23 +134,27 @@
 
   call put_var("controle","Control parameters",tab_cntrl)
-  
+
+  ! Write the controle array descriptor
+  call put_var(nid_restart,"controle_descriptor",&
+               "Description of control parameters",dscrpt_tab_cntrl)
+
   ! Write the mid-layer depths
   call put_var("soildepth","Soil mid-layer depth",mlayer)
-  
+
   ! Write longitudes
   call put_field("longitude","Longitudes of physics grid",lonfi)
-  
+
   ! Write latitudes
   call put_field("latitude","Latitudes of physics grid",latfi)
-  
+
   ! Write mesh areas
   call put_field("area","Mesh area",cell_area)
-  
+
   ! Write surface geopotential
   call put_field("phisfi","Geopotential at the surface",phisfi)
-  
+
   ! Write surface albedo
   !call put_field("albedodat","Albedo of bare ground",alb)
-  
+
   ! Subgrid topogaphy variables
   call put_field("ZMEA","Relief: mean relief",zmea)
@@ -123,8 +163,8 @@
   call put_field("ZGAM","Relief: gamma parameter",zgam)
   call put_field("ZTHE","Relief: theta parameter",zthe)
-  
+
   ! Close file
   call close_restartphy
-  
+
 end subroutine physdem0
 
@@ -133,5 +173,5 @@
                     emis,alb,q2,qsurf)
   ! write time-dependent variable to restart file
-  use iostart, only : open_restartphy, close_restartphy, & 
+  use iostart, only : open_restartphy, close_restartphy, &
                       put_var, put_field
   use tracer_h, only: noms
@@ -155,5 +195,5 @@
 
   integer :: iq
-  
+
   ! Open file
   call open_restartphy(filename)
@@ -162,5 +202,5 @@
   ! set time counter in file
   !call put_var("Time","Temps de simulation",time)
-  
+
   ! Surface temperature
   call put_field("tsurf","Surface temperature",tsurf)
@@ -168,8 +208,8 @@
   ! Soil inertia
   call put_field("inertiedat","Soil thermal inertia",inertiesoil)
-  
+
   ! Soil temperature
   call put_field("tsoil","Soil temperature",tsoil)
-  
+
   ! Emissivity of the surface
   call put_field("emis","Surface emissivity",emis)
@@ -177,8 +217,8 @@
   ! Albedo of the surface
   call put_field("albedodat","Albedo of bare ground",alb)
-  
+
   ! Planetary Boundary Layer
   call put_field("q2","pbl wind variance",q2)
-  
+
 ! cloud fraction and sea ice !AF24: removed