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