source: LMDZ6/trunk/libf/phylmd/ecrad/susrtm.F90 @ 4190

Last change on this file since 4190 was 3908, checked in by idelkadi, 3 years ago

Online implementation of the radiative transfer code ECRAD in the LMDZ model.

  • Inclusion of the ecrad directory containing the sources of the ECRAD code
    • interface routine : radiation_scheme.F90
  • Adaptation of compilation scripts :
    • compilation under CPP key CPP_ECRAD
    • compilation with option "-rad ecard" or "-ecard true"
    • The "-rad old/rtm/ecran" build option will need to replace the "-rrtm true" and "-ecrad true" options in the future.
  • Runing LMDZ simulations with ecrad, you need :
    • logical key iflag_rrtm = 2 in physiq.def
    • namelist_ecrad (DefLists?)
    • the directory "data" containing the configuration files is temporarily placed in ../libfphylmd/ecrad/
  • Compilation and execution are tested in the 1D case. The repository under svn would allow to continue the implementation work: tests, verification of the results, ...
File size: 17.8 KB
RevLine 
[3908]1SUBROUTINE SUSRTM
2
3!     Adapted from E.J. Mlawer, J. Delamere, Atmospheric & Environmental Research.
4!     by JJMorcrette, ECMWF
5!     Modified to add arrays relevant to mapping for g-point reduction,
6!     M.J. Iacono, Atmospheric & Environmental Research, Inc.
7!     JJMorcrette 20010610 Flexible configuration for number of g-points
8!     ------------------------------------------------------------------
9
10USE PARKIND1  ,ONLY : JPRB ,   JPIM
11USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
12
13USE YOESRTM  , ONLY : JPGPT, NGBSW, NGN
14USE YOESRTWN , ONLY : NG      , NSPA, NSPB   , NMPSRTM, &
15 & PREF    , PREFLOG , TREF   , &
16 & NGM     , WT      , NGC    , NGS
17! & NGM     , WT      , NGC    , NGS , NGN    , NGBSW
18! & WAVENUM1, WAVENUM2, DELWAVE, PREF, PREFLOG, TREF   , &
19
20!     ------------------------------------------------------------------
21
22IMPLICIT NONE
23
24INTEGER(KIND=JPIM) :: IGC56(14), IGC112(14) , IGC224(14)
25INTEGER(KIND=JPIM) :: IGS56(14), IGS112(14) , IGS224(14)
26
27INTEGER(KIND=JPIM) :: IGM56(224),IGM112(224), IGM224(224)
28
29INTEGER(KIND=JPIM) :: IGN56(56), IGN112(112), IGN224(224)
30INTEGER(KIND=JPIM) :: IGB56(56), IGB112(112), IGB224(224)
31
32REAL(KIND=JPRB) :: ZHOOK_HANDLE
33!-----------------------------------------------------------------------
34IF (LHOOK) CALL DR_HOOK('SUSRTM',0,ZHOOK_HANDLE)
35
36NG(:)     =(/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16 /)
37NSPA(:)   =(/  9, 9, 9, 9, 1, 9, 9, 1, 9, 1, 0, 1, 9, 1 /)
38NSPB(:)   =(/  1, 5, 1, 1, 1, 5, 1, 0, 1, 0, 0, 1, 5, 1 /)
39NMPSRTM(:)=(/  6, 6, 5, 5, 5, 5, 5, 4, 4, 3, 2, 2, 1, 6 /)
40
41!WAVENUM1( :) = (/&
42! & 2600._JPRB, 3250._JPRB, 4000._JPRB, 4650._JPRB, 5150._JPRB, 6150._JPRB, 7700._JPRB &
43! & , 8050._JPRB,12850._JPRB,16000._JPRB,22650._JPRB,29000._JPRB,38000._JPRB,  820._JPRB /) 
44!WAVENUM2( :) = (/&
45! & 3250._JPRB, 4000._JPRB, 4650._JPRB, 5150._JPRB, 6150._JPRB, 7700._JPRB, 8050._JPRB &
46! & ,12850._JPRB,16000._JPRB,22650._JPRB,29000._JPRB,38000._JPRB,50000._JPRB, 2600._JPRB /) 
47!DELWAVE( :) = (/&
48! & 650._JPRB,  750._JPRB,  650._JPRB,  500._JPRB, 1000._JPRB, 1550._JPRB,  350._JPRB &
49! & , 4800._JPRB, 3150._JPRB, 6650._JPRB, 6350._JPRB, 9000._JPRB,12000._JPRB, 1780._JPRB /) 
50
51!=====================================================================
52! Set arrays needed for the g-point reduction from 224 to
53! - either 112 for the high-resolution forecast model configuration
54! - or 56 for the EPS-type configuration 
55! in the 14 SW bands:
56
57! NB: This mapping from 224 to 112 points has been carefully selected to
58! minimize the effect on the resulting fluxes and cooling rates, and
59! caution should be used if the mapping is modified.
60!     The further reduction to 56 for EPS configuration is considered
61! acceptable, only because of the random perturbations introduced on
62! the total heating rates produced by the physical parametrization package.
63! While a reduction to 56 obviously speeds up the model, it as obviously
64! reduces the accuracy that could be expected from the radiation scheme.
65
66! JPGPT   The total number of new g-points (NGPT)
67! NGC     The number of new g-points in each band (14)
68! NGS     The cumulative sum of new g-points for each band (14)
69! NGM     The index of each new g-point relative to the original
70!         16 g-points for each band.
71! NGN     The number of original g-points that are combined to make
72!         each new g-point in each band.
73! NGB     The band index for each new g-point.
74! WT      RRTM weights for 16 g-points. (16)
75
76!-- ECMWF EPS model RRTM_SW configuration with 56 g-points
77IGC56(:) = (/ 3, 6, 4, 4, 5, 5, 1, 5, 4, 3, 3, 4, 3, 6 /)
78IGS56(:) = (/ 3, 9,13,17,22,27,28,33,37,40,43,47, 50, 56 /)
79
80IGM56(:) = (/ 1,1,1,1,2,2,2,2,3,3,3,3,3,3,3,3, &            ! Band 16
81            & 1,1,2,2,3,3,3,4,4,4,5,5,5,6,6,6, &            ! Band 17
82            & 1,1,2,2,3,3,3,3,4,4,4,4,4,4,4,4, &            ! Band 18
83            & 1,1,2,2,3,3,3,3,4,4,4,4,4,4,4,4, &            ! Band 19
84            & 1,1,2,2,3,3,4,4,5,5,5,5,5,5,5,5, &            ! Band 20
85            & 1,1,2,2,3,3,4,4,5,5,5,5,5,5,5,5, &            ! Band 21
86            & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &            ! Band 22
87            & 1,1,1,1,2,2,3,3,4,4,5,5,5,5,5,5, &            ! Band 23
88            & 1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4, &            ! Band 24
89            & 1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,3, &            ! Band 25
90            & 1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,3, &            ! Band 26
91            & 1,1,2,2,3,3,4,4,4,4,4,4,4,4,4,4, &            ! Band 27
92            & 1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,3, &            ! Band 28
93            & 1,1,2,2,3,3,3,3,4,4,4,4,5,5,6,6 /)            ! Band 29
94
95IGN56(:) = (/ 4,4,8, &                                     ! Band 16
96            & 2,2,3,3,3,3, &                               ! Band 17
97            & 2,2,4,8, &                                   ! Band 18
98            & 2,2,4,8, &                                   ! Band 19
99            & 2,2,2,2,8, &                                 ! Band 20
100            & 2,2,2,2,8, &                                 ! Band 21
101            & 16, &                                        ! Band 22
102            & 4,2,2,2,6, &                                 ! Band 23
103            & 4,4,4,4, &                                   ! Band 24
104            & 2,4,10, &                                    ! Band 25
105            & 2,4,10, &                                    ! Band 26
106            & 2,2,2,10, &                                  ! Band 27
107            & 2,4,10, &                                    ! Band 28
108            & 2,2,4,4,2,2 /)                               ! Band 29
109
110IGB56(:) = (/ 16,16,16, &                                  ! Band 16
111            & 17,17,17,17,17,17, &                         ! Band 17
112            & 18,18,18,18, &                               ! Band 18
113            & 19,19,19,19, &                               ! Band 19
114            & 20,20,20,20,20, &                            ! Band 20
115            & 21,21,21,21,21, &                            ! Band 21
116            & 22, &                                        ! Band 22
117            & 23,23,23,23,23, &                            ! Band 23
118            & 24,24,24,24, &                               ! Band 24
119            & 25,25,25, &                                  ! Band 25
120            & 26,26,26, &                                  ! Band 26
121            & 27,27,27,27, &                               ! Band 27
122            & 28,28,28, &                                  ! Band 28
123            & 29,29,29,29,29,29 /)                         ! Band 29
124
125!-------------------------------------------------------------------------------
126!-- ECMWF high-resolution model RRTM_SW configuration with 112 g-points
127! Use this NGC, NGS, NGM, and NGN for reduced (112) g-point set
128! (A related code change is required in modules parsrtm.F90 and yoesrtwn.F90)
129
130IGC112(:) = (/ 6,12, 8, 8,10,10, 2,10, 8, 6, 6, 8, 6,12 /)
131IGS112(:) = (/ 6,18,26,34,44,54,56,66,74,80,86,94,100,112 /)
132
133!NGM(:)
134IGM112(:) = (/ 1,1,2,2,3,3,4,4,5,5,5,5,6,6,6,6, &           ! Band 16
135             & 1,2,3,4,5,6,6,7,8,8,9,10,10,11,12,12, &      ! Band 17
136             & 1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8, &           ! Band 18
137             & 1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8, &           ! Band 19
138             & 1,2,3,4,5,6,7,8,9,9,10,10,10,10,10,10, &     ! Band 20
139             & 1,2,3,4,5,6,7,8,9,9,10,10,10,10,10,10, &     ! Band 21
140             & 1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2, &           ! Band 22
141             & 1,1,2,2,3,4,5,6,7,8,9,9,10,10,10,10, &       ! Band 23
142             & 1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8, &           ! Band 24
143             & 1,2,3,3,4,4,5,5,5,5,6,6,6,6,6,6, &           ! Band 25
144             & 1,2,3,3,4,4,5,5,5,5,6,6,6,6,6,6, &           ! Band 26
145             & 1,2,3,4,5,6,7,7,7,7,8,8,8,8,8,8, &           ! Band 27
146             & 1,2,3,3,4,4,5,5,5,5,6,6,6,6,6,6, &           ! Band 28
147             & 1,2,3,4,5,5,6,6,7,7,8,8,9,10,11,12 /)        ! Band 29
148!NGN(:)
149IGN112(:) = (/ 2,2,2,2,4,4, &                               ! Band 16
150             & 1,1,1,1,1,2,1,2,1,2,1,2, &                   ! Band 17
151             & 1,1,1,1,2,2,4,4, &                           ! Band 18
152             & 1,1,1,1,2,2,4,4, &                           ! Band 19
153             & 1,1,1,1,1,1,1,1,2,6, &                       ! Band 20
154             & 1,1,1,1,1,1,1,1,2,6, &                       ! Band 21
155             & 8,8, &                                       ! Band 22
156             & 2,2,1,1,1,1,1,1,2,4, &                       ! Band 23
157             & 2,2,2,2,2,2,2,2, &                           ! Band 24
158             & 1,1,2,2,4,6, &                               ! Band 25
159             & 1,1,2,2,4,6, &                               ! Band 26
160             & 1,1,1,1,1,1,4,6, &                           ! Band 27
161             & 1,1,2,2,4,6, &                               ! Band 28
162             & 1,1,1,1,2,2,2,2,1,1,1,1 /)                   ! Band 29
163!NGBSW(:)
164IGB112(:) = (/ 16,16,16,16,16,16, &                         ! Band 16
165             & 17,17,17,17,17,17,17,17,17,17,17,17, &       ! Band 17
166             & 18,18,18,18,18,18,18,18, &                   ! Band 18
167             & 19,19,19,19,19,19,19,19, &                   ! Band 19
168             & 20,20,20,20,20,20,20,20,20,20, &             ! Band 20
169             & 21,21,21,21,21,21,21,21,21,21, &             ! Band 21
170             & 22,22, &                                     ! Band 22
171             & 23,23,23,23,23,23,23,23,23,23, &             ! Band 23
172             & 24,24,24,24,24,24,24,24, &                   ! Band 24
173             & 25,25,25,25,25,25, &                         ! Band 25
174             & 26,26,26,26,26,26, &                         ! Band 26
175             & 27,27,27,27,27,27,27,27, &                   ! Band 27
176             & 28,28,28,28,28,28, &                         ! Band 28
177             & 29,29,29,29,29,29,29,29,29,29,29,29 /)       ! Band 29
178
179!-------------------------------------------------------------------------------
180!-- original RRTM_SW configuration with 224 (14*16 g-points)
181! Use this NGC, NGS, NGM, and NGN for full (224) g-point set
182! (A related code change is required in modules parsrtm.F90 and yoesrtwn.F90)
183IGC224(:) = (/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16 /)
184IGS224(:) = (/ 16,32,48,64,80,96,112,128,144,160,176,192,208,224 /)
185
186IGM224(:) = (/ 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 16
187             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 17
188             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 18
189             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 19
190             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 20
191             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 21
192             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 22
193             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 23
194             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 24
195             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 25
196             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 26
197             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 27
198             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, &    ! Band 28
199             & 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16 /)    ! Band 29
200
201IGN224(:) = (/ 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 16
202             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 17
203             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 18
204             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 19
205             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 20
206             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 21
207             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 22
208             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 23
209             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 24
210             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 25
211             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 26
212             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 27
213             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, &           ! Band 28
214             & 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1 /)           ! Band 29
215
216IGB224(:) = (/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16,16,16, &   ! Band 16
217             & 17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17, &   ! Band 17
218             & 18,18,18,18,18,18,18,18,18,18,18,18,18,18,18,18, &   ! Band 18
219             & 19,19,19,19,19,19,19,19,19,19,19,19,19,19,19,19, &   ! Band 19
220             & 20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20, &   ! Band 20
221             & 21,21,21,21,21,21,21,21,21,21,21,21,21,21,21,21, &   ! Band 21
222             & 22,22,22,22,22,22,22,22,22,22,22,22,22,22,22,22, &   ! Band 22
223             & 23,23,23,23,23,23,23,23,23,23,23,23,23,23,23,23, &   ! Band 23
224             & 24,24,24,24,24,24,24,24,24,24,24,24,24,24,24,24, &   ! Band 24
225             & 25,25,25,25,25,25,25,25,25,25,25,25,25,25,25,25, &   ! Band 25
226             & 26,26,26,26,26,26,26,26,26,26,26,26,26,26,26,26, &   ! Band 26
227             & 27,27,27,27,27,27,27,27,27,27,27,27,27,27,27,27, &   ! Band 27
228             & 28,28,28,28,28,28,28,28,28,28,28,28,28,28,28,28, &   ! Band 28
229             & 29,29,29,29,29,29,29,29,29,29,29,29,29,29,29,28 /)   ! Band 29
230
231!=============================================================================
232
233WT(:) =  (/ 0.1527534276_JPRB, 0.1491729617_JPRB, 0.1420961469_JPRB, &
234          & 0.1316886544_JPRB, 0.1181945205_JPRB, 0.1019300893_JPRB, &
235          & 0.0832767040_JPRB, 0.0626720116_JPRB, 0.0424925000_JPRB, &
236          & 0.0046269894_JPRB, 0.0038279891_JPRB, 0.0030260086_JPRB, &
237          & 0.0022199750_JPRB, 0.0014140010_JPRB, 0.0005330000_JPRB, &
238          & 0.0000750000_JPRB /)
239
240!=============================================================================
241
242! These pressures are chosen such that the ln of the first pressure
243! has only a few non-zero digits (i.e. ln(PREF(1)) = 6.96000) and
244!  each subsequent ln(pressure) differs from the previous one by 0.2.
245PREF = (/ &
246 & 1.05363E+03_JPRB,8.62642E+02_JPRB,7.06272E+02_JPRB,5.78246E+02_JPRB,4.73428E+02_JPRB, &
247 & 3.87610E+02_JPRB,3.17348E+02_JPRB,2.59823E+02_JPRB,2.12725E+02_JPRB,1.74164E+02_JPRB, &
248 & 1.42594E+02_JPRB,1.16746E+02_JPRB,9.55835E+01_JPRB,7.82571E+01_JPRB,6.40715E+01_JPRB, &
249 & 5.24573E+01_JPRB,4.29484E+01_JPRB,3.51632E+01_JPRB,2.87892E+01_JPRB,2.35706E+01_JPRB, &
250 & 1.92980E+01_JPRB,1.57998E+01_JPRB,1.29358E+01_JPRB,1.05910E+01_JPRB,8.67114E+00_JPRB, &
251 & 7.09933E+00_JPRB,5.81244E+00_JPRB,4.75882E+00_JPRB,3.89619E+00_JPRB,3.18993E+00_JPRB, &
252 & 2.61170E+00_JPRB,2.13828E+00_JPRB,1.75067E+00_JPRB,1.43333E+00_JPRB,1.17351E+00_JPRB, &
253 & 9.60789E-01_JPRB,7.86628E-01_JPRB,6.44036E-01_JPRB,5.27292E-01_JPRB,4.31710E-01_JPRB, &
254 & 3.53455E-01_JPRB,2.89384E-01_JPRB,2.36928E-01_JPRB,1.93980E-01_JPRB,1.58817E-01_JPRB, &
255 & 1.30029E-01_JPRB,1.06458E-01_JPRB,8.71608E-02_JPRB,7.13612E-02_JPRB,5.84256E-02_JPRB, &
256 & 4.78349E-02_JPRB,3.91639E-02_JPRB,3.20647E-02_JPRB,2.62523E-02_JPRB,2.14936E-02_JPRB, &
257 & 1.75975E-02_JPRB,1.44076E-02_JPRB,1.17959E-02_JPRB,9.65769E-03_JPRB /) 
258PREFLOG = (/ &
259 & 6.9600E+00_JPRB, 6.7600E+00_JPRB, 6.5600E+00_JPRB, 6.3600E+00_JPRB, 6.1600E+00_JPRB, &
260 & 5.9600E+00_JPRB, 5.7600E+00_JPRB, 5.5600E+00_JPRB, 5.3600E+00_JPRB, 5.1600E+00_JPRB, &
261 & 4.9600E+00_JPRB, 4.7600E+00_JPRB, 4.5600E+00_JPRB, 4.3600E+00_JPRB, 4.1600E+00_JPRB, &
262 & 3.9600E+00_JPRB, 3.7600E+00_JPRB, 3.5600E+00_JPRB, 3.3600E+00_JPRB, 3.1600E+00_JPRB, &
263 & 2.9600E+00_JPRB, 2.7600E+00_JPRB, 2.5600E+00_JPRB, 2.3600E+00_JPRB, 2.1600E+00_JPRB, &
264 & 1.9600E+00_JPRB, 1.7600E+00_JPRB, 1.5600E+00_JPRB, 1.3600E+00_JPRB, 1.1600E+00_JPRB, &
265 & 9.6000E-01_JPRB, 7.6000E-01_JPRB, 5.6000E-01_JPRB, 3.6000E-01_JPRB, 1.6000E-01_JPRB, &
266 & -4.0000E-02_JPRB,-2.4000E-01_JPRB,-4.4000E-01_JPRB,-6.4000E-01_JPRB,-8.4000E-01_JPRB, &
267 & -1.0400E+00_JPRB,-1.2400E+00_JPRB,-1.4400E+00_JPRB,-1.6400E+00_JPRB,-1.8400E+00_JPRB, &
268 & -2.0400E+00_JPRB,-2.2400E+00_JPRB,-2.4400E+00_JPRB,-2.6400E+00_JPRB,-2.8400E+00_JPRB, &
269 & -3.0400E+00_JPRB,-3.2400E+00_JPRB,-3.4400E+00_JPRB,-3.6400E+00_JPRB,-3.8400E+00_JPRB, &
270 & -4.0400E+00_JPRB,-4.2400E+00_JPRB,-4.4400E+00_JPRB,-4.6400E+00_JPRB /) 
271! These are the temperatures associated with the respective
272! pressures for the MLS standard atmosphere.
273TREF = (/ &
274 & 2.9420E+02_JPRB, 2.8799E+02_JPRB, 2.7894E+02_JPRB, 2.6925E+02_JPRB, 2.5983E+02_JPRB, &
275 & 2.5017E+02_JPRB, 2.4077E+02_JPRB, 2.3179E+02_JPRB, 2.2306E+02_JPRB, 2.1578E+02_JPRB, &
276 & 2.1570E+02_JPRB, 2.1570E+02_JPRB, 2.1570E+02_JPRB, 2.1706E+02_JPRB, 2.1858E+02_JPRB, &
277 & 2.2018E+02_JPRB, 2.2174E+02_JPRB, 2.2328E+02_JPRB, 2.2479E+02_JPRB, 2.2655E+02_JPRB, &
278 & 2.2834E+02_JPRB, 2.3113E+02_JPRB, 2.3401E+02_JPRB, 2.3703E+02_JPRB, 2.4022E+02_JPRB, &
279 & 2.4371E+02_JPRB, 2.4726E+02_JPRB, 2.5085E+02_JPRB, 2.5457E+02_JPRB, 2.5832E+02_JPRB, &
280 & 2.6216E+02_JPRB, 2.6606E+02_JPRB, 2.6999E+02_JPRB, 2.7340E+02_JPRB, 2.7536E+02_JPRB, &
281 & 2.7568E+02_JPRB, 2.7372E+02_JPRB, 2.7163E+02_JPRB, 2.6955E+02_JPRB, 2.6593E+02_JPRB, &
282 & 2.6211E+02_JPRB, 2.5828E+02_JPRB, 2.5360E+02_JPRB, 2.4854E+02_JPRB, 2.4348E+02_JPRB, &
283 & 2.3809E+02_JPRB, 2.3206E+02_JPRB, 2.2603E+02_JPRB, 2.2000E+02_JPRB, 2.1435E+02_JPRB, &
284 & 2.0887E+02_JPRB, 2.0340E+02_JPRB, 1.9792E+02_JPRB, 1.9290E+02_JPRB, 1.8809E+02_JPRB, &
285 & 1.8329E+02_JPRB, 1.7849E+02_JPRB, 1.7394E+02_JPRB, 1.7212E+02_JPRB /) 
286!     -----------------------------------------------------------------
287
288IF (JPGPT == 56) THEN
289
290!- 14
291  NGC(:)=IGC56(:)
292  NGS(:)=IGS56(:)
293!- 14*16=224
294  NGM(:)=IGM56(:)
295
296  NGN(1:56)=IGN56(1:56)
297  NGBSW(1:56)=IGB56(1:56)
298
299ELSEIF (JPGPT == 112) THEN
300!- 14
301  NGC(:)=IGC112(:)
302  NGS(:)=IGS112(:)
303!- 14*16=224
304  NGM(:)=IGM112(:)
305
306  NGN(1:112)=IGN112(1:112)
307  NGBSW(1:112)=IGB112(1:112)
308
309ELSEIF (JPGPT == 224) THEN
310!- 14
311  NGC(:)=IGC224(:)
312  NGS(:)=IGS224(:)
313!- 14*16=224
314  NGM(:)=IGM224(:)
315
316  NGN(1:224)=IGN224(1:224)
317  NGBSW(1:224)=IGB224(1:224)
318
319ENDIF
320
321!     -----------------------------------------------------------------
322IF (LHOOK) CALL DR_HOOK('SUSRTM',1,ZHOOK_HANDLE)
323END SUBROUTINE SUSRTM
324
Note: See TracBrowser for help on using the repository browser.