source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/watersat_grad.F90

Last change on this file was 222, checked in by ymipsl, 10 years ago

Creating temporary dynamico/lmdz/saturn branche

YM
File size: 1.0 KB
Line 
1subroutine watersat_grad(T,qsat,dqsat)
2
3  use watercommon_h, only: T_h2O_ice_liq, RLVTT, RCPD,T_coup
4  implicit none
5
6!==================================================================
7!     Purpose
8!     -------
9!     Compute the L/cp*d (q_sat)/d T
10!     for a given temperature (K)
11!
12!     Authors
13!     -------
14!     Robin Wordsworth (2010)
15!
16!==================================================================
17
18!   input
19  real T,qsat
20 
21!   output
22  real dqsat
23
24!  if (T.lt.T_coup) then ! solid / vapour !why use T_coup?????????? JL12
25  if (T.lt.T_h2O_ice_liq) then ! solid / vapour
26     dqsat = RLVTT/RCPD*qsat*(3.56654/T             &
27          +2484.896*LOG(10.)/T**2                   &
28          -0.00320991*LOG(10.))
29  else                 ! liquid / vapour
30     dqsat = RLVTT/RCPD*qsat*LOG(10.)*              &
31          (2948.964/T**2-5.028/LOG(10.)/T           &
32          +25.21935*2999.924/T**2*EXP(-2999.924/T)  &
33          +29810.16*0.0699382*EXP(-0.0699382*T))
34  end if
35
36  return
37end subroutine watersat_grad
38
Note: See TracBrowser for help on using the repository browser.