Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

susrtm.F90 @ 5407

Last change on this file since 5407 was 4773, checked in by idelkadi, 12 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

Rev	Line
[4773]	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