1 |
|
1 |
SUBROUTINE SURRTRF |
2 |
|
|
|
3 |
|
|
! Adapted from Eli J. Mlawer, Atmospheric & Environmental Research. |
4 |
|
|
! by JJMorcrette, ECMWF |
5 |
|
|
! ------------------------------------------------------------------ |
6 |
|
|
|
7 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
8 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
9 |
|
|
|
10 |
|
|
USE YOERRTRF , ONLY : PREF ,PREFLOG ,TREF |
11 |
|
|
|
12 |
|
|
! ------------------------------------------------------------------ |
13 |
|
|
|
14 |
|
|
! These pressures are chosen such that the ln of the first pressure |
15 |
|
|
! has only a few non-zero digits (i.e. ln(PREF(1)) = 6.96000) and |
16 |
|
|
! each subsequent ln(pressure) differs from the previous one by 0.2. |
17 |
|
|
|
18 |
|
|
IMPLICIT NONE |
19 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
20 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SURRTRF',0,ZHOOK_HANDLE) |
21 |
|
|
PREF( :) = (/& |
22 |
|
|
& 1.05363E+03_JPRB,8.62642E+02_JPRB,7.06272E+02_JPRB,5.78246E+02_JPRB,4.73428E+02_JPRB,& |
23 |
|
|
& 3.87610E+02_JPRB,3.17348E+02_JPRB,2.59823E+02_JPRB,2.12725E+02_JPRB,1.74164E+02_JPRB,& |
24 |
|
|
& 1.42594E+02_JPRB,1.16746E+02_JPRB,9.55835E+01_JPRB,7.82571E+01_JPRB,6.40715E+01_JPRB,& |
25 |
|
|
& 5.24573E+01_JPRB,4.29484E+01_JPRB,3.51632E+01_JPRB,2.87892E+01_JPRB,2.35706E+01_JPRB,& |
26 |
|
|
& 1.92980E+01_JPRB,1.57998E+01_JPRB,1.29358E+01_JPRB,1.05910E+01_JPRB,8.67114E+00_JPRB,& |
27 |
|
|
& 7.09933E+00_JPRB,5.81244E+00_JPRB,4.75882E+00_JPRB,3.89619E+00_JPRB,3.18993E+00_JPRB,& |
28 |
|
|
& 2.61170E+00_JPRB,2.13828E+00_JPRB,1.75067E+00_JPRB,1.43333E+00_JPRB,1.17351E+00_JPRB,& |
29 |
|
|
& 9.60789E-01_JPRB,7.86628E-01_JPRB,6.44036E-01_JPRB,5.27292E-01_JPRB,4.31710E-01_JPRB,& |
30 |
|
|
& 3.53455E-01_JPRB,2.89384E-01_JPRB,2.36928E-01_JPRB,1.93980E-01_JPRB,1.58817E-01_JPRB,& |
31 |
|
|
& 1.30029E-01_JPRB,1.06458E-01_JPRB,8.71608E-02_JPRB,7.13612E-02_JPRB,5.84256E-02_JPRB,& |
32 |
|
|
& 4.78349E-02_JPRB,3.91639E-02_JPRB,3.20647E-02_JPRB,2.62523E-02_JPRB,2.14936E-02_JPRB,& |
33 |
|
1 |
& 1.75975E-02_JPRB,1.44076E-02_JPRB,1.17959E-02_JPRB,9.65769E-03_JPRB/) |
34 |
|
|
|
35 |
|
|
PREFLOG( :) = (/& |
36 |
|
|
& 6.9600E+00_JPRB, 6.7600E+00_JPRB, 6.5600E+00_JPRB, 6.3600E+00_JPRB, 6.1600E+00_JPRB,& |
37 |
|
|
& 5.9600E+00_JPRB, 5.7600E+00_JPRB, 5.5600E+00_JPRB, 5.3600E+00_JPRB, 5.1600E+00_JPRB,& |
38 |
|
|
& 4.9600E+00_JPRB, 4.7600E+00_JPRB, 4.5600E+00_JPRB, 4.3600E+00_JPRB, 4.1600E+00_JPRB,& |
39 |
|
|
& 3.9600E+00_JPRB, 3.7600E+00_JPRB, 3.5600E+00_JPRB, 3.3600E+00_JPRB, 3.1600E+00_JPRB,& |
40 |
|
|
& 2.9600E+00_JPRB, 2.7600E+00_JPRB, 2.5600E+00_JPRB, 2.3600E+00_JPRB, 2.1600E+00_JPRB,& |
41 |
|
|
& 1.9600E+00_JPRB, 1.7600E+00_JPRB, 1.5600E+00_JPRB, 1.3600E+00_JPRB, 1.1600E+00_JPRB,& |
42 |
|
|
& 9.6000E-01_JPRB, 7.6000E-01_JPRB, 5.6000E-01_JPRB, 3.6000E-01_JPRB, 1.6000E-01_JPRB,& |
43 |
|
|
& -4.0000E-02_JPRB,-2.4000E-01_JPRB,-4.4000E-01_JPRB,-6.4000E-01_JPRB,-8.4000E-01_JPRB,& |
44 |
|
|
& -1.0400E+00_JPRB,-1.2400E+00_JPRB,-1.4400E+00_JPRB,-1.6400E+00_JPRB,-1.8400E+00_JPRB,& |
45 |
|
|
& -2.0400E+00_JPRB,-2.2400E+00_JPRB,-2.4400E+00_JPRB,-2.6400E+00_JPRB,-2.8400E+00_JPRB,& |
46 |
|
|
& -3.0400E+00_JPRB,-3.2400E+00_JPRB,-3.4400E+00_JPRB,-3.6400E+00_JPRB,-3.8400E+00_JPRB,& |
47 |
|
1 |
& -4.0400E+00_JPRB,-4.2400E+00_JPRB,-4.4400E+00_JPRB,-4.6400E+00_JPRB/) |
48 |
|
|
|
49 |
|
|
! These are the temperatures associated with the respective |
50 |
|
|
! pressures for the MLS standard atmosphere. |
51 |
|
|
TREF( :) = (/& |
52 |
|
|
& 2.9420E+02_JPRB, 2.8799E+02_JPRB, 2.7894E+02_JPRB, 2.6925E+02_JPRB, 2.5983E+02_JPRB,& |
53 |
|
|
& 2.5017E+02_JPRB, 2.4077E+02_JPRB, 2.3179E+02_JPRB, 2.2306E+02_JPRB, 2.1578E+02_JPRB,& |
54 |
|
|
& 2.1570E+02_JPRB, 2.1570E+02_JPRB, 2.1570E+02_JPRB, 2.1706E+02_JPRB, 2.1858E+02_JPRB,& |
55 |
|
|
& 2.2018E+02_JPRB, 2.2174E+02_JPRB, 2.2328E+02_JPRB, 2.2479E+02_JPRB, 2.2655E+02_JPRB,& |
56 |
|
|
& 2.2834E+02_JPRB, 2.3113E+02_JPRB, 2.3401E+02_JPRB, 2.3703E+02_JPRB, 2.4022E+02_JPRB,& |
57 |
|
|
& 2.4371E+02_JPRB, 2.4726E+02_JPRB, 2.5085E+02_JPRB, 2.5457E+02_JPRB, 2.5832E+02_JPRB,& |
58 |
|
|
& 2.6216E+02_JPRB, 2.6606E+02_JPRB, 2.6999E+02_JPRB, 2.7340E+02_JPRB, 2.7536E+02_JPRB,& |
59 |
|
|
& 2.7568E+02_JPRB, 2.7372E+02_JPRB, 2.7163E+02_JPRB, 2.6955E+02_JPRB, 2.6593E+02_JPRB,& |
60 |
|
|
& 2.6211E+02_JPRB, 2.5828E+02_JPRB, 2.5360E+02_JPRB, 2.4854E+02_JPRB, 2.4348E+02_JPRB,& |
61 |
|
|
& 2.3809E+02_JPRB, 2.3206E+02_JPRB, 2.2603E+02_JPRB, 2.2000E+02_JPRB, 2.1435E+02_JPRB,& |
62 |
|
|
& 2.0887E+02_JPRB, 2.0340E+02_JPRB, 1.9792E+02_JPRB, 1.9290E+02_JPRB, 1.8809E+02_JPRB,& |
63 |
|
1 |
& 1.8329E+02_JPRB, 1.7849E+02_JPRB, 1.7394E+02_JPRB, 1.7212E+02_JPRB/) |
64 |
|
|
|
65 |
|
|
! ----------------------------------------------------------------- |
66 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('SURRTRF',1,ZHOOK_HANDLE) |
67 |
|
1 |
END SUBROUTINE SURRTRF |