Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed

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

    r2949 r2950  
    185185     Evdw=1.0
    186186   ELSEIF (ok_vdw .EQ. 1) THEN  !--E(0) case
    187      AvdWi = AvdW/(k_B*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
     187     AvdWi = AvdW/(RKBOL*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
    188188     xvdW = LOG(1.+AvdWi)
    189189     EvdW = 1. + avdW1*xvdW + avdW3*xvdW**3
    190190   ELSEIF (ok_vdw .EQ. 2) THEN  !--E(infinity) case
    191      AvdWi = AvdW/(k_B*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
     191     AvdWi = AvdW/(RKBOL*t_seri(ilon,ilev))*(4.*radius(i)*radius(j))/(radius(i)+radius(j))**2.
    192192     xvdW = LOG(1.+AvdWi)
    193193     EvdW = 1. + SQRT(AvdWi/3.)/(1.+bvdW0*SQRT(AvdWi)) + bvdW1*xvdW + bvdW3*xvdW**3.
Note: See TracChangeset for help on using the changeset viewer.