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