Changeset 1480 for LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar

Timestamp:

Jan 31, 2011, 10:29:58 PM (14 years ago)

Author:

jghattas

Message:

Put back modification done in revision 1472. This modification was done
to simplify the understanding of the code but it did also slightly change
the results. With this commit we got back the same results as in the revision
before 1472.

File:

• LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/disvert.F90 (modified) (1 diff)

LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3dpar/disvert.F90

@@ -80 +80 @@
 
      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