--- trunk/dyn3d/disvert.f90	2014/02/05 17:51:07	78
+++ trunk/dyn3d/disvert.f	2014/03/12 21:16:36	90
@@ -6,24 +6,24 @@
 
   private llm
 
-  real ap(llm+1), pa ! in Pa
-  real bp(llm+1)
-  real presnivs(llm) ! pressions approximatives des milieux de couches, en Pa
-  real, parameter:: preff = 101325. ! in Pa
-  real nivsigs(llm), nivsig(llm+1)
+  real, save:: ap(llm+1), pa ! in Pa
+  real, save:: bp(llm+1)
+
+  real, save:: presnivs(llm)
+  ! pressions approximatives des milieux de couches, en Pa
 
-  save
+  real, parameter:: preff = 101325. ! in Pa
 
 contains
 
   SUBROUTINE disvert
 
-    ! From dyn3d/disvert.F, v 1.1.1.1 2004/05/19 12:53:05
+    ! From dyn3d/disvert.F, version 1.1.1.1, 2004/05/19 12:53:05
     ! Author: P. Le Van
 
     ! This procedure sets the vertical grid. It defines the host
-    ! variables "ap", "bp", "presnivs", "nivsigs" and "nivsig". "pa"
-    ! should be defined before this procedure is called.
+    ! variables "ap", "bp", "presnivs". "pa" should be defined before
+    ! this procedure is called.
 
     use jumble, only: new_unit
     use nr_util, only: pi, assert
@@ -36,7 +36,7 @@
     ! the interface between layers "l" and "l-1"
 
     real ds(llm)
-    ! "ds(l)" : épaisseur de la couche "l" dans la coordonnée "s"
+    ! "ds(l)" : \'epaisseur de la couche "l" dans la coordonn\'ee "s"
 
     INTEGER l, unit
     REAL alpha, x(llm)
@@ -48,7 +48,7 @@
     ! used only if vert_sampling == "param"
 
     ! These variables are used only in the case vert_sampling == "param":
-    real:: deltaz = 0.04 ! épaisseur de la première couche
+    real:: deltaz = 0.04 ! \'epaisseur de la premi\`ere couche
     real:: beta = 1.3 ! facteur d'accroissement en haut
     real:: k0 = 20. ! nombre de couches dans la transition surface
     real:: k1 = 1.2 ! nombre de couches dans la transition haute
@@ -59,9 +59,6 @@
 
     print *, "Call sequence information: disvert"
 
-    forall (l = 1: llm) nivsigs(l) = REAL(l)
-    forall (l = 1: llm + 1) nivsig(l) = REAL(l)
-
     print *, "Enter namelist 'disvert_nml'."
     read(unit=*, nml=disvert_nml)
     write(unit_nml, nml=disvert_nml)