Index: /LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/disvert.F90
===================================================================
--- /LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/disvert.F90	(revision 1479)
+++ /LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/disvert.F90	(revision 1480)
@@ -80,9 +80,9 @@
 
      DO l = 1, llm
-        x = pi * (l - 0.5) / (llm + 1)
+        x = 2*asin(1.) * (l - 0.5) / (llm + 1)
 
         IF (ok_strato) THEN
            dsig(l) =(dsigmin + 7. * SIN(x)**2) &
-                *(1. - tanh(2 * x / pi - 1.))**2 / 4.
+                *(0.5*(1.-tanh(1.*(x-asin(1.))/asin(1.))))**2
         ELSE
            dsig(l) = 1.0 + 7.0 * SIN(x)**2
Index: /LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/disvert.F90
===================================================================
--- /LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/disvert.F90	(revision 1479)
+++ /LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/disvert.F90	(revision 1480)
@@ -80,9 +80,9 @@
 
      DO l = 1, llm
-        x = pi * (l - 0.5) / (llm + 1)
+        x = 2*asin(1.) * (l - 0.5) / (llm + 1)
 
         IF (ok_strato) THEN
            dsig(l) =(dsigmin + 7. * SIN(x)**2) &
-                *(1. - tanh(2 * x / pi - 1.))**2 / 4.
+                *(0.5*(1.-tanh(1.*(x-asin(1.))/asin(1.))))**2
         ELSE
            dsig(l) = 1.0 + 7.0 * SIN(x)**2