Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

susrtm.F90 @ 5165

Last change on this file since 5165 was 4773, checked in by idelkadi, 11 months ago

Update of Ecrad in LMDZ The same organization of the Ecrad offline version is retained in order to facilitate the updating of Ecrad in LMDZ and the comparison between online and offline results. version 1.6.1 of Ecrad (https://github.com/lguez/ecrad.git)

Implementation of the double call of Ecrad in LMDZ

File size: 17.8 KB

Line
1	SUBROUTINE 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
10	USE PARKIND1 ,ONLY : JPRB , JPIM
11	USE YOMHOOK ,ONLY : LHOOK, DR_HOOK, JPHOOK
12
13	USE YOESRTM , ONLY : JPGPT, NGBSW, NGN
14	USE 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
22	IMPLICIT NONE
23
24	INTEGER(KIND=JPIM) :: IGC56(14), IGC112(14) , IGC224(14)
25	INTEGER(KIND=JPIM) :: IGS56(14), IGS112(14) , IGS224(14)
26
27	INTEGER(KIND=JPIM) :: IGM56(224),IGM112(224), IGM224(224)
28
29	INTEGER(KIND=JPIM) :: IGN56(56), IGN112(112), IGN224(224)
30	INTEGER(KIND=JPIM) :: IGB56(56), IGB112(112), IGB224(224)
31
32	REAL(KIND=JPHOOK) :: ZHOOK_HANDLE
33	!-----------------------------------------------------------------------
34	IF (LHOOK) CALL DR_HOOK('SUSRTM',0,ZHOOK_HANDLE)
35
36	NG(:) =(/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16 /)
37	NSPA(:) =(/ 9, 9, 9, 9, 1, 9, 9, 1, 9, 1, 0, 1, 9, 1 /)
38	NSPB(:) =(/ 1, 5, 1, 1, 1, 5, 1, 0, 1, 0, 0, 1, 5, 1 /)
39	NMPSRTM(:)=(/ 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
77	IGC56(:) = (/ 3, 6, 4, 4, 5, 5, 1, 5, 4, 3, 3, 4, 3, 6 /)
78	IGS56(:) = (/ 3, 9,13,17,22,27,28,33,37,40,43,47, 50, 56 /)
79
80	IGM56(:) = (/ 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
95	IGN56(:) = (/ 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
110	IGB56(:) = (/ 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
130	IGC112(:) = (/ 6,12, 8, 8,10,10, 2,10, 8, 6, 6, 8, 6,12 /)
131	IGS112(:) = (/ 6,18,26,34,44,54,56,66,74,80,86,94,100,112 /)
132
133	!NGM(:)
134	IGM112(:) = (/ 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(:)
149	IGN112(:) = (/ 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(:)
164	IGB112(:) = (/ 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)
183	IGC224(:) = (/ 16,16,16,16,16,16,16,16,16,16,16,16,16,16 /)
184	IGS224(:) = (/ 16,32,48,64,80,96,112,128,144,160,176,192,208,224 /)
185
186	IGM224(:) = (/ 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
201	IGN224(:) = (/ 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
216	IGB224(:) = (/ 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
233	WT(:) = (/ 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.
245	PREF = (/ &
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 /)
258	PREFLOG = (/ &
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.
273	TREF = (/ &
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
288	IF (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
299	ELSEIF (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
309	ELSEIF (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
319	ENDIF
320
321	! -----------------------------------------------------------------
322	IF (LHOOK) CALL DR_HOOK('SUSRTM',1,ZHOOK_HANDLE)
323	END SUBROUTINE SUSRTM
324

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format