Index: LMDZ4/trunk/libf/phylmd/radiation_AR4.F
===================================================================
--- LMDZ4/trunk/libf/phylmd/radiation_AR4.F	(revision 1107)
+++ LMDZ4/trunk/libf/phylmd/radiation_AR4.F	(revision 1279)
@@ -54,58 +54,62 @@
 C* ARGUMENTS:
 C
-      REAL*8 PSCT  ! constante solaire (valeur conseillee: 1370)
-cIM ctes ds clesphys.h   REAL*8 RCO2  ! concentration CO2 (IPCC: 353.E-06*44.011/28.97)
+      REAL(KIND=8) PSCT  ! constante solaire (valeur conseillee: 1370)
+cIM ctes ds clesphys.h   REAL(KIND=8) RCO2  ! concentration CO2 (IPCC: 353.E-06*44.011/28.97)
 #include "clesphys.h"
 C
-      REAL*8 PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
-      REAL*8 PDP(KDLON,KFLEV)    ! LAYER THICKNESS (PA)
-      REAL*8 PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
-C
-      REAL*8 PRMU0(KDLON)  ! COSINE OF ZENITHAL ANGLE
-      REAL*8 PFRAC(KDLON)  ! fraction de la journee
-C
-      REAL*8 PTAVE(KDLON,KFLEV)  ! LAYER TEMPERATURE (K)
-      REAL*8 PWV(KDLON,KFLEV)    ! SPECIFIC HUMIDITY (KG/KG)
-      REAL*8 PQS(KDLON,KFLEV)    ! SATURATED WATER VAPOUR (KG/KG)
-      REAL*8 POZON(KDLON,KFLEV)  ! OZONE CONCENTRATION (KG/KG)
-      REAL*8 PAER(KDLON,KFLEV,5) ! AEROSOLS' OPTICAL THICKNESS
-C
-      REAL*8 PALBD(KDLON,2)  ! albedo du sol (lumiere diffuse)
-      REAL*8 PALBP(KDLON,2)  ! albedo du sol (lumiere parallele)
-C
-      REAL*8 PCLDSW(KDLON,KFLEV)    ! CLOUD FRACTION
-      REAL*8 PTAU(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS
-      REAL*8 PCG(KDLON,2,KFLEV)     ! ASYMETRY FACTOR
-      REAL*8 POMEGA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
-C
-      REAL*8 PHEAT(KDLON,KFLEV) ! SHORTWAVE HEATING (K/DAY)
-      REAL*8 PHEAT0(KDLON,KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
-      REAL*8 PALBPLA(KDLON)     ! PLANETARY ALBEDO
-      REAL*8 PTOPSW(KDLON)      ! SHORTWAVE FLUX AT T.O.A.
-      REAL*8 PSOLSW(KDLON)      ! SHORTWAVE FLUX AT SURFACE
-      REAL*8 PTOPSW0(KDLON)     ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
-      REAL*8 PSOLSW0(KDLON)     ! SHORTWAVE FLUX AT SURFACE (CLEAR-SKY)
+      REAL(KIND=8) PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
+      REAL(KIND=8) PDP(KDLON,KFLEV)    ! LAYER THICKNESS (PA)
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
+C
+      REAL(KIND=8) PRMU0(KDLON)  ! COSINE OF ZENITHAL ANGLE
+      REAL(KIND=8) PFRAC(KDLON)  ! fraction de la journee
+C
+      REAL(KIND=8) PTAVE(KDLON,KFLEV)  ! LAYER TEMPERATURE (K)
+      REAL(KIND=8) PWV(KDLON,KFLEV)    ! SPECIFIC HUMIDITY (KG/KG)
+      REAL(KIND=8) PQS(KDLON,KFLEV)    ! SATURATED WATER VAPOUR (KG/KG)
+      REAL(KIND=8) POZON(KDLON,KFLEV)  ! OZONE CONCENTRATION (KG/KG)
+      REAL(KIND=8) PAER(KDLON,KFLEV,5) ! AEROSOLS' OPTICAL THICKNESS
+C
+      REAL(KIND=8) PALBD(KDLON,2)  ! albedo du sol (lumiere diffuse)
+      REAL(KIND=8) PALBP(KDLON,2)  ! albedo du sol (lumiere parallele)
+C
+      REAL(KIND=8) PCLDSW(KDLON,KFLEV)    ! CLOUD FRACTION
+      REAL(KIND=8) PTAU(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS
+      REAL(KIND=8) PCG(KDLON,2,KFLEV)     ! ASYMETRY FACTOR
+      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
+C
+      REAL(KIND=8) PHEAT(KDLON,KFLEV) ! SHORTWAVE HEATING (K/DAY)
+      REAL(KIND=8) PHEAT0(KDLON,KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
+      REAL(KIND=8) PALBPLA(KDLON)     ! PLANETARY ALBEDO
+      REAL(KIND=8) PTOPSW(KDLON)      ! SHORTWAVE FLUX AT T.O.A.
+      REAL(KIND=8) PSOLSW(KDLON)      ! SHORTWAVE FLUX AT SURFACE
+      REAL(KIND=8) PTOPSW0(KDLON)     ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
+      REAL(KIND=8) PSOLSW0(KDLON)     ! SHORTWAVE FLUX AT SURFACE (CLEAR-SKY)
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZOZ(KDLON,KFLEV)
-      REAL*8 ZAKI(KDLON,2)     
-      REAL*8 ZCLD(KDLON,KFLEV)
-      REAL*8 ZCLEAR(KDLON) 
-      REAL*8 ZDSIG(KDLON,KFLEV)
-      REAL*8 ZFACT(KDLON)
-      REAL*8 ZFD(KDLON,KFLEV+1)
-      REAL*8 ZFDOWN(KDLON,KFLEV+1)
-      REAL*8 ZFU(KDLON,KFLEV+1)
-      REAL*8 ZFUP(KDLON,KFLEV+1)
-      REAL*8 ZRMU(KDLON)
-      REAL*8 ZSEC(KDLON)
-      REAL*8 ZUD(KDLON,5,KFLEV+1)
-      REAL*8 ZCLDSW0(KDLON,KFLEV)
+      real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
+
+      REAL(KIND=8) ZOZ(KDLON,KFLEV)
+!     column-density of ozone in layer, in kilo-Dobsons
+
+      REAL(KIND=8) ZAKI(KDLON,2)     
+      REAL(KIND=8) ZCLD(KDLON,KFLEV)
+      REAL(KIND=8) ZCLEAR(KDLON) 
+      REAL(KIND=8) ZDSIG(KDLON,KFLEV)
+      REAL(KIND=8) ZFACT(KDLON)
+      REAL(KIND=8) ZFD(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFDOWN(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFU(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFUP(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRMU(KDLON)
+      REAL(KIND=8) ZSEC(KDLON)
+      REAL(KIND=8) ZUD(KDLON,5,KFLEV+1)
+      REAL(KIND=8) ZCLDSW0(KDLON,KFLEV)
 c
-      REAL*8 ZFSUP(KDLON,KFLEV+1)
-      REAL*8 ZFSDN(KDLON,KFLEV+1)
-      REAL*8 ZFSUP0(KDLON,KFLEV+1)
-      REAL*8 ZFSDN0(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFSUP(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFSDN(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFSUP0(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFSDN0(KDLON,KFLEV+1)
 C
       INTEGER inu, jl, jk, i, k, kpl1
@@ -122,23 +126,23 @@
 c$OMP THREADPRIVATE(itapsw)
 cjq-Introduced for aerosol forcings
-      real*8 flag_aer
+      real(kind=8) flag_aer
       logical ok_ade, ok_aie    ! use aerosol forcings or not?
-      real*8 tauae(kdlon,kflev,2)  ! aerosol optical properties
-      real*8 pizae(kdlon,kflev,2)  ! (see aeropt.F)
-      real*8 cgae(kdlon,kflev,2)   ! -"-
-      REAL*8 PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (pre-industrial value)
-      REAL*8 POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
-      REAL*8 PTOPSWAD(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
-      REAL*8 PSOLSWAD(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
-      REAL*8 PTOPSWAI(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND)
-      REAL*8 PSOLSWAI(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND)
+      real(kind=8) tauae(kdlon,kflev,2)  ! aerosol optical properties
+      real(kind=8) pizae(kdlon,kflev,2)  ! (see aeropt.F)
+      real(kind=8) cgae(kdlon,kflev,2)   ! -"-
+      REAL(KIND=8) PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (pre-industrial value)
+      REAL(KIND=8) POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
+      REAL(KIND=8) PTOPSWAD(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
+      REAL(KIND=8) PSOLSWAD(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
+      REAL(KIND=8) PTOPSWAI(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND)
+      REAL(KIND=8) PSOLSWAI(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND)
 cjq - Fluxes including aerosol effects
-      REAL*8,allocatable,save :: ZFSUPAD(:,:)
+      REAL(KIND=8),allocatable,save :: ZFSUPAD(:,:)
 c$OMP THREADPRIVATE(ZFSUPAD)
-      REAL*8,allocatable,save :: ZFSDNAD(:,:)
+      REAL(KIND=8),allocatable,save :: ZFSDNAD(:,:)
 c$OMP THREADPRIVATE(ZFSDNAD)
-      REAL*8,allocatable,save :: ZFSUPAI(:,:)
+      REAL(KIND=8),allocatable,save :: ZFSUPAI(:,:)
 c$OMP THREADPRIVATE(ZFSUPAI)
-      REAL*8,allocatable,save :: ZFSDNAI(:,:)
+      REAL(KIND=8),allocatable,save :: ZFSDNAI(:,:)
 c$OMP THREADPRIVATE(ZFSDNAI)
       logical initialized
@@ -151,4 +155,5 @@
 c$OMP THREADPRIVATE(initialized)
 cjq-end
+      REAL tmp_
       if(.not.initialized) then
         flag_aer=0.
@@ -158,9 +163,10 @@
         allocate(ZFSUPAI(KDLON,KFLEV+1))
         allocate(ZFSDNAI(KDLON,KFLEV+1))
-	ZFSUPAD(:,:)=0.
-        ZFSDNAD(:,:)=0.
-        ZFSUPAI(:,:)=0.
-        ZFSDNAI(:,:)=0.
-	
+        DO JK = 1 , KDLON*(KFLEV+1)
+          ZFSUPAD(JK,1) = 0.0     ! ZFSUPAD(:,:)=0.
+          ZFSDNAD(JK,1) = 0.0     ! ZFSDNAD(:,:)=0.
+          ZFSUPAI(JK,1) = 0.0     ! ZFSUPAI(:,:)=0.
+          ZFSDNAI(JK,1) = 0.0     ! ZFSDNAI(:,:)=0.
+        END DO
       endif
 !rv
@@ -175,10 +181,11 @@
       IF (MOD(itapsw,swpas).EQ.0) THEN
 c
+      tmp_ = 1./( dobson_u * 1e3 * RG)
+!cdir collapse
       DO JK = 1 , KFLEV
-      DO JL = 1, KDLON
-         ZCLDSW0(JL,JK) = 0.0
-         ZOZ(JL,JK) = POZON(JL,JK)*46.6968/RG
-     .               *PDP(JL,JK)*(101325.0/PPSOL(JL))
-      ENDDO
+        DO JL = 1, KDLON
+          ZCLDSW0(JL,JK) = 0.0
+          ZOZ(JL,JK) = POZON(JL,JK)*tmp_*PDP(JL,JK)
+        ENDDO
       ENDDO
 C
@@ -339,4 +346,6 @@
      S                PRMU,PSEC,PUD)
       USE dimphy
+      USE radiation_AR4_param, only :
+     S     ZPDH2O,ZPDUMG,ZPRH2O,ZPRUMG,RTDH2O,RTDUMG,RTH2O,RTUMG
       IMPLICIT none
 cym#include "dimensions.h"
@@ -349,59 +358,41 @@
 C* ARGUMENTS:
 C
-      REAL*8 PSCT
-cIM ctes ds clesphys.h   REAL*8 RCO2
+      REAL(KIND=8) PSCT
+cIM ctes ds clesphys.h   REAL(KIND=8) RCO2
 #include "clesphys.h"
-      REAL*8 PCLDSW(KDLON,KFLEV)
-      REAL*8 PPMB(KDLON,KFLEV+1)
-      REAL*8 PPSOL(KDLON)
-      REAL*8 PRMU0(KDLON)
-      REAL*8 PFRAC(KDLON)
-      REAL*8 PTAVE(KDLON,KFLEV)
-      REAL*8 PWV(KDLON,KFLEV)
-C
-      REAL*8 PAKI(KDLON,2)
-      REAL*8 PCLD(KDLON,KFLEV)
-      REAL*8 PCLEAR(KDLON)
-      REAL*8 PDSIG(KDLON,KFLEV)
-      REAL*8 PFACT(KDLON)
-      REAL*8 PRMU(KDLON)
-      REAL*8 PSEC(KDLON)
-      REAL*8 PUD(KDLON,5,KFLEV+1)
+      REAL(KIND=8) PCLDSW(KDLON,KFLEV)
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1)
+      REAL(KIND=8) PPSOL(KDLON)
+      REAL(KIND=8) PRMU0(KDLON)
+      REAL(KIND=8) PFRAC(KDLON)
+      REAL(KIND=8) PTAVE(KDLON,KFLEV)
+      REAL(KIND=8) PWV(KDLON,KFLEV)
+C
+      REAL(KIND=8) PAKI(KDLON,2)
+      REAL(KIND=8) PCLD(KDLON,KFLEV)
+      REAL(KIND=8) PCLEAR(KDLON)
+      REAL(KIND=8) PDSIG(KDLON,KFLEV)
+      REAL(KIND=8) PFACT(KDLON)
+      REAL(KIND=8) PRMU(KDLON)
+      REAL(KIND=8) PSEC(KDLON)
+      REAL(KIND=8) PUD(KDLON,5,KFLEV+1)
 C
 C* LOCAL VARIABLES:
 C
       INTEGER IIND(2)
-      REAL*8 ZC1J(KDLON,KFLEV+1)
-      REAL*8 ZCLEAR(KDLON)
-      REAL*8 ZCLOUD(KDLON)
-      REAL*8 ZN175(KDLON)
-      REAL*8 ZN190(KDLON)
-      REAL*8 ZO175(KDLON)
-      REAL*8 ZO190(KDLON)
-      REAL*8 ZSIGN(KDLON)
-      REAL*8 ZR(KDLON,2) 
-      REAL*8 ZSIGO(KDLON)
-      REAL*8 ZUD(KDLON,2)
-      REAL*8 ZRTH, ZRTU, ZWH2O, ZDSCO2, ZDSH2O, ZFPPW
+      REAL(KIND=8) ZC1J(KDLON,KFLEV+1)
+      REAL(KIND=8) ZCLEAR(KDLON)
+      REAL(KIND=8) ZCLOUD(KDLON)
+      REAL(KIND=8) ZN175(KDLON)
+      REAL(KIND=8) ZN190(KDLON)
+      REAL(KIND=8) ZO175(KDLON)
+      REAL(KIND=8) ZO190(KDLON)
+      REAL(KIND=8) ZSIGN(KDLON)
+      REAL(KIND=8) ZR(KDLON,2) 
+      REAL(KIND=8) ZSIGO(KDLON)
+      REAL(KIND=8) ZUD(KDLON,2)
+      REAL(KIND=8) ZRTH, ZRTU, ZWH2O, ZDSCO2, ZDSH2O, ZFPPW
       INTEGER jl, jk, jkp1, jkl, jklp1, ja
 C
-C* Prescribed Data:
-c
-      REAL*8 ZPDH2O,ZPDUMG
-      SAVE ZPDH2O,ZPDUMG
-c$OMP THREADPRIVATE(ZPDH2O,ZPDUMG)
-      REAL*8 ZPRH2O,ZPRUMG
-      SAVE ZPRH2O,ZPRUMG
-c$OMP THREADPRIVATE(ZPRH2O,ZPRUMG)
-      REAL*8 RTDH2O,RTDUMG
-      SAVE RTDH2O,RTDUMG
-c$OMP THREADPRIVATE(RTDH2O,RTDUMG)
-      REAL*8 RTH2O ,RTUMG
-      SAVE RTH2O ,RTUMG
-c$OMP THREADPRIVATE(RTH2O ,RTUMG)
-      DATA ZPDH2O,ZPDUMG / 0.8   , 0.75 /
-      DATA ZPRH2O,ZPRUMG / 30000., 30000. /
-      DATA RTDH2O,RTDUMG /  0.40  , 0.375 /
-      DATA RTH2O ,RTUMG  /  240.  , 240.  /
 C     ------------------------------------------------------------------
 C
@@ -534,4 +525,5 @@
      S  ,  PFD   , PFU)
       USE dimphy
+      USE radiation_AR4_param, only : RSUN, RRAY
       IMPLICIT none
 cym#include "dimensions.h"
@@ -573,25 +565,25 @@
       INTEGER KNU
 c-OB
-      real*8 flag_aer
-      real*8 tauae(kdlon,kflev,2)
-      real*8 pizae(kdlon,kflev,2)
-      real*8 cgae(kdlon,kflev,2)
-      REAL*8 PAER(KDLON,KFLEV,5)
-      REAL*8 PALBD(KDLON,2)
-      REAL*8 PALBP(KDLON,2)
-      REAL*8 PCG(KDLON,2,KFLEV)  
-      REAL*8 PCLD(KDLON,KFLEV)
-      REAL*8 PCLDSW(KDLON,KFLEV)
-      REAL*8 PCLEAR(KDLON)
-      REAL*8 PDSIG(KDLON,KFLEV)
-      REAL*8 POMEGA(KDLON,2,KFLEV)
-      REAL*8 POZ(KDLON,KFLEV)
-      REAL*8 PRMU(KDLON)
-      REAL*8 PSEC(KDLON)
-      REAL*8 PTAU(KDLON,2,KFLEV)
-      REAL*8 PUD(KDLON,5,KFLEV+1)
-C
-      REAL*8 PFD(KDLON,KFLEV+1)
-      REAL*8 PFU(KDLON,KFLEV+1)
+      real(kind=8) flag_aer
+      real(kind=8) tauae(kdlon,kflev,2)
+      real(kind=8) pizae(kdlon,kflev,2)
+      real(kind=8) cgae(kdlon,kflev,2)
+      REAL(KIND=8) PAER(KDLON,KFLEV,5)
+      REAL(KIND=8) PALBD(KDLON,2)
+      REAL(KIND=8) PALBP(KDLON,2)
+      REAL(KIND=8) PCG(KDLON,2,KFLEV)  
+      REAL(KIND=8) PCLD(KDLON,KFLEV)
+      REAL(KIND=8) PCLDSW(KDLON,KFLEV)
+      REAL(KIND=8) PCLEAR(KDLON)
+      REAL(KIND=8) PDSIG(KDLON,KFLEV)
+      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)
+      REAL(KIND=8) POZ(KDLON,KFLEV)
+      REAL(KIND=8) PRMU(KDLON)
+      REAL(KIND=8) PSEC(KDLON)
+      REAL(KIND=8) PTAU(KDLON,2,KFLEV)
+      REAL(KIND=8) PUD(KDLON,5,KFLEV+1)
+C
+      REAL(KIND=8) PFD(KDLON,KFLEV+1)
+      REAL(KIND=8) PFU(KDLON,KFLEV+1)
 C
 C* LOCAL VARIABLES:
@@ -599,42 +591,26 @@
       INTEGER IIND(4)
 C      
-      REAL*8 ZCGAZ(KDLON,KFLEV) 
-      REAL*8 ZDIFF(KDLON)
-      REAL*8 ZDIRF(KDLON)        
-      REAL*8 ZPIZAZ(KDLON,KFLEV)
-      REAL*8 ZRAYL(KDLON)
-      REAL*8 ZRAY1(KDLON,KFLEV+1)
-      REAL*8 ZRAY2(KDLON,KFLEV+1)
-      REAL*8 ZREFZ(KDLON,2,KFLEV+1)
-      REAL*8 ZRJ(KDLON,6,KFLEV+1)
-      REAL*8 ZRJ0(KDLON,6,KFLEV+1)
-      REAL*8 ZRK(KDLON,6,KFLEV+1)
-      REAL*8 ZRK0(KDLON,6,KFLEV+1)
-      REAL*8 ZRMUE(KDLON,KFLEV+1)
-      REAL*8 ZRMU0(KDLON,KFLEV+1)
-      REAL*8 ZR(KDLON,4)
-      REAL*8 ZTAUAZ(KDLON,KFLEV)
-      REAL*8 ZTRA1(KDLON,KFLEV+1)
-      REAL*8 ZTRA2(KDLON,KFLEV+1)
-      REAL*8 ZW(KDLON,4)
+      REAL(KIND=8) ZCGAZ(KDLON,KFLEV) 
+      REAL(KIND=8) ZDIFF(KDLON)
+      REAL(KIND=8) ZDIRF(KDLON)        
+      REAL(KIND=8) ZPIZAZ(KDLON,KFLEV)
+      REAL(KIND=8) ZRAYL(KDLON)
+      REAL(KIND=8) ZRAY1(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRAY2(KDLON,KFLEV+1)
+      REAL(KIND=8) ZREFZ(KDLON,2,KFLEV+1)
+      REAL(KIND=8) ZRJ(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRJ0(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRK(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRK0(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRMUE(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRMU0(KDLON,KFLEV+1)
+      REAL(KIND=8) ZR(KDLON,4)
+      REAL(KIND=8) ZTAUAZ(KDLON,KFLEV)
+      REAL(KIND=8) ZTRA1(KDLON,KFLEV+1)
+      REAL(KIND=8) ZTRA2(KDLON,KFLEV+1)
+      REAL(KIND=8) ZW(KDLON,4)
 C
       INTEGER jl, jk, k, jaj, ikm1, ikl
-c
-c Prescribed Data:
-c
-      REAL*8 RSUN(2)
-      SAVE RSUN
-c$OMP THREADPRIVATE(RSUN)
-      REAL*8 RRAY(2,6)
-      SAVE RRAY
-c$OMP THREADPRIVATE(RRAY)
-      DATA RSUN(1) / 0.441676 /
-      DATA RSUN(2) / 0.558324 /
-      DATA (RRAY(1,K),K=1,6) /
-     S .428937E-01, .890743E+00,-.288555E+01,
-     S .522744E+01,-.469173E+01, .161645E+01/
-      DATA (RRAY(2,K),K=1,6) /
-     S .697200E-02, .173297E-01,-.850903E-01,
-     S .248261E+00,-.302031E+00, .129662E+00/
+
 C     ------------------------------------------------------------------
 C
@@ -777,4 +753,5 @@
      S  ,  PFDOWN,PFUP                                            )
       USE dimphy
+      USE radiation_AR4_param, only : RSUN, RRAY
       IMPLICIT none
 cym#include "dimensions.h"
@@ -821,96 +798,81 @@
       INTEGER KNU
 c-OB
-      real*8 flag_aer
-      real*8 tauae(kdlon,kflev,2)
-      real*8 pizae(kdlon,kflev,2)
-      real*8 cgae(kdlon,kflev,2)
-      REAL*8 PAER(KDLON,KFLEV,5)
-      REAL*8 PAKI(KDLON,2)
-      REAL*8 PALBD(KDLON,2)
-      REAL*8 PALBP(KDLON,2)
-      REAL*8 PCG(KDLON,2,KFLEV)
-      REAL*8 PCLD(KDLON,KFLEV)
-      REAL*8 PCLDSW(KDLON,KFLEV)
-      REAL*8 PCLEAR(KDLON)
-      REAL*8 PDSIG(KDLON,KFLEV)
-      REAL*8 POMEGA(KDLON,2,KFLEV)
-      REAL*8 POZ(KDLON,KFLEV)
-      REAL*8 PQS(KDLON,KFLEV)
-      REAL*8 PRMU(KDLON)
-      REAL*8 PSEC(KDLON)
-      REAL*8 PTAU(KDLON,2,KFLEV)
-      REAL*8 PUD(KDLON,5,KFLEV+1)
-      REAL*8 PWV(KDLON,KFLEV)
-C
-      REAL*8 PFDOWN(KDLON,KFLEV+1)
-      REAL*8 PFUP(KDLON,KFLEV+1)
+      real(kind=8) flag_aer
+      real(kind=8) tauae(kdlon,kflev,2)
+      real(kind=8) pizae(kdlon,kflev,2)
+      real(kind=8) cgae(kdlon,kflev,2)
+      REAL(KIND=8) PAER(KDLON,KFLEV,5)
+      REAL(KIND=8) PAKI(KDLON,2)
+      REAL(KIND=8) PALBD(KDLON,2)
+      REAL(KIND=8) PALBP(KDLON,2)
+      REAL(KIND=8) PCG(KDLON,2,KFLEV)
+      REAL(KIND=8) PCLD(KDLON,KFLEV)
+      REAL(KIND=8) PCLDSW(KDLON,KFLEV)
+      REAL(KIND=8) PCLEAR(KDLON)
+      REAL(KIND=8) PDSIG(KDLON,KFLEV)
+      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)
+      REAL(KIND=8) POZ(KDLON,KFLEV)
+      REAL(KIND=8) PQS(KDLON,KFLEV)
+      REAL(KIND=8) PRMU(KDLON)
+      REAL(KIND=8) PSEC(KDLON)
+      REAL(KIND=8) PTAU(KDLON,2,KFLEV)
+      REAL(KIND=8) PUD(KDLON,5,KFLEV+1)
+      REAL(KIND=8) PWV(KDLON,KFLEV)
+C
+      REAL(KIND=8) PFDOWN(KDLON,KFLEV+1)
+      REAL(KIND=8) PFUP(KDLON,KFLEV+1)
 C
 C* LOCAL VARIABLES:
 C
       INTEGER IIND2(2), IIND3(3)
-      REAL*8 ZCGAZ(KDLON,KFLEV)
-      REAL*8 ZFD(KDLON,KFLEV+1)
-      REAL*8 ZFU(KDLON,KFLEV+1) 
-      REAL*8 ZG(KDLON)
-      REAL*8 ZGG(KDLON)
-      REAL*8 ZPIZAZ(KDLON,KFLEV)
-      REAL*8 ZRAYL(KDLON)
-      REAL*8 ZRAY1(KDLON,KFLEV+1)
-      REAL*8 ZRAY2(KDLON,KFLEV+1)
-      REAL*8 ZREF(KDLON)
-      REAL*8 ZREFZ(KDLON,2,KFLEV+1)
-      REAL*8 ZRE1(KDLON)
-      REAL*8 ZRE2(KDLON)
-      REAL*8 ZRJ(KDLON,6,KFLEV+1)
-      REAL*8 ZRJ0(KDLON,6,KFLEV+1)
-      REAL*8 ZRK(KDLON,6,KFLEV+1)
-      REAL*8 ZRK0(KDLON,6,KFLEV+1)
-      REAL*8 ZRL(KDLON,8)
-      REAL*8 ZRMUE(KDLON,KFLEV+1)
-      REAL*8 ZRMU0(KDLON,KFLEV+1)
-      REAL*8 ZRMUZ(KDLON)
-      REAL*8 ZRNEB(KDLON)
-      REAL*8 ZRUEF(KDLON,8)
-      REAL*8 ZR1(KDLON) 
-      REAL*8 ZR2(KDLON,2)
-      REAL*8 ZR3(KDLON,3)
-      REAL*8 ZR4(KDLON)
-      REAL*8 ZR21(KDLON)
-      REAL*8 ZR22(KDLON)
-      REAL*8 ZS(KDLON)
-      REAL*8 ZTAUAZ(KDLON,KFLEV)
-      REAL*8 ZTO1(KDLON)
-      REAL*8 ZTR(KDLON,2,KFLEV+1)
-      REAL*8 ZTRA1(KDLON,KFLEV+1)
-      REAL*8 ZTRA2(KDLON,KFLEV+1)
-      REAL*8 ZTR1(KDLON)
-      REAL*8 ZTR2(KDLON)
-      REAL*8 ZW(KDLON)   
-      REAL*8 ZW1(KDLON)
-      REAL*8 ZW2(KDLON,2)
-      REAL*8 ZW3(KDLON,3)
-      REAL*8 ZW4(KDLON)
-      REAL*8 ZW5(KDLON)
+      REAL(KIND=8) ZCGAZ(KDLON,KFLEV)
+      REAL(KIND=8) ZFD(KDLON,KFLEV+1)
+      REAL(KIND=8) ZFU(KDLON,KFLEV+1) 
+      REAL(KIND=8) ZG(KDLON)
+      REAL(KIND=8) ZGG(KDLON)
+      REAL(KIND=8) ZPIZAZ(KDLON,KFLEV)
+      REAL(KIND=8) ZRAYL(KDLON)
+      REAL(KIND=8) ZRAY1(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRAY2(KDLON,KFLEV+1)
+      REAL(KIND=8) ZREF(KDLON)
+      REAL(KIND=8) ZREFZ(KDLON,2,KFLEV+1)
+      REAL(KIND=8) ZRE1(KDLON)
+      REAL(KIND=8) ZRE2(KDLON)
+      REAL(KIND=8) ZRJ(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRJ0(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRK(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRK0(KDLON,6,KFLEV+1)
+      REAL(KIND=8) ZRL(KDLON,8)
+      REAL(KIND=8) ZRMUE(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRMU0(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRMUZ(KDLON)
+      REAL(KIND=8) ZRNEB(KDLON)
+      REAL(KIND=8) ZRUEF(KDLON,8)
+      REAL(KIND=8) ZR1(KDLON) 
+      REAL(KIND=8) ZR2(KDLON,2)
+      REAL(KIND=8) ZR3(KDLON,3)
+      REAL(KIND=8) ZR4(KDLON)
+      REAL(KIND=8) ZR21(KDLON)
+      REAL(KIND=8) ZR22(KDLON)
+      REAL(KIND=8) ZS(KDLON)
+      REAL(KIND=8) ZTAUAZ(KDLON,KFLEV)
+      REAL(KIND=8) ZTO1(KDLON)
+      REAL(KIND=8) ZTR(KDLON,2,KFLEV+1)
+      REAL(KIND=8) ZTRA1(KDLON,KFLEV+1)
+      REAL(KIND=8) ZTRA2(KDLON,KFLEV+1)
+      REAL(KIND=8) ZTR1(KDLON)
+      REAL(KIND=8) ZTR2(KDLON)
+      REAL(KIND=8) ZW(KDLON)   
+      REAL(KIND=8) ZW1(KDLON)
+      REAL(KIND=8) ZW2(KDLON,2)
+      REAL(KIND=8) ZW3(KDLON,3)
+      REAL(KIND=8) ZW4(KDLON)
+      REAL(KIND=8) ZW5(KDLON)
 C
       INTEGER jl, jk, k, jaj, ikm1, ikl, jn, jabs, jkm1
       INTEGER jref, jkl, jklp1, jajp, jkki, jkkp4, jn2j, iabs
-      REAL*8 ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
-C
-C* Prescribed Data:
-C
-      REAL*8 RSUN(2)
-      SAVE RSUN
-c$OMP THREADPRIVATE(RSUN)
-      REAL*8 RRAY(2,6)
-      SAVE RRAY
-c$OMP THREADPRIVATE(RRAY)
-      DATA RSUN(1) / 0.441676 /
-      DATA RSUN(2) / 0.558324 /
-      DATA (RRAY(1,K),K=1,6) /
-     S .428937E-01, .890743E+00,-.288555E+01,
-     S .522744E+01,-.469173E+01, .161645E+01/
-      DATA (RRAY(2,K),K=1,6) /
-     S .697200E-02, .173297E-01,-.850903E-01,
-     S .248261E+00,-.302031E+00, .129662E+00/
+      REAL(KIND=8) ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
+C
+
 C
 C     ------------------------------------------------------------------
@@ -1326,4 +1288,5 @@
      S  , PRK   , PRMU0 , PTAUAZ, PTRA1 , PTRA2                   )
       USE dimphy
+      USE radiation_AR4_param, only : TAUA, RPIZA, RCGA
       IMPLICIT none
 cym#include "dimensions.h"
@@ -1357,56 +1320,43 @@
       INTEGER KNU
 c-OB
-      real*8 flag_aer
-      real*8 tauae(kdlon,kflev,2)
-      real*8 pizae(kdlon,kflev,2)
-      real*8 cgae(kdlon,kflev,2)
-      REAL*8 PAER(KDLON,KFLEV,5)
-      REAL*8 PALBP(KDLON,2)
-      REAL*8 PDSIG(KDLON,KFLEV)
-      REAL*8 PRAYL(KDLON)
-      REAL*8 PSEC(KDLON)
-C
-      REAL*8 PCGAZ(KDLON,KFLEV)     
-      REAL*8 PPIZAZ(KDLON,KFLEV)
-      REAL*8 PRAY1(KDLON,KFLEV+1)
-      REAL*8 PRAY2(KDLON,KFLEV+1)
-      REAL*8 PREFZ(KDLON,2,KFLEV+1)
-      REAL*8 PRJ(KDLON,6,KFLEV+1)
-      REAL*8 PRK(KDLON,6,KFLEV+1)
-      REAL*8 PRMU0(KDLON,KFLEV+1)
-      REAL*8 PTAUAZ(KDLON,KFLEV)
-      REAL*8 PTRA1(KDLON,KFLEV+1)
-      REAL*8 PTRA2(KDLON,KFLEV+1)
+      real(kind=8) flag_aer
+      real(kind=8) tauae(kdlon,kflev,2)
+      real(kind=8) pizae(kdlon,kflev,2)
+      real(kind=8) cgae(kdlon,kflev,2)
+      REAL(KIND=8) PAER(KDLON,KFLEV,5)
+      REAL(KIND=8) PALBP(KDLON,2)
+      REAL(KIND=8) PDSIG(KDLON,KFLEV)
+      REAL(KIND=8) PRAYL(KDLON)
+      REAL(KIND=8) PSEC(KDLON)
+C
+      REAL(KIND=8) PCGAZ(KDLON,KFLEV)     
+      REAL(KIND=8) PPIZAZ(KDLON,KFLEV)
+      REAL(KIND=8) PRAY1(KDLON,KFLEV+1)
+      REAL(KIND=8) PRAY2(KDLON,KFLEV+1)
+      REAL(KIND=8) PREFZ(KDLON,2,KFLEV+1)
+      REAL(KIND=8) PRJ(KDLON,6,KFLEV+1)
+      REAL(KIND=8) PRK(KDLON,6,KFLEV+1)
+      REAL(KIND=8) PRMU0(KDLON,KFLEV+1)
+      REAL(KIND=8) PTAUAZ(KDLON,KFLEV)
+      REAL(KIND=8) PTRA1(KDLON,KFLEV+1)
+      REAL(KIND=8) PTRA2(KDLON,KFLEV+1)
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZC0I(KDLON,KFLEV+1)       
-      REAL*8 ZCLE0(KDLON,KFLEV)
-      REAL*8 ZCLEAR(KDLON)
-      REAL*8 ZR21(KDLON)
-      REAL*8 ZR23(KDLON)
-      REAL*8 ZSS0(KDLON)
-      REAL*8 ZSCAT(KDLON)
-      REAL*8 ZTR(KDLON,2,KFLEV+1)
+      REAL(KIND=8) ZC0I(KDLON,KFLEV+1)       
+      REAL(KIND=8) ZCLE0(KDLON,KFLEV)
+      REAL(KIND=8) ZCLEAR(KDLON)
+      REAL(KIND=8) ZR21(KDLON)
+      REAL(KIND=8) ZR23(KDLON)
+      REAL(KIND=8) ZSS0(KDLON)
+      REAL(KIND=8) ZSCAT(KDLON)
+      REAL(KIND=8) ZTR(KDLON,2,KFLEV+1)
 C
       INTEGER jl, jk, ja, jae, jkl, jklp1, jaj, jkm1, in
-      REAL*8 ZTRAY, ZGAR, ZRATIO, ZFF, ZFACOA, ZCORAE
-      REAL*8 ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZMU1, ZDEN1
-      REAL*8 ZBMU0, ZBMU1, ZRE11
-C
-C* Prescribed Data for Aerosols:
-C
-      REAL*8 TAUA(2,5), RPIZA(2,5), RCGA(2,5)
-      SAVE TAUA, RPIZA, RCGA
-c$OMP THREADPRIVATE(TAUA, RPIZA, RCGA)
-      DATA ((TAUA(IN,JA),JA=1,5),IN=1,2) /
-     S .730719, .912819, .725059, .745405, .682188 ,
-     S .730719, .912819, .725059, .745405, .682188 /
-      DATA ((RPIZA(IN,JA),JA=1,5),IN=1,2) /
-     S .872212, .982545, .623143, .944887, .997975 ,
-     S .872212, .982545, .623143, .944887, .997975 /
-      DATA ((RCGA (IN,JA),JA=1,5),IN=1,2) /
-     S .647596, .739002, .580845, .662657, .624246 ,
-     S .647596, .739002, .580845, .662657, .624246 /
+      REAL(KIND=8) ZTRAY, ZGAR, ZRATIO, ZFF, ZFACOA, ZCORAE
+      REAL(KIND=8) ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZMU1, ZDEN1
+      REAL(KIND=8) ZBMU0, ZBMU1, ZRE11
+C
+
 C     ------------------------------------------------------------------
 C
@@ -1416,4 +1366,5 @@
  100  CONTINUE
 C
+!cdir collapse
       DO 103 JK = 1 , KFLEV+1
       DO 102 JA = 1 , 6
@@ -1724,51 +1675,51 @@
 C
       INTEGER KNU
-      REAL*8 PALBD(KDLON,2)
-      REAL*8 PCG(KDLON,2,KFLEV)
-      REAL*8 PCLD(KDLON,KFLEV)
-      REAL*8 PDSIG(KDLON,KFLEV)
-      REAL*8 POMEGA(KDLON,2,KFLEV)
-      REAL*8 PRAYL(KDLON)
-      REAL*8 PSEC(KDLON)
-      REAL*8 PTAU(KDLON,2,KFLEV)
-C
-      REAL*8 PRAY1(KDLON,KFLEV+1)
-      REAL*8 PRAY2(KDLON,KFLEV+1)
-      REAL*8 PREFZ(KDLON,2,KFLEV+1)
-      REAL*8 PRJ(KDLON,6,KFLEV+1)
-      REAL*8 PRK(KDLON,6,KFLEV+1)
-      REAL*8 PRMUE(KDLON,KFLEV+1)
-      REAL*8 PCGAZ(KDLON,KFLEV)
-      REAL*8 PPIZAZ(KDLON,KFLEV)
-      REAL*8 PTAUAZ(KDLON,KFLEV)
-      REAL*8 PTRA1(KDLON,KFLEV+1)
-      REAL*8 PTRA2(KDLON,KFLEV+1)
+      REAL(KIND=8) PALBD(KDLON,2)
+      REAL(KIND=8) PCG(KDLON,2,KFLEV)
+      REAL(KIND=8) PCLD(KDLON,KFLEV)
+      REAL(KIND=8) PDSIG(KDLON,KFLEV)
+      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)
+      REAL(KIND=8) PRAYL(KDLON)
+      REAL(KIND=8) PSEC(KDLON)
+      REAL(KIND=8) PTAU(KDLON,2,KFLEV)
+C
+      REAL(KIND=8) PRAY1(KDLON,KFLEV+1)
+      REAL(KIND=8) PRAY2(KDLON,KFLEV+1)
+      REAL(KIND=8) PREFZ(KDLON,2,KFLEV+1)
+      REAL(KIND=8) PRJ(KDLON,6,KFLEV+1)
+      REAL(KIND=8) PRK(KDLON,6,KFLEV+1)
+      REAL(KIND=8) PRMUE(KDLON,KFLEV+1)
+      REAL(KIND=8) PCGAZ(KDLON,KFLEV)
+      REAL(KIND=8) PPIZAZ(KDLON,KFLEV)
+      REAL(KIND=8) PTAUAZ(KDLON,KFLEV)
+      REAL(KIND=8) PTRA1(KDLON,KFLEV+1)
+      REAL(KIND=8) PTRA2(KDLON,KFLEV+1)
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZC1I(KDLON,KFLEV+1)
-      REAL*8 ZCLEQ(KDLON,KFLEV)
-      REAL*8 ZCLEAR(KDLON)
-      REAL*8 ZCLOUD(KDLON)
-      REAL*8 ZGG(KDLON)
-      REAL*8 ZREF(KDLON)
-      REAL*8 ZRE1(KDLON)
-      REAL*8 ZRE2(KDLON)
-      REAL*8 ZRMUZ(KDLON)
-      REAL*8 ZRNEB(KDLON)
-      REAL*8 ZR21(KDLON)
-      REAL*8 ZR22(KDLON)
-      REAL*8 ZR23(KDLON)
-      REAL*8 ZSS1(KDLON)
-      REAL*8 ZTO1(KDLON)
-      REAL*8 ZTR(KDLON,2,KFLEV+1)
-      REAL*8 ZTR1(KDLON)
-      REAL*8 ZTR2(KDLON)
-      REAL*8 ZW(KDLON)
+      REAL(KIND=8) ZC1I(KDLON,KFLEV+1)
+      REAL(KIND=8) ZCLEQ(KDLON,KFLEV)
+      REAL(KIND=8) ZCLEAR(KDLON)
+      REAL(KIND=8) ZCLOUD(KDLON)
+      REAL(KIND=8) ZGG(KDLON)
+      REAL(KIND=8) ZREF(KDLON)
+      REAL(KIND=8) ZRE1(KDLON)
+      REAL(KIND=8) ZRE2(KDLON)
+      REAL(KIND=8) ZRMUZ(KDLON)
+      REAL(KIND=8) ZRNEB(KDLON)
+      REAL(KIND=8) ZR21(KDLON)
+      REAL(KIND=8) ZR22(KDLON)
+      REAL(KIND=8) ZR23(KDLON)
+      REAL(KIND=8) ZSS1(KDLON)
+      REAL(KIND=8) ZTO1(KDLON)
+      REAL(KIND=8) ZTR(KDLON,2,KFLEV+1)
+      REAL(KIND=8) ZTR1(KDLON)
+      REAL(KIND=8) ZTR2(KDLON)
+      REAL(KIND=8) ZW(KDLON)
 C
       INTEGER jk, jl, ja, jkl, jklp1, jkm1, jaj
-      REAL*8 ZFACOA, ZFACOC, ZCORAE, ZCORCD
-      REAL*8 ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZDEN1
-      REAL*8 ZMU1, ZRE11, ZBMU0, ZBMU1
+      REAL(KIND=8) ZFACOA, ZFACOC, ZCORAE, ZCORCD
+      REAL(KIND=8) ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZDEN1
+      REAL(KIND=8) ZMU1, ZRE11, ZBMU0, ZBMU1
 C
 C     ------------------------------------------------------------------
@@ -2077,25 +2028,26 @@
 C* ARGUMENTS:
 C
-      REAL*8 PGG(KDLON)   ! ASSYMETRY FACTOR
-      REAL*8 PREF(KDLON)  ! REFLECTIVITY OF THE UNDERLYING LAYER
-      REAL*8 PRMUZ(KDLON) ! COSINE OF SOLAR ZENITH ANGLE
-      REAL*8 PTO1(KDLON)  ! OPTICAL THICKNESS
-      REAL*8 PW(KDLON)    ! SINGLE SCATTERING ALBEDO
-      REAL*8 PRE1(KDLON)  ! LAYER REFLECTIVITY (NO UNDERLYING-LAYER REFLECTION)
-      REAL*8 PRE2(KDLON)  ! LAYER REFLECTIVITY
-      REAL*8 PTR1(KDLON)  ! LAYER TRANSMISSIVITY (NO UNDERLYING-LAYER REFLECTION)
-      REAL*8 PTR2(KDLON)  ! LAYER TRANSMISSIVITY
+      REAL(KIND=8) PGG(KDLON)   ! ASSYMETRY FACTOR
+      REAL(KIND=8) PREF(KDLON)  ! REFLECTIVITY OF THE UNDERLYING LAYER
+      REAL(KIND=8) PRMUZ(KDLON) ! COSINE OF SOLAR ZENITH ANGLE
+      REAL(KIND=8) PTO1(KDLON)  ! OPTICAL THICKNESS
+      REAL(KIND=8) PW(KDLON)    ! SINGLE SCATTERING ALBEDO
+      REAL(KIND=8) PRE1(KDLON)  ! LAYER REFLECTIVITY (NO UNDERLYING-LAYER REFLECTION)
+      REAL(KIND=8) PRE2(KDLON)  ! LAYER REFLECTIVITY
+      REAL(KIND=8) PTR1(KDLON)  ! LAYER TRANSMISSIVITY (NO UNDERLYING-LAYER REFLECTION)
+      REAL(KIND=8) PTR2(KDLON)  ! LAYER TRANSMISSIVITY
 C
 C* LOCAL VARIABLES:
 C
       INTEGER jl
-      REAL*8 ZFF, ZGP, ZTOP, ZWCP, ZDT, ZX1, ZWM
-      REAL*8 ZRM2, ZRK, ZX2, ZRP, ZALPHA, ZBETA, ZARG
-      REAL*8 ZEXMU0, ZARG2, ZEXKP, ZEXKM, ZXP2P, ZXM2P, ZAP2B, ZAM2B
-      REAL*8 ZA11, ZA12, ZA13, ZA21, ZA22, ZA23
-      REAL*8 ZDENA, ZC1A, ZC2A, ZRI0A, ZRI1A
-      REAL*8 ZRI0B, ZRI1B
-      REAL*8 ZB21, ZB22, ZB23, ZDENB, ZC1B, ZC2B
-      REAL*8 ZRI0C, ZRI1C, ZRI0D, ZRI1D
+      REAL(KIND=8) ZFF, ZGP, ZTOP, ZWCP, ZDT, ZX1, ZWM
+      REAL(KIND=8) ZRM2, ZRK, ZX2, ZRP, ZALPHA, ZBETA, ZARG
+      REAL(KIND=8) ZEXMU0, ZARG2, ZEXKP, ZEXKM, ZXP2P, ZXM2P, ZAP2B,
+     $     ZAM2B
+      REAL(KIND=8) ZA11, ZA12, ZA13, ZA21, ZA22, ZA23
+      REAL(KIND=8) ZDENA, ZC1A, ZC2A, ZRI0A, ZRI1A
+      REAL(KIND=8) ZRI0B, ZRI1B
+      REAL(KIND=8) ZB21, ZB22, ZB23, ZDENB, ZC1B, ZC2B
+      REAL(KIND=8) ZRI0C, ZRI1C, ZRI0D, ZRI1D
 C     ------------------------------------------------------------------
 C
@@ -2175,4 +2127,5 @@
       SUBROUTINE SWTT_LMDAR4 (KNU,KA,PU,PTR)
       USE dimphy
+      USE radiation_AR4_param, only : APAD, BPAD, D
       IMPLICIT none
 cym#include "dimensions.h"
@@ -2212,54 +2165,14 @@
       INTEGER KNU     ! INDEX OF THE SPECTRAL INTERVAL
       INTEGER KA      ! INDEX OF THE ABSORBER
-      REAL*8 PU(KDLON)  ! ABSORBER AMOUNT
-C
-      REAL*8 PTR(KDLON) ! TRANSMISSION FUNCTION
+      REAL(KIND=8) PU(KDLON)  ! ABSORBER AMOUNT
+C
+      REAL(KIND=8) PTR(KDLON) ! TRANSMISSION FUNCTION
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZR1(KDLON), ZR2(KDLON)
+      REAL(KIND=8) ZR1(KDLON), ZR2(KDLON)
       INTEGER jl, i,j
 C
-C* Prescribed Data:
-C
-      REAL*8 APAD(2,3,7), BPAD(2,3,7), D(2,3)
-      SAVE APAD, BPAD, D
-c$OMP THREADPRIVATE(APAD, BPAD, D)
-      DATA ((APAD(1,I,J),I=1,3),J=1,7) /
-     S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
-     S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01,
-     S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00,
-     S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 /
-      DATA ((APAD(2,I,J),I=1,3),J=1,7) /
-     S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
-     S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02,
-     S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00,
-     S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00,
-     S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00,
-     S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00,
-     S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /
-C
-      DATA ((BPAD(1,I,J),I=1,3),J=1,7) /
-     S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
-     S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01,
-     S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00,
-     S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 /
-      DATA ((BPAD(2,I,J),I=1,3),J=1,7) /
-     S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
-     S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02,
-     S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01,
-     S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00,
-     S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00,
-     S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00,
-     S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /
-c
-      DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 /
-      DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 /
+
 C
 C-----------------------------------------------------------------------
@@ -2293,4 +2206,5 @@
       SUBROUTINE SWTT1_LMDAR4(KNU,KABS,KIND, PU, PTR)
       USE dimphy
+      USE radiation_AR4_param, only : APAD, BPAD, D
       IMPLICIT none
 cym#include "dimensions.h"
@@ -2330,56 +2244,16 @@
       INTEGER KABS         ! NUMBER OF ABSORBERS
       INTEGER KIND(KABS)   ! INDICES OF THE ABSORBERS
-      REAL*8 PU(KDLON,KABS)  ! ABSORBER AMOUNT
-C
-      REAL*8 PTR(KDLON,KABS) ! TRANSMISSION FUNCTION
+      REAL(KIND=8) PU(KDLON,KABS)  ! ABSORBER AMOUNT
+C
+      REAL(KIND=8) PTR(KDLON,KABS) ! TRANSMISSION FUNCTION
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZR1(KDLON)
-      REAL*8 ZR2(KDLON)
-      REAL*8 ZU(KDLON)
+      REAL(KIND=8) ZR1(KDLON)
+      REAL(KIND=8) ZR2(KDLON)
+      REAL(KIND=8) ZU(KDLON)
       INTEGER jl, ja, i, j, ia
 C
-C* Prescribed Data:
-C
-      REAL*8 APAD(2,3,7), BPAD(2,3,7), D(2,3)
-      SAVE APAD, BPAD, D
-c$OMP THREADPRIVATE(APAD, BPAD, D)
-      DATA ((APAD(1,I,J),I=1,3),J=1,7) /
-     S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
-     S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01,
-     S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00,
-     S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 /
-      DATA ((APAD(2,I,J),I=1,3),J=1,7) /
-     S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
-     S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02,
-     S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00,
-     S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00,
-     S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00,
-     S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00,
-     S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /
-C
-      DATA ((BPAD(1,I,J),I=1,3),J=1,7) /
-     S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
-     S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01,
-     S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00,
-     S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02,
-     S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 /
-      DATA ((BPAD(2,I,J),I=1,3),J=1,7) /
-     S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
-     S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02,
-     S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01,
-     S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00,
-     S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00,
-     S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00,
-     S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /
-c
-      DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 /
-      DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 /
+
 C-----------------------------------------------------------------------
 C
@@ -2464,58 +2338,60 @@
 C-----------------------------------------------------------------------
 cIM ctes ds clesphys.h
-c     REAL*8 RCO2   ! CO2 CONCENTRATION (IPCC:353.E-06* 44.011/28.97)
-c     REAL*8 RCH4   ! CH4 CONCENTRATION (IPCC: 1.72E-06* 16.043/28.97)
-c     REAL*8 RN2O   ! N2O CONCENTRATION (IPCC: 310.E-09* 44.013/28.97)
-c     REAL*8 RCFC11 ! CFC11 CONCENTRATION (IPCC: 280.E-12* 137.3686/28.97)
-c     REAL*8 RCFC12 ! CFC12 CONCENTRATION (IPCC: 484.E-12* 120.9140/28.97)
+c     REAL(KIND=8) RCO2   ! CO2 CONCENTRATION (IPCC:353.E-06* 44.011/28.97)
+c     REAL(KIND=8) RCH4   ! CH4 CONCENTRATION (IPCC: 1.72E-06* 16.043/28.97)
+c     REAL(KIND=8) RN2O   ! N2O CONCENTRATION (IPCC: 310.E-09* 44.013/28.97)
+c     REAL(KIND=8) RCFC11 ! CFC11 CONCENTRATION (IPCC: 280.E-12* 137.3686/28.97)
+c     REAL(KIND=8) RCFC12 ! CFC12 CONCENTRATION (IPCC: 484.E-12* 120.9140/28.97)
 #include "clesphys.h"
-      REAL*8 PCLDLD(KDLON,KFLEV)  ! DOWNWARD EFFECTIVE CLOUD COVER
-      REAL*8 PCLDLU(KDLON,KFLEV)  ! UPWARD EFFECTIVE CLOUD COVER
-      REAL*8 PDP(KDLON,KFLEV)     ! LAYER PRESSURE THICKNESS (Pa)
-      REAL*8 PDT0(KDLON)          ! SURFACE TEMPERATURE DISCONTINUITY (K)
-      REAL*8 PEMIS(KDLON)         ! SURFACE EMISSIVITY
-      REAL*8 PPMB(KDLON,KFLEV+1)  ! HALF LEVEL PRESSURE (mb)
-      REAL*8 PPSOL(KDLON)         ! SURFACE PRESSURE (Pa)
-      REAL*8 POZON(KDLON,KFLEV)   ! O3 CONCENTRATION (kg/kg)
-      REAL*8 PTL(KDLON,KFLEV+1)   ! HALF LEVEL TEMPERATURE (K)
-      REAL*8 PAER(KDLON,KFLEV,5)  ! OPTICAL THICKNESS OF THE AEROSOLS
-      REAL*8 PTAVE(KDLON,KFLEV)   ! LAYER TEMPERATURE (K)
-      REAL*8 PVIEW(KDLON)         ! COSECANT OF VIEWING ANGLE
-      REAL*8 PWV(KDLON,KFLEV)     ! SPECIFIC HUMIDITY (kg/kg)
-C
-      REAL*8 PCOLR(KDLON,KFLEV)   ! LONG-WAVE TENDENCY (K/day)
-      REAL*8 PCOLR0(KDLON,KFLEV)  ! LONG-WAVE TENDENCY (K/day) clear-sky
-      REAL*8 PTOPLW(KDLON)        ! LONGWAVE FLUX AT T.O.A.
-      REAL*8 PSOLLW(KDLON)        ! LONGWAVE FLUX AT SURFACE
-      REAL*8 PTOPLW0(KDLON)       ! LONGWAVE FLUX AT T.O.A. (CLEAR-SKY)
-      REAL*8 PSOLLW0(KDLON)       ! LONGWAVE FLUX AT SURFACE (CLEAR-SKY)
+      REAL(KIND=8) PCLDLD(KDLON,KFLEV)  ! DOWNWARD EFFECTIVE CLOUD COVER
+      REAL(KIND=8) PCLDLU(KDLON,KFLEV)  ! UPWARD EFFECTIVE CLOUD COVER
+      REAL(KIND=8) PDP(KDLON,KFLEV)     ! LAYER PRESSURE THICKNESS (Pa)
+      REAL(KIND=8) PDT0(KDLON)          ! SURFACE TEMPERATURE DISCONTINUITY (K)
+      REAL(KIND=8) PEMIS(KDLON)         ! SURFACE EMISSIVITY
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1)  ! HALF LEVEL PRESSURE (mb)
+      REAL(KIND=8) PPSOL(KDLON)         ! SURFACE PRESSURE (Pa)
+      REAL(KIND=8) POZON(KDLON,KFLEV)   ! O3 mass fraction
+      REAL(KIND=8) PTL(KDLON,KFLEV+1)   ! HALF LEVEL TEMPERATURE (K)
+      REAL(KIND=8) PAER(KDLON,KFLEV,5)  ! OPTICAL THICKNESS OF THE AEROSOLS
+      REAL(KIND=8) PTAVE(KDLON,KFLEV)   ! LAYER TEMPERATURE (K)
+      REAL(KIND=8) PVIEW(KDLON)         ! COSECANT OF VIEWING ANGLE
+      REAL(KIND=8) PWV(KDLON,KFLEV)     ! SPECIFIC HUMIDITY (kg/kg)
+C
+      REAL(KIND=8) PCOLR(KDLON,KFLEV)   ! LONG-WAVE TENDENCY (K/day)
+      REAL(KIND=8) PCOLR0(KDLON,KFLEV)  ! LONG-WAVE TENDENCY (K/day) clear-sky
+      REAL(KIND=8) PTOPLW(KDLON)        ! LONGWAVE FLUX AT T.O.A.
+      REAL(KIND=8) PSOLLW(KDLON)        ! LONGWAVE FLUX AT SURFACE
+      REAL(KIND=8) PTOPLW0(KDLON)       ! LONGWAVE FLUX AT T.O.A. (CLEAR-SKY)
+      REAL(KIND=8) PSOLLW0(KDLON)       ! LONGWAVE FLUX AT SURFACE (CLEAR-SKY)
 c Rajout LF
-      real*8 psollwdown(kdlon)    ! LONGWAVE downwards flux at surface
+      real(kind=8) psollwdown(kdlon)    ! LONGWAVE downwards flux at surface
 c Rajout IM
-cIM   real*8 psollwdownclr(kdlon) ! LONGWAVE CS downwards flux at surface
-cIM   real*8 ptoplwdown(kdlon)    ! LONGWAVE downwards flux at T.O.A.
-cIM   real*8 ptoplwdownclr(kdlon) ! LONGWAVE CS downwards flux at T.O.A.
+cIM   real(kind=8) psollwdownclr(kdlon) ! LONGWAVE CS downwards flux at surface
+cIM   real(kind=8) ptoplwdown(kdlon)    ! LONGWAVE downwards flux at T.O.A.
+cIM   real(kind=8) ptoplwdownclr(kdlon) ! LONGWAVE CS downwards flux at T.O.A.
 cIM
-      REAL*8 plwup(KDLON,KFLEV+1)  ! LW up total sky
-      REAL*8 plwup0(KDLON,KFLEV+1) ! LW up clear sky
-      REAL*8 plwdn(KDLON,KFLEV+1)  ! LW down total sky
-      REAL*8 plwdn0(KDLON,KFLEV+1) ! LW down clear sky
+      REAL(KIND=8) plwup(KDLON,KFLEV+1)  ! LW up total sky
+      REAL(KIND=8) plwup0(KDLON,KFLEV+1) ! LW up clear sky
+      REAL(KIND=8) plwdn(KDLON,KFLEV+1)  ! LW down total sky
+      REAL(KIND=8) plwdn0(KDLON,KFLEV+1) ! LW down clear sky
 C-------------------------------------------------------------------------
-      REAL*8 ZABCU(KDLON,NUA,3*KFLEV+1)
-      REAL*8 ZOZ(KDLON,KFLEV)
-c
-cym      REAL*8 ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
-cym      REAL*8 ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
-cym      REAL*8 ZBINT(KDLON,KFLEV+1)            ! Intermediate variable
-cym      REAL*8 ZBSUI(KDLON)                    ! Intermediate variable
-cym      REAL*8,ZCTS(KDLON,KFLEV)               ! Intermediate variable
-cym      REAL*8 ZCNTRB(KDLON,KFLEV+1,KFLEV+1)   ! Intermediate variable
+      REAL(KIND=8) ZABCU(KDLON,NUA,3*KFLEV+1)
+
+      REAL(KIND=8) ZOZ(KDLON,KFLEV)
+!     equivalent pressure of ozone in a layer, in Pa
+
+cym      REAL(KIND=8) ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
+cym      REAL(KIND=8) ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
+cym      REAL(KIND=8) ZBINT(KDLON,KFLEV+1)            ! Intermediate variable
+cym      REAL(KIND=8) ZBSUI(KDLON)                    ! Intermediate variable
+cym      REAL(KIND=8) ZCTS(KDLON,KFLEV)               ! Intermediate variable
+cym      REAL(KIND=8) ZCNTRB(KDLON,KFLEV+1,KFLEV+1)   ! Intermediate variable
 cym      SAVE ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB
-      REAL*8,allocatable,save :: ZFLUX(:,:,:) ! RADIATIVE FLUXES (1:up; 2:down)
-      REAL*8,allocatable,save :: ZFLUC(:,:,:) ! CLEAR-SKY RADIATIVE FLUXES
-      REAL*8,allocatable,save :: ZBINT(:,:)            ! Intermediate variable
-      REAL*8,allocatable,save :: ZBSUI(:)                    ! Intermediate variable
-      REAL*8,allocatable,save :: ZCTS(:,:)               ! Intermediate variable
-      REAL*8,allocatable,save :: ZCNTRB(:,:,:)   ! Intermediate variable
+      REAL(KIND=8),allocatable,save :: ZFLUX(:,:,:) ! RADIATIVE FLUXES (1:up; 2:down)
+      REAL(KIND=8),allocatable,save :: ZFLUC(:,:,:) ! CLEAR-SKY RADIATIVE FLUXES
+      REAL(KIND=8),allocatable,save :: ZBINT(:,:)            ! Intermediate variable
+      REAL(KIND=8),allocatable,save :: ZBSUI(:)                    ! Intermediate variable
+      REAL(KIND=8),allocatable,save :: ZCTS(:,:)               ! Intermediate variable
+      REAL(KIND=8),allocatable,save :: ZCNTRB(:,:,:)   ! Intermediate variable
 c$OMP THREADPRIVATE(ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB)
 c
@@ -2550,10 +2426,9 @@
 C
       IF (MOD(itaplw0,lw0pas).EQ.0) THEN
-      DO k = 1, KFLEV  ! convertir ozone de kg/kg en pa/pa
-      DO i = 1, KDLON
-c convertir ozone de kg/kg en pa (modif MPL 100505)
-         ZOZ(i,k) = POZON(i,k)*PDP(i,k) * RMD/RMO3
-c        print *,'LW: ZOZ*10**6=',ZOZ(i,k)*1000000.
-      ENDDO
+c     Compute equivalent pressure of ozone from mass fraction:
+      DO k = 1, KFLEV
+         DO i = 1, KDLON
+            ZOZ(i,k) = POZON(i,k)*PDP(i,k)
+         ENDDO
       ENDDO
 cIM ctes ds clesphys.h   CALL LWU(RCO2,RCH4, RN2O, RCFC11, RCFC12,
@@ -2609,4 +2484,5 @@
      S               PABCU)
       USE dimphy
+      USE radiation_AR4_param, only : TREF, RT1, RAER, AT, BT, OCT
       IMPLICIT none
 cym#include "dimensions.h"
@@ -2647,42 +2523,42 @@
 C* ARGUMENTS:
 cIM ctes ds clesphys.h
-c     REAL*8 RCO2
-c     REAL*8 RCH4, RN2O, RCFC11, RCFC12
+c     REAL(KIND=8) RCO2
+c     REAL(KIND=8) RCH4, RN2O, RCFC11, RCFC12
 #include "clesphys.h"
-      REAL*8 PAER(KDLON,KFLEV,5)
-      REAL*8 PDP(KDLON,KFLEV)
-      REAL*8 PPMB(KDLON,KFLEV+1)
-      REAL*8 PPSOL(KDLON)
-      REAL*8 POZ(KDLON,KFLEV)
-      REAL*8 PTAVE(KDLON,KFLEV)
-      REAL*8 PVIEW(KDLON)
-      REAL*8 PWV(KDLON,KFLEV)
-C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
+      REAL(KIND=8) PAER(KDLON,KFLEV,5)
+      REAL(KIND=8) PDP(KDLON,KFLEV)
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1)
+      REAL(KIND=8) PPSOL(KDLON)
+      REAL(KIND=8) POZ(KDLON,KFLEV)
+      REAL(KIND=8) PTAVE(KDLON,KFLEV)
+      REAL(KIND=8) PVIEW(KDLON)
+      REAL(KIND=8) PWV(KDLON,KFLEV)
+C
+      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
 C
 C-----------------------------------------------------------------------
 C* LOCAL VARIABLES:
-      REAL*8 ZABLY(KDLON,NUA,3*KFLEV+1)
-      REAL*8 ZDUC(KDLON,3*KFLEV+1)
-      REAL*8 ZPHIO(KDLON)
-      REAL*8 ZPSC2(KDLON)
-      REAL*8 ZPSC3(KDLON)
-      REAL*8 ZPSH1(KDLON)
-      REAL*8 ZPSH2(KDLON)
-      REAL*8 ZPSH3(KDLON)
-      REAL*8 ZPSH4(KDLON)
-      REAL*8 ZPSH5(KDLON)
-      REAL*8 ZPSH6(KDLON)
-      REAL*8 ZPSIO(KDLON)
-      REAL*8 ZTCON(KDLON)
-      REAL*8 ZPHM6(KDLON)
-      REAL*8 ZPSM6(KDLON)
-      REAL*8 ZPHN6(KDLON)
-      REAL*8 ZPSN6(KDLON)
-      REAL*8 ZSSIG(KDLON,3*KFLEV+1)
-      REAL*8 ZTAVI(KDLON)
-      REAL*8 ZUAER(KDLON,Ninter)
-      REAL*8 ZXOZ(KDLON)
-      REAL*8 ZXWV(KDLON)
+      REAL(KIND=8) ZABLY(KDLON,NUA,3*KFLEV+1)
+      REAL(KIND=8) ZDUC(KDLON,3*KFLEV+1)
+      REAL(KIND=8) ZPHIO(KDLON)
+      REAL(KIND=8) ZPSC2(KDLON)
+      REAL(KIND=8) ZPSC3(KDLON)
+      REAL(KIND=8) ZPSH1(KDLON)
+      REAL(KIND=8) ZPSH2(KDLON)
+      REAL(KIND=8) ZPSH3(KDLON)
+      REAL(KIND=8) ZPSH4(KDLON)
+      REAL(KIND=8) ZPSH5(KDLON)
+      REAL(KIND=8) ZPSH6(KDLON)
+      REAL(KIND=8) ZPSIO(KDLON)
+      REAL(KIND=8) ZTCON(KDLON)
+      REAL(KIND=8) ZPHM6(KDLON)
+      REAL(KIND=8) ZPSM6(KDLON)
+      REAL(KIND=8) ZPHN6(KDLON)
+      REAL(KIND=8) ZPSN6(KDLON)
+      REAL(KIND=8) ZSSIG(KDLON,3*KFLEV+1)
+      REAL(KIND=8) ZTAVI(KDLON)
+      REAL(KIND=8) ZUAER(KDLON,Ninter)
+      REAL(KIND=8) ZXOZ(KDLON)
+      REAL(KIND=8) ZXWV(KDLON)
 C
       INTEGER jl, jk, jkj, jkjr, jkjp, ig1
@@ -2691,74 +2567,16 @@
       INTEGER jae1, jae2, jae3, jae, jjpn
       INTEGER ir, jc, jcp1
-      REAL*8 zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup
-      REAL*8 zfppw, ztx, ztx2, zzably
-      REAL*8 zcah1, zcbh1, zcah2, zcbh2, zcah3, zcbh3
-      REAL*8 zcah4, zcbh4, zcah5, zcbh5, zcah6, zcbh6
-      REAL*8 zcac8, zcbc8
-      REAL*8 zalup, zdiff
+      REAL(KIND=8) zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup
+      REAL(KIND=8) zfppw, ztx, ztx2, zzably
+      REAL(KIND=8) zcah1, zcbh1, zcah2, zcbh2, zcah3, zcbh3
+      REAL(KIND=8) zcah4, zcbh4, zcah5, zcbh5, zcah6, zcbh6
+      REAL(KIND=8) zcac8, zcbc8
+      REAL(KIND=8) zalup, zdiff
 c
-      REAL*8 PVGCO2, PVGH2O, PVGO3
-C
-      REAL*8 R10E  ! DECIMAL/NATURAL LOG.FACTOR
+      REAL(KIND=8) PVGCO2, PVGH2O, PVGO3
+C
+      REAL(KIND=8) R10E  ! DECIMAL/NATURAL LOG.FACTOR
       PARAMETER (R10E=0.4342945)
-c
-c Used Data Block:
-c
-      REAL*8 TREF
-      SAVE TREF
-c$OMP THREADPRIVATE(TREF)
-      REAL*8 RT1(2)
-      SAVE RT1
-c$OMP THREADPRIVATE(RT1)
-      REAL*8 RAER(5,5)
-      SAVE RAER
-c$OMP THREADPRIVATE(RAER)
-      REAL*8 AT(8,3), BT(8,3)
-      SAVE AT, BT
-c$OMP THREADPRIVATE(AT, BT)
-      REAL*8 OCT(4)
-      SAVE OCT
-c$OMP THREADPRIVATE(OCT)
-      DATA TREF /250.0/
-      DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 /
-      DATA RAER / .038520, .037196, .040532, .054934, .038520
-     1          , .12613 , .18313 , .10357 , .064106, .126130
-     2          , .012579, .013649, .018652, .025181, .012579
-     3          , .011890, .016142, .021105, .028908, .011890
-     4          , .013792, .026810, .052203, .066338, .013792 /
-      DATA (AT(1,IR),IR=1,3) /
-     S 0.298199E-02,-.394023E-03,0.319566E-04 /
-      DATA (BT(1,IR),IR=1,3) /
-     S-0.106432E-04,0.660324E-06,0.174356E-06 /
-      DATA (AT(2,IR),IR=1,3) /
-     S 0.143676E-01,0.366501E-02,-.160822E-02 /
-      DATA (BT(2,IR),IR=1,3) /
-     S-0.553979E-04,-.101701E-04,0.920868E-05 /
-      DATA (AT(3,IR),IR=1,3) /
-     S 0.197861E-01,0.315541E-02,-.174547E-02 /
-      DATA (BT(3,IR),IR=1,3) /
-     S-0.877012E-04,0.513302E-04,0.523138E-06 /
-      DATA (AT(4,IR),IR=1,3) /
-     S 0.289560E-01,-.208807E-02,-.121943E-02 /
-      DATA (BT(4,IR),IR=1,3) /
-     S-0.165960E-03,0.157704E-03,-.146427E-04 /
-      DATA (AT(5,IR),IR=1,3) /
-     S 0.103800E-01,0.436296E-02,-.161431E-02 /
-      DATA (BT(5,IR),IR=1,3) /
-     S -.276744E-04,-.327381E-04,0.127646E-04 /
-      DATA (AT(6,IR),IR=1,3) /
-     S 0.868859E-02,-.972752E-03,0.000000E-00 /
-      DATA (BT(6,IR),IR=1,3) /
-     S -.278412E-04,-.713940E-06,0.117469E-05 /
-      DATA (AT(7,IR),IR=1,3) /
-     S 0.250073E-03,0.455875E-03,0.109242E-03 /
-      DATA (BT(7,IR),IR=1,3) /
-     S 0.199846E-05,-.216313E-05,0.175991E-06 /
-      DATA (AT(8,IR),IR=1,3) /
-     S 0.307423E-01,0.110879E-02,-.322172E-03 /
-      DATA (BT(8,IR),IR=1,3) /
-     S-0.108482E-03,0.258096E-05,-.814575E-06 /
-c
-      DATA OCT /-.326E-03, -.102E-05, .137E-02, -.535E-05/
+
 C-----------------------------------------------------------------------
 c
@@ -3048,32 +2866,32 @@
       INTEGER KLIM
 C
-      REAL*8 PDP(KDLON,KFLEV)
-      REAL*8 PDT0(KDLON)
-      REAL*8 PEMIS(KDLON)
-      REAL*8 PPMB(KDLON,KFLEV+1)
-      REAL*8 PTL(KDLON,KFLEV+1)
-      REAL*8 PTAVE(KDLON,KFLEV)
-C
-      REAL*8 PFLUC(KDLON,2,KFLEV+1)
+      REAL(KIND=8) PDP(KDLON,KFLEV)
+      REAL(KIND=8) PDT0(KDLON)
+      REAL(KIND=8) PEMIS(KDLON)
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1)
+      REAL(KIND=8) PTL(KDLON,KFLEV+1)
+      REAL(KIND=8) PTAVE(KDLON,KFLEV)
+C
+      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1)
 C     
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1)
-      REAL*8 PBINT(KDLON,KFLEV+1)
-      REAL*8 PBSUI(KDLON)
-      REAL*8 PCTS(KDLON,KFLEV)
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1)
+      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1)
+      REAL(KIND=8) PBINT(KDLON,KFLEV+1)
+      REAL(KIND=8) PBSUI(KDLON)
+      REAL(KIND=8) PCTS(KDLON,KFLEV)
+      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1)
 C
 C-------------------------------------------------------------------------
 C
 C* LOCAL VARIABLES:
-      REAL*8 ZB(KDLON,Ninter,KFLEV+1)
-      REAL*8 ZBSUR(KDLON,Ninter)
-      REAL*8 ZBTOP(KDLON,Ninter)
-      REAL*8 ZDBSL(KDLON,Ninter,KFLEV*2)
-      REAL*8 ZGA(KDLON,8,2,KFLEV)
-      REAL*8 ZGB(KDLON,8,2,KFLEV)
-      REAL*8 ZGASUR(KDLON,8,2)
-      REAL*8 ZGBSUR(KDLON,8,2)
-      REAL*8 ZGATOP(KDLON,8,2)
-      REAL*8 ZGBTOP(KDLON,8,2)
+      REAL(KIND=8) ZB(KDLON,Ninter,KFLEV+1)
+      REAL(KIND=8) ZBSUR(KDLON,Ninter)
+      REAL(KIND=8) ZBTOP(KDLON,Ninter)
+      REAL(KIND=8) ZDBSL(KDLON,Ninter,KFLEV*2)
+      REAL(KIND=8) ZGA(KDLON,8,2,KFLEV)
+      REAL(KIND=8) ZGB(KDLON,8,2,KFLEV)
+      REAL(KIND=8) ZGASUR(KDLON,8,2)
+      REAL(KIND=8) ZGBSUR(KDLON,8,2)
+      REAL(KIND=8) ZGATOP(KDLON,8,2)
+      REAL(KIND=8) ZGBTOP(KDLON,8,2)
 C
       INTEGER nuaer, ntraer
@@ -3151,28 +2969,29 @@
 C* ARGUMENTS:
       INTEGER klim
-      REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
-      REAL*8 PBINT(KDLON,KFLEV+1)   ! HALF LEVEL PLANCK FUNCTION
-      REAL*8 PBSUIN(KDLON)          ! SURFACE PLANCK FUNCTION
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) !CLEAR-SKY ENERGY EXCHANGE
-      REAL*8 PCTS(KDLON,KFLEV)      ! CLEAR-SKY LAYER COOLING-TO-SPACE
+      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
+      REAL(KIND=8) PBINT(KDLON,KFLEV+1)   ! HALF LEVEL PLANCK FUNCTION
+      REAL(KIND=8) PBSUIN(KDLON)          ! SURFACE PLANCK FUNCTION
+      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) !CLEAR-SKY ENERGY EXCHANGE
+      REAL(KIND=8) PCTS(KDLON,KFLEV)      ! CLEAR-SKY LAYER COOLING-TO-SPACE
 c
-      REAL*8 PCLDLD(KDLON,KFLEV)
-      REAL*8 PCLDLU(KDLON,KFLEV)
-      REAL*8 PEMIS(KDLON)
-C
-      REAL*8 PFLUX(KDLON,2,KFLEV+1)
+      REAL(KIND=8) PCLDLD(KDLON,KFLEV)
+      REAL(KIND=8) PCLDLU(KDLON,KFLEV)
+      REAL(KIND=8) PEMIS(KDLON)
+C
+      REAL(KIND=8) PFLUX(KDLON,2,KFLEV+1)
 C-----------------------------------------------------------------------
 C* LOCAL VARIABLES:
       INTEGER IMX(KDLON), IMXP(KDLON)
 C
-      REAL*8 ZCLEAR(KDLON),ZCLOUD(KDLON),ZDNF(KDLON,KFLEV+1,KFLEV+1)
+      REAL(KIND=8) ZCLEAR(KDLON),ZCLOUD(KDLON),
+     $     ZDNF(KDLON,KFLEV+1,KFLEV+1)
      S  , ZFD(KDLON), ZFN10(KDLON), ZFU(KDLON)
      S  , ZUPF(KDLON,KFLEV+1,KFLEV+1)
-      REAL*8 ZCLM(KDLON,KFLEV+1,KFLEV+1)
+      REAL(KIND=8) ZCLM(KDLON,KFLEV+1,KFLEV+1)
 C
       INTEGER jk, jl, imaxc, imx1, imx2, jkj, jkp1, jkm1
       INTEGER jk1, jk2, jkc, jkcp1, jcloud
       INTEGER imxm1, imxp1
-      REAL*8 zcfrac
+      REAL(KIND=8) zcfrac
 C     ------------------------------------------------------------------
 C
@@ -3504,4 +3323,5 @@
      S  , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP)
       USE dimphy
+      USE radiation_AR4_param, only : TINTP, XP, GA, GB
       IMPLICIT none
 cym#include "dimensions.h"
@@ -3560,26 +3380,26 @@
 C ARGUMENTS:
 C
-      REAL*8 PDT0(KDLON)
-      REAL*8 PTAVE(KDLON,KFLEV)
-      REAL*8 PTL(KDLON,KFLEV+1)
-C
-      REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION
-      REAL*8 PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION
-      REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
-      REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
-      REAL*8 PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION
-      REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
-      REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
-      REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
-      REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
-      REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
+      REAL(KIND=8) PDT0(KDLON)
+      REAL(KIND=8) PTAVE(KDLON,KFLEV)
+      REAL(KIND=8) PTL(KDLON,KFLEV+1)
+C
+      REAL(KIND=8) PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION
+      REAL(KIND=8) PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION
+      REAL(KIND=8) PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
+      REAL(KIND=8) PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
+      REAL(KIND=8) PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION
+      REAL(KIND=8) PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
+      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
+      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
+      REAL(KIND=8) PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
+      REAL(KIND=8) PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
+      REAL(KIND=8) PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
+      REAL(KIND=8) PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
 C
 C-------------------------------------------------------------------------
 C*  LOCAL VARIABLES:
       INTEGER INDB(KDLON),INDS(KDLON)
-      REAL*8 ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
-      REAL*8 ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)
+      REAL(KIND=8) ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
+      REAL(KIND=8) ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)
 c
       INTEGER jk, jl, ic, jnu, jf, jg
@@ -3587,10 +3407,10 @@
       INTEGER k, j, ixtox, indto, ixtx, indt
       INTEGER indsu, indtp
-      REAL*8 zdsto1, zdstox, zdst1, zdstx
+      REAL(KIND=8) zdsto1, zdstox, zdst1, zdstx
 c
 C* Quelques parametres:
-      REAL*8 TSTAND
+      REAL(KIND=8) TSTAND
       PARAMETER (TSTAND=250.0)
-      REAL*8 TSTP
+      REAL(KIND=8) TSTP
       PARAMETER (TSTP=12.5)
       INTEGER MXIXT
@@ -3598,16 +3418,16 @@
 C
 C* Used Data Block:
-      REAL*8 TINTP(11)
-      SAVE TINTP
-c$OMP THREADPRIVATE(TINTP)
-      REAL*8 GA(11,16,3), GB(11,16,3)
-      SAVE GA, GB
-c$OMP THREADPRIVATE(GA, GB)
-      REAL*8 XP(6,6)
-      SAVE XP
-c$OMP THREADPRIVATE(XP)
+c     REAL*8 TINTP(11)
+c     SAVE TINTP
+cc$OMP THREADPRIVATE(TINTP)
+c     REAL*8 GA(11,16,3), GB(11,16,3)
+c     SAVE GA, GB
+cc$OMP THREADPRIVATE(GA, GB)
+c     REAL*8 XP(6,6)
+c     SAVE XP
+cc$OMP THREADPRIVATE(XP)
 c
-      DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
-     S             262.5, 275., 287.5, 300., 312.5 /
+c     DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
+c    S             262.5, 275., 287.5, 300., 312.5 /
 C-----------------------------------------------------------------------
 C-- WATER VAPOR -- INT.1 -- 0- 500 CM-1 -- FROM ABS225 ----------------
@@ -3622,132 +3442,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 1, 1,IC),IC=1,3) /
-     S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/
-      DATA (GB( 1, 1,IC),IC=1,3) /
-     S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/
-      DATA (GA( 1, 2,IC),IC=1,3) /
-     S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/
-      DATA (GB( 1, 2,IC),IC=1,3) /
-     S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/
+C     DATA (GA( 1, 1,IC),IC=1,3) /
+C    S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/
+C     DATA (GB( 1, 1,IC),IC=1,3) /
+C    S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/
+C     DATA (GA( 1, 2,IC),IC=1,3) /
+C    S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/
+C     DATA (GB( 1, 2,IC),IC=1,3) /
+C    S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 2, 1,IC),IC=1,3) /
-     S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/
-      DATA (GB( 2, 1,IC),IC=1,3) /
-     S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/
-      DATA (GA( 2, 2,IC),IC=1,3) /
-     S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/
-      DATA (GB( 2, 2,IC),IC=1,3) /
-     S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/
+C     DATA (GA( 2, 1,IC),IC=1,3) /
+C    S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/
+C     DATA (GB( 2, 1,IC),IC=1,3) /
+C    S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/
+C     DATA (GA( 2, 2,IC),IC=1,3) /
+C    S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/
+C     DATA (GB( 2, 2,IC),IC=1,3) /
+C    S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 3, 1,IC),IC=1,3) /
-     S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/
-      DATA (GB( 3, 1,IC),IC=1,3) /
-     S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/
-      DATA (GA( 3, 2,IC),IC=1,3) /
-     S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/
-      DATA (GB( 3, 2,IC),IC=1,3) /
-     S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/
+C     DATA (GA( 3, 1,IC),IC=1,3) /
+C    S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/
+C     DATA (GB( 3, 1,IC),IC=1,3) /
+C    S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/
+C     DATA (GA( 3, 2,IC),IC=1,3) /
+C    S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/
+C     DATA (GB( 3, 2,IC),IC=1,3) /
+C    S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 4, 1,IC),IC=1,3) /
-     S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/
-      DATA (GB( 4, 1,IC),IC=1,3) /
-     S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/
-      DATA (GA( 4, 2,IC),IC=1,3) /
-     S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/
-      DATA (GB( 4, 2,IC),IC=1,3) /
-     S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/
+C     DATA (GA( 4, 1,IC),IC=1,3) /
+C    S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/
+C     DATA (GB( 4, 1,IC),IC=1,3) /
+C    S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/
+C     DATA (GA( 4, 2,IC),IC=1,3) /
+C    S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/
+C     DATA (GB( 4, 2,IC),IC=1,3) /
+C    S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 5, 1,IC),IC=1,3) /
-     S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/
-      DATA (GB( 5, 1,IC),IC=1,3) /
-     S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/
-      DATA (GA( 5, 2,IC),IC=1,3) /
-     S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/
-      DATA (GB( 5, 2,IC),IC=1,3) /
-     S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/
+C     DATA (GA( 5, 1,IC),IC=1,3) /
+C    S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/
+C     DATA (GB( 5, 1,IC),IC=1,3) /
+C    S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/
+C     DATA (GA( 5, 2,IC),IC=1,3) /
+C    S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/
+C     DATA (GB( 5, 2,IC),IC=1,3) /
+C    S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 6, 1,IC),IC=1,3) /
-     S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/
-      DATA (GB( 6, 1,IC),IC=1,3) /
-     S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/
-      DATA (GA( 6, 2,IC),IC=1,3) /
-     S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/
-      DATA (GB( 6, 2,IC),IC=1,3) /
-     S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/
+C     DATA (GA( 6, 1,IC),IC=1,3) /
+C    S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/
+C     DATA (GB( 6, 1,IC),IC=1,3) /
+C    S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/
+C     DATA (GA( 6, 2,IC),IC=1,3) /
+C    S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/
+C     DATA (GB( 6, 2,IC),IC=1,3) /
+C    S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 7, 1,IC),IC=1,3) /
-     S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/
-      DATA (GB( 7, 1,IC),IC=1,3) /
-     S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/
-      DATA (GA( 7, 2,IC),IC=1,3) /
-     S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/
-      DATA (GB( 7, 2,IC),IC=1,3) /
-     S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/
+C     DATA (GA( 7, 1,IC),IC=1,3) /
+C    S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/
+C     DATA (GB( 7, 1,IC),IC=1,3) /
+C    S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/
+C     DATA (GA( 7, 2,IC),IC=1,3) /
+C    S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/
+C     DATA (GB( 7, 2,IC),IC=1,3) /
+C    S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 8, 1,IC),IC=1,3) /
-     S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/
-      DATA (GB( 8, 1,IC),IC=1,3) /
-     S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/
-      DATA (GA( 8, 2,IC),IC=1,3) /
-     S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/
-      DATA (GB( 8, 2,IC),IC=1,3) /
-     S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/
+C     DATA (GA( 8, 1,IC),IC=1,3) /
+C    S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/
+C     DATA (GB( 8, 1,IC),IC=1,3) /
+C    S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/
+C     DATA (GA( 8, 2,IC),IC=1,3) /
+C    S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/
+C     DATA (GB( 8, 2,IC),IC=1,3) /
+C    S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 9, 1,IC),IC=1,3) /
-     S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/
-      DATA (GB( 9, 1,IC),IC=1,3) /
-     S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/
-      DATA (GA( 9, 2,IC),IC=1,3) /
-     S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/
-      DATA (GB( 9, 2,IC),IC=1,3) /
-     S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/
+C     DATA (GA( 9, 1,IC),IC=1,3) /
+C    S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/
+C     DATA (GB( 9, 1,IC),IC=1,3) /
+C    S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/
+C     DATA (GA( 9, 2,IC),IC=1,3) /
+C    S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/
+C     DATA (GB( 9, 2,IC),IC=1,3) /
+C    S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA(10, 1,IC),IC=1,3) /
-     S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/
-      DATA (GB(10, 1,IC),IC=1,3) /
-     S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/
-      DATA (GA(10, 2,IC),IC=1,3) /
-     S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/
-      DATA (GB(10, 2,IC),IC=1,3) /
-     S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/
+C     DATA (GA(10, 1,IC),IC=1,3) /
+C    S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/
+C     DATA (GB(10, 1,IC),IC=1,3) /
+C    S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/
+C     DATA (GA(10, 2,IC),IC=1,3) /
+C    S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/
+C     DATA (GB(10, 2,IC),IC=1,3) /
+C    S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 1 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA(11, 1,IC),IC=1,3) /
-     S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/
-      DATA (GB(11, 1,IC),IC=1,3) /
-     S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/
-      DATA (GA(11, 2,IC),IC=1,3) /
-     S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/
-      DATA (GB(11, 2,IC),IC=1,3) /
-     S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/
+C     DATA (GA(11, 1,IC),IC=1,3) /
+C    S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/
+C     DATA (GB(11, 1,IC),IC=1,3) /
+C    S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/
+C     DATA (GA(11, 2,IC),IC=1,3) /
+C    S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/
+C     DATA (GB(11, 2,IC),IC=1,3) /
+C    S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/
 C
 C
@@ -3764,132 +3584,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 1, 3,IC),IC=1,3) /
-     S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/
-      DATA (GB( 1, 3,IC),IC=1,3) /
-     S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/
-      DATA (GA( 1, 4,IC),IC=1,3) /
-     S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/
-      DATA (GB( 1, 4,IC),IC=1,3) /
-     S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/
+C     DATA (GA( 1, 3,IC),IC=1,3) /
+C    S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/
+C     DATA (GB( 1, 3,IC),IC=1,3) /
+C    S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/
+C     DATA (GA( 1, 4,IC),IC=1,3) /
+C    S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/
+C     DATA (GB( 1, 4,IC),IC=1,3) /
+C    S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 2, 3,IC),IC=1,3) /
-     S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/
-      DATA (GB( 2, 3,IC),IC=1,3) /
-     S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/
-      DATA (GA( 2, 4,IC),IC=1,3) /
-     S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/
-      DATA (GB( 2, 4,IC),IC=1,3) /
-     S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/
+C     DATA (GA( 2, 3,IC),IC=1,3) /
+C    S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/
+C     DATA (GB( 2, 3,IC),IC=1,3) /
+C    S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/
+C     DATA (GA( 2, 4,IC),IC=1,3) /
+C    S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/
+C     DATA (GB( 2, 4,IC),IC=1,3) /
+C    S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 3, 3,IC),IC=1,3) /
-     S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/
-      DATA (GB( 3, 3,IC),IC=1,3) /
-     S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/
-      DATA (GA( 3, 4,IC),IC=1,3) /
-     S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/
-      DATA (GB( 3, 4,IC),IC=1,3) /
-     S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/
+C     DATA (GA( 3, 3,IC),IC=1,3) /
+C    S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/
+C     DATA (GB( 3, 3,IC),IC=1,3) /
+C    S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/
+C     DATA (GA( 3, 4,IC),IC=1,3) /
+C    S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/
+C     DATA (GB( 3, 4,IC),IC=1,3) /
+C    S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 4, 3,IC),IC=1,3) /
-     S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/
-      DATA (GB( 4, 3,IC),IC=1,3) /
-     S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/
-      DATA (GA( 4, 4,IC),IC=1,3) /
-     S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/
-      DATA (GB( 4, 4,IC),IC=1,3) /
-     S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/
+C     DATA (GA( 4, 3,IC),IC=1,3) /
+C    S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/
+C     DATA (GB( 4, 3,IC),IC=1,3) /
+C    S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/
+C     DATA (GA( 4, 4,IC),IC=1,3) /
+C    S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/
+C     DATA (GB( 4, 4,IC),IC=1,3) /
+C    S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 5, 3,IC),IC=1,3) /
-     S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/
-      DATA (GB( 5, 3,IC),IC=1,3) /
-     S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/
-      DATA (GA( 5, 4,IC),IC=1,3) /
-     S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/
-      DATA (GB( 5, 4,IC),IC=1,3) /
-     S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/
+C     DATA (GA( 5, 3,IC),IC=1,3) /
+C    S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/
+C     DATA (GB( 5, 3,IC),IC=1,3) /
+C    S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/
+C     DATA (GA( 5, 4,IC),IC=1,3) /
+C    S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/
+C     DATA (GB( 5, 4,IC),IC=1,3) /
+C    S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 6, 3,IC),IC=1,3) /
-     S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/
-      DATA (GB( 6, 3,IC),IC=1,3) /
-     S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/
-      DATA (GA( 6, 4,IC),IC=1,3) /
-     S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/
-      DATA (GB( 6, 4,IC),IC=1,3) /
-     S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/
+C     DATA (GA( 6, 3,IC),IC=1,3) /
+C    S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/
+C     DATA (GB( 6, 3,IC),IC=1,3) /
+C    S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/
+C     DATA (GA( 6, 4,IC),IC=1,3) /
+C    S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/
+C     DATA (GB( 6, 4,IC),IC=1,3) /
+C    S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 7, 3,IC),IC=1,3) /
-     S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/
-      DATA (GB( 7, 3,IC),IC=1,3) /
-     S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/
-      DATA (GA( 7, 4,IC),IC=1,3) /
-     S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/
-      DATA (GB( 7, 4,IC),IC=1,3) /
-     S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/
+C     DATA (GA( 7, 3,IC),IC=1,3) /
+C    S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/
+C     DATA (GB( 7, 3,IC),IC=1,3) /
+C    S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/
+C     DATA (GA( 7, 4,IC),IC=1,3) /
+C    S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/
+C     DATA (GB( 7, 4,IC),IC=1,3) /
+C    S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 8, 3,IC),IC=1,3) /
-     S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/
-      DATA (GB( 8, 3,IC),IC=1,3) /
-     S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/
-      DATA (GA( 8, 4,IC),IC=1,3) /
-     S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/
-      DATA (GB( 8, 4,IC),IC=1,3) /
-     S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/
+C     DATA (GA( 8, 3,IC),IC=1,3) /
+C    S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/
+C     DATA (GB( 8, 3,IC),IC=1,3) /
+C    S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/
+C     DATA (GA( 8, 4,IC),IC=1,3) /
+C    S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/
+C     DATA (GB( 8, 4,IC),IC=1,3) /
+C    S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 9, 3,IC),IC=1,3) /
-     S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/
-      DATA (GB( 9, 3,IC),IC=1,3) /
-     S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/
-      DATA (GA( 9, 4,IC),IC=1,3) /
-     S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/
-      DATA (GB( 9, 4,IC),IC=1,3) /
-     S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/
+C     DATA (GA( 9, 3,IC),IC=1,3) /
+C    S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/
+C     DATA (GB( 9, 3,IC),IC=1,3) /
+C    S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/
+C     DATA (GA( 9, 4,IC),IC=1,3) /
+C    S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/
+C     DATA (GB( 9, 4,IC),IC=1,3) /
+C    S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA(10, 3,IC),IC=1,3) /
-     S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/
-      DATA (GB(10, 3,IC),IC=1,3) /
-     S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/
-      DATA (GA(10, 4,IC),IC=1,3) /
-     S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/
-      DATA (GB(10, 4,IC),IC=1,3) /
-     S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/
+C     DATA (GA(10, 3,IC),IC=1,3) /
+C    S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/
+C     DATA (GB(10, 3,IC),IC=1,3) /
+C    S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/
+C     DATA (GA(10, 4,IC),IC=1,3) /
+C    S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/
+C     DATA (GB(10, 4,IC),IC=1,3) /
+C    S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA(11, 3,IC),IC=1,3) /
-     S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/
-      DATA (GB(11, 3,IC),IC=1,3) /
-     S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/
-      DATA (GA(11, 4,IC),IC=1,3) /
-     S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/
-      DATA (GB(11, 4,IC),IC=1,3) /
-     S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/
+C     DATA (GA(11, 3,IC),IC=1,3) /
+C    S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/
+C     DATA (GB(11, 3,IC),IC=1,3) /
+C    S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/
+C     DATA (GA(11, 4,IC),IC=1,3) /
+C    S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/
+C     DATA (GB(11, 4,IC),IC=1,3) /
+C    S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/
 C
 C
@@ -3910,132 +3730,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 1, 7,IC),IC=1,3) /
-     S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/
-      DATA (GB( 1, 7,IC),IC=1,3) /
-     S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/
-      DATA (GA( 1, 8,IC),IC=1,3) /
-     S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/
-      DATA (GB( 1, 8,IC),IC=1,3) /
-     S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/
+C     DATA (GA( 1, 7,IC),IC=1,3) /
+C    S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/
+C     DATA (GB( 1, 7,IC),IC=1,3) /
+C    S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/
+C     DATA (GA( 1, 8,IC),IC=1,3) /
+C    S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/
+C     DATA (GB( 1, 8,IC),IC=1,3) /
+C    S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 2, 7,IC),IC=1,3) /
-     S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/
-      DATA (GB( 2, 7,IC),IC=1,3) /
-     S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/
-      DATA (GA( 2, 8,IC),IC=1,3) /
-     S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/
-      DATA (GB( 2, 8,IC),IC=1,3) /
-     S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/
+C     DATA (GA( 2, 7,IC),IC=1,3) /
+C    S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/
+C     DATA (GB( 2, 7,IC),IC=1,3) /
+C    S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/
+C     DATA (GA( 2, 8,IC),IC=1,3) /
+C    S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/
+C     DATA (GB( 2, 8,IC),IC=1,3) /
+C    S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 3, 7,IC),IC=1,3) /
-     S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/
-      DATA (GB( 3, 7,IC),IC=1,3) /
-     S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/
-      DATA (GA( 3, 8,IC),IC=1,3) /
-     S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/
-      DATA (GB( 3, 8,IC),IC=1,3) /
-     S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/
+C     DATA (GA( 3, 7,IC),IC=1,3) /
+C    S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/
+C     DATA (GB( 3, 7,IC),IC=1,3) /
+C    S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/
+C     DATA (GA( 3, 8,IC),IC=1,3) /
+C    S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/
+C     DATA (GB( 3, 8,IC),IC=1,3) /
+C    S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 4, 7,IC),IC=1,3) /
-     S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/
-      DATA (GB( 4, 7,IC),IC=1,3) /
-     S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/
-      DATA (GA( 4, 8,IC),IC=1,3) /
-     S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/
-      DATA (GB( 4, 8,IC),IC=1,3) /
-     S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/
+C     DATA (GA( 4, 7,IC),IC=1,3) /
+C    S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/
+C     DATA (GB( 4, 7,IC),IC=1,3) /
+C    S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/
+C     DATA (GA( 4, 8,IC),IC=1,3) /
+C    S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/
+C     DATA (GB( 4, 8,IC),IC=1,3) /
+C    S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 5, 7,IC),IC=1,3) /
-     S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/
-      DATA (GB( 5, 7,IC),IC=1,3) /
-     S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/
-      DATA (GA( 5, 8,IC),IC=1,3) /
-     S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/
-      DATA (GB( 5, 8,IC),IC=1,3) /
-     S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/
+C     DATA (GA( 5, 7,IC),IC=1,3) /
+C    S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/
+C     DATA (GB( 5, 7,IC),IC=1,3) /
+C    S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/
+C     DATA (GA( 5, 8,IC),IC=1,3) /
+C    S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/
+C     DATA (GB( 5, 8,IC),IC=1,3) /
+C    S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 6, 7,IC),IC=1,3) /
-     S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/
-      DATA (GB( 6, 7,IC),IC=1,3) /
-     S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/
-      DATA (GA( 6, 8,IC),IC=1,3) /
-     S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/
-      DATA (GB( 6, 8,IC),IC=1,3) /
-     S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/
+C     DATA (GA( 6, 7,IC),IC=1,3) /
+C    S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/
+C     DATA (GB( 6, 7,IC),IC=1,3) /
+C    S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/
+C     DATA (GA( 6, 8,IC),IC=1,3) /
+C    S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/
+C     DATA (GB( 6, 8,IC),IC=1,3) /
+C    S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 7, 7,IC),IC=1,3) /
-     S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/
-      DATA (GB( 7, 7,IC),IC=1,3) /
-     S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/
-      DATA (GA( 7, 8,IC),IC=1,3) /
-     S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/
-      DATA (GB( 7, 8,IC),IC=1,3) /
-     S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/
+C     DATA (GA( 7, 7,IC),IC=1,3) /
+C    S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/
+C     DATA (GB( 7, 7,IC),IC=1,3) /
+C    S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/
+C     DATA (GA( 7, 8,IC),IC=1,3) /
+C    S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/
+C     DATA (GB( 7, 8,IC),IC=1,3) /
+C    S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 8, 7,IC),IC=1,3) /
-     S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/
-      DATA (GB( 8, 7,IC),IC=1,3) /
-     S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/
-      DATA (GA( 8, 8,IC),IC=1,3) /
-     S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/
-      DATA (GB( 8, 8,IC),IC=1,3) /
-     S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/
+C     DATA (GA( 8, 7,IC),IC=1,3) /
+C    S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/
+C     DATA (GB( 8, 7,IC),IC=1,3) /
+C    S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/
+C     DATA (GA( 8, 8,IC),IC=1,3) /
+C    S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/
+C     DATA (GB( 8, 8,IC),IC=1,3) /
+C    S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 9, 7,IC),IC=1,3) /
-     S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/
-      DATA (GB( 9, 7,IC),IC=1,3) /
-     S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/
-      DATA (GA( 9, 8,IC),IC=1,3) /
-     S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/
-      DATA (GB( 9, 8,IC),IC=1,3) /
-     S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/
+C     DATA (GA( 9, 7,IC),IC=1,3) /
+C    S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/
+C     DATA (GB( 9, 7,IC),IC=1,3) /
+C    S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/
+C     DATA (GA( 9, 8,IC),IC=1,3) /
+C    S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/
+C     DATA (GB( 9, 8,IC),IC=1,3) /
+C    S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA(10, 7,IC),IC=1,3) /
-     S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/
-      DATA (GB(10, 7,IC),IC=1,3) /
-     S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/
-      DATA (GA(10, 8,IC),IC=1,3) /
-     S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/
-      DATA (GB(10, 8,IC),IC=1,3) /
-     S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/
+C     DATA (GA(10, 7,IC),IC=1,3) /
+C    S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/
+C     DATA (GB(10, 7,IC),IC=1,3) /
+C    S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/
+C     DATA (GA(10, 8,IC),IC=1,3) /
+C    S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/
+C     DATA (GB(10, 8,IC),IC=1,3) /
+C    S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 3 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA(11, 7,IC),IC=1,3) /
-     S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/
-      DATA (GB(11, 7,IC),IC=1,3) /
-     S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/
-      DATA (GA(11, 8,IC),IC=1,3) /
-     S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/
-      DATA (GB(11, 8,IC),IC=1,3) /
-     S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/
+C     DATA (GA(11, 7,IC),IC=1,3) /
+C    S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/
+C     DATA (GB(11, 7,IC),IC=1,3) /
+C    S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/
+C     DATA (GA(11, 8,IC),IC=1,3) /
+C    S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/
+C     DATA (GB(11, 8,IC),IC=1,3) /
+C    S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/
 C
 C
@@ -4047,132 +3867,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 1, 9,IC),IC=1,3) /
-     S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/
-      DATA (GB( 1, 9,IC),IC=1,3) /
-     S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/
-      DATA (GA( 1,10,IC),IC=1,3) /
-     S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/
-      DATA (GB( 1,10,IC),IC=1,3) /
-     S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/
+C     DATA (GA( 1, 9,IC),IC=1,3) /
+C    S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/
+C     DATA (GB( 1, 9,IC),IC=1,3) /
+C    S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/
+C     DATA (GA( 1,10,IC),IC=1,3) /
+C    S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/
+C     DATA (GB( 1,10,IC),IC=1,3) /
+C    S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 2, 9,IC),IC=1,3) /
-     S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/
-      DATA (GB( 2, 9,IC),IC=1,3) /
-     S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/
-      DATA (GA( 2,10,IC),IC=1,3) /
-     S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/
-      DATA (GB( 2,10,IC),IC=1,3) /
-     S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/
+C     DATA (GA( 2, 9,IC),IC=1,3) /
+C    S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/
+C     DATA (GB( 2, 9,IC),IC=1,3) /
+C    S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/
+C     DATA (GA( 2,10,IC),IC=1,3) /
+C    S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/
+C     DATA (GB( 2,10,IC),IC=1,3) /
+C    S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 3, 9,IC),IC=1,3) /
-     S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/
-      DATA (GB( 3, 9,IC),IC=1,3) /
-     S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/
-      DATA (GA( 3,10,IC),IC=1,3) /
-     S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/
-      DATA (GB( 3,10,IC),IC=1,3) /
-     S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/
+C     DATA (GA( 3, 9,IC),IC=1,3) /
+C    S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/
+C     DATA (GB( 3, 9,IC),IC=1,3) /
+C    S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/
+C     DATA (GA( 3,10,IC),IC=1,3) /
+C    S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/
+C     DATA (GB( 3,10,IC),IC=1,3) /
+C    S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 4, 9,IC),IC=1,3) /
-     S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/
-      DATA (GB( 4, 9,IC),IC=1,3) /
-     S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/
-      DATA (GA( 4,10,IC),IC=1,3) /
-     S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/
-      DATA (GB( 4,10,IC),IC=1,3) /
-     S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/
+C     DATA (GA( 4, 9,IC),IC=1,3) /
+C    S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/
+C     DATA (GB( 4, 9,IC),IC=1,3) /
+C    S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/
+C     DATA (GA( 4,10,IC),IC=1,3) /
+C    S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/
+C     DATA (GB( 4,10,IC),IC=1,3) /
+C    S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 5, 9,IC),IC=1,3) /
-     S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/
-      DATA (GB( 5, 9,IC),IC=1,3) /
-     S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/
-      DATA (GA( 5,10,IC),IC=1,3) /
-     S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/
-      DATA (GB( 5,10,IC),IC=1,3) /
-     S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/
+C     DATA (GA( 5, 9,IC),IC=1,3) /
+C    S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/
+C     DATA (GB( 5, 9,IC),IC=1,3) /
+C    S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/
+C     DATA (GA( 5,10,IC),IC=1,3) /
+C    S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/
+C     DATA (GB( 5,10,IC),IC=1,3) /
+C    S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 6, 9,IC),IC=1,3) /
-     S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/
-      DATA (GB( 6, 9,IC),IC=1,3) /
-     S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/
-      DATA (GA( 6,10,IC),IC=1,3) /
-     S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/
-      DATA (GB( 6,10,IC),IC=1,3) /
-     S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/
+C     DATA (GA( 6, 9,IC),IC=1,3) /
+C    S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/
+C     DATA (GB( 6, 9,IC),IC=1,3) /
+C    S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/
+C     DATA (GA( 6,10,IC),IC=1,3) /
+C    S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/
+C     DATA (GB( 6,10,IC),IC=1,3) /
+C    S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 7, 9,IC),IC=1,3) /
-     S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/
-      DATA (GB( 7, 9,IC),IC=1,3) /
-     S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/
-      DATA (GA( 7,10,IC),IC=1,3) /
-     S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/
-      DATA (GB( 7,10,IC),IC=1,3) /
-     S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/
+C     DATA (GA( 7, 9,IC),IC=1,3) /
+C    S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/
+C     DATA (GB( 7, 9,IC),IC=1,3) /
+C    S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/
+C     DATA (GA( 7,10,IC),IC=1,3) /
+C    S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/
+C     DATA (GB( 7,10,IC),IC=1,3) /
+C    S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 8, 9,IC),IC=1,3) /
-     S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/
-      DATA (GB( 8, 9,IC),IC=1,3) /
-     S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/
-      DATA (GA( 8,10,IC),IC=1,3) /
-     S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/
-      DATA (GB( 8,10,IC),IC=1,3) /
-     S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/
+C     DATA (GA( 8, 9,IC),IC=1,3) /
+C    S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/
+C     DATA (GB( 8, 9,IC),IC=1,3) /
+C    S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/
+C     DATA (GA( 8,10,IC),IC=1,3) /
+C    S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/
+C     DATA (GB( 8,10,IC),IC=1,3) /
+C    S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA( 9, 9,IC),IC=1,3) /
-     S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/
-      DATA (GB( 9, 9,IC),IC=1,3) /
-     S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/
-      DATA (GA( 9,10,IC),IC=1,3) /
-     S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/
-      DATA (GB( 9,10,IC),IC=1,3) /
-     S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/
+C     DATA (GA( 9, 9,IC),IC=1,3) /
+C    S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/
+C     DATA (GB( 9, 9,IC),IC=1,3) /
+C    S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/
+C     DATA (GA( 9,10,IC),IC=1,3) /
+C    S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/
+C     DATA (GB( 9,10,IC),IC=1,3) /
+C    S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA(10, 9,IC),IC=1,3) /
-     S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/
-      DATA (GB(10, 9,IC),IC=1,3) /
-     S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/
-      DATA (GA(10,10,IC),IC=1,3) /
-     S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/
-      DATA (GB(10,10,IC),IC=1,3) /
-     S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/
+C     DATA (GA(10, 9,IC),IC=1,3) /
+C    S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/
+C     DATA (GB(10, 9,IC),IC=1,3) /
+C    S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/
+C     DATA (GA(10,10,IC),IC=1,3) /
+C    S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/
+C     DATA (GB(10,10,IC),IC=1,3) /
+C    S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
-      DATA (GA(11, 9,IC),IC=1,3) /
-     S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/
-      DATA (GB(11, 9,IC),IC=1,3) /
-     S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/
-      DATA (GA(11,10,IC),IC=1,3) /
-     S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/
-      DATA (GB(11,10,IC),IC=1,3) /
-     S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/
+C     DATA (GA(11, 9,IC),IC=1,3) /
+C    S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/
+C     DATA (GB(11, 9,IC),IC=1,3) /
+C    S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/
+C     DATA (GA(11,10,IC),IC=1,3) /
+C    S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/
+C     DATA (GB(11,10,IC),IC=1,3) /
+C    S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/
 C
 C
@@ -4187,132 +4007,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 1, 5,IC),IC=1,3) /
-     S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/
-      DATA (GB( 1, 5,IC),IC=1,3) /
-     S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/
-      DATA (GA( 1, 6,IC),IC=1,3) /
-     S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/
-      DATA (GB( 1, 6,IC),IC=1,3) /
-     S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/
+C     DATA (GA( 1, 5,IC),IC=1,3) /
+C    S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/
+C     DATA (GB( 1, 5,IC),IC=1,3) /
+C    S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/
+C     DATA (GA( 1, 6,IC),IC=1,3) /
+C    S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/
+C     DATA (GB( 1, 6,IC),IC=1,3) /
+C    S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 2, 5,IC),IC=1,3) /
-     S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/
-      DATA (GB( 2, 5,IC),IC=1,3) /
-     S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/
-      DATA (GA( 2, 6,IC),IC=1,3) /
-     S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/
-      DATA (GB( 2, 6,IC),IC=1,3) /
-     S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/
+C     DATA (GA( 2, 5,IC),IC=1,3) /
+C    S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/
+C     DATA (GB( 2, 5,IC),IC=1,3) /
+C    S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/
+C     DATA (GA( 2, 6,IC),IC=1,3) /
+C    S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/
+C     DATA (GB( 2, 6,IC),IC=1,3) /
+C    S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 3, 5,IC),IC=1,3) /
-     S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/
-      DATA (GB( 3, 5,IC),IC=1,3) /
-     S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/
-      DATA (GA( 3, 6,IC),IC=1,3) /
-     S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/
-      DATA (GB( 3, 6,IC),IC=1,3) /
-     S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/
+C     DATA (GA( 3, 5,IC),IC=1,3) /
+C    S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/
+C     DATA (GB( 3, 5,IC),IC=1,3) /
+C    S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/
+C     DATA (GA( 3, 6,IC),IC=1,3) /
+C    S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/
+C     DATA (GB( 3, 6,IC),IC=1,3) /
+C    S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 4, 5,IC),IC=1,3) /
-     S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/
-      DATA (GB( 4, 5,IC),IC=1,3) /
-     S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/
-      DATA (GA( 4, 6,IC),IC=1,3) /
-     S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/
-      DATA (GB( 4, 6,IC),IC=1,3) /
-     S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/
+C     DATA (GA( 4, 5,IC),IC=1,3) /
+C    S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/
+C     DATA (GB( 4, 5,IC),IC=1,3) /
+C    S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/
+C     DATA (GA( 4, 6,IC),IC=1,3) /
+C    S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/
+C     DATA (GB( 4, 6,IC),IC=1,3) /
+C    S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 5, 5,IC),IC=1,3) /
-     S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/
-      DATA (GB( 5, 5,IC),IC=1,3) /
-     S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/
-      DATA (GA( 5, 6,IC),IC=1,3) /
-     S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/
-      DATA (GB( 5, 6,IC),IC=1,3) /
-     S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/
+C     DATA (GA( 5, 5,IC),IC=1,3) /
+C    S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/
+C     DATA (GB( 5, 5,IC),IC=1,3) /
+C    S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/
+C     DATA (GA( 5, 6,IC),IC=1,3) /
+C    S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/
+C     DATA (GB( 5, 6,IC),IC=1,3) /
+C    S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 6, 5,IC),IC=1,3) /
-     S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/
-      DATA (GB( 6, 5,IC),IC=1,3) /
-     S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/
-      DATA (GA( 6, 6,IC),IC=1,3) /
-     S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/
-      DATA (GB( 6, 6,IC),IC=1,3) /
-     S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/
+C     DATA (GA( 6, 5,IC),IC=1,3) /
+C    S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/
+C     DATA (GB( 6, 5,IC),IC=1,3) /
+C    S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/
+C     DATA (GA( 6, 6,IC),IC=1,3) /
+C    S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/
+C     DATA (GB( 6, 6,IC),IC=1,3) /
+C    S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 7, 5,IC),IC=1,3) /
-     S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/
-      DATA (GB( 7, 5,IC),IC=1,3) /
-     S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/
-      DATA (GA( 7, 6,IC),IC=1,3) /
-     S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/
-      DATA (GB( 7, 6,IC),IC=1,3) /
-     S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/
+C     DATA (GA( 7, 5,IC),IC=1,3) /
+C    S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/
+C     DATA (GB( 7, 5,IC),IC=1,3) /
+C    S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/
+C     DATA (GA( 7, 6,IC),IC=1,3) /
+C    S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/
+C     DATA (GB( 7, 6,IC),IC=1,3) /
+C    S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 8, 5,IC),IC=1,3) /
-     S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/
-      DATA (GB( 8, 5,IC),IC=1,3) /
-     S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/
-      DATA (GA( 8, 6,IC),IC=1,3) /
-     S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/
-      DATA (GB( 8, 6,IC),IC=1,3) /
-     S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/
+C     DATA (GA( 8, 5,IC),IC=1,3) /
+C    S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/
+C     DATA (GB( 8, 5,IC),IC=1,3) /
+C    S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/
+C     DATA (GA( 8, 6,IC),IC=1,3) /
+C    S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/
+C     DATA (GB( 8, 6,IC),IC=1,3) /
+C    S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 9, 5,IC),IC=1,3) /
-     S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/
-      DATA (GB( 9, 5,IC),IC=1,3) /
-     S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/
-      DATA (GA( 9, 6,IC),IC=1,3) /
-     S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/
-      DATA (GB( 9, 6,IC),IC=1,3) /
-     S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/
+C     DATA (GA( 9, 5,IC),IC=1,3) /
+C    S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/
+C     DATA (GB( 9, 5,IC),IC=1,3) /
+C    S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/
+C     DATA (GA( 9, 6,IC),IC=1,3) /
+C    S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/
+C     DATA (GB( 9, 6,IC),IC=1,3) /
+C    S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA(10, 5,IC),IC=1,3) /
-     S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/
-      DATA (GB(10, 5,IC),IC=1,3) /
-     S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/
-      DATA (GA(10, 6,IC),IC=1,3) /
-     S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/
-      DATA (GB(10, 6,IC),IC=1,3) /
-     S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/
+C     DATA (GA(10, 5,IC),IC=1,3) /
+C    S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/
+C     DATA (GB(10, 5,IC),IC=1,3) /
+C    S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/
+C     DATA (GA(10, 6,IC),IC=1,3) /
+C    S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/
+C     DATA (GB(10, 6,IC),IC=1,3) /
+C    S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/
 C
 C----- INTERVAL = 5 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA(11, 5,IC),IC=1,3) /
-     S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/
-      DATA (GB(11, 5,IC),IC=1,3) /
-     S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/
-      DATA (GA(11, 6,IC),IC=1,3) /
-     S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/
-      DATA (GB(11, 6,IC),IC=1,3) /
-     S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/
+C     DATA (GA(11, 5,IC),IC=1,3) /
+C    S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/
+C     DATA (GB(11, 5,IC),IC=1,3) /
+C    S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/
+C     DATA (GA(11, 6,IC),IC=1,3) /
+C    S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/
+C     DATA (GB(11, 6,IC),IC=1,3) /
+C    S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/
 C
 C
@@ -4326,132 +4146,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 1,11,IC),IC=1,3) /
-     S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/
-      DATA (GB( 1,11,IC),IC=1,3) /
-     S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/
-      DATA (GA( 1,12,IC),IC=1,3) /
-     S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/
-      DATA (GB( 1,12,IC),IC=1,3) /
-     S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/
+C     DATA (GA( 1,11,IC),IC=1,3) /
+C    S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/
+C     DATA (GB( 1,11,IC),IC=1,3) /
+C    S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/
+C     DATA (GA( 1,12,IC),IC=1,3) /
+C    S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/
+C     DATA (GB( 1,12,IC),IC=1,3) /
+C    S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 2,11,IC),IC=1,3) /
-     S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/
-      DATA (GB( 2,11,IC),IC=1,3) /
-     S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/
-      DATA (GA( 2,12,IC),IC=1,3) /
-     S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/
-      DATA (GB( 2,12,IC),IC=1,3) /
-     S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/
+C     DATA (GA( 2,11,IC),IC=1,3) /
+C    S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/
+C     DATA (GB( 2,11,IC),IC=1,3) /
+C    S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/
+C     DATA (GA( 2,12,IC),IC=1,3) /
+C    S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/
+C     DATA (GB( 2,12,IC),IC=1,3) /
+C    S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 3,11,IC),IC=1,3) /
-     S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/
-      DATA (GB( 3,11,IC),IC=1,3) /
-     S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/
-      DATA (GA( 3,12,IC),IC=1,3) /
-     S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/
-      DATA (GB( 3,12,IC),IC=1,3) /
-     S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/
+C     DATA (GA( 3,11,IC),IC=1,3) /
+C    S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/
+C     DATA (GB( 3,11,IC),IC=1,3) /
+C    S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/
+C     DATA (GA( 3,12,IC),IC=1,3) /
+C    S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/
+C     DATA (GB( 3,12,IC),IC=1,3) /
+C    S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 4,11,IC),IC=1,3) /
-     S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/
-      DATA (GB( 4,11,IC),IC=1,3) /
-     S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/
-      DATA (GA( 4,12,IC),IC=1,3) /
-     S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/
-      DATA (GB( 4,12,IC),IC=1,3) /
-     S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/
+C     DATA (GA( 4,11,IC),IC=1,3) /
+C    S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/
+C     DATA (GB( 4,11,IC),IC=1,3) /
+C    S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/
+C     DATA (GA( 4,12,IC),IC=1,3) /
+C    S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/
+C     DATA (GB( 4,12,IC),IC=1,3) /
+C    S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 5,11,IC),IC=1,3) /
-     S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/
-      DATA (GB( 5,11,IC),IC=1,3) /
-     S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/
-      DATA (GA( 5,12,IC),IC=1,3) /
-     S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/
-      DATA (GB( 5,12,IC),IC=1,3) /
-     S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/
+C     DATA (GA( 5,11,IC),IC=1,3) /
+C    S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/
+C     DATA (GB( 5,11,IC),IC=1,3) /
+C    S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/
+C     DATA (GA( 5,12,IC),IC=1,3) /
+C    S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/
+C     DATA (GB( 5,12,IC),IC=1,3) /
+C    S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 6,11,IC),IC=1,3) /
-     S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/
-      DATA (GB( 6,11,IC),IC=1,3) /
-     S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/
-      DATA (GA( 6,12,IC),IC=1,3) /
-     S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/
-      DATA (GB( 6,12,IC),IC=1,3) /
-     S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/
+C     DATA (GA( 6,11,IC),IC=1,3) /
+C    S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/
+C     DATA (GB( 6,11,IC),IC=1,3) /
+C    S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/
+C     DATA (GA( 6,12,IC),IC=1,3) /
+C    S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/
+C     DATA (GB( 6,12,IC),IC=1,3) /
+C    S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 7,11,IC),IC=1,3) /
-     S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/
-      DATA (GB( 7,11,IC),IC=1,3) /
-     S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/
-      DATA (GA( 7,12,IC),IC=1,3) /
-     S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/
-      DATA (GB( 7,12,IC),IC=1,3) /
-     S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/
+C     DATA (GA( 7,11,IC),IC=1,3) /
+C    S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/
+C     DATA (GB( 7,11,IC),IC=1,3) /
+C    S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/
+C     DATA (GA( 7,12,IC),IC=1,3) /
+C    S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/
+C     DATA (GB( 7,12,IC),IC=1,3) /
+C    S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 8,11,IC),IC=1,3) /
-     S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/
-      DATA (GB( 8,11,IC),IC=1,3) /
-     S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/
-      DATA (GA( 8,12,IC),IC=1,3) /
-     S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/
-      DATA (GB( 8,12,IC),IC=1,3) /
-     S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/
+C     DATA (GA( 8,11,IC),IC=1,3) /
+C    S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/
+C     DATA (GB( 8,11,IC),IC=1,3) /
+C    S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/
+C     DATA (GA( 8,12,IC),IC=1,3) /
+C    S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/
+C     DATA (GB( 8,12,IC),IC=1,3) /
+C    S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA( 9,11,IC),IC=1,3) /
-     S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/
-      DATA (GB( 9,11,IC),IC=1,3) /
-     S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/
-      DATA (GA( 9,12,IC),IC=1,3) /
-     S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/
-      DATA (GB( 9,12,IC),IC=1,3) /
-     S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/
+C     DATA (GA( 9,11,IC),IC=1,3) /
+C    S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/
+C     DATA (GB( 9,11,IC),IC=1,3) /
+C    S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/
+C     DATA (GA( 9,12,IC),IC=1,3) /
+C    S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/
+C     DATA (GB( 9,12,IC),IC=1,3) /
+C    S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA(10,11,IC),IC=1,3) /
-     S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/
-      DATA (GB(10,11,IC),IC=1,3) /
-     S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/
-      DATA (GA(10,12,IC),IC=1,3) /
-     S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/
-      DATA (GB(10,12,IC),IC=1,3) /
-     S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/
+C     DATA (GA(10,11,IC),IC=1,3) /
+C    S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/
+C     DATA (GB(10,11,IC),IC=1,3) /
+C    S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/
+C     DATA (GA(10,12,IC),IC=1,3) /
+C    S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/
+C     DATA (GB(10,12,IC),IC=1,3) /
+C    S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 6 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
-      DATA (GA(11,11,IC),IC=1,3) /
-     S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/
-      DATA (GB(11,11,IC),IC=1,3) /
-     S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/
-      DATA (GA(11,12,IC),IC=1,3) /
-     S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/
-      DATA (GB(11,12,IC),IC=1,3) /
-     S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/
+C     DATA (GA(11,11,IC),IC=1,3) /
+C    S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/
+C     DATA (GB(11,11,IC),IC=1,3) /
+C    S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/
+C     DATA (GA(11,12,IC),IC=1,3) /
+C    S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/
+C     DATA (GB(11,12,IC),IC=1,3) /
+C    S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/
 C
 C
@@ -4471,132 +4291,132 @@
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 1,13,IC),IC=1,3) /
-     S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/
-      DATA (GB( 1,13,IC),IC=1,3) /
-     S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/
-      DATA (GA( 1,14,IC),IC=1,3) /
-     S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/
-      DATA (GB( 1,14,IC),IC=1,3) /
-     S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/
+C     DATA (GA( 1,13,IC),IC=1,3) /
+C    S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/
+C     DATA (GB( 1,13,IC),IC=1,3) /
+C    S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/
+C     DATA (GA( 1,14,IC),IC=1,3) /
+C    S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/
+C     DATA (GB( 1,14,IC),IC=1,3) /
+C    S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 2,13,IC),IC=1,3) /
-     S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/
-      DATA (GB( 2,13,IC),IC=1,3) /
-     S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/
-      DATA (GA( 2,14,IC),IC=1,3) /
-     S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/
-      DATA (GB( 2,14,IC),IC=1,3) /
-     S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/
+C     DATA (GA( 2,13,IC),IC=1,3) /
+C    S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/
+C     DATA (GB( 2,13,IC),IC=1,3) /
+C    S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/
+C     DATA (GA( 2,14,IC),IC=1,3) /
+C    S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/
+C     DATA (GB( 2,14,IC),IC=1,3) /
+C    S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 3,13,IC),IC=1,3) /
-     S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/
-      DATA (GB( 3,13,IC),IC=1,3) /
-     S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/
-      DATA (GA( 3,14,IC),IC=1,3) /
-     S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/
-      DATA (GB( 3,14,IC),IC=1,3) /
-     S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/
+C     DATA (GA( 3,13,IC),IC=1,3) /
+C    S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/
+C     DATA (GB( 3,13,IC),IC=1,3) /
+C    S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/
+C     DATA (GA( 3,14,IC),IC=1,3) /
+C    S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/
+C     DATA (GB( 3,14,IC),IC=1,3) /
+C    S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 4,13,IC),IC=1,3) /
-     S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/
-      DATA (GB( 4,13,IC),IC=1,3) /
-     S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/
-      DATA (GA( 4,14,IC),IC=1,3) /
-     S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/
-      DATA (GB( 4,14,IC),IC=1,3) /
-     S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/
+C     DATA (GA( 4,13,IC),IC=1,3) /
+C    S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/
+C     DATA (GB( 4,13,IC),IC=1,3) /
+C    S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/
+C     DATA (GA( 4,14,IC),IC=1,3) /
+C    S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/
+C     DATA (GB( 4,14,IC),IC=1,3) /
+C    S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 5,13,IC),IC=1,3) /
-     S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/
-      DATA (GB( 5,13,IC),IC=1,3) /
-     S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/
-      DATA (GA( 5,14,IC),IC=1,3) /
-     S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/
-      DATA (GB( 5,14,IC),IC=1,3) /
-     S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/
+C     DATA (GA( 5,13,IC),IC=1,3) /
+C    S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/
+C     DATA (GB( 5,13,IC),IC=1,3) /
+C    S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/
+C     DATA (GA( 5,14,IC),IC=1,3) /
+C    S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/
+C     DATA (GB( 5,14,IC),IC=1,3) /
+C    S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 6,13,IC),IC=1,3) /
-     S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/
-      DATA (GB( 6,13,IC),IC=1,3) /
-     S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/
-      DATA (GA( 6,14,IC),IC=1,3) /
-     S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/
-      DATA (GB( 6,14,IC),IC=1,3) /
-     S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/
+C     DATA (GA( 6,13,IC),IC=1,3) /
+C    S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/
+C     DATA (GB( 6,13,IC),IC=1,3) /
+C    S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/
+C     DATA (GA( 6,14,IC),IC=1,3) /
+C    S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/
+C     DATA (GB( 6,14,IC),IC=1,3) /
+C    S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 7,13,IC),IC=1,3) /
-     S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/
-      DATA (GB( 7,13,IC),IC=1,3) /
-     S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/
-      DATA (GA( 7,14,IC),IC=1,3) /
-     S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/
-      DATA (GB( 7,14,IC),IC=1,3) /
-     S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/
+C     DATA (GA( 7,13,IC),IC=1,3) /
+C    S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/
+C     DATA (GB( 7,13,IC),IC=1,3) /
+C    S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/
+C     DATA (GA( 7,14,IC),IC=1,3) /
+C    S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/
+C     DATA (GB( 7,14,IC),IC=1,3) /
+C    S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 8,13,IC),IC=1,3) /
-     S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/
-      DATA (GB( 8,13,IC),IC=1,3) /
-     S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/
-      DATA (GA( 8,14,IC),IC=1,3) /
-     S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/
-      DATA (GB( 8,14,IC),IC=1,3) /
-     S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/
+C     DATA (GA( 8,13,IC),IC=1,3) /
+C    S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/
+C     DATA (GB( 8,13,IC),IC=1,3) /
+C    S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/
+C     DATA (GA( 8,14,IC),IC=1,3) /
+C    S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/
+C     DATA (GB( 8,14,IC),IC=1,3) /
+C    S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA( 9,13,IC),IC=1,3) /
-     S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/
-      DATA (GB( 9,13,IC),IC=1,3) /
-     S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/
-      DATA (GA( 9,14,IC),IC=1,3) /
-     S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/
-      DATA (GB( 9,14,IC),IC=1,3) /
-     S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/
+C     DATA (GA( 9,13,IC),IC=1,3) /
+C    S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/
+C     DATA (GB( 9,13,IC),IC=1,3) /
+C    S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/
+C     DATA (GA( 9,14,IC),IC=1,3) /
+C    S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/
+C     DATA (GB( 9,14,IC),IC=1,3) /
+C    S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA(10,13,IC),IC=1,3) /
-     S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/
-      DATA (GB(10,13,IC),IC=1,3) /
-     S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/
-      DATA (GA(10,14,IC),IC=1,3) /
-     S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/
-      DATA (GB(10,14,IC),IC=1,3) /
-     S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/
+C     DATA (GA(10,13,IC),IC=1,3) /
+C    S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/
+C     DATA (GB(10,13,IC),IC=1,3) /
+C    S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/
+C     DATA (GA(10,14,IC),IC=1,3) /
+C    S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/
+C     DATA (GB(10,14,IC),IC=1,3) /
+C    S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/
 C
 C----- INTERVAL = 2 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
-      DATA (GA(11,13,IC),IC=1,3) /
-     S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/
-      DATA (GB(11,13,IC),IC=1,3) /
-     S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/
-      DATA (GA(11,14,IC),IC=1,3) /
-     S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/
-      DATA (GB(11,14,IC),IC=1,3) /
-     S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/
+C     DATA (GA(11,13,IC),IC=1,3) /
+C    S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/
+C     DATA (GB(11,13,IC),IC=1,3) /
+C    S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/
+C     DATA (GA(11,14,IC),IC=1,3) /
+C    S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/
+C     DATA (GB(11,14,IC),IC=1,3) /
+C    S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/
 C
 C
@@ -4618,147 +4438,148 @@
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 1,15,IC),IC=1,3) /
-     S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/
-      DATA (GB( 1,15,IC),IC=1,3) /
-     S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/
-      DATA (GA( 1,16,IC),IC=1,3) /
-     S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/
-      DATA (GB( 1,16,IC),IC=1,3) /
-     S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/
+C     DATA (GA( 1,15,IC),IC=1,3) /
+C    S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/
+C     DATA (GB( 1,15,IC),IC=1,3) /
+C    S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/
+C     DATA (GA( 1,16,IC),IC=1,3) /
+C    S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/
+C     DATA (GB( 1,16,IC),IC=1,3) /
+C    S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  200.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 2,15,IC),IC=1,3) /
-     S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/
-      DATA (GB( 2,15,IC),IC=1,3) /
-     S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/
-      DATA (GA( 2,16,IC),IC=1,3) /
-     S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/
-      DATA (GB( 2,16,IC),IC=1,3) /
-     S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/
+C     DATA (GA( 2,15,IC),IC=1,3) /
+C    S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/
+C     DATA (GB( 2,15,IC),IC=1,3) /
+C    S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/
+C     DATA (GA( 2,16,IC),IC=1,3) /
+C    S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/
+C     DATA (GB( 2,16,IC),IC=1,3) /
+C    S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  212.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 3,15,IC),IC=1,3) /
-     S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/
-      DATA (GB( 3,15,IC),IC=1,3) /
-     S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/
-      DATA (GA( 3,16,IC),IC=1,3) /
-     S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/
-      DATA (GB( 3,16,IC),IC=1,3) /
-     S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/
+C     DATA (GA( 3,15,IC),IC=1,3) /
+C    S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/
+C     DATA (GB( 3,15,IC),IC=1,3) /
+C    S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/
+C     DATA (GA( 3,16,IC),IC=1,3) /
+C    S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/
+C     DATA (GB( 3,16,IC),IC=1,3) /
+C    S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  225.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 4,15,IC),IC=1,3) /
-     S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/
-      DATA (GB( 4,15,IC),IC=1,3) /
-     S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/
-      DATA (GA( 4,16,IC),IC=1,3) /
-     S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/
-      DATA (GB( 4,16,IC),IC=1,3) /
-     S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/
+C     DATA (GA( 4,15,IC),IC=1,3) /
+C    S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/
+C     DATA (GB( 4,15,IC),IC=1,3) /
+C    S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/
+C     DATA (GA( 4,16,IC),IC=1,3) /
+C    S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/
+C     DATA (GB( 4,16,IC),IC=1,3) /
+C    S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  237.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 5,15,IC),IC=1,3) /
-     S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/
-      DATA (GB( 5,15,IC),IC=1,3) /
-     S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/
-      DATA (GA( 5,16,IC),IC=1,3) /
-     S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/
-      DATA (GB( 5,16,IC),IC=1,3) /
-     S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/
+C     DATA (GA( 5,15,IC),IC=1,3) /
+C    S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/
+C     DATA (GB( 5,15,IC),IC=1,3) /
+C    S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/
+C     DATA (GA( 5,16,IC),IC=1,3) /
+C    S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/
+C     DATA (GB( 5,16,IC),IC=1,3) /
+C    S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  250.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 6,15,IC),IC=1,3) /
-     S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/
-      DATA (GB( 6,15,IC),IC=1,3) /
-     S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/
-      DATA (GA( 6,16,IC),IC=1,3) /
-     S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/
-      DATA (GB( 6,16,IC),IC=1,3) /
-     S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/
+C     DATA (GA( 6,15,IC),IC=1,3) /
+C    S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/
+C     DATA (GB( 6,15,IC),IC=1,3) /
+C    S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/
+C     DATA (GA( 6,16,IC),IC=1,3) /
+C    S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/
+C     DATA (GB( 6,16,IC),IC=1,3) /
+C    S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  262.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 7,15,IC),IC=1,3) /
-     S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/
-      DATA (GB( 7,15,IC),IC=1,3) /
-     S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/
-      DATA (GA( 7,16,IC),IC=1,3) /
-     S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/
-      DATA (GB( 7,16,IC),IC=1,3) /
-     S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/
+C     DATA (GA( 7,15,IC),IC=1,3) /
+C    S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/
+C     DATA (GB( 7,15,IC),IC=1,3) /
+C    S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/
+C     DATA (GA( 7,16,IC),IC=1,3) /
+C    S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/
+C     DATA (GB( 7,16,IC),IC=1,3) /
+C    S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  275.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 8,15,IC),IC=1,3) /
-     S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/
-      DATA (GB( 8,15,IC),IC=1,3) /
-     S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/
-      DATA (GA( 8,16,IC),IC=1,3) /
-     S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/
-      DATA (GB( 8,16,IC),IC=1,3) /
-     S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/
+C     DATA (GA( 8,15,IC),IC=1,3) /
+C    S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/
+C     DATA (GB( 8,15,IC),IC=1,3) /
+C    S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/
+C     DATA (GA( 8,16,IC),IC=1,3) /
+C    S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/
+C     DATA (GB( 8,16,IC),IC=1,3) /
+C    S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  287.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA( 9,15,IC),IC=1,3) /
-     S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/
-      DATA (GB( 9,15,IC),IC=1,3) /
-     S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/
-      DATA (GA( 9,16,IC),IC=1,3) /
-     S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/
-      DATA (GB( 9,16,IC),IC=1,3) /
-     S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/
+C     DATA (GA( 9,15,IC),IC=1,3) /
+C    S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/
+C     DATA (GB( 9,15,IC),IC=1,3) /
+C    S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/
+C     DATA (GA( 9,16,IC),IC=1,3) /
+C    S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/
+C     DATA (GB( 9,16,IC),IC=1,3) /
+C    S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  300.0
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA(10,15,IC),IC=1,3) /
-     S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/
-      DATA (GB(10,15,IC),IC=1,3) /
-     S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/
-      DATA (GA(10,16,IC),IC=1,3) /
-     S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/
-      DATA (GB(10,16,IC),IC=1,3) /
-     S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/
+C     DATA (GA(10,15,IC),IC=1,3) /
+C    S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/
+C     DATA (GB(10,15,IC),IC=1,3) /
+C    S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/
+C     DATA (GA(10,16,IC),IC=1,3) /
+C    S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/
+C     DATA (GB(10,16,IC),IC=1,3) /
+C    S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/
 C
 C----- INTERVAL = 4 ----- T =  312.5
 C
 C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
-      DATA (GA(11,15,IC),IC=1,3) /
-     S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/
-      DATA (GB(11,15,IC),IC=1,3) /
-     S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/
-      DATA (GA(11,16,IC),IC=1,3) /
-     S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/
-      DATA (GB(11,16,IC),IC=1,3) /
-     S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/
+C     DATA (GA(11,15,IC),IC=1,3) /
+C    S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/
+C     DATA (GB(11,15,IC),IC=1,3) /
+C    S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/
+C     DATA (GA(11,16,IC),IC=1,3) /
+C    S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/
+C     DATA (GB(11,16,IC),IC=1,3) /
+C    S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/
+C
+C     ------------------------------------------------------------------
+C     DATA (( XP(  J,K),J=1,6),       K=1,6) /
+C    S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,
+C    S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,
+C    S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,
+C    S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,
+C    S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,
+C    S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,
+C    S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,
+C    S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,
+C    S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,
+C    S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,
+C    S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,
+C    S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /
 
-C     ------------------------------------------------------------------
-      DATA (( XP(  J,K),J=1,6),       K=1,6) /
-     S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,
-     S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,
-     S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,
-     S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,
-     S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,
-     S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,
-     S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,
-     S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,
-     S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,
-     S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,
-     S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,
-     S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /
 C
 C
@@ -4768,4 +4589,5 @@
  100  CONTINUE
 C
+!cdir collapse
       DO 102 JK = 1 , KFLEV+1
       DO 101 JL = 1, KDLON
@@ -4959,31 +4781,31 @@
       INTEGER KUAER,KTRAER, KLIM
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
-      REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
-      REAL*8 PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
-      REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
-      REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
-      REAL*8 PBTOP(KDLON,Ninter) ! T.O.A. SPECTRAL PLANCK FUNCTION
-      REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
-      REAL*8 PEMIS(KDLON) ! SURFACE EMISSIVITY
-      REAL*8 PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
-      REAL*8 PTAVE(KDLON,KFLEV) ! TEMPERATURE
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGASUR(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGBSUR(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGATOP(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGBTOP(KDLON,8,2) ! PADE APPROXIMANTS
-C
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
-      REAL*8 PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
-      REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
+      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
+      REAL(KIND=8) PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
+      REAL(KIND=8) PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
+      REAL(KIND=8) PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
+      REAL(KIND=8) PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
+      REAL(KIND=8) PBTOP(KDLON,Ninter) ! T.O.A. SPECTRAL PLANCK FUNCTION
+      REAL(KIND=8) PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
+      REAL(KIND=8) PEMIS(KDLON) ! SURFACE EMISSIVITY
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
+      REAL(KIND=8) PTAVE(KDLON,KFLEV) ! TEMPERATURE
+      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGASUR(KDLON,8,2) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGBSUR(KDLON,8,2) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGATOP(KDLON,8,2) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGBTOP(KDLON,8,2) ! PADE APPROXIMANTS
+C
+      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
+      REAL(KIND=8) PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
+      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
 C-----------------------------------------------------------------------
 C LOCAL VARIABLES:
-      REAL*8 ZADJD(KDLON,KFLEV+1)
-      REAL*8 ZADJU(KDLON,KFLEV+1)
-      REAL*8 ZDBDT(KDLON,Ninter,KFLEV)
-      REAL*8 ZDISD(KDLON,KFLEV+1)
-      REAL*8 ZDISU(KDLON,KFLEV+1)
+      REAL(KIND=8) ZADJD(KDLON,KFLEV+1)
+      REAL(KIND=8) ZADJU(KDLON,KFLEV+1)
+      REAL(KIND=8) ZDBDT(KDLON,Ninter,KFLEV)
+      REAL(KIND=8) ZDISD(KDLON,KFLEV+1)
+      REAL(KIND=8) ZDISU(KDLON,KFLEV+1)
 C
       INTEGER jk, jl
@@ -5076,43 +4898,43 @@
       INTEGER KUAER,KTRAER, KLIM
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
-      REAL*8 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
-      REAL*8 PADJU(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
-      REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
-      REAL*8 PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
-      REAL*8 PBSUR(KDLON,Ninter) ! SPECTRAL SURFACE PLANCK FUNCTION
-      REAL*8 PBSUI(KDLON) ! SURFACE PLANCK FUNCTION
-      REAL*8 PBTOP(KDLON,Ninter) ! SPECTRAL T.O.A. PLANCK FUNCTION
-      REAL*8 PDISD(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
-      REAL*8 PDISU(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
-      REAL*8 PEMIS(KDLON) ! SURFACE EMISSIVITY
-      REAL*8 PPMB(KDLON,KFLEV+1) ! PRESSURE MB
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
-      REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
-      REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
-      REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
-C
-      REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
-      REAL*8 PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
+      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
+      REAL(KIND=8) PADJD(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
+      REAL(KIND=8) PADJU(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
+      REAL(KIND=8) PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
+      REAL(KIND=8) PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
+      REAL(KIND=8) PBSUR(KDLON,Ninter) ! SPECTRAL SURFACE PLANCK FUNCTION
+      REAL(KIND=8) PBSUI(KDLON) ! SURFACE PLANCK FUNCTION
+      REAL(KIND=8) PBTOP(KDLON,Ninter) ! SPECTRAL T.O.A. PLANCK FUNCTION
+      REAL(KIND=8) PDISD(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
+      REAL(KIND=8) PDISU(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
+      REAL(KIND=8) PEMIS(KDLON) ! SURFACE EMISSIVITY
+      REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! PRESSURE MB
+      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
+      REAL(KIND=8) PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
+      REAL(KIND=8) PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
+      REAL(KIND=8) PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
+C
+      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
+      REAL(KIND=8) PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZBGND(KDLON)
-      REAL*8 ZFD(KDLON)
-      REAL*8  ZFN10(KDLON)
-      REAL*8 ZFU(KDLON)
-      REAL*8  ZTT(KDLON,NTRA)
-      REAL*8 ZTT1(KDLON,NTRA)
-      REAL*8 ZTT2(KDLON,NTRA)
-      REAL*8  ZUU(KDLON,NUA) 
-      REAL*8 ZCNSOL(KDLON)
-      REAL*8 ZCNTOP(KDLON)
+      REAL(KIND=8) ZBGND(KDLON)
+      REAL(KIND=8) ZFD(KDLON)
+      REAL(KIND=8)  ZFN10(KDLON)
+      REAL(KIND=8) ZFU(KDLON)
+      REAL(KIND=8)  ZTT(KDLON,NTRA)
+      REAL(KIND=8) ZTT1(KDLON,NTRA)
+      REAL(KIND=8) ZTT2(KDLON,NTRA)
+      REAL(KIND=8)  ZUU(KDLON,NUA) 
+      REAL(KIND=8) ZCNSOL(KDLON)
+      REAL(KIND=8) ZCNTOP(KDLON)
 C
       INTEGER jk, jl, ja
       INTEGER jstra, jstru
       INTEGER ind1, ind2, ind3, ind4, in, jlim
-      REAL*8 zctstr
+      REAL(KIND=8) zctstr
 C-----------------------------------------------------------------------
 C
@@ -5401,25 +5223,25 @@
       INTEGER KUAER,KTRAER
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
-      REAL*8 PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-C
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! ENERGY EXCHANGE MATRIX
-      REAL*8 PDISD(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
-      REAL*8 PDISU(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
+      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
+      REAL(KIND=8) PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT
+      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+C
+      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! ENERGY EXCHANGE MATRIX
+      REAL(KIND=8) PDISD(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
+      REAL(KIND=8) PDISU(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 ZGLAYD(KDLON)
-      REAL*8 ZGLAYU(KDLON)
-      REAL*8 ZTT(KDLON,NTRA)
-      REAL*8 ZTT1(KDLON,NTRA)
-      REAL*8 ZTT2(KDLON,NTRA)
+      REAL(KIND=8) ZGLAYD(KDLON)
+      REAL(KIND=8) ZGLAYU(KDLON)
+      REAL(KIND=8) ZTT(KDLON,NTRA)
+      REAL(KIND=8) ZTT1(KDLON,NTRA)
+      REAL(KIND=8) ZTT2(KDLON,NTRA)
 C
       INTEGER jl, jk, ja, ikp1, ikn, ikd1, jkj, ikd2
       INTEGER ikjp1, ikm1, ikj, jlk, iku1, ijkl, iku2
       INTEGER ind1, ind2, ind3, ind4, itt
-      REAL*8 zww, zdzxdg, zdzxmg
+      REAL(KIND=8) zww, zdzxdg, zdzxmg
 C
 C*         1.    INITIALIZATION
@@ -5619,4 +5441,5 @@
      S  , PADJD,PADJU,PCNTRB,PDBDT)
        USE dimphy
+      USE radiation_AR4_param, only : WG1
       IMPLICIT none
 cym#include "dimensions.h"
@@ -5656,33 +5479,28 @@
       INTEGER KUAER,KTRAER
 C
-      REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
-      REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
-      REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
-C
-      REAL*8 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
-      REAL*8 PADJU(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
-      REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
-      REAL*8 PDBDT(KDLON,Ninter,KFLEV) !  LAYER PLANCK FUNCTION GRADIENT
+      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
+      REAL(KIND=8) PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
+      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
+C
+      REAL(KIND=8) PADJD(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
+      REAL(KIND=8) PADJU(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
+      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
+      REAL(KIND=8) PDBDT(KDLON,Ninter,KFLEV) !  LAYER PLANCK FUNCTION GRADIENT
 C
 C* LOCAL ARRAYS:
 C
-      REAL*8 ZGLAYD(KDLON)
-      REAL*8 ZGLAYU(KDLON)
-      REAL*8 ZTT(KDLON,NTRA)
-      REAL*8 ZTT1(KDLON,NTRA)
-      REAL*8 ZTT2(KDLON,NTRA)
-      REAL*8 ZUU(KDLON,NUA)
+      REAL(KIND=8) ZGLAYD(KDLON)
+      REAL(KIND=8) ZGLAYU(KDLON)
+      REAL(KIND=8) ZTT(KDLON,NTRA)
+      REAL(KIND=8) ZTT1(KDLON,NTRA)
+      REAL(KIND=8) ZTT2(KDLON,NTRA)
+      REAL(KIND=8) ZUU(KDLON,NUA)
 C
       INTEGER jk, jl, ja, im12, ind, inu, ixu, jg
       INTEGER ixd, ibs, idd, imu, jk1, jk2, jnu
-      REAL*8 zwtr
+      REAL(KIND=8) zwtr
 c
-C* Data Block:
-c
-      REAL*8 WG1(2)
-      SAVE WG1
-c$OMP THREADPRIVATE(WG1)
-      DATA (WG1(jk),jk=1,2) /1.0, 1.0/
+
 C-----------------------------------------------------------------------
 C
@@ -5865,27 +5683,27 @@
 C
 C-----------------------------------------------------------------------
-      REAL*8 O1H, O2H
+      REAL(KIND=8) O1H, O2H
       PARAMETER (O1H=2230.)
       PARAMETER (O2H=100.)
-      REAL*8 RPIALF0
+      REAL(KIND=8) RPIALF0
       PARAMETER (RPIALF0=2.0)
 C
 C* ARGUMENTS:
 C
-      REAL*8 PUU(KDLON,NUA)
-      REAL*8 PTT(KDLON,NTRA)
-      REAL*8 PGA(KDLON,8,2)
-      REAL*8 PGB(KDLON,8,2)
+      REAL(KIND=8) PUU(KDLON,NUA)
+      REAL(KIND=8) PTT(KDLON,NTRA)
+      REAL(KIND=8) PGA(KDLON,8,2)
+      REAL(KIND=8) PGB(KDLON,8,2)
 C
 C* LOCAL VARIABLES:
 C
-      REAL*8 zz, zxd, zxn
-      REAL*8 zpu, zpu10, zpu11, zpu12, zpu13
-      REAL*8 zeu, zeu10, zeu11, zeu12, zeu13
-      REAL*8 zx, zy, zsq1, zsq2, zvxy, zuxy
-      REAL*8 zaercn, zto1, zto2, zxch4, zych4, zxn2o, zyn2o
-      REAL*8 zsqn21, zodn21, zsqh42, zodh42
-      REAL*8 zsqh41, zodh41, zsqn22, zodn22, zttf11, zttf12
-      REAL*8 zuu11, zuu12, za11, za12
+      REAL(KIND=8) zz, zxd, zxn
+      REAL(KIND=8) zpu, zpu10, zpu11, zpu12, zpu13
+      REAL(KIND=8) zeu, zeu10, zeu11, zeu12, zeu13
+      REAL(KIND=8) zx, zy, zsq1, zsq2, zvxy, zuxy
+      REAL(KIND=8) zaercn, zto1, zto2, zxch4, zych4, zxn2o, zyn2o
+      REAL(KIND=8) zsqn21, zodn21, zsqh42, zodh42
+      REAL(KIND=8) zsqh41, zodh41, zsqn22, zodn22, zttf11, zttf12
+      REAL(KIND=8) zuu11, zuu12, za11, za12
       INTEGER jl, ja
 C     ------------------------------------------------------------------
@@ -5897,4 +5715,5 @@
 C
 C
+!cdir collapse
       DO 130 JA = 1 , 8
       DO 120 JL = 1, KDLON
@@ -6041,29 +5860,29 @@
 C
 C-----------------------------------------------------------------------
-      REAL*8 O1H, O2H
+      REAL(KIND=8) O1H, O2H
       PARAMETER (O1H=2230.)
       PARAMETER (O2H=100.)
-      REAL*8 RPIALF0
+      REAL(KIND=8) RPIALF0
       PARAMETER (RPIALF0=2.0)
 C
 C* ARGUMENTS:
 C
-      REAL*8 PGA(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PGB(KDLON,8,2) ! PADE APPROXIMANTS
-      REAL*8 PUU1(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 1
-      REAL*8 PUU2(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 2
-      REAL*8 PTT(KDLON,NTRA) ! TRANSMISSION FUNCTIONS
+      REAL(KIND=8) PGA(KDLON,8,2) ! PADE APPROXIMANTS
+      REAL(KIND=8) PGB(KDLON,8,2) ! PADE APPROXIMANTS
+      REAL(KIND=8) PUU1(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 1
+      REAL(KIND=8) PUU2(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 2
+      REAL(KIND=8) PTT(KDLON,NTRA) ! TRANSMISSION FUNCTIONS
 C
 C* LOCAL VARIABLES:
 C
       INTEGER ja, jl
-      REAL*8 zz, zxd, zxn
-      REAL*8 zpu, zpu10, zpu11, zpu12, zpu13
-      REAL*8 zeu, zeu10, zeu11, zeu12, zeu13
-      REAL*8 zx, zy, zuxy, zsq1, zsq2, zvxy, zaercn, zto1, zto2
-      REAL*8 zxch4, zych4, zsqh41, zodh41
-      REAL*8 zxn2o, zyn2o, zsqn21, zodn21, zsqh42, zodh42
-      REAL*8 zsqn22, zodn22, za11, zttf11, za12, zttf12
-      REAL*8 zuu11, zuu12
+      REAL(KIND=8) zz, zxd, zxn
+      REAL(KIND=8) zpu, zpu10, zpu11, zpu12, zpu13
+      REAL(KIND=8) zeu, zeu10, zeu11, zeu12, zeu13
+      REAL(KIND=8) zx, zy, zuxy, zsq1, zsq2, zvxy, zaercn, zto1, zto2
+      REAL(KIND=8) zxch4, zych4, zsqh41, zodh41
+      REAL(KIND=8) zxn2o, zyn2o, zsqn21, zodn21, zsqh42, zodh42
+      REAL(KIND=8) zsqn22, zodn22, za11, zttf11, za12, zttf12
+      REAL(KIND=8) zuu11, zuu12
 C     ------------------------------------------------------------------
 C
@@ -6074,4 +5893,8 @@
 C
 C
+
+!CDIR ON_ADB(PUU1)
+!CDIR ON_ADB(PUU2)
+!CDIR COLLAPSE
       DO 130 JA = 1 , 8
       DO 120 JL = 1, KDLON
