| 1 | ! ------------------------------------------------------------------ |
|---|
| 2 | ! This COMDECK includes the Thermodynamical functions for the cy39 |
|---|
| 3 | ! ECMWF Physics package. |
|---|
| 4 | ! Consistent with YOMCST Basic physics constants, assuming the |
|---|
| 5 | ! partial pressure of water vapour is given by a first order |
|---|
| 6 | ! Taylor expansion of Qs(T) w.r.t. to Temperature, using constants |
|---|
| 7 | ! in YOETHF |
|---|
| 8 | ! ------------------------------------------------------------------ |
|---|
| 9 | REAL PTARG, PDELARG, P5ARG, PQSARG, PCOARG |
|---|
| 10 | REAL FOEEW, FOEDE, qsats, qsatl, dqsats, dqsatl |
|---|
| 11 | LOGICAL thermcep |
|---|
| 12 | PARAMETER (thermcep=.TRUE.) |
|---|
| 13 | ! |
|---|
| 14 | FOEEW ( PTARG,PDELARG ) = EXP ( & |
|---|
| 15 | & (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT) & |
|---|
| 16 | & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) ) |
|---|
| 17 | ! |
|---|
| 18 | FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG & |
|---|
| 19 | & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2 |
|---|
| 20 | ! |
|---|
| 21 | qsats(ptarg) = 100.0 * 0.622 * 10.0 & |
|---|
| 22 | & ** (2.07023 - 0.00320991 * ptarg & |
|---|
| 23 | & - 2484.896 / ptarg + 3.56654 * LOG10(ptarg)) |
|---|
| 24 | qsatl(ptarg) = 100.0 * 0.622 * 10.0 & |
|---|
| 25 | & ** (23.8319 - 2948.964 / ptarg & |
|---|
| 26 | & - 5.028 * LOG10(ptarg) & |
|---|
| 27 | & - 29810.16 * EXP( - 0.0699382 * ptarg) & |
|---|
| 28 | & + 25.21935 * EXP( - 2999.924 / ptarg)) |
|---|
| 29 | ! |
|---|
| 30 | dqsats(ptarg,pqsarg) = RLVTT/RCPD*pqsarg * (3.56654/ptarg & |
|---|
| 31 | & +2484.896*LOG(10.)/ptarg**2 & |
|---|
| 32 | & -0.00320991*LOG(10.)) |
|---|
| 33 | dqsatl(ptarg,pqsarg) = RLVTT/RCPD*pqsarg*LOG(10.)* & |
|---|
| 34 | & (2948.964/ptarg**2-5.028/LOG(10.)/ptarg & |
|---|
| 35 | & +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg) & |
|---|
| 36 | & +29810.16*0.0699382*EXP(-0.0699382*ptarg)) |
|---|
