--- trunk/dyn3d/disvert.f	2014/03/05 14:57:53	82
+++ trunk/dyn3d/disvert.f	2014/03/06 15:12:00	83
@@ -6,13 +6,13 @@
 
   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
 
@@ -22,8 +22,8 @@
     ! 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
@@ -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)