Index: trunk/LMDZ.COMMON/libf/evolution/atmosphere.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/atmosphere.F90	(revision 4134)
+++ trunk/LMDZ.COMMON/libf/evolution/atmosphere.F90	(revision 4138)
@@ -408,5 +408,5 @@
 
 !=======================================================================
-SUBROUTINE compute_alt_coord(pa,preff,ap,bp,aps,bps)
+SUBROUTINE compute_alt_coord(pa,preff,ap_l,bp_l,aps,bps)
 !-----------------------------------------------------------------------
 ! NAME
@@ -449,8 +449,8 @@
 ! ---------
 real(dp)                                               :: pa, preff
-real(dp), dimension(nlayer + 1), intent(out)           :: ap, bp
+real(dp), dimension(nlayer + 1), intent(out)           :: ap_l, bp_l
 real(dp), dimension(nlayer),     intent(out), optional :: aps, bps
 
-! LOCAL VAIRABLES
+! LOCAL VARIABLES
 ! ---------------
 integer(di)                     :: i
@@ -471,11 +471,11 @@
         call compute_hybrid_sig(sig(i),pa,preff,newsig)
         bp(i) = exp(1._dp - 1._dp/(newsig**2))
-        ap(i) = pa*(newsig - bp(i))
+        ap_l(i) = pa*(newsig - bp(i))
     end do
-    ap(nlayer + 1) = 0._dp
+    ap_l(nlayer + 1) = 0._dp
     bp(nlayer + 1) = 0._dp
 else
     call print_msg('> Defining sigma altitude coordinates',LVL_NFO)
-    ap(:) = 0._dp
+    ap_l(:) = 0._dp
     bp(1:nlayer) = sig(1:nlayer)
     bp(nlayer + 1) = 0._dp
@@ -487,5 +487,5 @@
 
 do i = 1,nlayer - 1
-    aps(i) =  0.5_dp*(ap(i) + ap(i + 1))
+    aps(i) =  0.5_dp*(ap_l(i) + ap_l(i + 1))
     bps(i) =  0.5_dp*(bp(i) + bp(i + 1))
 end do