Ignore:
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.

Location:
LMDZ5/branches/LMDZ5V2.0-dev/libf
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/LMDZ5V2.0-dev/libf/dyn3d/disvert.F90

    r1472 r1480  
    8080
    8181     DO l = 1, llm
    82         x = pi * (l - 0.5) / (llm + 1)
     82        x = 2*asin(1.) * (l - 0.5) / (llm + 1)
    8383
    8484        IF (ok_strato) THEN
    8585           dsig(l) =(dsigmin + 7. * SIN(x)**2) &
    86                 *(1. - tanh(2 * x / pi - 1.))**2 / 4.
     86                *(0.5*(1.-tanh(1.*(x-asin(1.))/asin(1.))))**2
    8787        ELSE
    8888           dsig(l) = 1.0 + 7.0 * SIN(x)**2

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

    r1472 r1480  
    8080
    8181     DO l = 1, llm
    82         x = pi * (l - 0.5) / (llm + 1)
     82        x = 2*asin(1.) * (l - 0.5) / (llm + 1)
    8383
    8484        IF (ok_strato) THEN
    8585           dsig(l) =(dsigmin + 7. * SIN(x)**2) &
    86                 *(1. - tanh(2 * x / pi - 1.))**2 / 4.
     86                *(0.5*(1.-tanh(1.*(x-asin(1.))/asin(1.))))**2
    8787        ELSE
    8888           dsig(l) = 1.0 + 7.0 * SIN(x)**2
Note: See TracChangeset for help on using the changeset viewer.