source: LMDZ6/branches/LMDZ-QUEST/libf/phylmd/rrtm/susw15.F90

Last change on this file was 3043, checked in by Ehouarn Millour, 7 years ago

Some fixes to enable compiling with the pgf compiler and updated arch files for pgf on Ciclad.
EM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 5.0 KB
Line 
1!OPTIONS XOPT(NOEVAL)
2SUBROUTINE SUSW15
3
4!**** *SUSW15*   - INITIALIZE COMMON YOMSW15
5!****              FROZEN VERSION (CYCLE 15) OF SUSW
6
7!     PURPOSE.
8!     --------
9!           INITIALIZE YOMSW15, THE COMMON THAT CONTAINS COEFFICIENTS
10!           NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES
11
12!**   INTERFACE.
13!     ----------
14!        *CALL* *SUSW15
15
16!        EXPLICIT ARGUMENTS :
17!        --------------------
18!        NONE
19
20!        IMPLICIT ARGUMENTS :
21!        --------------------
22!        COMMON YOMSW15
23
24!     METHOD.
25!     -------
26!        SEE DOCUMENTATION
27
28!     EXTERNALS.
29!     ----------
30
31!     REFERENCE.
32!     ----------
33!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
34
35!     AUTHOR.
36!     -------
37!        96-11: Ph. Dandin. Meteo-France
38!        ORIGINAL : 88-12-15 BY JEAN-JACQUES MORCRETTE *ECMWF*
39
40!     MODIFICATIONS.
41!     --------------
42!        M.Hamrud      01-Oct-2003 CY28 Cleaning
43
44!     ------------------------------------------------------------------
45
46USE PARKIND1  ,ONLY : JPIM     ,JPRB
47USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
48
49USE YOMCST   , ONLY : RG
50USE YOMSW15  , ONLY : APAD15   ,BPAD15   ,RRAY15   ,RSUN15   ,&
51 & RPDH115  ,RPDU115  ,RPNH15   ,RPNU15   ,RSWCE15  ,&
52 & RSWCP15  ,RTDH2O15 ,RTDUMG15 ,RTH2O15  ,RTUMG15  ,D15 
53
54!      ----------------------------------------------------------------
55
56!*       1.    SET VALUES.
57!              -----------
58
59IMPLICIT NONE
60REAL(KIND=JPRB),SAVE   :: ZPDH2O=0.8_JPRB
61REAL(KIND=JPRB),SAVE   :: ZPDUMG=0.75_JPRB
62REAL(KIND=JPRB),SAVE   :: ZPRH2O=30000._JPRB
63REAL(KIND=JPRB),SAVE   :: ZPRUMG=30000._JPRB
64
65!$OMP THREADPRIVATE(zpdh2o,zpdumg,zprh2o,zprumg)
66
67REAL(KIND=JPRB) :: ZH2O, ZUMG
68REAL(KIND=JPRB) :: ZHOOK_HANDLE
69
70IF (LHOOK) CALL DR_HOOK('SUSW15',0,ZHOOK_HANDLE)
71RPDH115=ZPDH2O+1.0_JPRB
72RPDU115=ZPDUMG+1.0_JPRB
73ZH2O=1.0_JPRB/( 10._JPRB* RG * RPDH115 )
74ZUMG=1.0_JPRB/( 10._JPRB* RG * RPDU115 )
75RPNU15 = ZUMG/(ZPRUMG**ZPDUMG)
76RPNH15 = ZH2O/(ZPRH2O**ZPDH2O)
77
78RSWCP15=0.002_JPRB*RSWCE15
79
80RTDH2O15 = 0.40_JPRB
81RTDUMG15 = 0.375_JPRB
82RTH2O15 = 240._JPRB
83RTUMG15 = 240._JPRB
84RSWCE15 = 0._JPRB
85RSWCP15 = 0._JPRB
86
87RSUN15(1) = 0.441676_JPRB
88D15(1,:) = (/ 0.00_JPRB, 0.00_JPRB, 0.00_JPRB /)
89!* DERIVED FROM HITRAN APRIL 1991
90!       H2O:  Pref=300 hPa, Tref=240K, Pdep=0.8
91!       O3 :  unchanged
92
93APAD15(1,1:3,1:7) = RESHAPE((/&
94 & 0.912418292E+05_JPRB, 0.000000000E-00_JPRB, 0.925887084E-04_JPRB,&
95 & 0.723613782E+05_JPRB, 0.000000000E-00_JPRB, 0.129353723E-01_JPRB,&
96 & 0.596037057E+04_JPRB, 0.000000000E-00_JPRB, 0.800821928E+00_JPRB,&
97 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.242715973E+02_JPRB,&
98 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.878331486E+02_JPRB,&
99 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.191559725E+02_JPRB,&
100 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.000000000E+00_JPRB /)&
101 & ,(/3,7/)) 
102
103BPAD15(1,1:3,1:7) = RESHAPE((/&
104 & 0.912418292E+05_JPRB, 0.000000000E-00_JPRB, 0.925887084E-04_JPRB,&
105 & 0.724555318E+05_JPRB, 0.000000000E-00_JPRB, 0.131812683E-01_JPRB,&
106 & 0.602593328E+04_JPRB, 0.000000000E-00_JPRB, 0.812706117E+00_JPRB,&
107 & 0.100000000E+01_JPRB, 0.000000000E-00_JPRB, 0.249863591E+02_JPRB,&
108 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.931071925E+02_JPRB,&
109 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.252233437E+02_JPRB,&
110 & 0.000000000E-00_JPRB, 0.000000000E-00_JPRB, 0.100000000E+01_JPRB /)&
111 & ,(/3,7/)) 
112
113RRAY15(1,1:6) = (/&
114 & .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
115 & .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/) 
116
117RSUN15(2) = 0.558324_JPRB
118
119D15(2,1:3) = (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.800000000_JPRB /)
120
121!* INTERVAL 2:  0.68 - 4.00 MICRONS
122!* DERIVED FROM HITRAN APRIL 1991
123!       H2O:  Pref=300 hPa, Tref=240K, Pdep=0.80
124!       UMG:  Pref=300 hPa, Tref=240K, Pdep=0.75 (CO2+O2+CH4+N2O+CO)
125!       O3 :  unchanged
126
127APAD15(2,1:3,1:7) = RESHAPE((/&
128 & 0.376655383E-08_JPRB, 0.739646016E-08_JPRB, 0.410177786E+03_JPRB,&
129 & 0.978576773E-04_JPRB, 0.131849595E-03_JPRB, 0.672595424E+02_JPRB,&
130 & 0.387714006E+00_JPRB, 0.437772681E+00_JPRB, 0.000000000E-00_JPRB,&
131 & 0.118461660E+03_JPRB, 0.151345118E+03_JPRB, 0.000000000E-00_JPRB,&
132 & 0.119079797E+04_JPRB, 0.233628890E+04_JPRB, 0.000000000E-00_JPRB,&
133 & 0.293353397E+03_JPRB, 0.797219934E+03_JPRB, 0.000000000E-00_JPRB,&
134 & 0.000000000E+00_JPRB, 0.000000000E+00_JPRB, 0.000000000E+00_JPRB /)&
135 & ,(/3,7/)) 
136
137BPAD15(2,1:3,1:7) = RESHAPE((/&
138 & 0.376655383E-08_JPRB, 0.739646016E-08_JPRB, 0.410177786E+03_JPRB,&
139 & 0.979023421E-04_JPRB, 0.131861712E-03_JPRB, 0.731185438E+02_JPRB,&
140 & 0.388611139E+00_JPRB, 0.437949001E+00_JPRB, 0.100000000E+01_JPRB,&
141 & 0.120291383E+03_JPRB, 0.151692730E+03_JPRB, 0.000000000E+00_JPRB,&
142 & 0.130531005E+04_JPRB, 0.237071130E+04_JPRB, 0.000000000E+00_JPRB,&
143 & 0.415049409E+03_JPRB, 0.867914360E+03_JPRB, 0.000000000E+00_JPRB,&
144 & 0.100000000E+01_JPRB, 0.100000000E+01_JPRB, 0.000000000E+00_JPRB /)&
145 & ,(/3,7/)) 
146
147RRAY15(2,1:6) = (/&
148 & .697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,&
149 & .248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/) 
150IF (LHOOK) CALL DR_HOOK('SUSW15',1,ZHOOK_HANDLE)
151END SUBROUTINE SUSW15
Note: See TracBrowser for help on using the repository browser.