| 1 | ! |
|---|
| 2 | ! $Header$ |
|---|
| 3 | ! |
|---|
| 4 | ! |
|---|
| 5 | ! ATTENTION!!!!: ce fichier include est compatible format fixe/format libre |
|---|
| 6 | ! veillez n'utiliser que des ! pour les commentaires |
|---|
| 7 | ! et bien positionner les & des lignes de continuation |
|---|
| 8 | ! (les placer en colonne 6 et en colonne 73) |
|---|
| 9 | ! |
|---|
| 10 | ! ------------------------------------------------------------------ |
|---|
| 11 | ! This COMDECK includes the Thermodynamical functions for the cy39 |
|---|
| 12 | ! ECMWF Physics package. |
|---|
| 13 | ! Consistent with YOMCST Basic physics constants, assuming the |
|---|
| 14 | ! partial pressure of water vapour is given by a first order |
|---|
| 15 | ! Taylor expansion of Qs(T) w.r.t. to Temperature, using constants |
|---|
| 16 | ! in YOETHF |
|---|
| 17 | ! ------------------------------------------------------------------ |
|---|
| 18 | REAL PTARG, PDELARG, P5ARG, PQSARG, PCOARG |
|---|
| 19 | REAL FOEEW, FOEDE, qsats, qsatl, dqsats, dqsatl |
|---|
| 20 | LOGICAL thermcep |
|---|
| 21 | PARAMETER (thermcep=.TRUE.) |
|---|
| 22 | ! |
|---|
| 23 | FOEEW ( PTARG,PDELARG ) = EXP ( & |
|---|
| 24 | & (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT) & |
|---|
| 25 | & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) ) |
|---|
| 26 | ! |
|---|
| 27 | FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG & |
|---|
| 28 | & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2 |
|---|
| 29 | ! |
|---|
| 30 | qsats(ptarg) = 100.0 * 0.622 * 10.0 & |
|---|
| 31 | & ** (2.07023 - 0.00320991 * ptarg & |
|---|
| 32 | & - 2484.896 / ptarg + 3.56654 * LOG10(ptarg)) |
|---|
| 33 | qsatl(ptarg) = 100.0 * 0.622 * 10.0 & |
|---|
| 34 | & ** (23.8319 - 2948.964 / ptarg & |
|---|
| 35 | & - 5.028 * LOG10(ptarg) & |
|---|
| 36 | & - 29810.16 * EXP( - 0.0699382 * ptarg) & |
|---|
| 37 | & + 25.21935 * EXP( - 2999.924 / ptarg)) |
|---|
| 38 | ! |
|---|
| 39 | dqsats(ptarg,pqsarg) = RLVTT/RCPD*pqsarg * (3.56654/ptarg & |
|---|
| 40 | & +2484.896*LOG(10.)/ptarg**2 & |
|---|
| 41 | & -0.00320991*LOG(10.)) |
|---|
| 42 | dqsatl(ptarg,pqsarg) = RLVTT/RCPD*pqsarg*LOG(10.)* & |
|---|
| 43 | & (2948.964/ptarg**2-5.028/LOG(10.)/ptarg & |
|---|
| 44 | & +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg) & |
|---|
| 45 | & +29810.16*0.0699382*EXP(-0.0699382*ptarg)) |
|---|
