Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

suolwn.F90 @ 4031

Last change on this file since 4031 was 3990, checked in by millour, 8 years ago
An updated LMDZ5 (based on LMDZ rev 2786) to play with. EM
File size: 39.3 KB

Rev	Line
[3990]	1	SUBROUTINE SUOLW
	2	!
	3	!**** SULW - INITIALIZE COMMON YOELW
	4	!
	5	! PURPOSE.
	6	! --------
	7	! INITIALIZE YOELW, THE COMMON THAT CONTAINS COEFFICIENTS
	8	! NEEDED TO RUN THE ECMWF LONGWAVE RADIATION ROUTINES
	9	!
	10	!** INTERFACE.
	11	! ----------
	12	! CALL *SULW
	13	!
	14	! EXPLICIT ARGUMENTS :
	15	! --------------------
	16	! NONE
	17	!
	18	! IMPLICIT ARGUMENTS :
	19	! --------------------
	20	! COMMON YOELW
	21	!
	22	! METHOD.
	23	! -------
	24	! SEE DOCUMENTATION
	25	!
	26	! EXTERNALS.
	27	! ----------
	28	!
	29	! REFERENCE.
	30	! ----------
	31	! ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
	32	!
	33	! AUTHOR.
	34	! -------
	35	! JEAN-JACQUES MORCRETTE ECMWF
	36	!
	37	! MODIFICATIONS.
	38	! --------------
	39	! ORIGINAL : 88-12-15
	40	! J.-J. Morcrette ECMWF 951215: Pseudo-Voigt Line Profiles
	41	! ------------------------------------------------------------------
	42	! IMPLICIT LOGICAL (L)
	43	!
	44	!#include "yomlun.h"
	45	!#include "yoerad.h"
	46	!#include "yoelw.h"
	47
	48	#include "tsmbkind.h"
	49
	50	USE YOEOLW , ONLY : MXIXT ,NISP ,NIPD ,NIPD2 ,&
	51	&NTR ,NTRA ,NUA ,NG1 ,NG1P1 ,&
	52	&AT ,BT ,GA ,GB ,&
	53	&O1H ,O2H ,OCT ,RPIALF0 ,&
	54	&RT1 ,TINTP ,TREF ,TSTAND ,TSTP ,&
	55	&WG1 ,XP ,RVGCO2 ,RVGH2O ,RVGO3
	56
	57
	58	IMPLICIT NONE
	59
	60
	61	! ----------------------------------------------------------------
	62	!
	63	!* 1. SET VALUES.
	64	! -----------
	65	!
	66	NG1=2
	67	NG1P1=NG1+1
	68	NISP=6
	69	NIPD=8
	70	NIPD2=2*NIPD
	71	NTR=11
	72	NTRA=15
	73	NUA=24
	74	!
	75	!
	76	!-- COEFFICIENTS FOR THE PSEUDO-VOIGT LINE PROFILE (PA)
	77	!
	78	RVGCO2= 60.
	79	RVGH2O= 30.
	80	RVGO3 =400.
	81	!
	82	!
	83	!-- REFERENCE TEMPERATURES FOR THE PLANCK FUNCTIONS ----- TSTAND
	84	! TEMPERATURE DEPENDENCE OF THE ABSORPTION ----- TREF
	85	TSTAND= 250.00
	86	TREF = 250.00
	87	!
	88
	89	!-- ROOTS AND WEIGHTS FOR THE 2-POINT GAUSSIAN QUADRATURE
	90	RT1 = (/ -0.577350269, +0.577350269 /)
	91	WG1 = (/ 1. , 1. /)
	92
	93	!-- TEMPERATURE STEP AND MAXIMUM INDEX FOR THE TRANSMISSIVITIES
	94	TSTP = 12.5
	95	MXIXT = 10
	96
	97	!
	98	!-- REFERENCE TEMPERATURE FOR THE TRANSMISSIVITIES
	99
	100	TINTP = (/ 187.5, 200., 212.5, 225., 237.5, &
	101	&250.,&
	102	&262.5, 275., 287.5, 300., 312.5 /)
	103	!
	104	!-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE PLANCK FUNCTIONS
	105	!
	106	XP( :, 1) = (/&
	107	&0.46430621E+02, 0.12928299E+03, 0.20732648E+03,&
	108	&0.31398411E+03, 0.18373177E+03,-0.11412303E+03/)
	109	XP( :, 2) = (/&
	110	&0.73604774E+02, 0.27887914E+03, 0.27076947E+03,&
	111	&-.57322111E+02,-0.64742459E+02, 0.87238280E+02/)
	112	XP( :, 3) = (/&
	113	&0.37050866E+02, 0.20498759E+03, 0.37558029E+03,&
	114	&0.17401171E+03,-0.13350302E+03,-0.37651795E+02/)
	115	XP( :, 4) = (/&
	116	&0.14930141E+02, 0.89161160E+02, 0.17793062E+03,&
	117	&0.93433860E+02,-0.70646020E+02,-0.26373150E+02/)
	118	XP( :, 5) = (/&
	119	&0.40386780E+02, 0.10855270E+03, 0.50755010E+02,&
	120	&-0.31496190E+02, 0.12791300E+00, 0.18017770E+01/)
	121	XP( :, 6) = (/&
	122	&0.90811926E+01, 0.75073923E+02, 0.24654438E+03,&
	123	&0.39332612E+03, 0.29385281E+03, 0.89107921E+02/)
	124
	125	!
	126	!-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE TEMPERATURE DEPENDENCE OF
	127	! THE ABSORPTION
	128	!
	129	!- from Zhong & Haigh
	130	!
	131	!-- H2O -- 0 - 350
	132	AT(1,:)= (/ 0.298199E-02,-.394023E-03,0.319566E-04 /)
	133	BT(1,:)= (/-0.106432E-04,0.660324E-06,0.174356E-06 /)
	134	!
	135	!-- H2O -- 500 - 800
	136	AT(2,:)= (/ 0.143676E-01,0.366501E-02,-.160822E-02 /)
	137	BT(2,:)= (/-0.553979E-04,-.101701E-04,0.920868E-05 /)
	138	!
	139	!-- H2O -- 800 - 970 + 1110 - 1250
	140	AT(3,:)= (/ 0.197861E-01,0.315541E-02,-.174547E-02 /)
	141	BT(3,:)= (/-0.877012E-04,0.513302E-04,0.523138E-06 /)
	142	!
	143	!-- H2O -- 970 - 1110
	144	AT(4,:)= (/ 0.289560E-01,-.208807E-02,-.121943E-02 /)
	145	BT(4,:)= (/-0.165960E-03,0.157704E-03,-.146427E-04 /)
	146	!
	147	!-- H2O -- 350 - 500
	148	AT(5,:)= (/ 0.103800E-01,0.436296E-02,-.161431E-02 /)
	149	BT(5,:)= (/ -.276744E-04,-.327381E-04,0.127646E-04 /)
	150	!
	151	!-- H2O -- 1250 - 1450 + 1880 - 2820
	152	AT(6,:)= (/ 0.868859E-02,-.972752E-03,0.000000E-00 /)
	153	BT(6,:)= (/ -.278412E-04,-.713940E-06,0.117469E-05 /)
	154	!
	155
	156	!-------CO2
	157
	158	AT(7,:)= (/&
	159	& 0.250073E-03,0.455875E-03,0.109242E-03 /)
	160	BT(7,:)= (/&
	161	& 0.199846E-05,-.216313E-05,0.175991E-06 /)
	162	AT(8,:)= (/&
	163	& 0.307423E-01,0.110879E-02,-.322172E-03 /)
	164	BT(8,:)= (/&
	165	&-0.108482E-03,0.258096E-05,-.814575E-06 /)
	166	!
	167	!-- COEFFICIENTS FOR OZONE ABSORPTION
	168	!
	169	O1H = 2230.
	170	O2H = 100.
	171	RPIALF0= 2.
	172	OCT(1)= -.326E-03
	173	OCT(2)= -.102E-05
	174	OCT(3)= .137E-02
	175	OCT(4)= -.535E-05
	176	!
	177	!-- COEFFICIENTS FOR H2O E-TYPE CONTINUUM ABSORPTION
	178	!
	179	! DATA RNTNU / 70.0, 26.77, 8.45, 5.78 /
	180	!
	181	!-----------------------------------------------------------------------
	182	! PADE APPROXIMANTS FOR THE TRANSMISSIVITIES
	183	! ==========================================
	184	!-----------------------------------------------------------------------
	185	!
	186	!
	187	!
	188	!
	189	!-----------------------------------------------------------------------
	190	! PADE APPROXIMANTS FOR THE TRANSMISSIVITIES
	191	!
	192	!- H2O 0-350 cm-1
	193	!
	194	!-- R.D. -- G = - 0.2 SLA
	195	!
	196	!
	197	!----- INTERVAL = 1 ----- T = 187.5
	198	!
	199	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	200	GA( 1, 1, :) = (/&
	201	& 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/)
	202	GB( 1, 1, :) = (/&
	203	& 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/)
	204	GA( 1, 2, :) = (/&
	205	& 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/)
	206	GB( 1, 2, :) = (/&
	207	& 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/)
	208	!
	209	!----- INTERVAL = 1 ----- T = 200.0
	210	!
	211	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	212	GA( 2, 1, :) = (/&
	213	& 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/)
	214	GB( 2, 1, :) = (/&
	215	& 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/)
	216	GA( 2, 2, :) = (/&
	217	& 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/)
	218	GB( 2, 2, :) = (/&
	219	& 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/)
	220	!
	221	!----- INTERVAL = 1 ----- T = 212.5
	222	!
	223	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	224	GA( 3, 1, :) = (/&
	225	& 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/)
	226	GB( 3, 1, :) = (/&
	227	& 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/)
	228	GA( 3, 2, :) = (/&
	229	& 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/)
	230	GB( 3, 2, :) = (/&
	231	& 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/)
	232	!
	233	!----- INTERVAL = 1 ----- T = 225.0
	234	!
	235	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	236	GA( 4, 1, :) = (/&
	237	& 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/)
	238	GB( 4, 1, :) = (/&
	239	& 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/)
	240	GA( 4, 2, :) = (/&
	241	& 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/)
	242	GB( 4, 2, :) = (/&
	243	& 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/)
	244	!
	245	!----- INTERVAL = 1 ----- T = 237.5
	246	!
	247	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	248	GA( 5, 1, :) = (/&
	249	& 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/)
	250	GB( 5, 1, :) = (/&
	251	& 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/)
	252	GA( 5, 2, :) = (/&
	253	& 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/)
	254	GB( 5, 2, :) = (/&
	255	& 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/)
	256	!
	257	!----- INTERVAL = 1 ----- T = 250.0
	258	!
	259	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	260	GA( 6, 1, :) = (/&
	261	& 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/)
	262	GB( 6, 1, :) = (/&
	263	& 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/)
	264	GA( 6, 2, :) = (/&
	265	& 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/)
	266	GB( 6, 2, :) = (/&
	267	& 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/)
	268	!
	269	!----- INTERVAL = 1 ----- T = 262.5
	270	!
	271	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	272	GA( 7, 1, :) = (/&
	273	& 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/)
	274	GB( 7, 1, :) = (/&
	275	& 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/)
	276	GA( 7, 2, :) = (/&
	277	& 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/)
	278	GB( 7, 2, :) = (/&
	279	& 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/)
	280	!
	281	!----- INTERVAL = 1 ----- T = 275.0
	282	!
	283	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	284	GA( 8, 1, :) = (/&
	285	& 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/)
	286	GB( 8, 1, :) = (/&
	287	& 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/)
	288	GA( 8, 2, :) = (/&
	289	& 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/)
	290	GB( 8, 2, :) = (/&
	291	& 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/)
	292	!
	293	!----- INTERVAL = 1 ----- T = 287.5
	294	!
	295	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	296	GA( 9, 1, :) = (/&
	297	& 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/)
	298	GB( 9, 1, :) = (/&
	299	& 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/)
	300	GA( 9, 2, :) = (/&
	301	& 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/)
	302	GB( 9, 2, :) = (/&
	303	& 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/)
	304	!
	305	!----- INTERVAL = 1 ----- T = 300.0
	306	!
	307	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	308	GA(10, 1, :) = (/&
	309	& 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/)
	310	GB(10, 1, :) = (/&
	311	& 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/)
	312	GA(10, 2, :) = (/&
	313	& 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/)
	314	GB(10, 2, :) = (/&
	315	& 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/)
	316	!
	317	!----- INTERVAL = 1 ----- T = 312.5
	318	!
	319	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	320	GA(11, 1, :) = (/&
	321	& 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/)
	322	GB(11, 1, :) = (/&
	323	& 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/)
	324	GA(11, 2, :) = (/&
	325	& 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/)
	326	GB(11, 2, :) = (/&
	327	& 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/)
	328	!
	329	!!
	330	!- H2O 500-800 cm-1
	331	!
	332	!--- R.D. --- G = 0.02 + 0.50 / ( 1 + 4.5 U )
	333	!
	334	!
	335	!----- INTERVAL = 2 ----- T = 187.5
	336	!
	337	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	338	GA( 1, 3, :) = (/&
	339	& 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/)
	340	GB( 1, 3, :) = (/&
	341	& 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/)
	342	GA( 1, 4, :) = (/&
	343	& 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/)
	344	GB( 1, 4, :) = (/&
	345	& 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/)
	346	!
	347	!----- INTERVAL = 2 ----- T = 200.0
	348	!
	349	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	350	GA( 2, 3, :) = (/&
	351	& 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/)
	352	GB( 2, 3, :) = (/&
	353	& 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/)
	354	GA( 2, 4, :) = (/&
	355	& 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/)
	356	GB( 2, 4, :) = (/&
	357	& 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/)
	358	!
	359	!----- INTERVAL = 2 ----- T = 212.5
	360	!
	361	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	362	GA( 3, 3, :) = (/&
	363	& 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/)
	364	GB( 3, 3, :) = (/&
	365	& 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/)
	366	GA( 3, 4, :) = (/&
	367	& 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/)
	368	GB( 3, 4, :) = (/&
	369	& 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/)
	370	!
	371	!----- INTERVAL = 2 ----- T = 225.0
	372	!
	373	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	374	GA( 4, 3, :) = (/&
	375	& 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/)
	376	GB( 4, 3, :) = (/&
	377	& 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/)
	378	GA( 4, 4, :) = (/&
	379	& 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/)
	380	GB( 4, 4, :) = (/&
	381	& 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/)
	382	!
	383	!----- INTERVAL = 2 ----- T = 237.5
	384	!
	385	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	386	GA( 5, 3, :) = (/&
	387	& 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/)
	388	GB( 5, 3, :) = (/&
	389	& 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/)
	390	GA( 5, 4, :) = (/&
	391	& 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/)
	392	GB( 5, 4, :) = (/&
	393	& 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/)
	394	!
	395	!----- INTERVAL = 2 ----- T = 250.0
	396	!
	397	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	398	GA( 6, 3, :) = (/&
	399	& 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/)
	400	GB( 6, 3, :) = (/&
	401	& 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/)
	402	GA( 6, 4, :) = (/&
	403	& 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/)
	404	GB( 6, 4, :) = (/&
	405	& 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/)
	406	!
	407	!----- INTERVAL = 2 ----- T = 262.5
	408	!
	409	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	410	GA( 7, 3, :) = (/&
	411	& 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/)
	412	GB( 7, 3, :) = (/&
	413	& 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/)
	414	GA( 7, 4, :) = (/&
	415	& 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/)
	416	GB( 7, 4, :) = (/&
	417	& 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/)
	418	!
	419	!----- INTERVAL = 2 ----- T = 275.0
	420	!
	421	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	422	GA( 8, 3, :) = (/&
	423	& 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/)
	424	GB( 8, 3, :) = (/&
	425	& 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/)
	426	GA( 8, 4, :) = (/&
	427	& 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/)
	428	GB( 8, 4, :) = (/&
	429	& 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/)
	430	!
	431	!----- INTERVAL = 2 ----- T = 287.5
	432	!
	433	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	434	GA( 9, 3, :) = (/&
	435	& 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/)
	436	GB( 9, 3, :) = (/&
	437	& 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/)
	438	GA( 9, 4, :) = (/&
	439	& 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/)
	440	GB( 9, 4, :) = (/&
	441	& 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/)
	442	!
	443	!----- INTERVAL = 2 ----- T = 300.0
	444	!
	445	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	446	GA(10, 3, :) = (/&
	447	& 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/)
	448	GB(10, 3, :) = (/&
	449	& 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/)
	450	GA(10, 4, :) = (/&
	451	& 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/)
	452	GB(10, 4, :) = (/&
	453	& 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/)
	454	!
	455	!----- INTERVAL = 2 ----- T = 312.5
	456	!
	457	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	458	GA(11, 3, :) = (/&
	459	& 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/)
	460	GB(11, 3, :) = (/&
	461	& 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/)
	462	GA(11, 4, :) = (/&
	463	& 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/)
	464	GB(11, 4, :) = (/&
	465	& 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/)
	466	!
	467	!
	468	!
	469	!
	470	!- H2O 350-500 cm-1
	471	!
	472	!-- G = - 0.2*SLA, 0.0 +0.5/(1+0.5U)
	473	!
	474	!----- INTERVAL = 5 ----- T = 187.5
	475	!
	476	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	477	GA( 1, 5, :) = (/&
	478	& 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/)
	479	GB( 1, 5, :) = (/&
	480	& 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/)
	481	GA( 1, 6, :) = (/&
	482	& 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/)
	483	GB( 1, 6, :) = (/&
	484	& 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/)
	485	!
	486	!----- INTERVAL = 5 ----- T = 200.0
	487	!
	488	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	489	GA( 2, 5, :) = (/&
	490	& 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/)
	491	GB( 2, 5, :) = (/&
	492	& 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/)
	493	GA( 2, 6, :) = (/&
	494	& 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/)
	495	GB( 2, 6, :) = (/&
	496	& 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/)
	497	!
	498	!----- INTERVAL = 5 ----- T = 212.5
	499	!
	500	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	501	GA( 3, 5, :) = (/&
	502	& 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/)
	503	GB( 3, 5, :) = (/&
	504	& 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/)
	505	GA( 3, 6, :) = (/&
	506	& 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/)
	507	GB( 3, 6, :) = (/&
	508	& 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/)
	509	!
	510	!----- INTERVAL = 5 ----- T = 225.0
	511	!
	512	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	513	GA( 4, 5, :) = (/&
	514	& 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/)
	515	GB( 4, 5, :) = (/&
	516	& 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/)
	517	GA( 4, 6, :) = (/&
	518	& 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/)
	519	GB( 4, 6, :) = (/&
	520	& 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/)
	521	!
	522	!----- INTERVAL = 5 ----- T = 237.5
	523	!
	524	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	525	GA( 5, 5, :) = (/&
	526	& 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/)
	527	GB( 5, 5, :) = (/&
	528	& 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/)
	529	GA( 5, 6, :) = (/&
	530	& 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/)
	531	GB( 5, 6, :) = (/&
	532	& 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/)
	533	!
	534	!----- INTERVAL = 5 ----- T = 250.0
	535	!
	536	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	537	GA( 6, 5, :) = (/&
	538	& 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/)
	539	GB( 6, 5, :) = (/&
	540	& 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/)
	541	GA( 6, 6, :) = (/&
	542	& 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/)
	543	GB( 6, 6, :) = (/&
	544	& 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/)
	545	!
	546	!----- INTERVAL = 5 ----- T = 262.5
	547	!
	548	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	549	GA( 7, 5, :) = (/&
	550	& 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/)
	551	GB( 7, 5, :) = (/&
	552	& 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/)
	553	GA( 7, 6, :) = (/&
	554	& 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/)
	555	GB( 7, 6, :) = (/&
	556	& 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/)
	557	!
	558	!----- INTERVAL = 5 ----- T = 275.0
	559	!
	560	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	561	GA( 8, 5, :) = (/&
	562	& 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/)
	563	GB( 8, 5, :) = (/&
	564	& 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/)
	565	GA( 8, 6, :) = (/&
	566	& 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/)
	567	GB( 8, 6, :) = (/&
	568	& 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/)
	569	!
	570	!----- INTERVAL = 5 ----- T = 287.5
	571	!
	572	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	573	GA( 9, 5, :) = (/&
	574	& 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/)
	575	GB( 9, 5, :) = (/&
	576	& 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/)
	577	GA( 9, 6, :) = (/&
	578	& 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/)
	579	GB( 9, 6, :) = (/&
	580	& 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/)
	581	!
	582	!----- INTERVAL = 5 ----- T = 300.0
	583	!
	584	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	585	GA(10, 5, :) = (/&
	586	& 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/)
	587	GB(10, 5, :) = (/&
	588	& 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/)
	589	GA(10, 6, :) = (/&
	590	& 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/)
	591	GB(10, 6, :) = (/&
	592	& 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/)
	593	!
	594	!----- INTERVAL = 5 ----- T = 312.5
	595	!
	596	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	597	GA(11, 5, :) = (/&
	598	& 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/)
	599	GB(11, 5, :) = (/&
	600	& 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/)
	601	GA(11, 6, :) = (/&
	602	& 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/)
	603	GB(11, 6, :) = (/&
	604	& 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/)
	605
	606
	607
	608	!
	609	!
	610	!- H2O 800-970 + 1110-1250 cm-1
	611	!
	612	!--- G = 3.875E-03 ---------------
	613	!
	614	!----- INTERVAL = 3 ----- T = 187.5
	615	!
	616	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	617	GA( 1, 7, :) = (/&
	618	& 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/)
	619	GB( 1, 7, :) = (/&
	620	& 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/)
	621	GA( 1, 8, :) = (/&
	622	& 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/)
	623	GB( 1, 8, :) = (/&
	624	& 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/)
	625	!
	626	!----- INTERVAL = 3 ----- T = 200.0
	627	!
	628	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	629	GA( 2, 7, :) = (/&
	630	& 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/)
	631	GB( 2, 7, :) = (/&
	632	& 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/)
	633	GA( 2, 8, :) = (/&
	634	& 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/)
	635	GB( 2, 8, :) = (/&
	636	& 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/)
	637	!
	638	!----- INTERVAL = 3 ----- T = 212.5
	639	!
	640	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	641	GA( 3, 7, :) = (/&
	642	& 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/)
	643	GB( 3, 7, :) = (/&
	644	& 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/)
	645	GA( 3, 8, :) = (/&
	646	& 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/)
	647	GB( 3, 8, :) = (/&
	648	& 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/)
	649	!
	650	!----- INTERVAL = 3 ----- T = 225.0
	651	!
	652	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	653	GA( 4, 7, :) = (/&
	654	& 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/)
	655	GB( 4, 7, :) = (/&
	656	& 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/)
	657	GA( 4, 8, :) = (/&
	658	& 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/)
	659	GB( 4, 8, :) = (/&
	660	& 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/)
	661	!
	662	!----- INTERVAL = 3 ----- T = 237.5
	663	!
	664	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	665	GA( 5, 7, :) = (/&
	666	& 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/)
	667	GB( 5, 7, :) = (/&
	668	& 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/)
	669	GA( 5, 8, :) = (/&
	670	& 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/)
	671	GB( 5, 8, :) = (/&
	672	& 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/)
	673	!
	674	!----- INTERVAL = 3 ----- T = 250.0
	675	!
	676	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	677	GA( 6, 7, :) = (/&
	678	& 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/)
	679	GB( 6, 7, :) = (/&
	680	& 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/)
	681	GA( 6, 8, :) = (/&
	682	& 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/)
	683	GB( 6, 8, :) = (/&
	684	& 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/)
	685	!
	686	!----- INTERVAL = 3 ----- T = 262.5
	687	!
	688	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	689	GA( 7, 7, :) = (/&
	690	& 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/)
	691	GB( 7, 7, :) = (/&
	692	& 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/)
	693	GA( 7, 8, :) = (/&
	694	& 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/)
	695	GB( 7, 8, :) = (/&
	696	& 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/)
	697	!
	698	!----- INTERVAL = 3 ----- T = 275.0
	699	!
	700	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	701	GA( 8, 7, :) = (/&
	702	& 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/)
	703	GB( 8, 7, :) = (/&
	704	& 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/)
	705	GA( 8, 8, :) = (/&
	706	& 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/)
	707	GB( 8, 8, :) = (/&
	708	& 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/)
	709	!
	710	!----- INTERVAL = 3 ----- T = 287.5
	711	!
	712	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	713	GA( 9, 7, :) = (/&
	714	& 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/)
	715	GB( 9, 7, :) = (/&
	716	& 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/)
	717	GA( 9, 8, :) = (/&
	718	& 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/)
	719	GB( 9, 8, :) = (/&
	720	& 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/)
	721	!
	722	!----- INTERVAL = 3 ----- T = 300.0
	723	!
	724	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	725	GA(10, 7, :) = (/&
	726	& 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/)
	727	GB(10, 7, :) = (/&
	728	& 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/)
	729	GA(10, 8, :) = (/&
	730	& 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/)
	731	GB(10, 8, :) = (/&
	732	& 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/)
	733	!
	734	!----- INTERVAL = 3 ----- T = 312.5
	735	!
	736	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	737	GA(11, 7, :) = (/&
	738	& 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/)
	739	GB(11, 7, :) = (/&
	740	& 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/)
	741	GA(11, 8, :) = (/&
	742	& 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/)
	743	GB(11, 8, :) = (/&
	744	& 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/)
	745	!
	746	!
	747	!
	748	!- H2O 970 - 1110 cm-1
	749	!
	750	!!-- G = 3.6E-03
	751	!
	752	!----- INTERVAL = 4 ----- T = 187.5
	753	!
	754	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	755	GA( 1, 9, :) = (/&
	756	& 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/)
	757	GB( 1, 9, :) = (/&
	758	& 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/)
	759	GA( 1,10, :) = (/&
	760	& 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/)
	761	GB( 1,10, :) = (/&
	762	& 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/)
	763	!
	764	!----- INTERVAL = 4 ----- T = 200.0
	765	!
	766	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	767	GA( 2, 9, :) = (/&
	768	& 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/)
	769	GB( 2, 9, :) = (/&
	770	& 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/)
	771	GA( 2,10, :) = (/&
	772	& 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/)
	773	GB( 2,10, :) = (/&
	774	& 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/)
	775	!
	776	!----- INTERVAL = 4 ----- T = 212.5
	777	!
	778	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	779	GA( 3, 9, :) = (/&
	780	& 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/)
	781	GB( 3, 9, :) = (/&
	782	& 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/)
	783	GA( 3,10, :) = (/&
	784	& 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/)
	785	GB( 3,10, :) = (/&
	786	& 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/)
	787	!
	788	!----- INTERVAL = 4 ----- T = 225.0
	789	!
	790	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	791	GA( 4, 9, :) = (/&
	792	& 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/)
	793	GB( 4, 9, :) = (/&
	794	& 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/)
	795	GA( 4,10, :) = (/&
	796	& 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/)
	797	GB( 4,10, :) = (/&
	798	& 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/)
	799	!
	800	!----- INTERVAL = 4 ----- T = 237.5
	801	!
	802	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	803	GA( 5, 9, :) = (/&
	804	& 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/)
	805	GB( 5, 9, :) = (/&
	806	& 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/)
	807	GA( 5,10, :) = (/&
	808	& 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/)
	809	GB( 5,10, :) = (/&
	810	& 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/)
	811	!
	812	!----- INTERVAL = 4 ----- T = 250.0
	813	!
	814	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	815	GA( 6, 9, :) = (/&
	816	& 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/)
	817	GB( 6, 9, :) = (/&
	818	& 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/)
	819	GA( 6,10, :) = (/&
	820	& 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/)
	821	GB( 6,10, :) = (/&
	822	& 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/)
	823	!
	824	!----- INTERVAL = 4 ----- T = 262.5
	825	!
	826	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	827	GA( 7, 9, :) = (/&
	828	& 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/)
	829	GB( 7, 9, :) = (/&
	830	& 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/)
	831	GA( 7,10, :) = (/&
	832	& 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/)
	833	GB( 7,10, :) = (/&
	834	& 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/)
	835	!
	836	!----- INTERVAL = 4 ----- T = 275.0
	837	!
	838	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	839	GA( 8, 9, :) = (/&
	840	& 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/)
	841	GB( 8, 9, :) = (/&
	842	& 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/)
	843	GA( 8,10, :) = (/&
	844	& 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/)
	845	GB( 8,10, :) = (/&
	846	& 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/)
	847	!
	848	!----- INTERVAL = 4 ----- T = 287.5
	849	!
	850	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	851	GA( 9, 9, :) = (/&
	852	& 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/)
	853	GB( 9, 9, :) = (/&
	854	& 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/)
	855	GA( 9,10, :) = (/&
	856	& 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/)
	857	GB( 9,10, :) = (/&
	858	& 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/)
	859	!
	860	!----- INTERVAL = 4 ----- T = 300.0
	861	!
	862	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	863	GA(10, 9, :) = (/&
	864	& 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/)
	865	GB(10, 9, :) = (/&
	866	& 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/)
	867	GA(10,10, :) = (/&
	868	& 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/)
	869	GB(10,10, :) = (/&
	870	& 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/)
	871	!
	872	!----- INTERVAL = 4 ----- T = 312.5
	873	!
	874	!-- INDICES FOR PADE APPROXIMATION 1 28 37 45
	875	GA(11, 9, :) = (/&
	876	& 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/)
	877	GB(11, 9, :) = (/&
	878	& 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/)
	879	GA(11,10, :) = (/&
	880	& 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/)
	881	GB(11,10, :) = (/&
	882	& 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/)
	883	!
	884	!
	885	!
	886	!- H2O 1250-1450 + 1880-2820 cm-1
	887	!
	888	!--- G = 0.0
	889	!
	890	!
	891	!----- INTERVAL = 6 ----- T = 187.5
	892	!
	893	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	894	GA( 1,11, :) = (/&
	895	& 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/)
	896	GB( 1,11, :) = (/&
	897	& 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/)
	898	GA( 1,12, :) = (/&
	899	& 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/)
	900	GB( 1,12, :) = (/&
	901	& 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/)
	902	!
	903	!----- INTERVAL = 6 ----- T = 200.0
	904	!
	905	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	906	GA( 2,11, :) = (/&
	907	& 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/)
	908	GB( 2,11, :) = (/&
	909	& 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/)
	910	GA( 2,12, :) = (/&
	911	& 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/)
	912	GB( 2,12, :) = (/&
	913	& 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/)
	914	!
	915	!----- INTERVAL = 6 ----- T = 212.5
	916	!
	917	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	918	GA( 3,11, :) = (/&
	919	& 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/)
	920	GB( 3,11, :) = (/&
	921	& 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/)
	922	GA( 3,12, :) = (/&
	923	& 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/)
	924	GB( 3,12, :) = (/&
	925	& 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/)
	926	!
	927	!----- INTERVAL = 6 ----- T = 225.0
	928	!
	929	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	930	GA( 4,11, :) = (/&
	931	& 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/)
	932	GB( 4,11, :) = (/&
	933	& 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/)
	934	GA( 4,12, :) = (/&
	935	& 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/)
	936	GB( 4,12, :) = (/&
	937	& 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/)
	938	!
	939	!----- INTERVAL = 6 ----- T = 237.5
	940	!
	941	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	942	GA( 5,11, :) = (/&
	943	& 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/)
	944	GB( 5,11, :) = (/&
	945	& 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/)
	946	GA( 5,12, :) = (/&
	947	& 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/)
	948	GB( 5,12, :) = (/&
	949	& 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/)
	950	!
	951	!----- INTERVAL = 6 ----- T = 250.0
	952	!
	953	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	954	GA( 6,11, :) = (/&
	955	& 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/)
	956	GB( 6,11, :) = (/&
	957	& 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/)
	958	GA( 6,12, :) = (/&
	959	&-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/)
	960	GB( 6,12, :) = (/&
	961	&-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/)
	962	!
	963	!----- INTERVAL = 6 ----- T = 262.5
	964	!
	965	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	966	GA( 7,11, :) = (/&
	967	&-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/)
	968	GB( 7,11, :) = (/&
	969	&-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/)
	970	GA( 7,12, :) = (/&
	971	& 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/)
	972	GB( 7,12, :) = (/&
	973	& 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/)
	974	!
	975	!----- INTERVAL = 6 ----- T = 275.0
	976	!
	977	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	978	GA( 8,11, :) = (/&
	979	&-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/)
	980	GB( 8,11, :) = (/&
	981	&-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/)
	982	GA( 8,12, :) = (/&
	983	& 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/)
	984	GB( 8,12, :) = (/&
	985	& 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/)
	986	!
	987	!----- INTERVAL = 6 ----- T = 287.5
	988	!
	989	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	990	GA( 9,11, :) = (/&
	991	& 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/)
	992	GB( 9,11, :) = (/&
	993	& 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/)
	994	GA( 9,12, :) = (/&
	995	& 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/)
	996	GB( 9,12, :) = (/&
	997	& 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/)
	998	!
	999	!----- INTERVAL = 6 ----- T = 300.0
	1000	!
	1001	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	1002	GA(10,11, :) = (/&
	1003	& 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/)
	1004	GB(10,11, :) = (/&
	1005	& 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/)
	1006	GA(10,12, :) = (/&
	1007	& 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/)
	1008	GB(10,12, :) = (/&
	1009	& 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/)
	1010	!
	1011	!----- INTERVAL = 6 ----- T = 312.5
	1012	!
	1013	!-- INDICES FOR PADE APPROXIMATION 1 35 40 45
	1014	GA(11,11, :) = (/&
	1015	& 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/)
	1016	GB(11,11, :) = (/&
	1017	& 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/)
	1018	GA(11,12, :) = (/&
	1019	& 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/)
	1020	GB(11,12, :) = (/&
	1021	& 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/)
	1022	!
	1023	!
	1024
	1025	!
	1026	!
	1027	!
	1028	!-- CO2 -- INT.2 -- 500-800 CM-1 --- FROM ABS225 ----------------------
	1029	!
	1030	!
	1031	!
	1032	!-- FIU = 0.8 + MAX(0.35,(7-IU)*0.9) , X/T, 9
	1033	!
	1034	!----- INTERVAL = 2 ----- T = 187.5
	1035	!
	1036	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1037	GA( 1,13, :) = (/&
	1038	& 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/)
	1039	GB( 1,13, :) = (/&
	1040	& 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/)
	1041	GA( 1,14, :) = (/&
	1042	& 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/)
	1043	GB( 1,14, :) = (/&
	1044	& 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/)
	1045	!
	1046	!----- INTERVAL = 2 ----- T = 200.0
	1047	!
	1048	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1049	GA( 2,13, :) = (/&
	1050	& 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/)
	1051	GB( 2,13, :) = (/&
	1052	& 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/)
	1053	GA( 2,14, :) = (/&
	1054	& 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/)
	1055	GB( 2,14, :) = (/&
	1056	& 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/)
	1057	!
	1058	!----- INTERVAL = 2 ----- T = 212.5
	1059	!
	1060	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1061	GA( 3,13, :) = (/&
	1062	& 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/)
	1063	GB( 3,13, :) = (/&
	1064	& 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/)
	1065	GA( 3,14, :) = (/&
	1066	& 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/)
	1067	GB( 3,14, :) = (/&
	1068	& 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/)
	1069	!
	1070	!----- INTERVAL = 2 ----- T = 225.0
	1071	!
	1072	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1073	GA( 4,13, :) = (/&
	1074	& 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/)
	1075	GB( 4,13, :) = (/&
	1076	& 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/)
	1077	GA( 4,14, :) = (/&
	1078	& 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/)
	1079	GB( 4,14, :) = (/&
	1080	& 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/)
	1081	!
	1082	!----- INTERVAL = 2 ----- T = 237.5
	1083	!
	1084	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1085	GA( 5,13, :) = (/&
	1086	& 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/)
	1087	GB( 5,13, :) = (/&
	1088	& 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/)
	1089	GA( 5,14, :) = (/&
	1090	& 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/)
	1091	GB( 5,14, :) = (/&
	1092	& 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/)
	1093	!
	1094	!----- INTERVAL = 2 ----- T = 250.0
	1095	!
	1096	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1097	GA( 6,13, :) = (/&
	1098	& 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/)
	1099	GB( 6,13, :) = (/&
	1100	& 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/)
	1101	GA( 6,14, :) = (/&
	1102	& 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/)
	1103	GB( 6,14, :) = (/&
	1104	& 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/)
	1105	!
	1106	!----- INTERVAL = 2 ----- T = 262.5
	1107	!
	1108	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1109	GA( 7,13, :) = (/&
	1110	& 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/)
	1111	GB( 7,13, :) = (/&
	1112	& 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/)
	1113	GA( 7,14, :) = (/&
	1114	& 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/)
	1115	GB( 7,14, :) = (/&
	1116	& 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/)
	1117	!
	1118	!----- INTERVAL = 2 ----- T = 275.0
	1119	!
	1120	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1121	GA( 8,13, :) = (/&
	1122	& 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/)
	1123	GB( 8,13, :) = (/&
	1124	& 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/)
	1125	GA( 8,14, :) = (/&
	1126	& 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/)
	1127	GB( 8,14, :) = (/&
	1128	& 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/)
	1129	!
	1130	!----- INTERVAL = 2 ----- T = 287.5
	1131	!
	1132	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1133	GA( 9,13, :) = (/&
	1134	& 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/)
	1135	GB( 9,13, :) = (/&
	1136	& 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/)
	1137	GA( 9,14, :) = (/&
	1138	& 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/)
	1139	GB( 9,14, :) = (/&
	1140	& 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/)
	1141	!
	1142	!----- INTERVAL = 2 ----- T = 300.0
	1143	!
	1144	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1145	GA(10,13, :) = (/&
	1146	& 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/)
	1147	GB(10,13, :) = (/&
	1148	& 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/)
	1149	GA(10,14, :) = (/&
	1150	& 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/)
	1151	GB(10,14, :) = (/&
	1152	& 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/)
	1153	!
	1154	!----- INTERVAL = 2 ----- T = 312.5
	1155	!
	1156	!-- INDICES FOR PADE APPROXIMATION 1 30 38 45
	1157	GA(11,13, :) = (/&
	1158	& 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/)
	1159	GB(11,13, :) = (/&
	1160	& 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/)
	1161	GA(11,14, :) = (/&
	1162	& 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/)
	1163	GB(11,14, :) = (/&
	1164	& 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/)
	1165	!
	1166	!
	1167	!
	1168	!
	1169	!
	1170	!
	1171	!
	1172	!
	1173	!
	1174	!
	1175	!-- CARBON DIOXIDE LINES IN THE WINDOW REGION (800-1250 CM-1)
	1176	!
	1177	!
	1178	!-- G = 0.0
	1179	!
	1180	!
	1181	!----- INTERVAL = 4 ----- T = 187.5
	1182	!
	1183	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1184	GA( 1,15, :) = (/&
	1185	& 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/)
	1186	GB( 1,15, :) = (/&
	1187	& 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/)
	1188	GA( 1,16, :) = (/&
	1189	& 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/)
	1190	GB( 1,16, :) = (/&
	1191	& 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/)
	1192	!
	1193	!----- INTERVAL = 4 ----- T = 200.0
	1194	!
	1195	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1196	GA( 2,15, :) = (/&
	1197	& 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/)
	1198	GB( 2,15, :) = (/&
	1199	& 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/)
	1200	GA( 2,16, :) = (/&
	1201	& 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/)
	1202	GB( 2,16, :) = (/&
	1203	& 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/)
	1204	!
	1205	!----- INTERVAL = 4 ----- T = 212.5
	1206	!
	1207	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1208	GA( 3,15, :) = (/&
	1209	& 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/)
	1210	GB( 3,15, :) = (/&
	1211	& 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/)
	1212	GA( 3,16, :) = (/&
	1213	& 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/)
	1214	GB( 3,16, :) = (/&
	1215	& 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/)
	1216	!
	1217	!----- INTERVAL = 4 ----- T = 225.0
	1218	!
	1219	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1220	GA( 4,15, :) = (/&
	1221	& 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/)
	1222	GB( 4,15, :) = (/&
	1223	& 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/)
	1224	GA( 4,16, :) = (/&
	1225	& 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/)
	1226	GB( 4,16, :) = (/&
	1227	& 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/)
	1228	!
	1229	!----- INTERVAL = 4 ----- T = 237.5
	1230	!
	1231	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1232	GA( 5,15, :) = (/&
	1233	& 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/)
	1234	GB( 5,15, :) = (/&
	1235	& 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/)
	1236	GA( 5,16, :) = (/&
	1237	& 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/)
	1238	GB( 5,16, :) = (/&
	1239	& 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/)
	1240	!
	1241	!----- INTERVAL = 4 ----- T = 250.0
	1242	!
	1243	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1244	GA( 6,15, :) = (/&
	1245	& 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/)
	1246	GB( 6,15, :) = (/&
	1247	& 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/)
	1248	GA( 6,16, :) = (/&
	1249	& 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/)
	1250	GB( 6,16, :) = (/&
	1251	& 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/)
	1252	!
	1253	!----- INTERVAL = 4 ----- T = 262.5
	1254	!
	1255	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1256	GA( 7,15, :) = (/&
	1257	& 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/)
	1258	GB( 7,15, :) = (/&
	1259	& 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/)
	1260	GA( 7,16, :) = (/&
	1261	& 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/)
	1262	GB( 7,16, :) = (/&
	1263	& 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/)
	1264	!
	1265	!----- INTERVAL = 4 ----- T = 275.0
	1266	!
	1267	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1268	GA( 8,15, :) = (/&
	1269	& 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/)
	1270	GB( 8,15, :) = (/&
	1271	& 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/)
	1272	GA( 8,16, :) = (/&
	1273	& 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/)
	1274	GB( 8,16, :) = (/&
	1275	& 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/)
	1276	!
	1277	!----- INTERVAL = 4 ----- T = 287.5
	1278	!
	1279	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1280	GA( 9,15, :) = (/&
	1281	& 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/)
	1282	GB( 9,15, :) = (/&
	1283	& 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/)
	1284	GA( 9,16, :) = (/&
	1285	& 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/)
	1286	GB( 9,16, :) = (/&
	1287	& 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/)
	1288	!
	1289	!----- INTERVAL = 4 ----- T = 300.0
	1290	!
	1291	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1292	GA(10,15, :) = (/&
	1293	& 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/)
	1294	GB(10,15, :) = (/&
	1295	& 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/)
	1296	GA(10,16, :) = (/&
	1297	& 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/)
	1298	GB(10,16, :) = (/&
	1299	& 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/)
	1300	!
	1301	!----- INTERVAL = 4 ----- T = 312.5
	1302	!
	1303	!-- INDICES FOR PADE APPROXIMATION 1 15 29 45
	1304	GA(11,15, :) = (/&
	1305	& 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/)
	1306	GB(11,15, :) = (/&
	1307	& 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/)
	1308	GA(11,16, :) = (/&
	1309	& 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/)
	1310	GB(11,16, :) = (/&
	1311	& 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/)
	1312	!
	1313	RETURN
	1314	END SUBROUTINE SUOLW

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

Download in other formats:

Original Format