Changeset 2950 for LMDZ5/trunk/libf/phylmd/StratAer

Timestamp:

Jul 14, 2017, 5:20:51 PM (7 years ago)

Author:

oboucher

Message:

Correcting a small mistake from last commit.
k_B should be RKBOL

File:

• LMDZ5/trunk/libf/phylmd/StratAer/coagulate.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/StratAer/coagulate.F90

@@ -185 +185 @@
      Evdw=1.0
    ELSEIF (ok_vdw .EQ. 1) THEN  !--E(0) case
-     AvdWi = AvdW/(k_B*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
+     AvdWi = AvdW/(RKBOL*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
      xvdW = LOG(1.+AvdWi)
      EvdW = 1. + avdW1*xvdW + avdW3*xvdW**3
    ELSEIF (ok_vdw .EQ. 2) THEN  !--E(infinity) case
-     AvdWi = AvdW/(k_B*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
+     AvdWi = AvdW/(RKBOL*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
      xvdW = LOG(1.+AvdWi)
      EvdW = 1. + SQRT(AvdWi/3.)/(1.+bvdW0*SQRT(AvdWi)) + bvdW1*xvdW + bvdW3*xvdW**3.