Index: trunk/LMDZ.MARS/libf/dyn3d/comvert_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/comvert_mod.F90	(revision 1594)
+++ trunk/LMDZ.MARS/libf/dyn3d/comvert_mod.F90	(revision 1621)
@@ -10,4 +10,5 @@
 ! Mars Ce qui suit vient de gcm
 REAL sig(llm+1),ds(llm),aps(llm),bps(llm),pseudoalt(llm)
+REAL scaleheight ! atmospheric scale height (km)
 
 END MODULE comvert_mod
Index: trunk/LMDZ.MARS/libf/dyn3d/disvert.F
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/disvert.F	(revision 1594)
+++ trunk/LMDZ.MARS/libf/dyn3d/disvert.F	(revision 1621)
@@ -5,5 +5,5 @@
 
       USE comvert_mod, ONLY: ap,bp,sig,nivsigs,nivsig,pa,preff,
-     .                  aps,bps,presnivs,pseudoalt
+     .                  aps,bps,presnivs,pseudoalt,scaleheight
       USE comconst_mod, ONLY: kappa,pi
       USE logic_mod, ONLY: hybrid
@@ -25,5 +25,5 @@
       INTEGER l,ll
       REAL snorm
-      REAL alpha,beta,gama,delta,deltaz,h,quoi,quand
+      REAL alpha,beta,gama,delta,deltaz,quoi,quand
       REAL zsig(llm)
       INTEGER np,ierr
@@ -70,5 +70,5 @@
          PRINT*,'WARNING Lecture de esasig.def'
          PRINT*,'*****************************'
-         READ(99,*) h
+         READ(99,*) scaleheight
          READ(99,*) dz0
          READ(99,*) dz1
@@ -76,6 +76,6 @@
          CLOSE(99)
 
-         dz0=dz0/h
-         dz1=dz1/h
+         dz0=dz0/scaleheight
+         dz1=dz1/scaleheight
 
          sig1=(1.-dz1)/tanh(.5*(llm-1)/nhaut)
@@ -120,5 +120,5 @@
          PRINT*,'****************************'
 
-         READ(99,*) h
+         READ(99,*) scaleheight
          do l=1,llm
             read(99,*) zsig(l)
@@ -128,5 +128,6 @@
          sig(1) =1
          do l=2,llm
-           sig(l) = 0.5 * ( exp(-zsig(l)/h) + exp(-zsig(l-1)/h) )
+           sig(l) = 0.5 * ( exp(-zsig(l)/scaleheight) + 
+     &                      exp(-zsig(l-1)/scaleheight) )
          end do
          sig(llm+1) =0
