Context Navigation

suswn.F90 @ 2172

Last change on this file since 2172 was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM

Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence

Property copyright set to
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

File size: 43.6 KB

Rev	Line
[1989]	1	SUBROUTINE SUSWN (KTSW, KSW)
	2
	3	!**** SUSW - INITIALIZE COMMON YOESW
	4
	5	! PURPOSE.
	6	! --------
	7	! INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS
	8	! NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES
	9
	10	!** INTERFACE.
	11	! ----------
	12	! CALL *SUSW
	13
	14	! EXPLICIT ARGUMENTS :
	15	! --------------------
	16	! NONE
	17
	18	! IMPLICIT ARGUMENTS :
	19	! --------------------
	20	! COMMON YOESW
	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	! 97-04-16 JJ Morcrette 2 and 4 interval spectral resolution
	41	! 00-10-24 JJ Morcrette sea-ice albedo revisited
	42	! 00-12-14 JJ Morcrette
	43	! and Ph.Dubuisson B.Bonnel 6 spectral interval resolution
	44	! 01-04-17 Ph.Dubuisson, B.Bonnel, JJ.Morcrette 6 sp.int.resolu'n
	45	! 01-06-28 B.Bonnel, JJ.Morcrette, Ph.Dubuisson Rayleigh (2/4/6)
	46	! 01-11-05 Ph.Dubuisson, JJMorcrette (new 2 intervals for TL/AD)
	47	! 03-04-24 JJ Morcrette SRTM with 14 spectral intervals
	48	! M.Hamrud 01-Oct-2003 CY28 Cleaning
	49
	50	! ------------------------------------------------------------------
	51
	52	USE PARKIND1 ,ONLY : JPIM ,JPRB
	53	USE YOMHOOK ,ONLY : LHOOK, DR_HOOK
	54
	55	USE YOMCST , ONLY : RG
	56	USE YOERAD , ONLY : NTSW ,LONEWSW ,LRAYL
	57	USE YOESW , ONLY : APAD ,BPAD ,RRAY ,RSUN ,&
	58	& RPDH1 ,RPDU1 ,RPNH ,RPNU ,RSWCE ,&
	59	& RSWCP ,RTDH2O ,RTDUMG ,RTH2O ,RTUMG ,&
	60	& D ,RROMA ,RROMB ,RRASY ,RHSRA ,&
	61	& RHSRB ,RHSRC ,RHSRD ,RHSRE ,RHSRF ,&
	62	& RHSRTA ,RHSRTB ,REXPO3 ,&
	63	& RADJUST ,NEXPO3
	64
	65	IMPLICIT NONE
	66
	67	INTEGER(KIND=JPIM),INTENT(IN) :: KTSW
	68	INTEGER(KIND=JPIM),INTENT(IN) :: KSW
	69	INTEGER(KIND=JPIM) :: IMPSRTM(14)
	70
	71	! ----------------------------------------------------------------
	72	REAL(KIND=JPRB) :: ZAPAD2(2,3,7) , ZBPAD2(2,3,7) , ZD2(2,3)&
	73	& , ZRAY2(2,6) , ZSUN2(2) , ZSWCE2(2) , ZSWCP2(2)&
	74	& , ZROMA2(2) , ZROMB2(2) , ZRASY2(2)&
	75	& , ZRA2(2) , ZRB2(2) , ZRC2(2)&
	76	& , ZRD2(2) , ZRE2(2) , ZRF2(2)
	77
	78	REAL(KIND=JPRB) :: ZAPAD4(4,3,7) , ZBPAD4(4,3,7) , ZD4(4,3)&
	79	& , ZRAY4(4,6) , ZSUN4(4) , ZSWCE4(4) , ZSWCP4(4)&
	80	& , ZROMA4(4) , ZROMB4(4) , ZRASY4(4)&
	81	& , ZRA4(4) , ZRB4(4) , ZRC4(4)&
	82	& , ZRD4(4) , ZRE4(4) , ZRF4(4)
	83
	84	REAL(KIND=JPRB) :: ZAPAD6(6,3,7) , ZBPAD6(6,3,7) , ZD6(6,3)&
	85	& , ZRAY6(6,6) , ZSUN6(6) , ZSWCE6(6) , ZSWCP6(6)&
	86	& , ZROMA6(6) , ZROMB6(6) , ZRASY6(6)&
	87	& , ZRA6(6) , ZRB6(6) , ZRC6(6)&
	88	& , ZRD6(6) , ZRE6(6) , ZRF6(6)&
	89	& , ZRBY6(6,6)
	90
	91
	92	INTEGER(KIND=JPIM) :: JC3, JC6, JI, JJ, JW, I_K, ISW6
	93
	94	REAL(KIND=JPRB) :: ZH2O, ZPDHNIS, ZPDH2O, ZPDUNIS,&
	95	& ZPDUMG, ZPRHNIS, ZPRH2O, &
	96	& ZPRUNIS, ZPRUMG, ZRTO1, ZRTO2, ZTDHNIS, &
	97	& ZTDUNIS, ZTHNIS, &
	98	& ZTUNIS, ZUMG, ZADJUS2, ZADJUS4, ZADJUS6
	99
	100	REAL(KIND=JPRB) :: ZTDH2IS, ZTDU2IS, ZTH2IS, ZTU2IS, ZPDH2IS, ZPDU2IS, ZPRH2IS, ZPRU2IS
	101	REAL(KIND=JPRB) :: ZTDH4IS, ZTDU4IS, ZTH4IS, ZTU4IS, ZPDH4IS, ZPDU4IS, ZPRH4IS, ZPRU4IS
	102	REAL(KIND=JPRB) :: ZTDH6IS, ZTDU6IS, ZTH6IS, ZTU6IS, ZPDH6IS, ZPDU6IS, ZPRH6IS, ZPRU6IS
	103	REAL(KIND=JPRB) :: ZHOOK_HANDLE
	104
	105	#include "abor1.intfb.h"
	106
	107	! ----------------------------------------------------------------
	108
	109	!* 1. CLEAR-SKY ABSORPTION COEFFICIENTS FOR N SPECTRAL INTERVALS
	110	! --------------------------------------------------------
	111
	112	!* spectrally independent
	113
	114	!-- New coefficients
	115
	116	IF (LHOOK) CALL DR_HOOK('SUSWN',0,ZHOOK_HANDLE)
	117	ZTDHNIS = 0.450_JPRB
	118	ZTDUNIS = 0.375_JPRB
	119	ZTHNIS = 273._JPRB
	120	ZTUNIS = 273._JPRB
	121	ZPDHNIS = 0.90_JPRB
	122	ZPDUNIS = 0.75_JPRB
	123	ZPRHNIS = 101300._JPRB
	124	ZPRUNIS = 101300._JPRB
	125	ISW6=6
	126
	127	! ----------------------------------------------------------------
	128
	129	!* 1.1 COEFFICIENTS FOR TWO SPECTRAL INTERVALS
	130	! ---------------------------------------
	131
	132	!* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
	133	! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
	134	! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
	135	! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
	136
	137	ZADJUS2 = 1.0_JPRB
	138	ZTDH2IS = 0.450_JPRB
	139	ZTDU2IS = 0.375_JPRB
	140	ZTH2IS = 273._JPRB
	141	ZTU2IS = 273._JPRB
	142	ZPDH2IS = 0.90_JPRB
	143	ZPDU2IS = 0.75_JPRB
	144	ZPRH2IS = 101300._JPRB
	145	ZPRU2IS = 101300._JPRB
	146
	147	!* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
	148
	149	ZSUN2(1) = 0.459760_JPRB
	150
	151	!-- New coefficients
	152	ZD2(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
	153
	154	ZAPAD2(1, 1, :) = (/&
	155	& 0.1095121E+03_JPRB,&
	156	& 0.8906647E+03_JPRB,&
	157	& 0.1030601E+03_JPRB,&
	158	& 0.000000000E-00_JPRB,&
	159	& 0.000000000E-00_JPRB,&
	160	& 0.000000000E-00_JPRB,&
	161	& 0.000000000E-00_JPRB/)
	162	ZAPAD2(1, 2, :) = (/&
	163	& 0.5124541E+03_JPRB,&
	164	& 0.3478800E+04_JPRB,&
	165	& 0.1244272E+04_JPRB,&
	166	& 0.000000000E-00_JPRB,&
	167	& 0.000000000E-00_JPRB,&
	168	& 0.000000000E-00_JPRB,&
	169	& 0.000000000E-00_JPRB/)
	170	ZAPAD2(1, 3, :) = (/&
	171	& 0.4766857E-04_JPRB,&
	172	& 0.4019806E-02_JPRB,&
	173	& 0.1936247E+00_JPRB,&
	174	& 0.1886522E+02_JPRB,&
	175	& 0.6498215E+02_JPRB,&
	176	& 0.1799387E+02_JPRB,&
	177	& 0.000000000E+00_JPRB/)
	178
	179	ZBPAD2(1, 1, :) = (/&
	180	& 0.1095121E+03_JPRB,&
	181	& 0.8919103E+03_JPRB,&
	182	& 0.1142218E+03_JPRB,&
	183	& 0.100000000E+01_JPRB,&
	184	& 0.000000000E-00_JPRB,&
	185	& 0.000000000E-00_JPRB,&
	186	& 0.000000000E-00_JPRB/)
	187	ZBPAD2(1, 2, :) = (/&
	188	& 0.5124541E+03_JPRB,&
	189	& 0.3499516E+04_JPRB,&
	190	& 0.1278186E+04_JPRB,&
	191	& 0.100000000E+01_JPRB,&
	192	& 0.000000000E-00_JPRB,&
	193	& 0.000000000E-00_JPRB,&
	194	& 0.000000000E-00_JPRB/)
	195	ZBPAD2(1, 3, :) = (/&
	196	& 0.4766857E-04_JPRB,&
	197	& 0.4153343E-02_JPRB,&
	198	& 0.1925810E+00_JPRB,&
	199	& 0.1938938E+02_JPRB,&
	200	& 0.6918191E+02_JPRB,&
	201	& 0.2278795E+02_JPRB,&
	202	& 0.100000000E+01_JPRB/)
	203
	204	ZRAY2(1,:)= (/&
	205	& .389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,&
	206	& .659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/)
	207
	208	!-- Original coefficients of ECMWF model (before 27 June 2000)
	209
	210	!ZTDH2IS = 0.450_JPRB
	211	!ZTDU2IS = 0.375_JPRB
	212	!ZTH2IS = 296._JPRB
	213	!ZTU2IS = 296._JPRB
	214	!ZPDH2IS = 0.90_JPRB
	215	!ZPDU2IS = 0.75_JPRB
	216	!ZPRH2IS = 100000._JPRB
	217	!ZPRU2IS = 100000._JPRB
	218
	219	!ZD2(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
	220
	221	!ZAPAD2(1, 1, :) = (/&
	222	! &0.184678379E+06_JPRB,&
	223	! &0.553080884E+05_JPRB,&
	224	! &0.248143712E+04_JPRB,&
	225	! &0.000000000E-00_JPRB,&
	226	! &0.000000000E-00_JPRB,&
	227	! &0.000000000E-00_JPRB,&
	228	! &0.000000000E-00_JPRB/)
	229	!ZAPAD2(1, 2, :) = (/&
	230	! &0.715303869E+01_JPRB,&
	231	! &0.219386847E+03_JPRB,&
	232	! &0.830001089E+03_JPRB,&
	233	! &0.000000000E-00_JPRB,&
	234	! &0.000000000E-00_JPRB,&
	235	! &0.000000000E-00_JPRB,&
	236	! &0.000000000E-00_JPRB/)
	237	!ZAPAD2(1, 3, :) = (/&
	238	! &0.925887084E-04_JPRB,&
	239	! &0.129353723E-01_JPRB,&
	240	! &0.800821928E+00_JPRB,&
	241	! &0.242715973E+02_JPRB,&
	242	! &0.878331486E+02_JPRB,&
	243	! &0.191559725E+02_JPRB,&
	244	! &0.000000000E+00_JPRB/)
	245
	246	!ZBPAD2(1, 1, :) = (/&
	247	! &0.184678379E+06_JPRB,&
	248	! &0.555188347E+05_JPRB,&
	249	! &0.253257443E+04_JPRB,&
	250	! &0.100000000E+01_JPRB,&
	251	! &0.000000000E-00_JPRB,&
	252	! &0.000000000E-00_JPRB,&
	253	! &0.000000000E-00_JPRB/)
	254	!ZBPAD2(1, 2, :) = (/&
	255	! &0.715303869E+01_JPRB,&
	256	! &0.219441875E+03_JPRB,&
	257	! &0.831119997E+03_JPRB,&
	258	! &0.100000000E+01_JPRB,&
	259	! &0.000000000E-00_JPRB,&
	260	! &0.000000000E-00_JPRB,&
	261	! &0.000000000E-00_JPRB/)
	262	!ZBPAD2(1, 3, :) = (/&
	263	! &0.925887084E-04_JPRB,&
	264	! &0.131812683E-01_JPRB,&
	265	! &0.812706117E+00_JPRB,&
	266	! &0.249863591E+02_JPRB,&
	267	! &0.931071925E+02_JPRB,&
	268	! &0.252233437E+02_JPRB,&
	269	! &0.100000000E+01_JPRB/)
	270
	271	ZRAY2(1,:)= (/&
	272	& .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
	273	& .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
	274
	275	ZSWCE2(1) = 0.0_JPRB
	276	ZSWCP2(1) = 0.0_JPRB
	277
	278	! ----------------------------------------------------------------
	279
	280	!* Near-Infrared (0.69 - 4.0 Microns) is not sub-divided
	281
	282	! ----------------------------------------------------------------
	283
	284	!* 0.68 - 4.00 Micron
	285
	286	ZSUN2(2) = 0.540240_JPRB
	287
	288	!-- New coefficients
	289
	290	ZD2(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB/)
	291
	292	ZAPAD2(2, 1, :) = (/&
	293	& 0.1345066E-03_JPRB,&
	294	& 0.9276888E+00_JPRB,&
	295	& 0.3449105E+03_JPRB,&
	296	& 0.7143025E+04_JPRB,&
	297	& 0.6531248E+04_JPRB,&
	298	& 0.2210815E+03_JPRB,&
	299	& 0.0000000E+00_JPRB/)
	300	ZAPAD2(2, 2, :) = (/&
	301	& 0.7776703E-05_JPRB,&
	302	& 0.6851592E-02_JPRB,&
	303	& 0.5533947E+00_JPRB,&
	304	& 0.2538585E+01_JPRB,&
	305	& 0.000000000E+00_JPRB,&
	306	& 0.000000000E+00_JPRB,&
	307	& 0.000000000E+00_JPRB/)
	308	ZAPAD2(2, 3, :) = (/&
	309	& 0.932012720E+02_JPRB,&
	310	& 0.437866983E+03_JPRB,&
	311	& 0.000000000E-00_JPRB,&
	312	& 0.000000000E-00_JPRB,&
	313	& 0.000000000E-00_JPRB,&
	314	& 0.000000000E-00_JPRB,&
	315	& 0.000000000E+00_JPRB/)
	316
	317	ZBPAD2(2, 1, :) = (/&
	318	& 0.1345066E-03_JPRB,&
	319	& 0.9297728E+00_JPRB,&
	320	& 0.3523821E+03_JPRB,&
	321	& 0.7931590E+04_JPRB,&
	322	& 0.8979086E+04_JPRB,&
	323	& 0.5050864E+03_JPRB,&
	324	& 0.1000000E+01_JPRB/)
	325	ZBPAD2(2, 2, :) = (/&
	326	& 0.7776703E-05_JPRB,&
	327	& 0.6976301E-02_JPRB,&
	328	& 0.6074521E+00_JPRB,&
	329	& 0.3517521E+01_JPRB,&
	330	& 0.100000000E+01_JPRB,&
	331	& 0.000000000E+00_JPRB,&
	332	& 0.000000000E+00_JPRB/)
	333	ZBPAD2(2, 3, :) = (/&
	334	& 0.932012720E+02_JPRB,&
	335	& 0.439078599E+03_JPRB,&
	336	& 0.100000000E+01_JPRB,&
	337	& 0.000000000E+00_JPRB,&
	338	& 0.000000000E+00_JPRB,&
	339	& 0.000000000E+00_JPRB,&
	340	& 0.000000000E+00_JPRB/)
	341
	342	ZRAY2(2,:)= (/&
	343	& .104635E-01_JPRB, .102535E+00_JPRB,-.536771E+00_JPRB,&
	344	& .123640E+01_JPRB,-.128930E+01_JPRB, .502426E+00_JPRB/)
	345
	346	!-- Original coefficients of ECMWF model (before 27 June 2000)
	347
	348	!ZD2(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB/)
	349
	350	!ZAPAD2(2, 1, :) = (/&
	351	! &0.577160808E-06_JPRB,&
	352	! &0.681760991E-02_JPRB,&
	353	! &0.502416370E+01_JPRB,&
	354	! &0.209771270E+03_JPRB,&
	355	! &0.644493139E+03_JPRB,&
	356	! &0.120654163E+03_JPRB,&
	357	! &0.000000000E+00_JPRB/)
	358	!ZAPAD2(2, 2, :) = (/&
	359	! &0.422874202E-09_JPRB,&
	360	! &0.257499477E-05_JPRB,&
	361	! &0.349570458E-02_JPRB,&
	362	! &0.917758876E+00_JPRB,&
	363	! &0.327669989E+02_JPRB,&
	364	! &0.874155115E+02_JPRB,&
	365	! &0.000000000E+00_JPRB/)
	366	!ZAPAD2(2, 3, :) = (/&
	367	! &0.932012720E+02_JPRB,&
	368	! &0.437866983E+03_JPRB,&
	369	! &0.000000000E-00_JPRB,&
	370	! &0.000000000E-00_JPRB,&
	371	! &0.000000000E-00_JPRB,&
	372	! &0.000000000E-00_JPRB,&
	373	! &0.000000000E+00_JPRB/)
	374
	375	!ZBPAD2(2, 1, :) = (/&
	376	! &0.577160808E-06_JPRB,&
	377	! &0.682595589E-02_JPRB,&
	378	! &0.508797283E+01_JPRB,&
	379	! &0.223901000E+03_JPRB,&
	380	! &0.770640908E+03_JPRB,&
	381	! &0.186812480E+03_JPRB,&
	382	! &0.100000000E+01_JPRB/)
	383	!ZBPAD2(2, 2, :) = (/&
	384	! &0.422874202E-09_JPRB,&
	385	! &0.257567242E-05_JPRB,&
	386	! &0.349912408E-02_JPRB,&
	387	! &0.920770481E+00_JPRB,&
	388	! &0.330990960E+02_JPRB,&
	389	! &0.903249595E+02_JPRB,&
	390	! &0.100000000E+01_JPRB/)
	391	!ZBPAD2(2, 3, :) = (/&
	392	! &0.932012720E+02_JPRB,&
	393	! &0.439078599E+03_JPRB,&
	394	! &0.100000000E+01_JPRB,&
	395	! &0.000000000E+00_JPRB,&
	396	! &0.000000000E+00_JPRB,&
	397	! &0.000000000E+00_JPRB,&
	398	! &0.000000000E+00_JPRB/)
	399
	400	!ZRAY2(2,:)= (/&
	401	! &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,&
	402	! &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/)
	403
	404	ZSWCE2(2) = 0.0_JPRB
	405	ZSWCP2(2) = 0.0_JPRB
	406
	407	!=====================================================================
	408
	409	!* 1.2 COEFFICIENTS FOR FOUR SPECTRAL INTERVALS
	410	! ----------------------------------------
	411
	412	!* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
	413	! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
	414	! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
	415	! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
	416
	417	ZADJUS4 = 1.0_JPRB
	418
	419	!* 1st spectral interval: U.V. and Visible (0.25 - 0.69 Micron)
	420
	421	ZSUN4(1) = 0.459760_JPRB
	422
	423	!-- New coefficients
	424
	425	!ZD4(1,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
	426
	427	!ZAPAD4(1, 1, :) = (/&
	428	! &0.1095121E+03_JPRB,&
	429	! &0.8906647E+03_JPRB,&
	430	! &0.1030601E+03_JPRB,&
	431	! &0.0000000E-00_JPRB,&
	432	! &0.0000000E-00_JPRB,&
	433	! &0.0000000E-00_JPRB,&
	434	! &0.0000000E-00_JPRB/)
	435	!ZAPAD4(1, 2, :) = (/&
	436	! &0.5124541E+03_JPRB,&
	437	! &0.3478800E+04_JPRB,&
	438	! &0.1244272E+04_JPRB,&
	439	! &0.0000000E-00_JPRB,&
	440	! &0.0000000E-00_JPRB,&
	441	! &0.0000000E-00_JPRB,&
	442	! &0.0000000E-00_JPRB/)
	443	!ZAPAD4(1, 3, :) = (/&
	444	! &0.4766857E-04_JPRB,&
	445	! &0.4019806E-02_JPRB,&
	446	! &0.1936247E+00_JPRB,&
	447	! &0.1886522E+02_JPRB,&
	448	! &0.6498215E+02_JPRB,&
	449	! &0.1799387E+02_JPRB,&
	450	! &0.0000000E+00_JPRB/)
	451
	452	!ZBPAD4(1, 1, :) = (/&
	453	! &0.1095121E+03_JPRB,&
	454	! &0.8919103E+03_JPRB,&
	455	! &0.1142218E+03_JPRB,&
	456	! &0.1000000E+01_JPRB,&
	457	! &0.0000000E-00_JPRB,&
	458	! &0.0000000E-00_JPRB,&
	459	! &0.0000000E-00_JPRB/)
	460	!ZBPAD4(1, 2, :) = (/&
	461	! &0.5124541E+03_JPRB,&
	462	! &0.3499516E+04_JPRB,&
	463	! &0.1278186E+04_JPRB,&
	464	! &0.1000000E+01_JPRB,&
	465	! &0.0000000E-00_JPRB,&
	466	! &0.0000000E-00_JPRB,&
	467	! &0.0000000E-00_JPRB/)
	468	!ZBPAD4(1, 3, :) = (/&
	469	! &0.4766857E-04_JPRB,&
	470	! &0.4153343E-02_JPRB,&
	471	! &0.1925810E+00_JPRB,&
	472	! &0.1938938E+02_JPRB,&
	473	! &0.6918191E+02_JPRB,&
	474	! &0.2278795E+02_JPRB,&
	475	! &0.1000000E+01_JPRB/)
	476
	477	!ZRAY4(1,:)= (/&
	478	! &.389107E-01_JPRB, .101608E+01_JPRB,-.354921E+01_JPRB,&
	479	! &.659510E+01_JPRB,-.592423E+01_JPRB, .202348E+01_JPRB/)
	480
	481	!-- Original coefficients in ECMWF scheme
	482
	483	ZTDH4IS = 0.450_JPRB
	484	ZTDU4IS = 0.375_JPRB
	485	ZTH4IS = 296._JPRB
	486	ZTU4IS = 296._JPRB
	487	ZPDH4IS = 0.90_JPRB
	488	ZPDU4IS = 0.75_JPRB
	489	ZPRH4IS = 100000._JPRB
	490	ZPRU4IS = 100000._JPRB
	491
	492	ZD4(1,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
	493
	494	ZAPAD4(1, 1, :) = (/&
	495	& 0.184678379E+06_JPRB,&
	496	& 0.553080884E+05_JPRB,&
	497	& 0.248143712E+04_JPRB,&
	498	& 0.000000000E-00_JPRB,&
	499	& 0.000000000E-00_JPRB,&
	500	& 0.000000000E-00_JPRB,&
	501	& 0.000000000E-00_JPRB/)
	502	ZAPAD4(1, 2, :) = (/&
	503	& 0.715303869E+01_JPRB,&
	504	& 0.219386847E+03_JPRB,&
	505	& 0.830001089E+03_JPRB,&
	506	& 0.000000000E-00_JPRB,&
	507	& 0.000000000E-00_JPRB,&
	508	& 0.000000000E-00_JPRB,&
	509	& 0.000000000E-00_JPRB/)
	510	ZAPAD4(1, 3, :) = (/&
	511	& 0.925887084E-04_JPRB,&
	512	& 0.129353723E-01_JPRB,&
	513	& 0.800821928E+00_JPRB,&
	514	& 0.242715973E+02_JPRB,&
	515	& 0.878331486E+02_JPRB,&
	516	& 0.191559725E+02_JPRB,&
	517	& 0.000000000E+00_JPRB/)
	518
	519	ZBPAD4(1, 1, :) = (/&
	520	& 0.184678379E+06_JPRB,&
	521	& 0.555188347E+05_JPRB,&
	522	& 0.253257443E+04_JPRB,&
	523	& 0.100000000E+01_JPRB,&
	524	& 0.000000000E-00_JPRB,&
	525	& 0.000000000E-00_JPRB,&
	526	& 0.000000000E-00_JPRB/)
	527	ZBPAD4(1, 2, :) = (/&
	528	& 0.715303869E+01_JPRB,&
	529	& 0.219441875E+03_JPRB,&
	530	& 0.831119997E+03_JPRB,&
	531	& 0.100000000E+01_JPRB,&
	532	& 0.000000000E-00_JPRB,&
	533	& 0.000000000E-00_JPRB,&
	534	& 0.000000000E-00_JPRB/)
	535	ZBPAD4(1, 3, :) = (/&
	536	& 0.925887084E-04_JPRB,&
	537	& 0.131812683E-01_JPRB,&
	538	& 0.812706117E+00_JPRB,&
	539	& 0.249863591E+02_JPRB,&
	540	& 0.931071925E+02_JPRB,&
	541	& 0.252233437E+02_JPRB,&
	542	& 0.100000000E+01_JPRB/)
	543
	544	ZRAY4(1,:)= (/&
	545	& .428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
	546	& .522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
	547
	548	ZSWCE4(1) = 0.0_JPRB
	549	ZSWCP4(1) = 0.0_JPRB
	550
	551	! ----------------------------------------------------------------
	552
	553	!* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
	554
	555	! ----------------------------------------------------------------
	556
	557	!* 0.69 - 1.19 Micron
	558
	559	ZSUN4(2) = 0.326158_JPRB
	560
	561	!-- New coefficients
	562
	563	!ZD4(2,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /)
	564
	565	!ZAPAD4(2, 1, :) = (/&
	566	! &0.1335726E+02_JPRB,&
	567	! &0.2939136E+04_JPRB,&
	568	! &0.4010585E+05_JPRB,&
	569	! &0.7195030E+05_JPRB,&
	570	! &0.1648338E+05_JPRB,&
	571	! &0.3373738E+03_JPRB,&
	572	! &0.0000000E+00_JPRB/)
	573	!ZAPAD4(2, 2, :) = (/&
	574	! &0.2001271E-01_JPRB,&
	575	! &0.2480831E+01_JPRB,&
	576	! &0.3444162E+02_JPRB,&
	577	! &0.4788946E+02_JPRB,&
	578	! &0.0000000E+00_JPRB,&
	579	! &0.0000000E+00_JPRB,&
	580	! &0.0000000E+00_JPRB/)
	581	!ZAPAD4(2, 3, :) = (/&
	582	! &0.000000000E+00_JPRB,&
	583	! &0.000000000E+00_JPRB,&
	584	! &0.000000000E+00_JPRB,&
	585	! &0.000000000E+00_JPRB,&
	586	! &0.000000000E+00_JPRB,&
	587	! &0.000000000E+00_JPRB,&
	588	! &0.000000000E+00_JPRB/)
	589
	590	!ZBPAD4(2, 1, :) = (/&
	591	! &0.1335726E+02_JPRB,&
	592	! &0.2942327E+04_JPRB,&
	593	! &0.4077237E+05_JPRB,&
	594	! &0.7749017E+05_JPRB,&
	595	! &0.2123132E+05_JPRB,&
	596	! &0.6659687E+03_JPRB,&
	597	! &0.1000000E+01_JPRB/)
	598	!ZBPAD4(2, 2, :) = (/&
	599	! &0.2001271E-01_JPRB,&
	600	! &0.2549067E+01_JPRB,&
	601	! &0.3752433E+02_JPRB,&
	602	! &0.6276637E+02_JPRB,&
	603	! &0.1000000E+01_JPRB,&
	604	! &0.0000000E+00_JPRB,&
	605	! &0.0000000E+00_JPRB/)
	606	!ZBPAD4(2, 3, :) = (/&
	607	! &1.000000000E+00_JPRB,&
	608	! &0.000000000E+00_JPRB,&
	609	! &0.000000000E+00_JPRB,&
	610	! &0.000000000E+00_JPRB,&
	611	! &0.000000000E+00_JPRB,&
	612	! &0.000000000E+00_JPRB,&
	613	! &0.000000000E+00_JPRB/)
	614
	615	!DATA (ZRAY4(2,K),K=1,6) /&
	616	! & .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,&
	617	! &-.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/
	618
	619	!-- Original coefficients in the ECMWF model
	620
	621	ZD4(2,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /)
	622
	623	ZAPAD4(2, 1, :) = (/&
	624	& 0.690730834E-02_JPRB,&
	625	& 0.151704275E+01_JPRB,&
	626	& 0.751477543E+02_JPRB,&
	627	& 0.759770236E+03_JPRB,&
	628	& 0.109800326E+04_JPRB,&
	629	& 0.148407574E+03_JPRB,&
	630	& 0.000000000E+00_JPRB/)
	631	ZAPAD4(2, 2, :) = (/&
	632	& 0.863790752E-03_JPRB,&
	633	& 0.448762291E+00_JPRB,&
	634	& 0.332530367E+02_JPRB,&
	635	& 0.190914146E+03_JPRB,&
	636	& 0.000000000E+00_JPRB,&
	637	& 0.000000000E+00_JPRB,&
	638	& 0.000000000E+00_JPRB/)
	639	ZAPAD4(2, 3, :) = (/&
	640	& 0.000000000E+00_JPRB,&
	641	& 0.000000000E+00_JPRB,&
	642	& 0.000000000E+00_JPRB,&
	643	& 0.000000000E+00_JPRB,&
	644	& 0.000000000E+00_JPRB,&
	645	& 0.000000000E+00_JPRB,&
	646	& 0.000000000E+00_JPRB/)
	647
	648	ZBPAD4(2, 1, :) = (/&
	649	& 0.690730834E-02_JPRB,&
	650	& 0.151954406E+01_JPRB,&
	651	& 0.756512527E+02_JPRB,&
	652	& 0.779384997E+03_JPRB,&
	653	& 0.121113108E+04_JPRB,&
	654	& 0.207678436E+03_JPRB,&
	655	& 0.100000000E+01_JPRB/)
	656	ZBPAD4(2, 2, :) = (/&
	657	& 0.863790752E-03_JPRB,&
	658	& 0.448948107E+00_JPRB,&
	659	& 0.333186750E+02_JPRB,&
	660	& 0.192727216E+03_JPRB,&
	661	& 0.100000000E+01_JPRB,&
	662	& 0.000000000E+00_JPRB,&
	663	& 0.000000000E+00_JPRB/)
	664	ZBPAD4(2, 3, :) = (/&
	665	& 1.000000000E+00_JPRB,&
	666	& 0.000000000E+00_JPRB,&
	667	& 0.000000000E+00_JPRB,&
	668	& 0.000000000E+00_JPRB,&
	669	& 0.000000000E+00_JPRB,&
	670	& 0.000000000E+00_JPRB,&
	671	& 0.000000000E+00_JPRB/)
	672
	673	DATA (ZRAY4(2,I_K),I_K=1,6) /&
	674	& .164261E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	675	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/
	676
	677	ZSWCE4(2) = 0.0_JPRB
	678	ZSWCP4(2) = 0.0_JPRB
	679
	680	! ----------------------------------------------------------------
	681
	682	!* 1.19 - 2.38 Microns
	683
	684	ZSUN4(3) = 0.180608_JPRB
	685
	686	!-- New coefficients
	687
	688	!ZD4(3,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
	689
	690	!ZAPAD4(3, 1, :) = (/&
	691	! &0.3325841E-03_JPRB,&
	692	! &0.6194496E+00_JPRB,&
	693	! &0.1497138E+03_JPRB,&
	694	! &0.2314864E+04_JPRB,&
	695	! &0.2380109E+04_JPRB,&
	696	! &0.9553823E+02_JPRB,&
	697	! &0.0000000E+00_JPRB/)
	698	!ZAPAD4(3, 2, :) = (/&
	699	! &0.4552471E-03_JPRB,&
	700	! &0.4084154E+00_JPRB,&
	701	! &0.6114905E+01_JPRB,&
	702	! &0.7102540E+01_JPRB,&
	703	! &0.0000000E+00_JPRB,&
	704	! &0.0000000E+00_JPRB,&
	705	! &0.0000000E+00_JPRB/)
	706	!ZAPAD4(3, 3, :) = (/&
	707	! &0.000000000E+00_JPRB,&
	708	! &0.000000000E+00_JPRB,&
	709	! &0.000000000E+00_JPRB,&
	710	! &0.000000000E+00_JPRB,&
	711	! &0.000000000E+00_JPRB,&
	712	! &0.000000000E+00_JPRB,&
	713	! &0.000000000E+00_JPRB/)
	714
	715	!ZBPAD4(3, 1, :) = (/&
	716	! &0.3325841E-03_JPRB,&
	717	! &0.6231947E+00_JPRB,&
	718	! &0.1553098E+03_JPRB,&
	719	! &0.2822458E+04_JPRB,&
	720	! &0.3885194E+04_JPRB,&
	721	! &0.2700235E+03_JPRB,&
	722	! &0.1000000E+01_JPRB/)
	723	!ZBPAD4(3, 2, :) = (/&
	724	! &0.4552471E-03_JPRB,&
	725	! &0.4088242E+00_JPRB,&
	726	! &0.6411905E+01_JPRB,&
	727	! &0.9444439E+01_JPRB,&
	728	! &0.1000000E+01_JPRB,&
	729	! &0.0000000E+00_JPRB,&
	730	! &0.0000000E+00_JPRB/)
	731	!ZBPAD4(3, 3, :) = (/&
	732	! &1.000000000E+00_JPRB,&
	733	! &0.000000000E+00_JPRB,&
	734	! &0.000000000E+00_JPRB,&
	735	! &0.000000000E+00_JPRB,&
	736	! &0.000000000E+00_JPRB,&
	737	! &0.000000000E+00_JPRB,&
	738	! &0.000000000E+00_JPRB/)
	739
	740	!ZRAY4(3,:)= (/&
	741	! &.288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,&
	742	! &.540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/)
	743
	744	ZD4(3,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 1.000000000_JPRB /)
	745
	746	ZAPAD4(3, 1, :) = (/&
	747	& 0.837531303E-05_JPRB,&
	748	& 0.173886341E-01_JPRB,&
	749	& 0.518852799E+01_JPRB,&
	750	& 0.159078416E+03_JPRB,&
	751	& 0.493273523E+03_JPRB,&
	752	& 0.102567293E+03_JPRB,&
	753	& 0.000000000E+00_JPRB/)
	754	ZAPAD4(3, 2, :) = (/&
	755	& 0.657978575E-02_JPRB,&
	756	& 0.752617872E+00_JPRB,&
	757	& 0.158209734E+02_JPRB,&
	758	& 0.410274915E+02_JPRB,&
	759	& 0.000000000E+00_JPRB,&
	760	& 0.000000000E+00_JPRB,&
	761	& 0.000000000E+00_JPRB/)
	762	ZAPAD4(3, 3, :) = (/&
	763	& 0.000000000E+00_JPRB,&
	764	& 0.000000000E+00_JPRB,&
	765	& 0.000000000E+00_JPRB,&
	766	& 0.000000000E+00_JPRB,&
	767	& 0.000000000E+00_JPRB,&
	768	& 0.000000000E+00_JPRB,&
	769	& 0.000000000E+00_JPRB/)
	770
	771	ZBPAD4(3, 1, :) = (/&
	772	& 0.837531303E-05_JPRB,&
	773	& 0.174882536E-01_JPRB,&
	774	& 0.534536580E+01_JPRB,&
	775	& 0.180351767E+03_JPRB,&
	776	& 0.673126838E+03_JPRB,&
	777	& 0.182718543E+03_JPRB,&
	778	& 0.100000000E+01_JPRB/)
	779	ZBPAD4(3, 2, :) = (/&
	780	& 0.657978575E-02_JPRB,&
	781	& 0.753752065E+00_JPRB,&
	782	& 0.159286262E+02_JPRB,&
	783	& 0.424278450E+02_JPRB,&
	784	& 0.100000000E+01_JPRB,&
	785	& 0.000000000E+00_JPRB,&
	786	& 0.000000000E+00_JPRB/)
	787	ZBPAD4(3, 3, :) = (/&
	788	& 1.000000000E+00_JPRB,&
	789	& 0.000000000E+00_JPRB,&
	790	& 0.000000000E+00_JPRB,&
	791	& 0.000000000E+00_JPRB,&
	792	& 0.000000000E+00_JPRB,&
	793	& 0.000000000E+00_JPRB,&
	794	& 0.000000000E+00_JPRB/)
	795
	796	ZRAY4(3,:)= (/&
	797	& .180438E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	798	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	799
	800	ZSWCE4(3) = 0.0_JPRB
	801	ZSWCP4(3) = 0.0_JPRB
	802
	803	! ----------------------------------------------------------------
	804
	805	!* 2.38 - 4.00 Microns
	806
	807	ZSUN4(4) = 0.033474_JPRB
	808
	809	!-- New coefficients
	810
	811	!ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
	812
	813	!ZAPAD4(4, 1, :) = (/&
	814	! &0.2122889E-06_JPRB,&
	815	! &0.9030576E-03_JPRB,&
	816	! &0.2431282E+00_JPRB,&
	817	! &0.4901345E+01_JPRB,&
	818	! &0.3996347E+01_JPRB,&
	819	! &0.3910227E+01_JPRB,&
	820	! &0.0000000E+00_JPRB/)
	821	!ZAPAD4(4, 2, :) = (/&
	822	! &0.1215163E-03_JPRB,&
	823	! &0.1222574E+00_JPRB,&
	824	! &0.9382420E+01_JPRB,&
	825	! &0.6875727E+02_JPRB,&
	826	! &0.2746421E+02_JPRB,&
	827	! &0.0000000E+00_JPRB,&
	828	! &0.0000000E+00_JPRB/)
	829	!ZAPAD4(4, 3, :) = (/&
	830	! &0.263068898E+02_JPRB,&
	831	! &0.146425875E+03_JPRB,&
	832	! &0.860137809E+02_JPRB,&
	833	! &0.000000000E+00_JPRB,&
	834	! &0.000000000E+00_JPRB,&
	835	! &0.000000000E+00_JPRB,&
	836	! &0.000000000E+00_JPRB/)
	837
	838	!ZBPAD4(4, 1, :) = (/&
	839	! &0.2122889E-06_JPRB,&
	840	! &0.9379083E-03_JPRB,&
	841	! &0.2957335E+00_JPRB,&
	842	! &0.8747190E+01_JPRB,&
	843	! &0.1015794E+02_JPRB,&
	844	! &0.1361277E+02_JPRB,&
	845	! &0.1000000E+01_JPRB/)
	846	!ZBPAD4(4, 2, :) = (/&
	847	! &0.1215163E-03_JPRB,&
	848	! &0.1255648E+00_JPRB,&
	849	! &0.1060119E+02_JPRB,&
	850	! &0.8414439E+02_JPRB,&
	851	! &0.4299438E+02_JPRB,&
	852	! &0.1000000E+01_JPRB,&
	853	! &0.0000000E+00_JPRB/)
	854	!ZBPAD4(4, 3, :) = (/&
	855	! &0.263068898E+02_JPRB,&
	856	! &0.152569217E+03_JPRB,&
	857	! &0.976791971E+02_JPRB,&
	858	! &0.100000000E+01_JPRB,&
	859	! &0.000000000E+00_JPRB,&
	860	! &0.000000000E+00_JPRB,&
	861	! &0.000000000E+00_JPRB/)
	862
	863	!ZRAY4(4,:)= (/&
	864	! &.493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,&
	865	! &.741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/)
	866
	867	ZD4(4,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
	868
	869	ZAPAD4(4, 1, :) = (/&
	870	& 0.122118185E-06_JPRB,&
	871	& 0.154042531E-02_JPRB,&
	872	& 0.141152193E+01_JPRB,&
	873	& 0.685368761E+02_JPRB,&
	874	& 0.216522281E+03_JPRB,&
	875	& 0.421228746E+02_JPRB,&
	876	& 0.000000000E+00_JPRB/)
	877	ZAPAD4(4, 2, :) = (/&
	878	& 0.364233560E-10_JPRB,&
	879	& 0.217340835E-06_JPRB,&
	880	& 0.292623386E-03_JPRB,&
	881	& 0.797100631E-01_JPRB,&
	882	& 0.319103672E+01_JPRB,&
	883	& 0.110530283E+02_JPRB,&
	884	& 0.000000000E+00_JPRB/)
	885	ZAPAD4(4, 3, :) = (/&
	886	& 0.263068898E+02_JPRB,&
	887	& 0.146425875E+03_JPRB,&
	888	& 0.860137809E+02_JPRB,&
	889	& 0.000000000E+00_JPRB,&
	890	& 0.000000000E+00_JPRB,&
	891	& 0.000000000E+00_JPRB,&
	892	& 0.000000000E+00_JPRB/)
	893
	894	ZBPAD4(4, 1, :) = (/&
	895	& 0.122118185E-06_JPRB,&
	896	& 0.156287582E-02_JPRB,&
	897	& 0.156994562E+01_JPRB,&
	898	& 0.102304103E+03_JPRB,&
	899	& 0.475479878E+03_JPRB,&
	900	& 0.188714799E+03_JPRB,&
	901	& 0.100000000E+01_JPRB/)
	902	ZBPAD4(4, 2, :) = (/&
	903	& 0.364233560E-10_JPRB,&
	904	& 0.218265806E-06_JPRB,&
	905	& 0.297085917E-03_JPRB,&
	906	& 0.834253403E-01_JPRB,&
	907	& 0.358290002E+01_JPRB,&
	908	& 0.139206164E+02_JPRB,&
	909	& 0.100000000E+01_JPRB/)
	910	ZBPAD4(4, 3, :) = (/&
	911	& 0.263068898E+02_JPRB,&
	912	& 0.152569217E+03_JPRB,&
	913	& 0.976791971E+02_JPRB,&
	914	& 0.100000000E+01_JPRB,&
	915	& 0.000000000E+00_JPRB,&
	916	& 0.000000000E+00_JPRB,&
	917	& 0.000000000E+00_JPRB/)
	918
	919	ZRAY4(4,:)= (/&
	920	& .136020E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	921	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	922
	923	ZSWCE4(4) = 0.0_JPRB
	924	ZSWCP4(4) = 0.0_JPRB
	925
	926	!=====================================================================
	927
	928	!* 1.3 COEFFICIENTS FOR SIX SPECTRAL INTERVALS
	929	! ---------------------------------------
	930
	931	!* DERIVED FROM HITRAN APRIL 1992 with LOWTRAN P AND T SCALING
	932	! H2O: Pref=1000hPa, Tref=296K, Pdep=0.9
	933	! UMG: Pref=1000hPa, Tref=296K, Pdep=0.75 (CO2+N2O+CO+CH4+O2)
	934	! O3 : unchanged in interval 1, from HITRAN 92 in interval 4
	935
	936	ZADJUS6 = 0.9905709_JPRB
	937	ZADJUS6 = 1.0_JPRB
	938
	939	ZTDH6IS = 0.450_JPRB
	940	ZTDU6IS = 0.375_JPRB
	941	ZTH6IS = 273._JPRB
	942	ZTU6IS = 273._JPRB
	943	ZPDH6IS = 0.90_JPRB
	944	ZPDU6IS = 0.75_JPRB
	945	ZPRH6IS = 101300._JPRB
	946	ZPRU6IS = 101300._JPRB
	947
	948	!* 1st spectral interval: U.V. (0.18 - 0.25 Micron)
	949
	950	ZSUN6(1) = 0.001917_JPRB
	951
	952	ZD6(1,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
	953
	954	ZAPAD6(1, 1, :) = (/&
	955	& 0.000000000E-00_JPRB,&
	956	& 0.000000000E-00_JPRB,&
	957	& 0.000000000E-00_JPRB,&
	958	& 0.000000000E-00_JPRB,&
	959	& 0.000000000E-00_JPRB,&
	960	& 0.000000000E-00_JPRB,&
	961	& 0.000000000E-00_JPRB/)
	962	ZAPAD6(1, 2, :) = (/&
	963	& 0.000000000E-00_JPRB,&
	964	& 0.000000000E-00_JPRB,&
	965	& 0.000000000E-00_JPRB,&
	966	& 0.000000000E-00_JPRB,&
	967	& 0.000000000E-00_JPRB,&
	968	& 0.000000000E-00_JPRB,&
	969	& 0.000000000E-00_JPRB/)
	970	ZAPAD6(1, 3, :) = (/&
	971	& 0.000000000E-00_JPRB,&
	972	& 0.000000000E-00_JPRB,&
	973	& 0.000000000E-00_JPRB,&
	974	& 0.000000000E-00_JPRB,&
	975	& 0.000000000E-00_JPRB,&
	976	& 0.000000000E-00_JPRB,&
	977	& 0.000000000E-00_JPRB/)
	978
	979	ZBPAD6(1, 1, :) = (/&
	980	& 0.100000000E+01_JPRB,&
	981	& 0.000000000E-00_JPRB,&
	982	& 0.000000000E-00_JPRB,&
	983	& 0.000000000E-00_JPRB,&
	984	& 0.000000000E-00_JPRB,&
	985	& 0.000000000E-00_JPRB,&
	986	& 0.000000000E-00_JPRB/)
	987	ZBPAD6(1, 2, :) = (/&
	988	& 0.100000000E+01_JPRB,&
	989	& 0.000000000E-00_JPRB,&
	990	& 0.000000000E-00_JPRB,&
	991	& 0.000000000E-00_JPRB,&
	992	& 0.000000000E-00_JPRB,&
	993	& 0.000000000E-00_JPRB,&
	994	& 0.000000000E-00_JPRB/)
	995	ZBPAD6(1, 3, :) = (/&
	996	& 0.100000000E+01_JPRB,&
	997	& 0.000000000E-00_JPRB,&
	998	& 0.000000000E-00_JPRB,&
	999	& 0.000000000E-00_JPRB,&
	1000	& 0.000000000E-00_JPRB,&
	1001	& 0.000000000E-00_JPRB,&
	1002	& 0.000000000E-00_JPRB/)
	1003
	1004	ZRAY6(1,:)= (/&
	1005	& .334018E+00_JPRB,-.298776E+01_JPRB, .613911E+02_JPRB,&
	1006	& -.164331E+03_JPRB, .197946E+03_JPRB,-.652206E+02_JPRB/)
	1007	ZRBY6(1,:)= (/&
	1008	& .43959002E+01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	1009	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	1010
	1011	ZSWCE6(1) = 0.0_JPRB
	1012	ZSWCP6(1) = 0.0_JPRB
	1013
	1014	NEXPO3(1) = 7
	1015	REXPO3(1, 1, :) =(/&
	1016	& 0.051395E+00_JPRB, 0.048250E+00_JPRB, 0.112339E+00_JPRB,&
	1017	& 0.101426E+00_JPRB, 0.007700E+00_JPRB, 0.441320E+00_JPRB,&
	1018	& 0.237571E+00_JPRB /)
	1019	REXPO3(1, 2, :) =(/&
	1020	& 0.100022E+02_JPRB, 0.851159E+02_JPRB, 0.346737E+03_JPRB,&
	1021	& 0.158501E+02_JPRB, 0.724223E+01_JPRB, 0.177828E+03_JPRB,&
	1022	& 0.467708E+02_JPRB /)
	1023
	1024	!* 2nd spectral interval: U.V. (0.25 - 0.44 Micron)
	1025
	1026	ZSUN6(2) = 0.135708_JPRB
	1027
	1028	ZD6(2,:)= (/ 1.000000000_JPRB, 1.000000000_JPRB, 0.000000000_JPRB /)
	1029
	1030	ZAPAD6(2, 1, :) = (/&
	1031	& 0.000000000E-00_JPRB,&
	1032	& 0.000000000E-00_JPRB,&
	1033	& 0.000000000E-00_JPRB,&
	1034	& 0.000000000E-00_JPRB,&
	1035	& 0.000000000E-00_JPRB,&
	1036	& 0.000000000E-00_JPRB,&
	1037	& 0.000000000E-00_JPRB/)
	1038	ZAPAD6(2, 2, :) = (/&
	1039	& 0.000000000E-00_JPRB,&
	1040	& 0.000000000E-00_JPRB,&
	1041	& 0.000000000E-00_JPRB,&
	1042	& 0.000000000E-00_JPRB,&
	1043	& 0.000000000E-00_JPRB,&
	1044	& 0.000000000E-00_JPRB,&
	1045	& 0.000000000E-00_JPRB/)
	1046	ZAPAD6(2, 3, :) = (/&
	1047	& 0.000000000E-00_JPRB,&
	1048	& 0.000000000E-00_JPRB,&
	1049	& 0.000000000E-00_JPRB,&
	1050	& 0.000000000E-00_JPRB,&
	1051	& 0.000000000E-00_JPRB,&
	1052	& 0.000000000E-00_JPRB,&
	1053	& 0.000000000E-00_JPRB/)
	1054
	1055	ZBPAD6(2, 1, :) = (/&
	1056	& 0.100000000E+01_JPRB,&
	1057	& 0.000000000E-00_JPRB,&
	1058	& 0.000000000E-00_JPRB,&
	1059	& 0.000000000E-00_JPRB,&
	1060	& 0.000000000E-00_JPRB,&
	1061	& 0.000000000E-00_JPRB,&
	1062	& 0.000000000E-00_JPRB/)
	1063	ZBPAD6(2, 2, :) = (/&
	1064	& 0.100000000E+01_JPRB,&
	1065	& 0.000000000E-00_JPRB,&
	1066	& 0.000000000E-00_JPRB,&
	1067	& 0.000000000E-00_JPRB,&
	1068	& 0.000000000E-00_JPRB,&
	1069	& 0.000000000E-00_JPRB,&
	1070	& 0.000000000E-00_JPRB/)
	1071	ZBPAD6(2, 3, :) = (/&
	1072	& 0.100000000E+01_JPRB,&
	1073	& 0.000000000E-00_JPRB,&
	1074	& 0.000000000E-00_JPRB,&
	1075	& 0.000000000E-00_JPRB,&
	1076	& 0.000000000E-00_JPRB,&
	1077	& 0.000000000E-00_JPRB,&
	1078	& 0.000000000E-00_JPRB/)
	1079
	1080	ZRAY6(2,:)= (/&
	1081	& -.895408E-03_JPRB, .283130E+01_JPRB,-.867917E+01_JPRB,&
	1082	& .148672E+02_JPRB,-.127282E+02_JPRB, .422212E+01_JPRB/)
	1083	ZRBY6(2,:)= (/&
	1084	& .55503070E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	1085	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	1086
	1087	ZSWCE6(2) = 0.0_JPRB
	1088	ZSWCP6(2) = 0.0_JPRB
	1089
	1090	NEXPO3(2) = 7
	1091	REXPO3(2, 1, :) =(/&
	1092	& 0.043801E+00_JPRB, 0.078893E+00_JPRB, 0.036839E+00_JPRB,&
	1093	& 0.022503E+00_JPRB, 0.042333E+00_JPRB, 0.037870E+00_JPRB,&
	1094	& 0.737762E+00_JPRB /)
	1095	REXPO3(2, 2, :) =(/&
	1096	& 0.234249E+01_JPRB, 0.125170E+00_JPRB, 0.549527E+02_JPRB,&
	1097	& 0.257041E+03_JPRB, 0.476838E+00_JPRB, 0.911993E+01_JPRB,&
	1098	& 0.000000E+00_JPRB /)
	1099
	1100	!* 3rd spectral interval: Visible (0.44 - 0.69 Micron)
	1101
	1102	ZSUN6(3) = 0.322135_JPRB
	1103
	1104	ZD6(3,:)= (/ 0.800000000_JPRB, 0.900000000_JPRB, 0.000000000_JPRB /)
	1105
	1106	ZAPAD6(3, 1, :) = (/&
	1107	& 0.1762097E+03_JPRB,&
	1108	& 0.1641762E+03_JPRB,&
	1109	& 0.8687919E+02_JPRB,&
	1110	& 0.0000000E-00_JPRB,&
	1111	& 0.0000000E-00_JPRB,&
	1112	& 0.0000000E-00_JPRB,&
	1113	& 0.0000000E-00_JPRB/)
	1114	ZAPAD6(3, 2, :) = (/&
	1115	& 0.5581224E+00_JPRB,&
	1116	& 0.1748430E+03_JPRB,&
	1117	& 0.1134123E+04_JPRB,&
	1118	& 0.3490429E+03_JPRB,&
	1119	& 0.0000000E-00_JPRB,&
	1120	& 0.0000000E-00_JPRB,&
	1121	& 0.0000000E-00_JPRB/)
	1122	ZAPAD6(3, 3, :) = (/&
	1123	& 0.000000000E-00_JPRB,&
	1124	& 0.000000000E-00_JPRB,&
	1125	& 0.000000000E-00_JPRB,&
	1126	& 0.000000000E-00_JPRB,&
	1127	& 0.000000000E-00_JPRB,&
	1128	& 0.000000000E-00_JPRB,&
	1129	& 0.000000000E-00_JPRB/)
	1130
	1131	ZBPAD6(3, 1, :) = (/&
	1132	& 0.1762097E+03_JPRB,&
	1133	& 0.1663950E+03_JPRB,&
	1134	& 0.8939724E+02_JPRB,&
	1135	& 0.1000000E+01_JPRB,&
	1136	& 0.0000000E-00_JPRB,&
	1137	& 0.0000000E-00_JPRB,&
	1138	& 0.0000000E-00_JPRB/)
	1139	ZBPAD6(3, 2, :) = (/&
	1140	& 0.5581224E+00_JPRB,&
	1141	& 0.1749251E+03_JPRB,&
	1142	& 0.1159910E+04_JPRB,&
	1143	& 0.3893268E+03_JPRB,&
	1144	& 0.1000000E+01_JPRB,&
	1145	& 0.0000000E-00_JPRB,&
	1146	& 0.0000000E-00_JPRB/)
	1147	ZBPAD6(3, 3, :) = (/&
	1148	& 0.100000000E+01_JPRB,&
	1149	& 0.000000000E-00_JPRB,&
	1150	& 0.000000000E-00_JPRB,&
	1151	& 0.000000000E-00_JPRB,&
	1152	& 0.000000000E-00_JPRB,&
	1153	& 0.000000000E-00_JPRB,&
	1154	& 0.000000000E-00_JPRB/)
	1155
	1156	ZRAY6(3,:)= (/&
	1157	& .474852E-01_JPRB, .605110E+00_JPRB,-.248070E+01_JPRB,&
	1158	& .485580E+01_JPRB,-.444043E+01_JPRB, .152369E+01_JPRB/)
	1159	ZRBY6(3,:)= (/&
	1160	& .10528199E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	1161	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	1162
	1163	ZSWCE6(3) = 0.0_JPRB
	1164	ZSWCP6(3) = 0.0_JPRB
	1165
	1166	NEXPO3(3) = 6
	1167	REXPO3(3, 1, :) =(/&
	1168	& 0.063442E+00_JPRB, 0.058550E+00_JPRB, 0.237534E+00_JPRB,&
	1169	& 0.412292E+00_JPRB, 0.126141E+00_JPRB, 0.102041E+00_JPRB,&
	1170	& 0.000000E+00_JPRB /)
	1171	REXPO3(3, 2, :) =(/&
	1172	& 0.125170E+00_JPRB, 0.119209E-01_JPRB, 0.119209E+00_JPRB,&
	1173	& 0.417233E-01_JPRB, 0.894070E-01_JPRB, 0.000000E+00_JPRB,&
	1174	& 0.000000E+00_JPRB /)
	1175
	1176	! ----------------------------------------------------------------
	1177
	1178	!* Near-Infrared (0.69 - 4.0 Microns) is sub-divided into:
	1179
	1180	! ----------------------------------------------------------------
	1181
	1182	!* 0.69 - 1.19 Micron
	1183
	1184	!* UMG is O2 only
	1185
	1186	ZSUN6(4) = 0.326158_JPRB
	1187
	1188	ZD6(4,:)= (/ 0.000000000_JPRB, 0.900000000_JPRB, 1.000000000_JPRB /)
	1189
	1190	ZAPAD6(4, 1, :) = (/&
	1191	& 0.1335726E+02_JPRB,&
	1192	& 0.2939136E+04_JPRB,&
	1193	& 0.4010585E+05_JPRB,&
	1194	& 0.7195030E+05_JPRB,&
	1195	& 0.1648338E+05_JPRB,&
	1196	& 0.3373738E+03_JPRB,&
	1197	& 0.0000000E+00_JPRB/)
	1198	ZAPAD6(4, 2, :) = (/&
	1199	& 0.2001271E-01_JPRB,&
	1200	& 0.2480831E+01_JPRB,&
	1201	& 0.3444162E+02_JPRB,&
	1202	& 0.4788946E+02_JPRB,&
	1203	& 0.0000000E+00_JPRB,&
	1204	& 0.0000000E+00_JPRB,&
	1205	& 0.0000000E+00_JPRB/)
	1206	ZAPAD6(4, 3, :) = (/&
	1207	& 0.000000000E+00_JPRB,&
	1208	& 0.000000000E+00_JPRB,&
	1209	& 0.000000000E+00_JPRB,&
	1210	& 0.000000000E+00_JPRB,&
	1211	& 0.000000000E+00_JPRB,&
	1212	& 0.000000000E+00_JPRB,&
	1213	& 0.000000000E+00_JPRB/)
	1214
	1215	ZBPAD6(4, 1, :) = (/&
	1216	& 0.1335726E+02_JPRB,&
	1217	& 0.2942327E+04_JPRB,&
	1218	& 0.4077237E+05_JPRB,&
	1219	& 0.7749017E+05_JPRB,&
	1220	& 0.2123132E+05_JPRB,&
	1221	& 0.6659687E+03_JPRB,&
	1222	& 0.1000000E+01_JPRB/)
	1223	ZBPAD6(4, 2, :) = (/&
	1224	& 0.2001271E-01_JPRB,&
	1225	& 0.2549067E+01_JPRB,&
	1226	& 0.3752433E+02_JPRB,&
	1227	& 0.6276637E+02_JPRB,&
	1228	& 0.1000000E+01_JPRB,&
	1229	& 0.0000000E+00_JPRB,&
	1230	& 0.0000000E+00_JPRB/)
	1231	ZBPAD6(4, 3, :) = (/&
	1232	& 1.000000000E+00_JPRB,&
	1233	& 0.000000000E+00_JPRB,&
	1234	& 0.000000000E+00_JPRB,&
	1235	& 0.000000000E+00_JPRB,&
	1236	& 0.000000000E+00_JPRB,&
	1237	& 0.000000000E+00_JPRB,&
	1238	& 0.000000000E+00_JPRB/)
	1239
	1240	ZRAY6(4,:)= (/&
	1241	& .519710E-01_JPRB,-.150667E+00_JPRB, .154170E+01_JPRB,&
	1242	& -.518901E+01_JPRB, .668011E+01_JPRB,-.294843E+01_JPRB/)
	1243	ZRBY6(4,:)= (/&
	1244	& .16436996E-01_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	1245	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	1246
	1247	ZSWCE6(4) = 0.0_JPRB
	1248	ZSWCP6(4) = 0.0_JPRB
	1249
	1250	NEXPO3(4) = 4
	1251	REXPO3(4, 1, :) =(/&
	1252	& 0.000074E+00_JPRB, 0.320194E+00_JPRB, 0.082915E+00_JPRB,&
	1253	& 0.596816E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1254	& 0.000000E+00_JPRB /)
	1255	REXPO3(4, 2, :) =(/&
	1256	& 0.232458E+00_JPRB, 0.119209E-01_JPRB, 0.178814E-01_JPRB,&
	1257	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1258	& 0.000000E+00_JPRB /)
	1259
	1260	! ----------------------------------------------------------------
	1261
	1262	!* 1.19 - 2.38 Microns
	1263
	1264	!* UMG is CO2 only
	1265
	1266	ZSUN6(5) = 0.180608_JPRB
	1267
	1268	ZD6(5,:)= (/ 0.000000000_JPRB, 0.800000000_JPRB, 1.000000000_JPRB /)
	1269
	1270	ZAPAD6(5, 1, :) = (/&
	1271	& 0.3325841E-03_JPRB,&
	1272	& 0.6194496E+00_JPRB,&
	1273	& 0.1497138E+03_JPRB,&
	1274	& 0.2314864E+04_JPRB,&
	1275	& 0.2380109E+04_JPRB,&
	1276	& 0.9553823E+02_JPRB,&
	1277	& 0.0000000E+00_JPRB/)
	1278	ZAPAD6(5, 2, :) = (/&
	1279	& 0.4552471E-03_JPRB,&
	1280	& 0.4084154E+00_JPRB,&
	1281	& 0.6114905E+01_JPRB,&
	1282	& 0.7102540E+01_JPRB,&
	1283	& 0.0000000E+00_JPRB,&
	1284	& 0.0000000E+00_JPRB,&
	1285	& 0.0000000E+00_JPRB/)
	1286	ZAPAD6(5, 3, :) = (/&
	1287	& 0.000000000E+00_JPRB,&
	1288	& 0.000000000E+00_JPRB,&
	1289	& 0.000000000E+00_JPRB,&
	1290	& 0.000000000E+00_JPRB,&
	1291	& 0.000000000E+00_JPRB,&
	1292	& 0.000000000E+00_JPRB,&
	1293	& 0.000000000E+00_JPRB/)
	1294
	1295	ZBPAD6(5, 1, :) = (/&
	1296	& 0.3325841E-03_JPRB,&
	1297	& 0.6231947E+00_JPRB,&
	1298	& 0.1553098E+03_JPRB,&
	1299	& 0.2822458E+04_JPRB,&
	1300	& 0.3885194E+04_JPRB,&
	1301	& 0.2700235E+03_JPRB,&
	1302	& 0.1000000E+01_JPRB/)
	1303	ZBPAD6(5, 2, :) = (/&
	1304	& 0.4552471E-03_JPRB,&
	1305	& 0.4088242E+00_JPRB,&
	1306	& 0.6411905E+01_JPRB,&
	1307	& 0.9444439E+01_JPRB,&
	1308	& 0.1000000E+01_JPRB,&
	1309	& 0.0000000E+00_JPRB,&
	1310	& 0.0000000E+00_JPRB/)
	1311	ZBPAD6(5, 3, :) = (/&
	1312	& 1.000000000E+00_JPRB,&
	1313	& 0.000000000E+00_JPRB,&
	1314	& 0.000000000E+00_JPRB,&
	1315	& 0.000000000E+00_JPRB,&
	1316	& 0.000000000E+00_JPRB,&
	1317	& 0.000000000E+00_JPRB,&
	1318	& 0.000000000E+00_JPRB/)
	1319
	1320	ZRAY6(5,:)= (/&
	1321	& .288474E-02_JPRB, .133984E-02_JPRB,-.152812E-01_JPRB,&
	1322	& .540845E-01_JPRB,-.747444E-01_JPRB, .364285E-01_JPRB/)
	1323	ZRBY6(5,:)= (/&
	1324	& .18073079E-02_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	1325	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	1326
	1327	ZSWCE6(5) = 0.0_JPRB
	1328	ZSWCP6(5) = 0.0_JPRB
	1329
	1330	NEXPO3(5) = 0
	1331	REXPO3(5, 1, :) =(/&
	1332	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1333	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1334	& 0.000000E+00_JPRB /)
	1335	REXPO3(5, 2, :) =(/&
	1336	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1337	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1338	& 0.000000E+00_JPRB /)
	1339
	1340	! ----------------------------------------------------------------
	1341
	1342	!* 2.38 - 4.00 Microns
	1343
	1344	ZSUN6(6) = 0.033474_JPRB
	1345
	1346	ZD6(6,:)= (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.000000000_JPRB /)
	1347
	1348	ZAPAD6(6, 1, :) = (/&
	1349	& 0.2122889E-06_JPRB,&
	1350	& 0.9030576E-03_JPRB,&
	1351	& 0.2431282E+00_JPRB,&
	1352	& 0.4901345E+01_JPRB,&
	1353	& 0.3996347E+01_JPRB,&
	1354	& 0.3910227E+01_JPRB,&
	1355	& 0.0000000E+00_JPRB/)
	1356	ZAPAD6(6, 2, :) = (/&
	1357	& 0.1215163E-03_JPRB,&
	1358	& 0.1222574E+00_JPRB,&
	1359	& 0.9382420E+01_JPRB,&
	1360	& 0.6875727E+02_JPRB,&
	1361	& 0.2746421E+02_JPRB,&
	1362	& 0.0000000E+00_JPRB,&
	1363	& 0.0000000E+00_JPRB/)
	1364	ZAPAD6(6, 3, :) = (/&
	1365	& 0.263068898E+02_JPRB,&
	1366	& 0.146425875E+03_JPRB,&
	1367	& 0.860137809E+02_JPRB,&
	1368	& 0.000000000E+00_JPRB,&
	1369	& 0.000000000E+00_JPRB,&
	1370	& 0.000000000E+00_JPRB,&
	1371	& 0.000000000E+00_JPRB/)
	1372
	1373	ZBPAD6(6, 1, :) = (/&
	1374	& 0.2122889E-06_JPRB,&
	1375	& 0.9379083E-03_JPRB,&
	1376	& 0.2957335E+00_JPRB,&
	1377	& 0.8747190E+01_JPRB,&
	1378	& 0.1015794E+02_JPRB,&
	1379	& 0.1361277E+02_JPRB,&
	1380	& 0.1000000E+01_JPRB/)
	1381	ZBPAD6(6, 2, :) = (/&
	1382	& 0.1215163E-03_JPRB,&
	1383	& 0.1255648E+00_JPRB,&
	1384	& 0.1060119E+02_JPRB,&
	1385	& 0.8414439E+02_JPRB,&
	1386	& 0.4299438E+02_JPRB,&
	1387	& 0.1000000E+01_JPRB,&
	1388	& 0.0000000E+00_JPRB/)
	1389	ZBPAD6(6, 3, :) = (/&
	1390	& 0.263068898E+02_JPRB,&
	1391	& 0.152569217E+03_JPRB,&
	1392	& 0.976791971E+02_JPRB,&
	1393	& 0.100000000E+01_JPRB,&
	1394	& 0.000000000E+00_JPRB,&
	1395	& 0.000000000E+00_JPRB,&
	1396	& 0.000000000E+00_JPRB/)
	1397
	1398	ZRAY6(6,:)= (/&
	1399	& .493545E-03_JPRB, .355696E-03_JPRB,-.284543E-02_JPRB,&
	1400	& .741248E-02_JPRB,-.773691E-02_JPRB, .370679E-02_JPRB/)
	1401	ZRBY6(6,:)= (/&
	1402	& .13618247E-03_JPRB, .000000E+00_JPRB, .000000E+00_JPRB,&
	1403	& .000000E+00_JPRB, .000000E+00_JPRB, .000000E+00_JPRB/)
	1404
	1405	ZSWCE6(6) = 0.0_JPRB
	1406	ZSWCP6(6) = 0.0_JPRB
	1407
	1408	NEXPO3(6) = 0
	1409	REXPO3(6, 1, :) =(/&
	1410	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1411	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1412	& 0.000000E+00_JPRB /)
	1413	REXPO3(6, 2, :) =(/&
	1414	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1415	& 0.000000E+00_JPRB, 0.000000E+00_JPRB, 0.000000E+00_JPRB,&
	1416	& 0.000000E+00_JPRB /)
	1417
	1418	!=====================================================================
	1419	!* 2.4 OPTICAL PARAMETERS FOR RAIN DROPS
	1420	! Savijarvi et al. (1996)
	1421
	1422	ZRTO1 = 0.003_JPRB
	1423	ZRTO2 = -0.22_JPRB
	1424	! CAUTION JUST TEMPORARY PARAMETERS
	1425	ZROMA2(1:2)= (/ 0.00008_JPRB , 0.264_JPRB /)
	1426	ZROMB2(1:2)= (/ 0.23_JPRB , 0.09_JPRB /)
	1427	ZRASY2(1:2)= (/ 0.88_JPRB , 0.94_JPRB /)
	1428
	1429	ZRA2(1:2)= (/ 1.5_JPRB , 1.5_JPRB /)
	1430	ZRB2(1:2)= (/ 0.50_JPRB , 1.13_JPRB /)
	1431	ZRC2(1:2)= (/ 5.58E-7_JPRB , 8.55E-4_JPRB/)
	1432	ZRD2(1:2)= (/ 1.25E-7_JPRB , 1.28E-3_JPRB/)
	1433	ZRE2(1:2)= (/ 0.841_JPRB , 0.786_JPRB /)
	1434	ZRF2(1:2)= (/ 2.08E-3_JPRB , 5.32E-3_JPRB/)
	1435
	1436	ZROMA4(1:4)= (/ 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB , 0.465_JPRB /)
	1437	ZROMB4(1:4)= (/ 0.23_JPRB , 0.22_JPRB , 0.09_JPRB , 0.001_JPRB /)
	1438	ZRASY4(1:4)= (/ 0.88_JPRB , 0.89_JPRB , 0.94_JPRB , 0.97_JPRB /)
	1439
	1440	ZRA4(1:4)= (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /)
	1441	ZRB4(1:4)= (/ 0.50_JPRB , 0.78_JPRB , 1.13_JPRB , 2.00_JPRB /)
	1442	ZRC4(1:4)= (/ 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /)
	1443	ZRD4(1:4)= (/ 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /)
	1444	ZRE4(1:4)= (/ 0.841_JPRB , 0.821_JPRB , 0.786_JPRB , 0.820_JPRB /)
	1445	ZRF4(1:4)= (/ 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /)
	1446
	1447	ZROMA6(1:6)= (/ 0.00008_JPRB , 0.00008_JPRB , 0.00008_JPRB , 0.0105_JPRB , 0.264_JPRB , 0.465_JPRB /)
	1448	ZROMB6(1:6)= (/ 0.23_JPRB , 0.23_JPRB , 0.23_JPRB , 0.22_JPRB , 0.09_JPRB , 0.001_JPRB /)
	1449	ZRASY6(1:6)= (/ 0.88_JPRB , 0.88_JPRB , 0.88_JPRB , 0.89_JPRB , 0.94_JPRB , 0.97_JPRB /)
	1450
	1451	ZRA6(1:6)= (/ 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB , 1.5_JPRB /)
	1452	ZRB6(1:6)= (/ 0.50_JPRB , 0.50_JPRB , 0.50_JPRB , 0.78_JPRB , 1.13_JPRB , 2.00_JPRB /)
	1453	ZRC6(1:6)= (/ 5.58E-7_JPRB , 5.58E-7_JPRB , 5.58E-7_JPRB , 2.18E-5_JPRB, 8.55E-4_JPRB, 1.94E-1_JPRB /)
	1454	ZRD6(1:6)= (/ 1.25E-7_JPRB , 1.25E-7_JPRB , 1.25E-7_JPRB , 2.25E-5_JPRB, 1.28E-3_JPRB, 8.04E-3_JPRB /)
	1455	ZRE6(1:6)= (/ 0.841_JPRB , 0.841_JPRB , 0.841_JPRB , 0.821_JPRB , 0.786_JPRB , 0.820_JPRB /)
	1456	ZRF6(1:6)= (/ 2.08E-3_JPRB , 2.08E-3_JPRB , 2.08E-3_JPRB , 3.06E-3_JPRB, 5.32E-3_JPRB, 5.59E-3_JPRB /)
	1457
	1458	!=====================================================================
	1459	! 3. FOURTEEN SPECTRAL INTERVALS FOR RRTM_SW
	1460
	1461	! SW : 14 spectral intervals : mapping to SW6 spectral intervals
	1462	! 0.185-0.25-0.44-0.69-1.19-2.38-4.00
	1463	! 3.846 - 3.077 6
	1464	! 3.077 - 2.500 6
	1465	! 2.500 - 2.150 5
	1466	! 2.150 - 1.942 5
	1467	! 1.942 - 1.626 5
	1468	! 1.626 - 1.299 5
	1469	! 1.299 - 1.242 5
	1470	! 1.242 - 0.7782 4
	1471	! 0.7782- 0.6250 4
	1472	! 0.6250- 0.4415 3
	1473	! 0.4415- 0.3448 2
	1474	! 0.3448- 0.2632 2
	1475	! 0.2632- 0.2000 1
	1476	! 12.195 - 3.846 6
	1477
	1478	IMPSRTM(1:14)= (/ 6,6,5,5,5,5,5,4,4,3,2,2,1,6 /)
	1479	!=====================================================================
	1480
	1481	!* 4. SET VALUES.
	1482	! -----------
	1483
	1484	IF (LONEWSW) THEN
	1485	ZPDH2O = ZPDHNIS
	1486	ZPDUMG = ZPDUNIS
	1487	ZPRH2O = ZPRHNIS
	1488	ZPRUMG = ZPRUNIS
	1489	RTDH2O = ZTDHNIS
	1490	RTDUMG = ZTDUNIS
	1491	RTH2O = ZTHNIS
	1492	RTUMG = ZTUNIS
	1493	ENDIF
	1494
	1495	RPDH1=ZPDH2O+1.0_JPRB
	1496	RPDU1=ZPDUMG+1.0_JPRB
	1497	ZH2O=1.0_JPRB/( 10._JPRB* RG * RPDH1 )
	1498	ZUMG=1.0_JPRB/( 10._JPRB* RG * RPDU1 )
	1499	RPNU = ZUMG/(ZPRUMG**ZPDUMG)
	1500	RPNH = ZH2O/(ZPRH2O**ZPDH2O)
	1501
	1502	IF(.NOT.ALLOCATED(RSUN)) ALLOCATE (RSUN(NTSW))
	1503
	1504	DO JW=1,ISW6
	1505	RSUN(JW)=0.0_JPRB
	1506	RSWCE(JW)=0.0_JPRB
	1507	RSWCP(JW)=0.0_JPRB
	1508
	1509	DO JC3=1,3
	1510	D(JW,JC3)=0.0_JPRB
	1511	ENDDO
	1512	DO JC6=1,6
	1513	RRAY(JW,JC6)=0.0_JPRB
	1514	ENDDO
	1515	DO JI=1,3
	1516	DO JJ=1,7
	1517	APAD(JW,JI,JJ)=0.0_JPRB
	1518	BPAD(JW,JI,JJ)=0.0_JPRB
	1519	ENDDO
	1520	ENDDO
	1521	ENDDO
	1522
	1523	IF (KSW == 2) THEN
	1524	ZPDH2O = ZPDH2IS
	1525	ZPDUMG = ZPDU2IS
	1526	ZPRH2O = ZPRH2IS
	1527	ZPRUMG = ZPRU2IS
	1528	RTDH2O = ZTDH2IS
	1529	RTDUMG = ZTDU2IS
	1530	RTH2O = ZTH2IS
	1531	RTUMG = ZTU2IS
	1532
	1533	RADJUST=ZADJUS2
	1534	RHSRTA=ZRTO1
	1535	RHSRTB=ZRTO2
	1536	DO JW=1,KSW
	1537	RSUN (JW)=ZSUN2(JW)
	1538
	1539	RSWCE(JW)=ZSWCE2(JW)
	1540	RSWCP(JW)=ZSWCP2(JW)
	1541
	1542	RROMA(JW)=ZROMA2(JW)
	1543	RROMB(JW)=ZROMB2(JW)
	1544	RRASY(JW)=ZRASY2(JW)
	1545	RHSRA(JW)=ZRA2(JW)
	1546	RHSRB(JW)=ZRB2(JW)
	1547	RHSRC(JW)=ZRC2(JW)
	1548	RHSRD(JW)=ZRD2(JW)
	1549	RHSRE(JW)=ZRE2(JW)
	1550	RHSRF(JW)=ZRF2(JW)
	1551
	1552	DO JC3=1,3
	1553	D(JW,JC3)=ZD2(JW,JC3)
	1554	ENDDO
	1555	DO JC6=1,6
	1556	RRAY(JW,JC6)=ZRAY2(JW,JC6)
	1557	ENDDO
	1558	DO JI=1,3
	1559	DO JJ=1,7
	1560	APAD(JW,JI,JJ)=ZAPAD2(JW,JI,JJ)
	1561	BPAD(JW,JI,JJ)=ZBPAD2(JW,JI,JJ)
	1562	ENDDO
	1563	ENDDO
	1564	ENDDO
	1565
	1566	ELSEIF (KSW == 4) THEN
	1567	ZPDH2O = ZPDH4IS
	1568	ZPDUMG = ZPDU4IS
	1569	ZPRH2O = ZPRH4IS
	1570	ZPRUMG = ZPRU4IS
	1571	RTDH2O = ZTDH4IS
	1572	RTDUMG = ZTDU4IS
	1573	RTH2O = ZTH4IS
	1574	RTUMG = ZTU4IS
	1575
	1576	RADJUST=ZADJUS4
	1577	RHSRTA=ZRTO1
	1578	RHSRTB=ZRTO2
	1579	DO JW=1,KSW
	1580	RSUN (JW)=ZSUN4(JW)
	1581
	1582	RSWCE(JW)=ZSWCE4(JW)
	1583	RSWCP(JW)=ZSWCP4(JW)
	1584
	1585	RROMA(JW)=ZROMA4(JW)
	1586	RROMB(JW)=ZROMB4(JW)
	1587	RRASY(JW)=ZRASY4(JW)
	1588	RHSRA(JW)=ZRA4(JW)
	1589	RHSRB(JW)=ZRB4(JW)
	1590	RHSRC(JW)=ZRC4(JW)
	1591	RHSRD(JW)=ZRD4(JW)
	1592	RHSRE(JW)=ZRE4(JW)
	1593	RHSRF(JW)=ZRF4(JW)
	1594
	1595	DO JC3=1,3
	1596	D(JW,JC3)=ZD4(JW,JC3)
	1597	ENDDO
	1598	DO JC6=1,6
	1599	RRAY(JW,JC6)=ZRAY4(JW,JC6)
	1600	ENDDO
	1601	DO JI=1,3
	1602	DO JJ=1,7
	1603	APAD(JW,JI,JJ)=ZAPAD4(JW,JI,JJ)
	1604	BPAD(JW,JI,JJ)=ZBPAD4(JW,JI,JJ)
	1605	ENDDO
	1606	ENDDO
	1607	ENDDO
	1608
	1609	ELSEIF (KSW == 6) THEN
	1610	ZPDH2O = ZPDH6IS
	1611	ZPDUMG = ZPDU6IS
	1612	ZPRH2O = ZPRH6IS
	1613	ZPRUMG = ZPRU6IS
	1614	RTDH2O = ZTDH6IS
	1615	RTDUMG = ZTDU6IS
	1616	RTH2O = ZTH6IS
	1617	RTUMG = ZTU6IS
	1618
	1619	RADJUST=ZADJUS6
	1620	RHSRTA=ZRTO1
	1621	RHSRTB=ZRTO2
	1622	DO JW=1,KSW
	1623	RSUN (JW)=ZSUN6(JW)
	1624
	1625	RSWCE(JW)=ZSWCE6(JW)
	1626	RSWCP(JW)=ZSWCP6(JW)
	1627
	1628	RROMA(JW)=ZROMA6(JW)
	1629	RROMB(JW)=ZROMB6(JW)
	1630	RRASY(JW)=ZRASY6(JW)
	1631	RHSRA(JW)=ZRA6(JW)
	1632	RHSRB(JW)=ZRB6(JW)
	1633	RHSRC(JW)=ZRC6(JW)
	1634	RHSRD(JW)=ZRD6(JW)
	1635	RHSRE(JW)=ZRE6(JW)
	1636	RHSRF(JW)=ZRF6(JW)
	1637
	1638	DO JC3=1,3
	1639	D(JW,JC3)=ZD6(JW,JC3)
	1640	ENDDO
	1641	IF (LRAYL) THEN
	1642	DO JC6=1,6
	1643	RRAY(JW,JC6)=ZRAY6(JW,JC6)
	1644	ENDDO
	1645	ELSE
	1646	DO JC6=1,6
	1647	RRAY(JW,JC6)=ZRBY6(JW,JC6)
	1648	ENDDO
	1649	ENDIF
	1650	DO JI=1,3
	1651	DO JJ=1,7
	1652	APAD(JW,JI,JJ)=ZAPAD6(JW,JI,JJ)
	1653	BPAD(JW,JI,JJ)=ZBPAD6(JW,JI,JJ)
	1654	ENDDO
	1655	ENDDO
	1656	ENDDO
	1657
	1658	ELSEIF (KTSW == 14 .AND. KSW == 14) THEN
	1659	PRINT *,'SUSWN: 14-SPECTRAL INTERVALS --> RRTM_SW'
	1660
	1661	ELSE
	1662	PRINT *,'SUSWN: NTSW, NSW: ',KTSW,KSW
	1663	CALL ABOR1('SUSWN: WRONG NUMBER OF SW INTERVALS')
	1664	ENDIF
	1665
	1666	! ----------------------------------------------------------------
	1667
	1668	!* 5. PRINT VALUES.
	1669	! ------------
	1670
	1671	! PRINT 9700,RPDH1,RPDU1,RPNU,RPNH
	1672	!9700 FORMAT(1X,'RPDH1,RPDU1,RPNU,RPNH ',4F10.5)
	1673	! DO JW=1,KSW
	1674	! PRINT *,' SW INTERVAL ',JW
	1675	! PRINT 9701,JW,RSUN(JW),RSWCE(JW),RSWCP(JW)
	1676	!9701 FORMAT(1X,I3,F10.7,2F10.5)
	1677	! PRINT 9702,JW,(D(JW,JC3),JC3=1,3)
	1678	!9702 FORMAT(1X,I3,3F10.5)
	1679	! PRINT 9703,JW,(RRAY(JW,JC6),JC6=1,6)
	1680	!9703 FORMAT(1X,I3,6E15.8)
	1681	! PRINT 9704,JW,(JABS,(APAD(JW,JABS,JJ),JJ=1,7),JABS=1,3)
	1682	!9704 FORMAT(1X,I3,(I3,7E15.8))
	1683	! PRINT 9706,JW,RROMA(JW),RROMB(JW),RRASY(JW),RHSRTA,RHSRTB
	1684	!9706 FORMAT(1X,I3,' Rain Optics from Rain Rate: ',F8.5,F6.3,3F6.3)
	1685	! PRINT 9707,JW,RHSRA(JW),RHSRB(JW),RHSRC(JW),RHSRD(JW)
	1686	! & ,RHSRE(JW),RHSRF(JW)
	1687	!9707 FORMAT(1X,I3,' Rain Optics as function of Reff: ',6E11.4)
	1688	! ENDDO
	1689
	1690	! ----------------------------------------------------------------
	1691
	1692	IF (LHOOK) CALL DR_HOOK('SUSWN',1,ZHOOK_HANDLE)
	1693	END SUBROUTINE SUSWN

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

Download in other formats:

Original Format