Changeset 1674 for trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/wsa_new.F90

Timestamp:

Mar 10, 2017, 3:55:34 PM (8 years ago)

Author:

slebonnois

Message:

SL: correction of a bug in microphysics

File:

• trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/wsa_new.F90 (modified) (3 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/wsa_new.F90

@@ -894 +894 @@
       REAL :: PWVSAS_GV
       REAL :: ADOT,BDOT,CDOT,DDOT
-      REAL :: MMHGPA
+      REAL :: RGAScal,MMHGPA
       REAL :: K1,K2
       REAL :: A,B,C,Dd,CP,L,F,ALFA
@@ -904 +904 @@
               DDOT=0.618943d-3)
       PARAMETER(                   &
+!     Gas constant (cal/(deg mole)):
+           RGAScal=1.98726,        &
 !     Natural logarith of conversion factor between atm. and Pa:     
            MMHGPA=11.52608845,     &
@@ -945 +947 @@
       ALFA=ALH2O(W,KHI,KLO)
 !
-      A=ADOT+(CP-K1*ALFA)/RGAS
-      B=BDOT+(L-K1*CP+K2*ALFA)/RGAS
-      C=CDOT+(CP+(F-L)/K1)/RGAS
-      Dd=DDOT-ALFA/(2.0d0*RGAS)
+      A=ADOT+(CP-K1*ALFA)/RGAScal
+      B=BDOT+(L-K1*CP+K2*ALFA)/RGAScal
+      C=CDOT+(CP+(F-L)/K1)/RGAScal
+      Dd=DDOT-ALFA/(2.0d0*RGAScal)
 !
 !     WRITE(*,*) 'TAIR= ',T,'  WSA= ',W