Changeset 3072

Timestamp:

Nov 15, 2017, 10:54:46 AM (7 years ago)

Author:

jyg

Message:

bug correction in alpale_th for iflag_clos_bl=3

File:

• LMDZ6/trunk/libf/phylmd/alpale_th.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/alpale_th.F90

@@ -62 +62 @@
  REAL x
  umexp(x) = max(sign(1.,x-1.e-3),0.)*(1.-exp(-x))/max(x,1.e-3) + &
-            (1.-max(sign(1.,x-1.e-3),0.))*(1.-0.5*x*(1.-x/3.*(1.-0.25*x)))
+            (1.-max(sign(1.,x-1.e-3),0.))*(-0.5*x*(1.-x/3.*(1.-0.25*x)))
+!!!            (1.-max(sign(1.,x-1.e-3),0.))*(1.-0.5*x*(1.-x/3.*(1.-0.25*x)))  !!! bug (jyg)
 !
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++