Ignore:
Timestamp:
Dec 10, 2009, 10:02:56 AM (15 years ago)
Author:
Laurent Fairhead
Message:

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

Location:
LMDZ4/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk

  • LMDZ4/trunk/libf/phylmd/radiation_AR4.F

    r1107 r1279  
    5454C* ARGUMENTS:
    5555C
    56       REAL*8 PSCT  ! constante solaire (valeur conseillee: 1370)
    57 cIM ctes ds clesphys.h   REAL*8 RCO2  ! concentration CO2 (IPCC: 353.E-06*44.011/28.97)
     56      REAL(KIND=8) PSCT  ! constante solaire (valeur conseillee: 1370)
     57cIM ctes ds clesphys.h   REAL(KIND=8) RCO2  ! concentration CO2 (IPCC: 353.E-06*44.011/28.97)
    5858#include "clesphys.h"
    5959C
    60       REAL*8 PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
    61       REAL*8 PDP(KDLON,KFLEV)    ! LAYER THICKNESS (PA)
    62       REAL*8 PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
    63 C
    64       REAL*8 PRMU0(KDLON)  ! COSINE OF ZENITHAL ANGLE
    65       REAL*8 PFRAC(KDLON)  ! fraction de la journee
    66 C
    67       REAL*8 PTAVE(KDLON,KFLEV)  ! LAYER TEMPERATURE (K)
    68       REAL*8 PWV(KDLON,KFLEV)    ! SPECIFIC HUMIDITY (KG/KG)
    69       REAL*8 PQS(KDLON,KFLEV)    ! SATURATED WATER VAPOUR (KG/KG)
    70       REAL*8 POZON(KDLON,KFLEV)  ! OZONE CONCENTRATION (KG/KG)
    71       REAL*8 PAER(KDLON,KFLEV,5) ! AEROSOLS' OPTICAL THICKNESS
    72 C
    73       REAL*8 PALBD(KDLON,2)  ! albedo du sol (lumiere diffuse)
    74       REAL*8 PALBP(KDLON,2)  ! albedo du sol (lumiere parallele)
    75 C
    76       REAL*8 PCLDSW(KDLON,KFLEV)    ! CLOUD FRACTION
    77       REAL*8 PTAU(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS
    78       REAL*8 PCG(KDLON,2,KFLEV)     ! ASYMETRY FACTOR
    79       REAL*8 POMEGA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
    80 C
    81       REAL*8 PHEAT(KDLON,KFLEV) ! SHORTWAVE HEATING (K/DAY)
    82       REAL*8 PHEAT0(KDLON,KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
    83       REAL*8 PALBPLA(KDLON)     ! PLANETARY ALBEDO
    84       REAL*8 PTOPSW(KDLON)      ! SHORTWAVE FLUX AT T.O.A.
    85       REAL*8 PSOLSW(KDLON)      ! SHORTWAVE FLUX AT SURFACE
    86       REAL*8 PTOPSW0(KDLON)     ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
    87       REAL*8 PSOLSW0(KDLON)     ! SHORTWAVE FLUX AT SURFACE (CLEAR-SKY)
     60      REAL(KIND=8) PPSOL(KDLON)        ! SURFACE PRESSURE (PA)
     61      REAL(KIND=8) PDP(KDLON,KFLEV)    ! LAYER THICKNESS (PA)
     62      REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
     63C
     64      REAL(KIND=8) PRMU0(KDLON)  ! COSINE OF ZENITHAL ANGLE
     65      REAL(KIND=8) PFRAC(KDLON)  ! fraction de la journee
     66C
     67      REAL(KIND=8) PTAVE(KDLON,KFLEV)  ! LAYER TEMPERATURE (K)
     68      REAL(KIND=8) PWV(KDLON,KFLEV)    ! SPECIFIC HUMIDITY (KG/KG)
     69      REAL(KIND=8) PQS(KDLON,KFLEV)    ! SATURATED WATER VAPOUR (KG/KG)
     70      REAL(KIND=8) POZON(KDLON,KFLEV)  ! OZONE CONCENTRATION (KG/KG)
     71      REAL(KIND=8) PAER(KDLON,KFLEV,5) ! AEROSOLS' OPTICAL THICKNESS
     72C
     73      REAL(KIND=8) PALBD(KDLON,2)  ! albedo du sol (lumiere diffuse)
     74      REAL(KIND=8) PALBP(KDLON,2)  ! albedo du sol (lumiere parallele)
     75C
     76      REAL(KIND=8) PCLDSW(KDLON,KFLEV)    ! CLOUD FRACTION
     77      REAL(KIND=8) PTAU(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS
     78      REAL(KIND=8) PCG(KDLON,2,KFLEV)     ! ASYMETRY FACTOR
     79      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
     80C
     81      REAL(KIND=8) PHEAT(KDLON,KFLEV) ! SHORTWAVE HEATING (K/DAY)
     82      REAL(KIND=8) PHEAT0(KDLON,KFLEV)! SHORTWAVE HEATING (K/DAY) clear-sky
     83      REAL(KIND=8) PALBPLA(KDLON)     ! PLANETARY ALBEDO
     84      REAL(KIND=8) PTOPSW(KDLON)      ! SHORTWAVE FLUX AT T.O.A.
     85      REAL(KIND=8) PSOLSW(KDLON)      ! SHORTWAVE FLUX AT SURFACE
     86      REAL(KIND=8) PTOPSW0(KDLON)     ! SHORTWAVE FLUX AT T.O.A. (CLEAR-SKY)
     87      REAL(KIND=8) PSOLSW0(KDLON)     ! SHORTWAVE FLUX AT SURFACE (CLEAR-SKY)
    8888C
    8989C* LOCAL VARIABLES:
    9090C
    91       REAL*8 ZOZ(KDLON,KFLEV)
    92       REAL*8 ZAKI(KDLON,2)     
    93       REAL*8 ZCLD(KDLON,KFLEV)
    94       REAL*8 ZCLEAR(KDLON)
    95       REAL*8 ZDSIG(KDLON,KFLEV)
    96       REAL*8 ZFACT(KDLON)
    97       REAL*8 ZFD(KDLON,KFLEV+1)
    98       REAL*8 ZFDOWN(KDLON,KFLEV+1)
    99       REAL*8 ZFU(KDLON,KFLEV+1)
    100       REAL*8 ZFUP(KDLON,KFLEV+1)
    101       REAL*8 ZRMU(KDLON)
    102       REAL*8 ZSEC(KDLON)
    103       REAL*8 ZUD(KDLON,5,KFLEV+1)
    104       REAL*8 ZCLDSW0(KDLON,KFLEV)
     91      real, parameter:: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
     92
     93      REAL(KIND=8) ZOZ(KDLON,KFLEV)
     94!     column-density of ozone in layer, in kilo-Dobsons
     95
     96      REAL(KIND=8) ZAKI(KDLON,2)     
     97      REAL(KIND=8) ZCLD(KDLON,KFLEV)
     98      REAL(KIND=8) ZCLEAR(KDLON)
     99      REAL(KIND=8) ZDSIG(KDLON,KFLEV)
     100      REAL(KIND=8) ZFACT(KDLON)
     101      REAL(KIND=8) ZFD(KDLON,KFLEV+1)
     102      REAL(KIND=8) ZFDOWN(KDLON,KFLEV+1)
     103      REAL(KIND=8) ZFU(KDLON,KFLEV+1)
     104      REAL(KIND=8) ZFUP(KDLON,KFLEV+1)
     105      REAL(KIND=8) ZRMU(KDLON)
     106      REAL(KIND=8) ZSEC(KDLON)
     107      REAL(KIND=8) ZUD(KDLON,5,KFLEV+1)
     108      REAL(KIND=8) ZCLDSW0(KDLON,KFLEV)
    105109c
    106       REAL*8 ZFSUP(KDLON,KFLEV+1)
    107       REAL*8 ZFSDN(KDLON,KFLEV+1)
    108       REAL*8 ZFSUP0(KDLON,KFLEV+1)
    109       REAL*8 ZFSDN0(KDLON,KFLEV+1)
     110      REAL(KIND=8) ZFSUP(KDLON,KFLEV+1)
     111      REAL(KIND=8) ZFSDN(KDLON,KFLEV+1)
     112      REAL(KIND=8) ZFSUP0(KDLON,KFLEV+1)
     113      REAL(KIND=8) ZFSDN0(KDLON,KFLEV+1)
    110114C
    111115      INTEGER inu, jl, jk, i, k, kpl1
    … …  
    122126c$OMP THREADPRIVATE(itapsw)
    123127cjq-Introduced for aerosol forcings
    124       real*8 flag_aer
     128      real(kind=8) flag_aer
    125129      logical ok_ade, ok_aie    ! use aerosol forcings or not?
    126       real*8 tauae(kdlon,kflev,2)  ! aerosol optical properties
    127       real*8 pizae(kdlon,kflev,2)  ! (see aeropt.F)
    128       real*8 cgae(kdlon,kflev,2)   ! -"-
    129       REAL*8 PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (pre-industrial value)
    130       REAL*8 POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
    131       REAL*8 PTOPSWAD(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
    132       REAL*8 PSOLSWAD(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
    133       REAL*8 PTOPSWAI(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND)
    134       REAL*8 PSOLSWAI(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND)
     130      real(kind=8) tauae(kdlon,kflev,2)  ! aerosol optical properties
     131      real(kind=8) pizae(kdlon,kflev,2)  ! (see aeropt.F)
     132      real(kind=8) cgae(kdlon,kflev,2)   ! -"-
     133      REAL(KIND=8) PTAUA(KDLON,2,KFLEV)    ! CLOUD OPTICAL THICKNESS (pre-industrial value)
     134      REAL(KIND=8) POMEGAA(KDLON,2,KFLEV)  ! SINGLE SCATTERING ALBEDO
     135      REAL(KIND=8) PTOPSWAD(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL DIR)
     136      REAL(KIND=8) PSOLSWAD(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL DIR)
     137      REAL(KIND=8) PTOPSWAI(KDLON)     ! SHORTWAVE FLUX AT T.O.A.(+AEROSOL IND)
     138      REAL(KIND=8) PSOLSWAI(KDLON)     ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND)
    135139cjq - Fluxes including aerosol effects
    136       REAL*8,allocatable,save :: ZFSUPAD(:,:)
     140      REAL(KIND=8),allocatable,save :: ZFSUPAD(:,:)
    137141c$OMP THREADPRIVATE(ZFSUPAD)
    138       REAL*8,allocatable,save :: ZFSDNAD(:,:)
     142      REAL(KIND=8),allocatable,save :: ZFSDNAD(:,:)
    139143c$OMP THREADPRIVATE(ZFSDNAD)
    140       REAL*8,allocatable,save :: ZFSUPAI(:,:)
     144      REAL(KIND=8),allocatable,save :: ZFSUPAI(:,:)
    141145c$OMP THREADPRIVATE(ZFSUPAI)
    142       REAL*8,allocatable,save :: ZFSDNAI(:,:)
     146      REAL(KIND=8),allocatable,save :: ZFSDNAI(:,:)
    143147c$OMP THREADPRIVATE(ZFSDNAI)
    144148      logical initialized
    … …  
    151155c$OMP THREADPRIVATE(initialized)
    152156cjq-end
     157      REAL tmp_
    153158      if(.not.initialized) then
    154159        flag_aer=0.
    … …  
    158163        allocate(ZFSUPAI(KDLON,KFLEV+1))
    159164        allocate(ZFSDNAI(KDLON,KFLEV+1))
    160         ZFSUPAD(:,:)=0.
    161         ZFSDNAD(:,:)=0.
    162         ZFSUPAI(:,:)=0.
    163         ZFSDNAI(:,:)=0.
    164        
     165        DO JK = 1 , KDLON*(KFLEV+1)
     166          ZFSUPAD(JK,1) = 0.0     ! ZFSUPAD(:,:)=0.
     167          ZFSDNAD(JK,1) = 0.0     ! ZFSDNAD(:,:)=0.
     168          ZFSUPAI(JK,1) = 0.0     ! ZFSUPAI(:,:)=0.
     169          ZFSDNAI(JK,1) = 0.0     ! ZFSDNAI(:,:)=0.
     170        END DO
    165171      endif
    166172!rv
    … …  
    175181      IF (MOD(itapsw,swpas).EQ.0) THEN
    176182c
     183      tmp_ = 1./( dobson_u * 1e3 * RG)
     184!cdir collapse
    177185      DO JK = 1 , KFLEV
    178       DO JL = 1, KDLON
    179          ZCLDSW0(JL,JK) = 0.0
    180          ZOZ(JL,JK) = POZON(JL,JK)*46.6968/RG
    181      .               *PDP(JL,JK)*(101325.0/PPSOL(JL))
    182       ENDDO
     186        DO JL = 1, KDLON
     187          ZCLDSW0(JL,JK) = 0.0
     188          ZOZ(JL,JK) = POZON(JL,JK)*tmp_*PDP(JL,JK)
     189        ENDDO
    183190      ENDDO
    184191C
    … …  
    339346     S                PRMU,PSEC,PUD)
    340347      USE dimphy
     348      USE radiation_AR4_param, only :
     349     S     ZPDH2O,ZPDUMG,ZPRH2O,ZPRUMG,RTDH2O,RTDUMG,RTH2O,RTUMG
    341350      IMPLICIT none
    342351cym#include "dimensions.h"
    … …  
    349358C* ARGUMENTS:
    350359C
    351       REAL*8 PSCT
    352 cIM ctes ds clesphys.h   REAL*8 RCO2
     360      REAL(KIND=8) PSCT
     361cIM ctes ds clesphys.h   REAL(KIND=8) RCO2
    353362#include "clesphys.h"
    354       REAL*8 PCLDSW(KDLON,KFLEV)
    355       REAL*8 PPMB(KDLON,KFLEV+1)
    356       REAL*8 PPSOL(KDLON)
    357       REAL*8 PRMU0(KDLON)
    358       REAL*8 PFRAC(KDLON)
    359       REAL*8 PTAVE(KDLON,KFLEV)
    360       REAL*8 PWV(KDLON,KFLEV)
    361 C
    362       REAL*8 PAKI(KDLON,2)
    363       REAL*8 PCLD(KDLON,KFLEV)
    364       REAL*8 PCLEAR(KDLON)
    365       REAL*8 PDSIG(KDLON,KFLEV)
    366       REAL*8 PFACT(KDLON)
    367       REAL*8 PRMU(KDLON)
    368       REAL*8 PSEC(KDLON)
    369       REAL*8 PUD(KDLON,5,KFLEV+1)
     363      REAL(KIND=8) PCLDSW(KDLON,KFLEV)
     364      REAL(KIND=8) PPMB(KDLON,KFLEV+1)
     365      REAL(KIND=8) PPSOL(KDLON)
     366      REAL(KIND=8) PRMU0(KDLON)
     367      REAL(KIND=8) PFRAC(KDLON)
     368      REAL(KIND=8) PTAVE(KDLON,KFLEV)
     369      REAL(KIND=8) PWV(KDLON,KFLEV)
     370C
     371      REAL(KIND=8) PAKI(KDLON,2)
     372      REAL(KIND=8) PCLD(KDLON,KFLEV)
     373      REAL(KIND=8) PCLEAR(KDLON)
     374      REAL(KIND=8) PDSIG(KDLON,KFLEV)
     375      REAL(KIND=8) PFACT(KDLON)
     376      REAL(KIND=8) PRMU(KDLON)
     377      REAL(KIND=8) PSEC(KDLON)
     378      REAL(KIND=8) PUD(KDLON,5,KFLEV+1)
    370379C
    371380C* LOCAL VARIABLES:
    372381C
    373382      INTEGER IIND(2)
    374       REAL*8 ZC1J(KDLON,KFLEV+1)
    375       REAL*8 ZCLEAR(KDLON)
    376       REAL*8 ZCLOUD(KDLON)
    377       REAL*8 ZN175(KDLON)
    378       REAL*8 ZN190(KDLON)
    379       REAL*8 ZO175(KDLON)
    380       REAL*8 ZO190(KDLON)
    381       REAL*8 ZSIGN(KDLON)
    382       REAL*8 ZR(KDLON,2)
    383       REAL*8 ZSIGO(KDLON)
    384       REAL*8 ZUD(KDLON,2)
    385       REAL*8 ZRTH, ZRTU, ZWH2O, ZDSCO2, ZDSH2O, ZFPPW
     383      REAL(KIND=8) ZC1J(KDLON,KFLEV+1)
     384      REAL(KIND=8) ZCLEAR(KDLON)
     385      REAL(KIND=8) ZCLOUD(KDLON)
     386      REAL(KIND=8) ZN175(KDLON)
     387      REAL(KIND=8) ZN190(KDLON)
     388      REAL(KIND=8) ZO175(KDLON)
     389      REAL(KIND=8) ZO190(KDLON)
     390      REAL(KIND=8) ZSIGN(KDLON)
     391      REAL(KIND=8) ZR(KDLON,2)
     392      REAL(KIND=8) ZSIGO(KDLON)
     393      REAL(KIND=8) ZUD(KDLON,2)
     394      REAL(KIND=8) ZRTH, ZRTU, ZWH2O, ZDSCO2, ZDSH2O, ZFPPW
    386395      INTEGER jl, jk, jkp1, jkl, jklp1, ja
    387396C
    388 C* Prescribed Data:
    389 c
    390       REAL*8 ZPDH2O,ZPDUMG
    391       SAVE ZPDH2O,ZPDUMG
    392 c$OMP THREADPRIVATE(ZPDH2O,ZPDUMG)
    393       REAL*8 ZPRH2O,ZPRUMG
    394       SAVE ZPRH2O,ZPRUMG
    395 c$OMP THREADPRIVATE(ZPRH2O,ZPRUMG)
    396       REAL*8 RTDH2O,RTDUMG
    397       SAVE RTDH2O,RTDUMG
    398 c$OMP THREADPRIVATE(RTDH2O,RTDUMG)
    399       REAL*8 RTH2O ,RTUMG
    400       SAVE RTH2O ,RTUMG
    401 c$OMP THREADPRIVATE(RTH2O ,RTUMG)
    402       DATA ZPDH2O,ZPDUMG / 0.8   , 0.75 /
    403       DATA ZPRH2O,ZPRUMG / 30000., 30000. /
    404       DATA RTDH2O,RTDUMG /  0.40  , 0.375 /
    405       DATA RTH2O ,RTUMG  /  240.  , 240.  /
    406397C     ------------------------------------------------------------------
    407398C
    … …  
    534525     S  ,  PFD   , PFU)
    535526      USE dimphy
     527      USE radiation_AR4_param, only : RSUN, RRAY
    536528      IMPLICIT none
    537529cym#include "dimensions.h"
    … …  
    573565      INTEGER KNU
    574566c-OB
    575       real*8 flag_aer
    576       real*8 tauae(kdlon,kflev,2)
    577       real*8 pizae(kdlon,kflev,2)
    578       real*8 cgae(kdlon,kflev,2)
    579       REAL*8 PAER(KDLON,KFLEV,5)
    580       REAL*8 PALBD(KDLON,2)
    581       REAL*8 PALBP(KDLON,2)
    582       REAL*8 PCG(KDLON,2,KFLEV) 
    583       REAL*8 PCLD(KDLON,KFLEV)
    584       REAL*8 PCLDSW(KDLON,KFLEV)
    585       REAL*8 PCLEAR(KDLON)
    586       REAL*8 PDSIG(KDLON,KFLEV)
    587       REAL*8 POMEGA(KDLON,2,KFLEV)
    588       REAL*8 POZ(KDLON,KFLEV)
    589       REAL*8 PRMU(KDLON)
    590       REAL*8 PSEC(KDLON)
    591       REAL*8 PTAU(KDLON,2,KFLEV)
    592       REAL*8 PUD(KDLON,5,KFLEV+1)
    593 C
    594       REAL*8 PFD(KDLON,KFLEV+1)
    595       REAL*8 PFU(KDLON,KFLEV+1)
     567      real(kind=8) flag_aer
     568      real(kind=8) tauae(kdlon,kflev,2)
     569      real(kind=8) pizae(kdlon,kflev,2)
     570      real(kind=8) cgae(kdlon,kflev,2)
     571      REAL(KIND=8) PAER(KDLON,KFLEV,5)
     572      REAL(KIND=8) PALBD(KDLON,2)
     573      REAL(KIND=8) PALBP(KDLON,2)
     574      REAL(KIND=8) PCG(KDLON,2,KFLEV) 
     575      REAL(KIND=8) PCLD(KDLON,KFLEV)
     576      REAL(KIND=8) PCLDSW(KDLON,KFLEV)
     577      REAL(KIND=8) PCLEAR(KDLON)
     578      REAL(KIND=8) PDSIG(KDLON,KFLEV)
     579      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)
     580      REAL(KIND=8) POZ(KDLON,KFLEV)
     581      REAL(KIND=8) PRMU(KDLON)
     582      REAL(KIND=8) PSEC(KDLON)
     583      REAL(KIND=8) PTAU(KDLON,2,KFLEV)
     584      REAL(KIND=8) PUD(KDLON,5,KFLEV+1)
     585C
     586      REAL(KIND=8) PFD(KDLON,KFLEV+1)
     587      REAL(KIND=8) PFU(KDLON,KFLEV+1)
    596588C
    597589C* LOCAL VARIABLES:
    … …  
    599591      INTEGER IIND(4)
    600592C     
    601       REAL*8 ZCGAZ(KDLON,KFLEV)
    602       REAL*8 ZDIFF(KDLON)
    603       REAL*8 ZDIRF(KDLON)       
    604       REAL*8 ZPIZAZ(KDLON,KFLEV)
    605       REAL*8 ZRAYL(KDLON)
    606       REAL*8 ZRAY1(KDLON,KFLEV+1)
    607       REAL*8 ZRAY2(KDLON,KFLEV+1)
    608       REAL*8 ZREFZ(KDLON,2,KFLEV+1)
    609       REAL*8 ZRJ(KDLON,6,KFLEV+1)
    610       REAL*8 ZRJ0(KDLON,6,KFLEV+1)
    611       REAL*8 ZRK(KDLON,6,KFLEV+1)
    612       REAL*8 ZRK0(KDLON,6,KFLEV+1)
    613       REAL*8 ZRMUE(KDLON,KFLEV+1)
    614       REAL*8 ZRMU0(KDLON,KFLEV+1)
    615       REAL*8 ZR(KDLON,4)
    616       REAL*8 ZTAUAZ(KDLON,KFLEV)
    617       REAL*8 ZTRA1(KDLON,KFLEV+1)
    618       REAL*8 ZTRA2(KDLON,KFLEV+1)
    619       REAL*8 ZW(KDLON,4)
     593      REAL(KIND=8) ZCGAZ(KDLON,KFLEV)
     594      REAL(KIND=8) ZDIFF(KDLON)
     595      REAL(KIND=8) ZDIRF(KDLON)       
     596      REAL(KIND=8) ZPIZAZ(KDLON,KFLEV)
     597      REAL(KIND=8) ZRAYL(KDLON)
     598      REAL(KIND=8) ZRAY1(KDLON,KFLEV+1)
     599      REAL(KIND=8) ZRAY2(KDLON,KFLEV+1)
     600      REAL(KIND=8) ZREFZ(KDLON,2,KFLEV+1)
     601      REAL(KIND=8) ZRJ(KDLON,6,KFLEV+1)
     602      REAL(KIND=8) ZRJ0(KDLON,6,KFLEV+1)
     603      REAL(KIND=8) ZRK(KDLON,6,KFLEV+1)
     604      REAL(KIND=8) ZRK0(KDLON,6,KFLEV+1)
     605      REAL(KIND=8) ZRMUE(KDLON,KFLEV+1)
     606      REAL(KIND=8) ZRMU0(KDLON,KFLEV+1)
     607      REAL(KIND=8) ZR(KDLON,4)
     608      REAL(KIND=8) ZTAUAZ(KDLON,KFLEV)
     609      REAL(KIND=8) ZTRA1(KDLON,KFLEV+1)
     610      REAL(KIND=8) ZTRA2(KDLON,KFLEV+1)
     611      REAL(KIND=8) ZW(KDLON,4)
    620612C
    621613      INTEGER jl, jk, k, jaj, ikm1, ikl
    622 c
    623 c Prescribed Data:
    624 c
    625       REAL*8 RSUN(2)
    626       SAVE RSUN
    627 c$OMP THREADPRIVATE(RSUN)
    628       REAL*8 RRAY(2,6)
    629       SAVE RRAY
    630 c$OMP THREADPRIVATE(RRAY)
    631       DATA RSUN(1) / 0.441676 /
    632       DATA RSUN(2) / 0.558324 /
    633       DATA (RRAY(1,K),K=1,6) /
    634      S .428937E-01, .890743E+00,-.288555E+01,
    635      S .522744E+01,-.469173E+01, .161645E+01/
    636       DATA (RRAY(2,K),K=1,6) /
    637      S .697200E-02, .173297E-01,-.850903E-01,
    638      S .248261E+00,-.302031E+00, .129662E+00/
     614
    639615C     ------------------------------------------------------------------
    640616C
    … …  
    777753     S  ,  PFDOWN,PFUP                                            )
    778754      USE dimphy
     755      USE radiation_AR4_param, only : RSUN, RRAY
    779756      IMPLICIT none
    780757cym#include "dimensions.h"
    … …  
    821798      INTEGER KNU
    822799c-OB
    823       real*8 flag_aer
    824       real*8 tauae(kdlon,kflev,2)
    825       real*8 pizae(kdlon,kflev,2)
    826       real*8 cgae(kdlon,kflev,2)
    827       REAL*8 PAER(KDLON,KFLEV,5)
    828       REAL*8 PAKI(KDLON,2)
    829       REAL*8 PALBD(KDLON,2)
    830       REAL*8 PALBP(KDLON,2)
    831       REAL*8 PCG(KDLON,2,KFLEV)
    832       REAL*8 PCLD(KDLON,KFLEV)
    833       REAL*8 PCLDSW(KDLON,KFLEV)
    834       REAL*8 PCLEAR(KDLON)
    835       REAL*8 PDSIG(KDLON,KFLEV)
    836       REAL*8 POMEGA(KDLON,2,KFLEV)
    837       REAL*8 POZ(KDLON,KFLEV)
    838       REAL*8 PQS(KDLON,KFLEV)
    839       REAL*8 PRMU(KDLON)
    840       REAL*8 PSEC(KDLON)
    841       REAL*8 PTAU(KDLON,2,KFLEV)
    842       REAL*8 PUD(KDLON,5,KFLEV+1)
    843       REAL*8 PWV(KDLON,KFLEV)
    844 C
    845       REAL*8 PFDOWN(KDLON,KFLEV+1)
    846       REAL*8 PFUP(KDLON,KFLEV+1)
     800      real(kind=8) flag_aer
     801      real(kind=8) tauae(kdlon,kflev,2)
     802      real(kind=8) pizae(kdlon,kflev,2)
     803      real(kind=8) cgae(kdlon,kflev,2)
     804      REAL(KIND=8) PAER(KDLON,KFLEV,5)
     805      REAL(KIND=8) PAKI(KDLON,2)
     806      REAL(KIND=8) PALBD(KDLON,2)
     807      REAL(KIND=8) PALBP(KDLON,2)
     808      REAL(KIND=8) PCG(KDLON,2,KFLEV)
     809      REAL(KIND=8) PCLD(KDLON,KFLEV)
     810      REAL(KIND=8) PCLDSW(KDLON,KFLEV)
     811      REAL(KIND=8) PCLEAR(KDLON)
     812      REAL(KIND=8) PDSIG(KDLON,KFLEV)
     813      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)
     814      REAL(KIND=8) POZ(KDLON,KFLEV)
     815      REAL(KIND=8) PQS(KDLON,KFLEV)
     816      REAL(KIND=8) PRMU(KDLON)
     817      REAL(KIND=8) PSEC(KDLON)
     818      REAL(KIND=8) PTAU(KDLON,2,KFLEV)
     819      REAL(KIND=8) PUD(KDLON,5,KFLEV+1)
     820      REAL(KIND=8) PWV(KDLON,KFLEV)
     821C
     822      REAL(KIND=8) PFDOWN(KDLON,KFLEV+1)
     823      REAL(KIND=8) PFUP(KDLON,KFLEV+1)
    847824C
    848825C* LOCAL VARIABLES:
    849826C
    850827      INTEGER IIND2(2), IIND3(3)
    851       REAL*8 ZCGAZ(KDLON,KFLEV)
    852       REAL*8 ZFD(KDLON,KFLEV+1)
    853       REAL*8 ZFU(KDLON,KFLEV+1)
    854       REAL*8 ZG(KDLON)
    855       REAL*8 ZGG(KDLON)
    856       REAL*8 ZPIZAZ(KDLON,KFLEV)
    857       REAL*8 ZRAYL(KDLON)
    858       REAL*8 ZRAY1(KDLON,KFLEV+1)
    859       REAL*8 ZRAY2(KDLON,KFLEV+1)
    860       REAL*8 ZREF(KDLON)
    861       REAL*8 ZREFZ(KDLON,2,KFLEV+1)
    862       REAL*8 ZRE1(KDLON)
    863       REAL*8 ZRE2(KDLON)
    864       REAL*8 ZRJ(KDLON,6,KFLEV+1)
    865       REAL*8 ZRJ0(KDLON,6,KFLEV+1)
    866       REAL*8 ZRK(KDLON,6,KFLEV+1)
    867       REAL*8 ZRK0(KDLON,6,KFLEV+1)
    868       REAL*8 ZRL(KDLON,8)
    869       REAL*8 ZRMUE(KDLON,KFLEV+1)
    870       REAL*8 ZRMU0(KDLON,KFLEV+1)
    871       REAL*8 ZRMUZ(KDLON)
    872       REAL*8 ZRNEB(KDLON)
    873       REAL*8 ZRUEF(KDLON,8)
    874       REAL*8 ZR1(KDLON)
    875       REAL*8 ZR2(KDLON,2)
    876       REAL*8 ZR3(KDLON,3)
    877       REAL*8 ZR4(KDLON)
    878       REAL*8 ZR21(KDLON)
    879       REAL*8 ZR22(KDLON)
    880       REAL*8 ZS(KDLON)
    881       REAL*8 ZTAUAZ(KDLON,KFLEV)
    882       REAL*8 ZTO1(KDLON)
    883       REAL*8 ZTR(KDLON,2,KFLEV+1)
    884       REAL*8 ZTRA1(KDLON,KFLEV+1)
    885       REAL*8 ZTRA2(KDLON,KFLEV+1)
    886       REAL*8 ZTR1(KDLON)
    887       REAL*8 ZTR2(KDLON)
    888       REAL*8 ZW(KDLON)   
    889       REAL*8 ZW1(KDLON)
    890       REAL*8 ZW2(KDLON,2)
    891       REAL*8 ZW3(KDLON,3)
    892       REAL*8 ZW4(KDLON)
    893       REAL*8 ZW5(KDLON)
     828      REAL(KIND=8) ZCGAZ(KDLON,KFLEV)
     829      REAL(KIND=8) ZFD(KDLON,KFLEV+1)
     830      REAL(KIND=8) ZFU(KDLON,KFLEV+1)
     831      REAL(KIND=8) ZG(KDLON)
     832      REAL(KIND=8) ZGG(KDLON)
     833      REAL(KIND=8) ZPIZAZ(KDLON,KFLEV)
     834      REAL(KIND=8) ZRAYL(KDLON)
     835      REAL(KIND=8) ZRAY1(KDLON,KFLEV+1)
     836      REAL(KIND=8) ZRAY2(KDLON,KFLEV+1)
     837      REAL(KIND=8) ZREF(KDLON)
     838      REAL(KIND=8) ZREFZ(KDLON,2,KFLEV+1)
     839      REAL(KIND=8) ZRE1(KDLON)
     840      REAL(KIND=8) ZRE2(KDLON)
     841      REAL(KIND=8) ZRJ(KDLON,6,KFLEV+1)
     842      REAL(KIND=8) ZRJ0(KDLON,6,KFLEV+1)
     843      REAL(KIND=8) ZRK(KDLON,6,KFLEV+1)
     844      REAL(KIND=8) ZRK0(KDLON,6,KFLEV+1)
     845      REAL(KIND=8) ZRL(KDLON,8)
     846      REAL(KIND=8) ZRMUE(KDLON,KFLEV+1)
     847      REAL(KIND=8) ZRMU0(KDLON,KFLEV+1)
     848      REAL(KIND=8) ZRMUZ(KDLON)
     849      REAL(KIND=8) ZRNEB(KDLON)
     850      REAL(KIND=8) ZRUEF(KDLON,8)
     851      REAL(KIND=8) ZR1(KDLON)
     852      REAL(KIND=8) ZR2(KDLON,2)
     853      REAL(KIND=8) ZR3(KDLON,3)
     854      REAL(KIND=8) ZR4(KDLON)
     855      REAL(KIND=8) ZR21(KDLON)
     856      REAL(KIND=8) ZR22(KDLON)
     857      REAL(KIND=8) ZS(KDLON)
     858      REAL(KIND=8) ZTAUAZ(KDLON,KFLEV)
     859      REAL(KIND=8) ZTO1(KDLON)
     860      REAL(KIND=8) ZTR(KDLON,2,KFLEV+1)
     861      REAL(KIND=8) ZTRA1(KDLON,KFLEV+1)
     862      REAL(KIND=8) ZTRA2(KDLON,KFLEV+1)
     863      REAL(KIND=8) ZTR1(KDLON)
     864      REAL(KIND=8) ZTR2(KDLON)
     865      REAL(KIND=8) ZW(KDLON)   
     866      REAL(KIND=8) ZW1(KDLON)
     867      REAL(KIND=8) ZW2(KDLON,2)
     868      REAL(KIND=8) ZW3(KDLON,3)
     869      REAL(KIND=8) ZW4(KDLON)
     870      REAL(KIND=8) ZW5(KDLON)
    894871C
    895872      INTEGER jl, jk, k, jaj, ikm1, ikl, jn, jabs, jkm1
    896873      INTEGER jref, jkl, jklp1, jajp, jkki, jkkp4, jn2j, iabs
    897       REAL*8 ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
    898 C
    899 C* Prescribed Data:
    900 C
    901       REAL*8 RSUN(2)
    902       SAVE RSUN
    903 c$OMP THREADPRIVATE(RSUN)
    904       REAL*8 RRAY(2,6)
    905       SAVE RRAY
    906 c$OMP THREADPRIVATE(RRAY)
    907       DATA RSUN(1) / 0.441676 /
    908       DATA RSUN(2) / 0.558324 /
    909       DATA (RRAY(1,K),K=1,6) /
    910      S .428937E-01, .890743E+00,-.288555E+01,
    911      S .522744E+01,-.469173E+01, .161645E+01/
    912       DATA (RRAY(2,K),K=1,6) /
    913      S .697200E-02, .173297E-01,-.850903E-01,
    914      S .248261E+00,-.302031E+00, .129662E+00/
     874      REAL(KIND=8) ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
     875C
     876
    915877C
    916878C     ------------------------------------------------------------------
    … …  
    13261288     S  , PRK   , PRMU0 , PTAUAZ, PTRA1 , PTRA2                   )
    13271289      USE dimphy
     1290      USE radiation_AR4_param, only : TAUA, RPIZA, RCGA
    13281291      IMPLICIT none
    13291292cym#include "dimensions.h"
    … …  
    13571320      INTEGER KNU
    13581321c-OB
    1359       real*8 flag_aer
    1360       real*8 tauae(kdlon,kflev,2)
    1361       real*8 pizae(kdlon,kflev,2)
    1362       real*8 cgae(kdlon,kflev,2)
    1363       REAL*8 PAER(KDLON,KFLEV,5)
    1364       REAL*8 PALBP(KDLON,2)
    1365       REAL*8 PDSIG(KDLON,KFLEV)
    1366       REAL*8 PRAYL(KDLON)
    1367       REAL*8 PSEC(KDLON)
    1368 C
    1369       REAL*8 PCGAZ(KDLON,KFLEV)     
    1370       REAL*8 PPIZAZ(KDLON,KFLEV)
    1371       REAL*8 PRAY1(KDLON,KFLEV+1)
    1372       REAL*8 PRAY2(KDLON,KFLEV+1)
    1373       REAL*8 PREFZ(KDLON,2,KFLEV+1)
    1374       REAL*8 PRJ(KDLON,6,KFLEV+1)
    1375       REAL*8 PRK(KDLON,6,KFLEV+1)
    1376       REAL*8 PRMU0(KDLON,KFLEV+1)
    1377       REAL*8 PTAUAZ(KDLON,KFLEV)
    1378       REAL*8 PTRA1(KDLON,KFLEV+1)
    1379       REAL*8 PTRA2(KDLON,KFLEV+1)
     1322      real(kind=8) flag_aer
     1323      real(kind=8) tauae(kdlon,kflev,2)
     1324      real(kind=8) pizae(kdlon,kflev,2)
     1325      real(kind=8) cgae(kdlon,kflev,2)
     1326      REAL(KIND=8) PAER(KDLON,KFLEV,5)
     1327      REAL(KIND=8) PALBP(KDLON,2)
     1328      REAL(KIND=8) PDSIG(KDLON,KFLEV)
     1329      REAL(KIND=8) PRAYL(KDLON)
     1330      REAL(KIND=8) PSEC(KDLON)
     1331C
     1332      REAL(KIND=8) PCGAZ(KDLON,KFLEV)     
     1333      REAL(KIND=8) PPIZAZ(KDLON,KFLEV)
     1334      REAL(KIND=8) PRAY1(KDLON,KFLEV+1)
     1335      REAL(KIND=8) PRAY2(KDLON,KFLEV+1)
     1336      REAL(KIND=8) PREFZ(KDLON,2,KFLEV+1)
     1337      REAL(KIND=8) PRJ(KDLON,6,KFLEV+1)
     1338      REAL(KIND=8) PRK(KDLON,6,KFLEV+1)
     1339      REAL(KIND=8) PRMU0(KDLON,KFLEV+1)
     1340      REAL(KIND=8) PTAUAZ(KDLON,KFLEV)
     1341      REAL(KIND=8) PTRA1(KDLON,KFLEV+1)
     1342      REAL(KIND=8) PTRA2(KDLON,KFLEV+1)
    13801343C
    13811344C* LOCAL VARIABLES:
    13821345C
    1383       REAL*8 ZC0I(KDLON,KFLEV+1)       
    1384       REAL*8 ZCLE0(KDLON,KFLEV)
    1385       REAL*8 ZCLEAR(KDLON)
    1386       REAL*8 ZR21(KDLON)
    1387       REAL*8 ZR23(KDLON)
    1388       REAL*8 ZSS0(KDLON)
    1389       REAL*8 ZSCAT(KDLON)
    1390       REAL*8 ZTR(KDLON,2,KFLEV+1)
     1346      REAL(KIND=8) ZC0I(KDLON,KFLEV+1)       
     1347      REAL(KIND=8) ZCLE0(KDLON,KFLEV)
     1348      REAL(KIND=8) ZCLEAR(KDLON)
     1349      REAL(KIND=8) ZR21(KDLON)
     1350      REAL(KIND=8) ZR23(KDLON)
     1351      REAL(KIND=8) ZSS0(KDLON)
     1352      REAL(KIND=8) ZSCAT(KDLON)
     1353      REAL(KIND=8) ZTR(KDLON,2,KFLEV+1)
    13911354C
    13921355      INTEGER jl, jk, ja, jae, jkl, jklp1, jaj, jkm1, in
    1393       REAL*8 ZTRAY, ZGAR, ZRATIO, ZFF, ZFACOA, ZCORAE
    1394       REAL*8 ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZMU1, ZDEN1
    1395       REAL*8 ZBMU0, ZBMU1, ZRE11
    1396 C
    1397 C* Prescribed Data for Aerosols:
    1398 C
    1399       REAL*8 TAUA(2,5), RPIZA(2,5), RCGA(2,5)
    1400       SAVE TAUA, RPIZA, RCGA
    1401 c$OMP THREADPRIVATE(TAUA, RPIZA, RCGA)
    1402       DATA ((TAUA(IN,JA),JA=1,5),IN=1,2) /
    1403      S .730719, .912819, .725059, .745405, .682188 ,
    1404      S .730719, .912819, .725059, .745405, .682188 /
    1405       DATA ((RPIZA(IN,JA),JA=1,5),IN=1,2) /
    1406      S .872212, .982545, .623143, .944887, .997975 ,
    1407      S .872212, .982545, .623143, .944887, .997975 /
    1408       DATA ((RCGA (IN,JA),JA=1,5),IN=1,2) /
    1409      S .647596, .739002, .580845, .662657, .624246 ,
    1410      S .647596, .739002, .580845, .662657, .624246 /
     1356      REAL(KIND=8) ZTRAY, ZGAR, ZRATIO, ZFF, ZFACOA, ZCORAE
     1357      REAL(KIND=8) ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZMU1, ZDEN1
     1358      REAL(KIND=8) ZBMU0, ZBMU1, ZRE11
     1359C
     1360
    14111361C     ------------------------------------------------------------------
    14121362C
    … …  
    14161366 100  CONTINUE
    14171367C
     1368!cdir collapse
    14181369      DO 103 JK = 1 , KFLEV+1
    14191370      DO 102 JA = 1 , 6
    … …  
    17241675C
    17251676      INTEGER KNU
    1726       REAL*8 PALBD(KDLON,2)
    1727       REAL*8 PCG(KDLON,2,KFLEV)
    1728       REAL*8 PCLD(KDLON,KFLEV)
    1729       REAL*8 PDSIG(KDLON,KFLEV)
    1730       REAL*8 POMEGA(KDLON,2,KFLEV)
    1731       REAL*8 PRAYL(KDLON)
    1732       REAL*8 PSEC(KDLON)
    1733       REAL*8 PTAU(KDLON,2,KFLEV)
    1734 C
    1735       REAL*8 PRAY1(KDLON,KFLEV+1)
    1736       REAL*8 PRAY2(KDLON,KFLEV+1)
    1737       REAL*8 PREFZ(KDLON,2,KFLEV+1)
    1738       REAL*8 PRJ(KDLON,6,KFLEV+1)
    1739       REAL*8 PRK(KDLON,6,KFLEV+1)
    1740       REAL*8 PRMUE(KDLON,KFLEV+1)
    1741       REAL*8 PCGAZ(KDLON,KFLEV)
    1742       REAL*8 PPIZAZ(KDLON,KFLEV)
    1743       REAL*8 PTAUAZ(KDLON,KFLEV)
    1744       REAL*8 PTRA1(KDLON,KFLEV+1)
    1745       REAL*8 PTRA2(KDLON,KFLEV+1)
     1677      REAL(KIND=8) PALBD(KDLON,2)
     1678      REAL(KIND=8) PCG(KDLON,2,KFLEV)
     1679      REAL(KIND=8) PCLD(KDLON,KFLEV)
     1680      REAL(KIND=8) PDSIG(KDLON,KFLEV)
     1681      REAL(KIND=8) POMEGA(KDLON,2,KFLEV)
     1682      REAL(KIND=8) PRAYL(KDLON)
     1683      REAL(KIND=8) PSEC(KDLON)
     1684      REAL(KIND=8) PTAU(KDLON,2,KFLEV)
     1685C
     1686      REAL(KIND=8) PRAY1(KDLON,KFLEV+1)
     1687      REAL(KIND=8) PRAY2(KDLON,KFLEV+1)
     1688      REAL(KIND=8) PREFZ(KDLON,2,KFLEV+1)
     1689      REAL(KIND=8) PRJ(KDLON,6,KFLEV+1)
     1690      REAL(KIND=8) PRK(KDLON,6,KFLEV+1)
     1691      REAL(KIND=8) PRMUE(KDLON,KFLEV+1)
     1692      REAL(KIND=8) PCGAZ(KDLON,KFLEV)
     1693      REAL(KIND=8) PPIZAZ(KDLON,KFLEV)
     1694      REAL(KIND=8) PTAUAZ(KDLON,KFLEV)
     1695      REAL(KIND=8) PTRA1(KDLON,KFLEV+1)
     1696      REAL(KIND=8) PTRA2(KDLON,KFLEV+1)
    17461697C
    17471698C* LOCAL VARIABLES:
    17481699C
    1749       REAL*8 ZC1I(KDLON,KFLEV+1)
    1750       REAL*8 ZCLEQ(KDLON,KFLEV)
    1751       REAL*8 ZCLEAR(KDLON)
    1752       REAL*8 ZCLOUD(KDLON)
    1753       REAL*8 ZGG(KDLON)
    1754       REAL*8 ZREF(KDLON)
    1755       REAL*8 ZRE1(KDLON)
    1756       REAL*8 ZRE2(KDLON)
    1757       REAL*8 ZRMUZ(KDLON)
    1758       REAL*8 ZRNEB(KDLON)
    1759       REAL*8 ZR21(KDLON)
    1760       REAL*8 ZR22(KDLON)
    1761       REAL*8 ZR23(KDLON)
    1762       REAL*8 ZSS1(KDLON)
    1763       REAL*8 ZTO1(KDLON)
    1764       REAL*8 ZTR(KDLON,2,KFLEV+1)
    1765       REAL*8 ZTR1(KDLON)
    1766       REAL*8 ZTR2(KDLON)
    1767       REAL*8 ZW(KDLON)
     1700      REAL(KIND=8) ZC1I(KDLON,KFLEV+1)
     1701      REAL(KIND=8) ZCLEQ(KDLON,KFLEV)
     1702      REAL(KIND=8) ZCLEAR(KDLON)
     1703      REAL(KIND=8) ZCLOUD(KDLON)
     1704      REAL(KIND=8) ZGG(KDLON)
     1705      REAL(KIND=8) ZREF(KDLON)
     1706      REAL(KIND=8) ZRE1(KDLON)
     1707      REAL(KIND=8) ZRE2(KDLON)
     1708      REAL(KIND=8) ZRMUZ(KDLON)
     1709      REAL(KIND=8) ZRNEB(KDLON)
     1710      REAL(KIND=8) ZR21(KDLON)
     1711      REAL(KIND=8) ZR22(KDLON)
     1712      REAL(KIND=8) ZR23(KDLON)
     1713      REAL(KIND=8) ZSS1(KDLON)
     1714      REAL(KIND=8) ZTO1(KDLON)
     1715      REAL(KIND=8) ZTR(KDLON,2,KFLEV+1)
     1716      REAL(KIND=8) ZTR1(KDLON)
     1717      REAL(KIND=8) ZTR2(KDLON)
     1718      REAL(KIND=8) ZW(KDLON)
    17681719C
    17691720      INTEGER jk, jl, ja, jkl, jklp1, jkm1, jaj
    1770       REAL*8 ZFACOA, ZFACOC, ZCORAE, ZCORCD
    1771       REAL*8 ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZDEN1
    1772       REAL*8 ZMU1, ZRE11, ZBMU0, ZBMU1
     1721      REAL(KIND=8) ZFACOA, ZFACOC, ZCORAE, ZCORCD
     1722      REAL(KIND=8) ZMUE, ZGAP, ZWW, ZTO, ZDEN, ZDEN1
     1723      REAL(KIND=8) ZMU1, ZRE11, ZBMU0, ZBMU1
    17731724C
    17741725C     ------------------------------------------------------------------
    … …  
    20772028C* ARGUMENTS:
    20782029C
    2079       REAL*8 PGG(KDLON)   ! ASSYMETRY FACTOR
    2080       REAL*8 PREF(KDLON)  ! REFLECTIVITY OF THE UNDERLYING LAYER
    2081       REAL*8 PRMUZ(KDLON) ! COSINE OF SOLAR ZENITH ANGLE
    2082       REAL*8 PTO1(KDLON)  ! OPTICAL THICKNESS
    2083       REAL*8 PW(KDLON)    ! SINGLE SCATTERING ALBEDO
    2084       REAL*8 PRE1(KDLON)  ! LAYER REFLECTIVITY (NO UNDERLYING-LAYER REFLECTION)
    2085       REAL*8 PRE2(KDLON)  ! LAYER REFLECTIVITY
    2086       REAL*8 PTR1(KDLON)  ! LAYER TRANSMISSIVITY (NO UNDERLYING-LAYER REFLECTION)
    2087       REAL*8 PTR2(KDLON)  ! LAYER TRANSMISSIVITY
     2030      REAL(KIND=8) PGG(KDLON)   ! ASSYMETRY FACTOR
     2031      REAL(KIND=8) PREF(KDLON)  ! REFLECTIVITY OF THE UNDERLYING LAYER
     2032      REAL(KIND=8) PRMUZ(KDLON) ! COSINE OF SOLAR ZENITH ANGLE
     2033      REAL(KIND=8) PTO1(KDLON)  ! OPTICAL THICKNESS
     2034      REAL(KIND=8) PW(KDLON)    ! SINGLE SCATTERING ALBEDO
     2035      REAL(KIND=8) PRE1(KDLON)  ! LAYER REFLECTIVITY (NO UNDERLYING-LAYER REFLECTION)
     2036      REAL(KIND=8) PRE2(KDLON)  ! LAYER REFLECTIVITY
     2037      REAL(KIND=8) PTR1(KDLON)  ! LAYER TRANSMISSIVITY (NO UNDERLYING-LAYER REFLECTION)
     2038      REAL(KIND=8) PTR2(KDLON)  ! LAYER TRANSMISSIVITY
    20882039C
    20892040C* LOCAL VARIABLES:
    20902041C
    20912042      INTEGER jl
    2092       REAL*8 ZFF, ZGP, ZTOP, ZWCP, ZDT, ZX1, ZWM
    2093       REAL*8 ZRM2, ZRK, ZX2, ZRP, ZALPHA, ZBETA, ZARG
    2094       REAL*8 ZEXMU0, ZARG2, ZEXKP, ZEXKM, ZXP2P, ZXM2P, ZAP2B, ZAM2B
    2095       REAL*8 ZA11, ZA12, ZA13, ZA21, ZA22, ZA23
    2096       REAL*8 ZDENA, ZC1A, ZC2A, ZRI0A, ZRI1A
    2097       REAL*8 ZRI0B, ZRI1B
    2098       REAL*8 ZB21, ZB22, ZB23, ZDENB, ZC1B, ZC2B
    2099       REAL*8 ZRI0C, ZRI1C, ZRI0D, ZRI1D
     2043      REAL(KIND=8) ZFF, ZGP, ZTOP, ZWCP, ZDT, ZX1, ZWM
     2044      REAL(KIND=8) ZRM2, ZRK, ZX2, ZRP, ZALPHA, ZBETA, ZARG
     2045      REAL(KIND=8) ZEXMU0, ZARG2, ZEXKP, ZEXKM, ZXP2P, ZXM2P, ZAP2B,
     2046     $     ZAM2B
     2047      REAL(KIND=8) ZA11, ZA12, ZA13, ZA21, ZA22, ZA23
     2048      REAL(KIND=8) ZDENA, ZC1A, ZC2A, ZRI0A, ZRI1A
     2049      REAL(KIND=8) ZRI0B, ZRI1B
     2050      REAL(KIND=8) ZB21, ZB22, ZB23, ZDENB, ZC1B, ZC2B
     2051      REAL(KIND=8) ZRI0C, ZRI1C, ZRI0D, ZRI1D
    21002052C     ------------------------------------------------------------------
    21012053C
    … …  
    21752127      SUBROUTINE SWTT_LMDAR4 (KNU,KA,PU,PTR)
    21762128      USE dimphy
     2129      USE radiation_AR4_param, only : APAD, BPAD, D
    21772130      IMPLICIT none
    21782131cym#include "dimensions.h"
    … …  
    22122165      INTEGER KNU     ! INDEX OF THE SPECTRAL INTERVAL
    22132166      INTEGER KA      ! INDEX OF THE ABSORBER
    2214       REAL*8 PU(KDLON)  ! ABSORBER AMOUNT
    2215 C
    2216       REAL*8 PTR(KDLON) ! TRANSMISSION FUNCTION
     2167      REAL(KIND=8) PU(KDLON)  ! ABSORBER AMOUNT
     2168C
     2169      REAL(KIND=8) PTR(KDLON) ! TRANSMISSION FUNCTION
    22172170C
    22182171C* LOCAL VARIABLES:
    22192172C
    2220       REAL*8 ZR1(KDLON), ZR2(KDLON)
     2173      REAL(KIND=8) ZR1(KDLON), ZR2(KDLON)
    22212174      INTEGER jl, i,j
    22222175C
    2223 C* Prescribed Data:
    2224 C
    2225       REAL*8 APAD(2,3,7), BPAD(2,3,7), D(2,3)
    2226       SAVE APAD, BPAD, D
    2227 c$OMP THREADPRIVATE(APAD, BPAD, D)
    2228       DATA ((APAD(1,I,J),I=1,3),J=1,7) /
    2229      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2230      S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01,
    2231      S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00,
    2232      S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02,
    2233      S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02,
    2234      S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02,
    2235      S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 /
    2236       DATA ((APAD(2,I,J),I=1,3),J=1,7) /
    2237      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2238      S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02,
    2239      S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00,
    2240      S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00,
    2241      S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00,
    2242      S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00,
    2243      S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /
    2244 C
    2245       DATA ((BPAD(1,I,J),I=1,3),J=1,7) /
    2246      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2247      S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01,
    2248      S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00,
    2249      S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02,
    2250      S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02,
    2251      S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02,
    2252      S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 /
    2253       DATA ((BPAD(2,I,J),I=1,3),J=1,7) /
    2254      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2255      S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02,
    2256      S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01,
    2257      S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00,
    2258      S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00,
    2259      S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00,
    2260      S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /
    2261 c
    2262       DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 /
    2263       DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 /
     2176
    22642177C
    22652178C-----------------------------------------------------------------------
    … …  
    22932206      SUBROUTINE SWTT1_LMDAR4(KNU,KABS,KIND, PU, PTR)
    22942207      USE dimphy
     2208      USE radiation_AR4_param, only : APAD, BPAD, D
    22952209      IMPLICIT none
    22962210cym#include "dimensions.h"
    … …  
    23302244      INTEGER KABS         ! NUMBER OF ABSORBERS
    23312245      INTEGER KIND(KABS)   ! INDICES OF THE ABSORBERS
    2332       REAL*8 PU(KDLON,KABS)  ! ABSORBER AMOUNT
    2333 C
    2334       REAL*8 PTR(KDLON,KABS) ! TRANSMISSION FUNCTION
     2246      REAL(KIND=8) PU(KDLON,KABS)  ! ABSORBER AMOUNT
     2247C
     2248      REAL(KIND=8) PTR(KDLON,KABS) ! TRANSMISSION FUNCTION
    23352249C
    23362250C* LOCAL VARIABLES:
    23372251C
    2338       REAL*8 ZR1(KDLON)
    2339       REAL*8 ZR2(KDLON)
    2340       REAL*8 ZU(KDLON)
     2252      REAL(KIND=8) ZR1(KDLON)
     2253      REAL(KIND=8) ZR2(KDLON)
     2254      REAL(KIND=8) ZU(KDLON)
    23412255      INTEGER jl, ja, i, j, ia
    23422256C
    2343 C* Prescribed Data:
    2344 C
    2345       REAL*8 APAD(2,3,7), BPAD(2,3,7), D(2,3)
    2346       SAVE APAD, BPAD, D
    2347 c$OMP THREADPRIVATE(APAD, BPAD, D)
    2348       DATA ((APAD(1,I,J),I=1,3),J=1,7) /
    2349      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2350      S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01,
    2351      S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00,
    2352      S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02,
    2353      S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02,
    2354      S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02,
    2355      S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 /
    2356       DATA ((APAD(2,I,J),I=1,3),J=1,7) /
    2357      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2358      S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02,
    2359      S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00,
    2360      S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00,
    2361      S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00,
    2362      S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00,
    2363      S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /
    2364 C
    2365       DATA ((BPAD(1,I,J),I=1,3),J=1,7) /
    2366      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2367      S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01,
    2368      S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00,
    2369      S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02,
    2370      S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02,
    2371      S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02,
    2372      S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 /
    2373       DATA ((BPAD(2,I,J),I=1,3),J=1,7) /
    2374      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2375      S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02,
    2376      S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01,
    2377      S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00,
    2378      S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00,
    2379      S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00,
    2380      S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /
    2381 c
    2382       DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 /
    2383       DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 /
     2257
    23842258C-----------------------------------------------------------------------
    23852259C
    … …  
    24642338C-----------------------------------------------------------------------
    24652339cIM ctes ds clesphys.h
    2466 c     REAL*8 RCO2   ! CO2 CONCENTRATION (IPCC:353.E-06* 44.011/28.97)
    2467 c     REAL*8 RCH4   ! CH4 CONCENTRATION (IPCC: 1.72E-06* 16.043/28.97)
    2468 c     REAL*8 RN2O   ! N2O CONCENTRATION (IPCC: 310.E-09* 44.013/28.97)
    2469 c     REAL*8 RCFC11 ! CFC11 CONCENTRATION (IPCC: 280.E-12* 137.3686/28.97)
    2470 c     REAL*8 RCFC12 ! CFC12 CONCENTRATION (IPCC: 484.E-12* 120.9140/28.97)
     2340c     REAL(KIND=8) RCO2   ! CO2 CONCENTRATION (IPCC:353.E-06* 44.011/28.97)
     2341c     REAL(KIND=8) RCH4   ! CH4 CONCENTRATION (IPCC: 1.72E-06* 16.043/28.97)
     2342c     REAL(KIND=8) RN2O   ! N2O CONCENTRATION (IPCC: 310.E-09* 44.013/28.97)
     2343c     REAL(KIND=8) RCFC11 ! CFC11 CONCENTRATION (IPCC: 280.E-12* 137.3686/28.97)
     2344c     REAL(KIND=8) RCFC12 ! CFC12 CONCENTRATION (IPCC: 484.E-12* 120.9140/28.97)
    24712345#include "clesphys.h"
    2472       REAL*8 PCLDLD(KDLON,KFLEV)  ! DOWNWARD EFFECTIVE CLOUD COVER
    2473       REAL*8 PCLDLU(KDLON,KFLEV)  ! UPWARD EFFECTIVE CLOUD COVER
    2474       REAL*8 PDP(KDLON,KFLEV)     ! LAYER PRESSURE THICKNESS (Pa)
    2475       REAL*8 PDT0(KDLON)          ! SURFACE TEMPERATURE DISCONTINUITY (K)
    2476       REAL*8 PEMIS(KDLON)         ! SURFACE EMISSIVITY
    2477       REAL*8 PPMB(KDLON,KFLEV+1)  ! HALF LEVEL PRESSURE (mb)
    2478       REAL*8 PPSOL(KDLON)         ! SURFACE PRESSURE (Pa)
    2479       REAL*8 POZON(KDLON,KFLEV)   ! O3 CONCENTRATION (kg/kg)
    2480       REAL*8 PTL(KDLON,KFLEV+1)   ! HALF LEVEL TEMPERATURE (K)
    2481       REAL*8 PAER(KDLON,KFLEV,5)  ! OPTICAL THICKNESS OF THE AEROSOLS
    2482       REAL*8 PTAVE(KDLON,KFLEV)   ! LAYER TEMPERATURE (K)
    2483       REAL*8 PVIEW(KDLON)         ! COSECANT OF VIEWING ANGLE
    2484       REAL*8 PWV(KDLON,KFLEV)     ! SPECIFIC HUMIDITY (kg/kg)
    2485 C
    2486       REAL*8 PCOLR(KDLON,KFLEV)   ! LONG-WAVE TENDENCY (K/day)
    2487       REAL*8 PCOLR0(KDLON,KFLEV)  ! LONG-WAVE TENDENCY (K/day) clear-sky
    2488       REAL*8 PTOPLW(KDLON)        ! LONGWAVE FLUX AT T.O.A.
    2489       REAL*8 PSOLLW(KDLON)        ! LONGWAVE FLUX AT SURFACE
    2490       REAL*8 PTOPLW0(KDLON)       ! LONGWAVE FLUX AT T.O.A. (CLEAR-SKY)
    2491       REAL*8 PSOLLW0(KDLON)       ! LONGWAVE FLUX AT SURFACE (CLEAR-SKY)
     2346      REAL(KIND=8) PCLDLD(KDLON,KFLEV)  ! DOWNWARD EFFECTIVE CLOUD COVER
     2347      REAL(KIND=8) PCLDLU(KDLON,KFLEV)  ! UPWARD EFFECTIVE CLOUD COVER
     2348      REAL(KIND=8) PDP(KDLON,KFLEV)     ! LAYER PRESSURE THICKNESS (Pa)
     2349      REAL(KIND=8) PDT0(KDLON)          ! SURFACE TEMPERATURE DISCONTINUITY (K)
     2350      REAL(KIND=8) PEMIS(KDLON)         ! SURFACE EMISSIVITY
     2351      REAL(KIND=8) PPMB(KDLON,KFLEV+1)  ! HALF LEVEL PRESSURE (mb)
     2352      REAL(KIND=8) PPSOL(KDLON)         ! SURFACE PRESSURE (Pa)
     2353      REAL(KIND=8) POZON(KDLON,KFLEV)   ! O3 mass fraction
     2354      REAL(KIND=8) PTL(KDLON,KFLEV+1)   ! HALF LEVEL TEMPERATURE (K)
     2355      REAL(KIND=8) PAER(KDLON,KFLEV,5)  ! OPTICAL THICKNESS OF THE AEROSOLS
     2356      REAL(KIND=8) PTAVE(KDLON,KFLEV)   ! LAYER TEMPERATURE (K)
     2357      REAL(KIND=8) PVIEW(KDLON)         ! COSECANT OF VIEWING ANGLE
     2358      REAL(KIND=8) PWV(KDLON,KFLEV)     ! SPECIFIC HUMIDITY (kg/kg)
     2359C
     2360      REAL(KIND=8) PCOLR(KDLON,KFLEV)   ! LONG-WAVE TENDENCY (K/day)
     2361      REAL(KIND=8) PCOLR0(KDLON,KFLEV)  ! LONG-WAVE TENDENCY (K/day) clear-sky
     2362      REAL(KIND=8) PTOPLW(KDLON)        ! LONGWAVE FLUX AT T.O.A.
     2363      REAL(KIND=8) PSOLLW(KDLON)        ! LONGWAVE FLUX AT SURFACE
     2364      REAL(KIND=8) PTOPLW0(KDLON)       ! LONGWAVE FLUX AT T.O.A. (CLEAR-SKY)
     2365      REAL(KIND=8) PSOLLW0(KDLON)       ! LONGWAVE FLUX AT SURFACE (CLEAR-SKY)
    24922366c Rajout LF
    2493       real*8 psollwdown(kdlon)    ! LONGWAVE downwards flux at surface
     2367      real(kind=8) psollwdown(kdlon)    ! LONGWAVE downwards flux at surface
    24942368c Rajout IM
    2495 cIM   real*8 psollwdownclr(kdlon) ! LONGWAVE CS downwards flux at surface
    2496 cIM   real*8 ptoplwdown(kdlon)    ! LONGWAVE downwards flux at T.O.A.
    2497 cIM   real*8 ptoplwdownclr(kdlon) ! LONGWAVE CS downwards flux at T.O.A.
     2369cIM   real(kind=8) psollwdownclr(kdlon) ! LONGWAVE CS downwards flux at surface
     2370cIM   real(kind=8) ptoplwdown(kdlon)    ! LONGWAVE downwards flux at T.O.A.
     2371cIM   real(kind=8) ptoplwdownclr(kdlon) ! LONGWAVE CS downwards flux at T.O.A.
    24982372cIM
    2499       REAL*8 plwup(KDLON,KFLEV+1)  ! LW up total sky
    2500       REAL*8 plwup0(KDLON,KFLEV+1) ! LW up clear sky
    2501       REAL*8 plwdn(KDLON,KFLEV+1)  ! LW down total sky
    2502       REAL*8 plwdn0(KDLON,KFLEV+1) ! LW down clear sky
     2373      REAL(KIND=8) plwup(KDLON,KFLEV+1)  ! LW up total sky
     2374      REAL(KIND=8) plwup0(KDLON,KFLEV+1) ! LW up clear sky
     2375      REAL(KIND=8) plwdn(KDLON,KFLEV+1)  ! LW down total sky
     2376      REAL(KIND=8) plwdn0(KDLON,KFLEV+1) ! LW down clear sky
    25032377C-------------------------------------------------------------------------
    2504       REAL*8 ZABCU(KDLON,NUA,3*KFLEV+1)
    2505       REAL*8 ZOZ(KDLON,KFLEV)
    2506 c
    2507 cym      REAL*8 ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
    2508 cym      REAL*8 ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
    2509 cym      REAL*8 ZBINT(KDLON,KFLEV+1)            ! Intermediate variable
    2510 cym      REAL*8 ZBSUI(KDLON)                    ! Intermediate variable
    2511 cym      REAL*8,ZCTS(KDLON,KFLEV)               ! Intermediate variable
    2512 cym      REAL*8 ZCNTRB(KDLON,KFLEV+1,KFLEV+1)   ! Intermediate variable
     2378      REAL(KIND=8) ZABCU(KDLON,NUA,3*KFLEV+1)
     2379
     2380      REAL(KIND=8) ZOZ(KDLON,KFLEV)
     2381!     equivalent pressure of ozone in a layer, in Pa
     2382
     2383cym      REAL(KIND=8) ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down)
     2384cym      REAL(KIND=8) ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
     2385cym      REAL(KIND=8) ZBINT(KDLON,KFLEV+1)            ! Intermediate variable
     2386cym      REAL(KIND=8) ZBSUI(KDLON)                    ! Intermediate variable
     2387cym      REAL(KIND=8) ZCTS(KDLON,KFLEV)               ! Intermediate variable
     2388cym      REAL(KIND=8) ZCNTRB(KDLON,KFLEV+1,KFLEV+1)   ! Intermediate variable
    25132389cym      SAVE ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB
    2514       REAL*8,allocatable,save :: ZFLUX(:,:,:) ! RADIATIVE FLUXES (1:up; 2:down)
    2515       REAL*8,allocatable,save :: ZFLUC(:,:,:) ! CLEAR-SKY RADIATIVE FLUXES
    2516       REAL*8,allocatable,save :: ZBINT(:,:)            ! Intermediate variable
    2517       REAL*8,allocatable,save :: ZBSUI(:)                    ! Intermediate variable
    2518       REAL*8,allocatable,save :: ZCTS(:,:)               ! Intermediate variable
    2519       REAL*8,allocatable,save :: ZCNTRB(:,:,:)   ! Intermediate variable
     2390      REAL(KIND=8),allocatable,save :: ZFLUX(:,:,:) ! RADIATIVE FLUXES (1:up; 2:down)
     2391      REAL(KIND=8),allocatable,save :: ZFLUC(:,:,:) ! CLEAR-SKY RADIATIVE FLUXES
     2392      REAL(KIND=8),allocatable,save :: ZBINT(:,:)            ! Intermediate variable
     2393      REAL(KIND=8),allocatable,save :: ZBSUI(:)                    ! Intermediate variable
     2394      REAL(KIND=8),allocatable,save :: ZCTS(:,:)               ! Intermediate variable
     2395      REAL(KIND=8),allocatable,save :: ZCNTRB(:,:,:)   ! Intermediate variable
    25202396c$OMP THREADPRIVATE(ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB)
    25212397c
    … …  
    25502426C
    25512427      IF (MOD(itaplw0,lw0pas).EQ.0) THEN
    2552       DO k = 1, KFLEV  ! convertir ozone de kg/kg en pa/pa
    2553       DO i = 1, KDLON
    2554 c convertir ozone de kg/kg en pa (modif MPL 100505)
    2555          ZOZ(i,k) = POZON(i,k)*PDP(i,k) * RMD/RMO3
    2556 c        print *,'LW: ZOZ*10**6=',ZOZ(i,k)*1000000.
    2557       ENDDO
     2428c     Compute equivalent pressure of ozone from mass fraction:
     2429      DO k = 1, KFLEV
     2430         DO i = 1, KDLON
     2431            ZOZ(i,k) = POZON(i,k)*PDP(i,k)
     2432         ENDDO
    25582433      ENDDO
    25592434cIM ctes ds clesphys.h   CALL LWU(RCO2,RCH4, RN2O, RCFC11, RCFC12,
    … …  
    26092484     S               PABCU)
    26102485      USE dimphy
     2486      USE radiation_AR4_param, only : TREF, RT1, RAER, AT, BT, OCT
    26112487      IMPLICIT none
    26122488cym#include "dimensions.h"
    … …  
    26472523C* ARGUMENTS:
    26482524cIM ctes ds clesphys.h
    2649 c     REAL*8 RCO2
    2650 c     REAL*8 RCH4, RN2O, RCFC11, RCFC12
     2525c     REAL(KIND=8) RCO2
     2526c     REAL(KIND=8) RCH4, RN2O, RCFC11, RCFC12
    26512527#include "clesphys.h"
    2652       REAL*8 PAER(KDLON,KFLEV,5)
    2653       REAL*8 PDP(KDLON,KFLEV)
    2654       REAL*8 PPMB(KDLON,KFLEV+1)
    2655       REAL*8 PPSOL(KDLON)
    2656       REAL*8 POZ(KDLON,KFLEV)
    2657       REAL*8 PTAVE(KDLON,KFLEV)
    2658       REAL*8 PVIEW(KDLON)
    2659       REAL*8 PWV(KDLON,KFLEV)
    2660 C
    2661       REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
     2528      REAL(KIND=8) PAER(KDLON,KFLEV,5)
     2529      REAL(KIND=8) PDP(KDLON,KFLEV)
     2530      REAL(KIND=8) PPMB(KDLON,KFLEV+1)
     2531      REAL(KIND=8) PPSOL(KDLON)
     2532      REAL(KIND=8) POZ(KDLON,KFLEV)
     2533      REAL(KIND=8) PTAVE(KDLON,KFLEV)
     2534      REAL(KIND=8) PVIEW(KDLON)
     2535      REAL(KIND=8) PWV(KDLON,KFLEV)
     2536C
     2537      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
    26622538C
    26632539C-----------------------------------------------------------------------
    26642540C* LOCAL VARIABLES:
    2665       REAL*8 ZABLY(KDLON,NUA,3*KFLEV+1)
    2666       REAL*8 ZDUC(KDLON,3*KFLEV+1)
    2667       REAL*8 ZPHIO(KDLON)
    2668       REAL*8 ZPSC2(KDLON)
    2669       REAL*8 ZPSC3(KDLON)
    2670       REAL*8 ZPSH1(KDLON)
    2671       REAL*8 ZPSH2(KDLON)
    2672       REAL*8 ZPSH3(KDLON)
    2673       REAL*8 ZPSH4(KDLON)
    2674       REAL*8 ZPSH5(KDLON)
    2675       REAL*8 ZPSH6(KDLON)
    2676       REAL*8 ZPSIO(KDLON)
    2677       REAL*8 ZTCON(KDLON)
    2678       REAL*8 ZPHM6(KDLON)
    2679       REAL*8 ZPSM6(KDLON)
    2680       REAL*8 ZPHN6(KDLON)
    2681       REAL*8 ZPSN6(KDLON)
    2682       REAL*8 ZSSIG(KDLON,3*KFLEV+1)
    2683       REAL*8 ZTAVI(KDLON)
    2684       REAL*8 ZUAER(KDLON,Ninter)
    2685       REAL*8 ZXOZ(KDLON)
    2686       REAL*8 ZXWV(KDLON)
     2541      REAL(KIND=8) ZABLY(KDLON,NUA,3*KFLEV+1)
     2542      REAL(KIND=8) ZDUC(KDLON,3*KFLEV+1)
     2543      REAL(KIND=8) ZPHIO(KDLON)
     2544      REAL(KIND=8) ZPSC2(KDLON)
     2545      REAL(KIND=8) ZPSC3(KDLON)
     2546      REAL(KIND=8) ZPSH1(KDLON)
     2547      REAL(KIND=8) ZPSH2(KDLON)
     2548      REAL(KIND=8) ZPSH3(KDLON)
     2549      REAL(KIND=8) ZPSH4(KDLON)
     2550      REAL(KIND=8) ZPSH5(KDLON)
     2551      REAL(KIND=8) ZPSH6(KDLON)
     2552      REAL(KIND=8) ZPSIO(KDLON)
     2553      REAL(KIND=8) ZTCON(KDLON)
     2554      REAL(KIND=8) ZPHM6(KDLON)
     2555      REAL(KIND=8) ZPSM6(KDLON)
     2556      REAL(KIND=8) ZPHN6(KDLON)
     2557      REAL(KIND=8) ZPSN6(KDLON)
     2558      REAL(KIND=8) ZSSIG(KDLON,3*KFLEV+1)
     2559      REAL(KIND=8) ZTAVI(KDLON)
     2560      REAL(KIND=8) ZUAER(KDLON,Ninter)
     2561      REAL(KIND=8) ZXOZ(KDLON)
     2562      REAL(KIND=8) ZXWV(KDLON)
    26872563C
    26882564      INTEGER jl, jk, jkj, jkjr, jkjp, ig1
    … …  
    26912567      INTEGER jae1, jae2, jae3, jae, jjpn
    26922568      INTEGER ir, jc, jcp1
    2693       REAL*8 zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup
    2694       REAL*8 zfppw, ztx, ztx2, zzably
    2695       REAL*8 zcah1, zcbh1, zcah2, zcbh2, zcah3, zcbh3
    2696       REAL*8 zcah4, zcbh4, zcah5, zcbh5, zcah6, zcbh6
    2697       REAL*8 zcac8, zcbc8
    2698       REAL*8 zalup, zdiff
     2569      REAL(KIND=8) zdpm, zupm, zupmh2o, zupmco2, zupmo3, zu6, zup
     2570      REAL(KIND=8) zfppw, ztx, ztx2, zzably
     2571      REAL(KIND=8) zcah1, zcbh1, zcah2, zcbh2, zcah3, zcbh3
     2572      REAL(KIND=8) zcah4, zcbh4, zcah5, zcbh5, zcah6, zcbh6
     2573      REAL(KIND=8) zcac8, zcbc8
     2574      REAL(KIND=8) zalup, zdiff
    26992575c
    2700       REAL*8 PVGCO2, PVGH2O, PVGO3
    2701 C
    2702       REAL*8 R10E  ! DECIMAL/NATURAL LOG.FACTOR
     2576      REAL(KIND=8) PVGCO2, PVGH2O, PVGO3
     2577C
     2578      REAL(KIND=8) R10E  ! DECIMAL/NATURAL LOG.FACTOR
    27032579      PARAMETER (R10E=0.4342945)
    2704 c
    2705 c Used Data Block:
    2706 c
    2707       REAL*8 TREF
    2708       SAVE TREF
    2709 c$OMP THREADPRIVATE(TREF)
    2710       REAL*8 RT1(2)
    2711       SAVE RT1
    2712 c$OMP THREADPRIVATE(RT1)
    2713       REAL*8 RAER(5,5)
    2714       SAVE RAER
    2715 c$OMP THREADPRIVATE(RAER)
    2716       REAL*8 AT(8,3), BT(8,3)
    2717       SAVE AT, BT
    2718 c$OMP THREADPRIVATE(AT, BT)
    2719       REAL*8 OCT(4)
    2720       SAVE OCT
    2721 c$OMP THREADPRIVATE(OCT)
    2722       DATA TREF /250.0/
    2723       DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 /
    2724       DATA RAER / .038520, .037196, .040532, .054934, .038520
    2725      1          , .12613 , .18313 , .10357 , .064106, .126130
    2726      2          , .012579, .013649, .018652, .025181, .012579
    2727      3          , .011890, .016142, .021105, .028908, .011890
    2728      4          , .013792, .026810, .052203, .066338, .013792 /
    2729       DATA (AT(1,IR),IR=1,3) /
    2730      S 0.298199E-02,-.394023E-03,0.319566E-04 /
    2731       DATA (BT(1,IR),IR=1,3) /
    2732      S-0.106432E-04,0.660324E-06,0.174356E-06 /
    2733       DATA (AT(2,IR),IR=1,3) /
    2734      S 0.143676E-01,0.366501E-02,-.160822E-02 /
    2735       DATA (BT(2,IR),IR=1,3) /
    2736      S-0.553979E-04,-.101701E-04,0.920868E-05 /
    2737       DATA (AT(3,IR),IR=1,3) /
    2738      S 0.197861E-01,0.315541E-02,-.174547E-02 /
    2739       DATA (BT(3,IR),IR=1,3) /
    2740      S-0.877012E-04,0.513302E-04,0.523138E-06 /
    2741       DATA (AT(4,IR),IR=1,3) /
    2742      S 0.289560E-01,-.208807E-02,-.121943E-02 /
    2743       DATA (BT(4,IR),IR=1,3) /
    2744      S-0.165960E-03,0.157704E-03,-.146427E-04 /
    2745       DATA (AT(5,IR),IR=1,3) /
    2746      S 0.103800E-01,0.436296E-02,-.161431E-02 /
    2747       DATA (BT(5,IR),IR=1,3) /
    2748      S -.276744E-04,-.327381E-04,0.127646E-04 /
    2749       DATA (AT(6,IR),IR=1,3) /
    2750      S 0.868859E-02,-.972752E-03,0.000000E-00 /
    2751       DATA (BT(6,IR),IR=1,3) /
    2752      S -.278412E-04,-.713940E-06,0.117469E-05 /
    2753       DATA (AT(7,IR),IR=1,3) /
    2754      S 0.250073E-03,0.455875E-03,0.109242E-03 /
    2755       DATA (BT(7,IR),IR=1,3) /
    2756      S 0.199846E-05,-.216313E-05,0.175991E-06 /
    2757       DATA (AT(8,IR),IR=1,3) /
    2758      S 0.307423E-01,0.110879E-02,-.322172E-03 /
    2759       DATA (BT(8,IR),IR=1,3) /
    2760      S-0.108482E-03,0.258096E-05,-.814575E-06 /
    2761 c
    2762       DATA OCT /-.326E-03, -.102E-05, .137E-02, -.535E-05/
     2580
    27632581C-----------------------------------------------------------------------
    27642582c
    … …  
    30482866      INTEGER KLIM
    30492867C
    3050       REAL*8 PDP(KDLON,KFLEV)
    3051       REAL*8 PDT0(KDLON)
    3052       REAL*8 PEMIS(KDLON)
    3053       REAL*8 PPMB(KDLON,KFLEV+1)
    3054       REAL*8 PTL(KDLON,KFLEV+1)
    3055       REAL*8 PTAVE(KDLON,KFLEV)
    3056 C
    3057       REAL*8 PFLUC(KDLON,2,KFLEV+1)
     2868      REAL(KIND=8) PDP(KDLON,KFLEV)
     2869      REAL(KIND=8) PDT0(KDLON)
     2870      REAL(KIND=8) PEMIS(KDLON)
     2871      REAL(KIND=8) PPMB(KDLON,KFLEV+1)
     2872      REAL(KIND=8) PTL(KDLON,KFLEV+1)
     2873      REAL(KIND=8) PTAVE(KDLON,KFLEV)
     2874C
     2875      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1)
    30582876C     
    3059       REAL*8 PABCU(KDLON,NUA,3*KFLEV+1)
    3060       REAL*8 PBINT(KDLON,KFLEV+1)
    3061       REAL*8 PBSUI(KDLON)
    3062       REAL*8 PCTS(KDLON,KFLEV)
    3063       REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1)
     2877      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1)
     2878      REAL(KIND=8) PBINT(KDLON,KFLEV+1)
     2879      REAL(KIND=8) PBSUI(KDLON)
     2880      REAL(KIND=8) PCTS(KDLON,KFLEV)
     2881      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1)
    30642882C
    30652883C-------------------------------------------------------------------------
    30662884C
    30672885C* LOCAL VARIABLES:
    3068       REAL*8 ZB(KDLON,Ninter,KFLEV+1)
    3069       REAL*8 ZBSUR(KDLON,Ninter)
    3070       REAL*8 ZBTOP(KDLON,Ninter)
    3071       REAL*8 ZDBSL(KDLON,Ninter,KFLEV*2)
    3072       REAL*8 ZGA(KDLON,8,2,KFLEV)
    3073       REAL*8 ZGB(KDLON,8,2,KFLEV)
    3074       REAL*8 ZGASUR(KDLON,8,2)
    3075       REAL*8 ZGBSUR(KDLON,8,2)
    3076       REAL*8 ZGATOP(KDLON,8,2)
    3077       REAL*8 ZGBTOP(KDLON,8,2)
     2886      REAL(KIND=8) ZB(KDLON,Ninter,KFLEV+1)
     2887      REAL(KIND=8) ZBSUR(KDLON,Ninter)
     2888      REAL(KIND=8) ZBTOP(KDLON,Ninter)
     2889      REAL(KIND=8) ZDBSL(KDLON,Ninter,KFLEV*2)
     2890      REAL(KIND=8) ZGA(KDLON,8,2,KFLEV)
     2891      REAL(KIND=8) ZGB(KDLON,8,2,KFLEV)
     2892      REAL(KIND=8) ZGASUR(KDLON,8,2)
     2893      REAL(KIND=8) ZGBSUR(KDLON,8,2)
     2894      REAL(KIND=8) ZGATOP(KDLON,8,2)
     2895      REAL(KIND=8) ZGBTOP(KDLON,8,2)
    30782896C
    30792897      INTEGER nuaer, ntraer
    … …  
    31512969C* ARGUMENTS:
    31522970      INTEGER klim
    3153       REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
    3154       REAL*8 PBINT(KDLON,KFLEV+1)   ! HALF LEVEL PLANCK FUNCTION
    3155       REAL*8 PBSUIN(KDLON)          ! SURFACE PLANCK FUNCTION
    3156       REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) !CLEAR-SKY ENERGY EXCHANGE
    3157       REAL*8 PCTS(KDLON,KFLEV)      ! CLEAR-SKY LAYER COOLING-TO-SPACE
     2971      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
     2972      REAL(KIND=8) PBINT(KDLON,KFLEV+1)   ! HALF LEVEL PLANCK FUNCTION
     2973      REAL(KIND=8) PBSUIN(KDLON)          ! SURFACE PLANCK FUNCTION
     2974      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) !CLEAR-SKY ENERGY EXCHANGE
     2975      REAL(KIND=8) PCTS(KDLON,KFLEV)      ! CLEAR-SKY LAYER COOLING-TO-SPACE
    31582976c
    3159       REAL*8 PCLDLD(KDLON,KFLEV)
    3160       REAL*8 PCLDLU(KDLON,KFLEV)
    3161       REAL*8 PEMIS(KDLON)
    3162 C
    3163       REAL*8 PFLUX(KDLON,2,KFLEV+1)
     2977      REAL(KIND=8) PCLDLD(KDLON,KFLEV)
     2978      REAL(KIND=8) PCLDLU(KDLON,KFLEV)
     2979      REAL(KIND=8) PEMIS(KDLON)
     2980C
     2981      REAL(KIND=8) PFLUX(KDLON,2,KFLEV+1)
    31642982C-----------------------------------------------------------------------
    31652983C* LOCAL VARIABLES:
    31662984      INTEGER IMX(KDLON), IMXP(KDLON)
    31672985C
    3168       REAL*8 ZCLEAR(KDLON),ZCLOUD(KDLON),ZDNF(KDLON,KFLEV+1,KFLEV+1)
     2986      REAL(KIND=8) ZCLEAR(KDLON),ZCLOUD(KDLON),
     2987     $     ZDNF(KDLON,KFLEV+1,KFLEV+1)
    31692988     S  , ZFD(KDLON), ZFN10(KDLON), ZFU(KDLON)
    31702989     S  , ZUPF(KDLON,KFLEV+1,KFLEV+1)
    3171       REAL*8 ZCLM(KDLON,KFLEV+1,KFLEV+1)
     2990      REAL(KIND=8) ZCLM(KDLON,KFLEV+1,KFLEV+1)
    31722991C
    31732992      INTEGER jk, jl, imaxc, imx1, imx2, jkj, jkp1, jkm1
    31742993      INTEGER jk1, jk2, jkc, jkcp1, jcloud
    31752994      INTEGER imxm1, imxp1
    3176       REAL*8 zcfrac
     2995      REAL(KIND=8) zcfrac
    31772996C     ------------------------------------------------------------------
    31782997C
    … …  
    35043323     S  , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP)
    35053324      USE dimphy
     3325      USE radiation_AR4_param, only : TINTP, XP, GA, GB
    35063326      IMPLICIT none
    35073327cym#include "dimensions.h"
    … …  
    35603380C ARGUMENTS:
    35613381C
    3562       REAL*8 PDT0(KDLON)
    3563       REAL*8 PTAVE(KDLON,KFLEV)
    3564       REAL*8 PTL(KDLON,KFLEV+1)
    3565 C
    3566       REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION
    3567       REAL*8 PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION
    3568       REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
    3569       REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
    3570       REAL*8 PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION
    3571       REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
    3572       REAL*8 PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
    3573       REAL*8 PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
    3574       REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
    3575       REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
    3576       REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
    3577       REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
     3382      REAL(KIND=8) PDT0(KDLON)
     3383      REAL(KIND=8) PTAVE(KDLON,KFLEV)
     3384      REAL(KIND=8) PTL(KDLON,KFLEV+1)
     3385C
     3386      REAL(KIND=8) PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF LEVEL PLANCK FUNCTION
     3387      REAL(KIND=8) PBINT(KDLON,KFLEV+1) ! HALF LEVEL PLANCK FUNCTION
     3388      REAL(KIND=8) PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
     3389      REAL(KIND=8) PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
     3390      REAL(KIND=8) PBTOP(KDLON,Ninter) ! TOP SPECTRAL PLANCK FUNCTION
     3391      REAL(KIND=8) PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
     3392      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
     3393      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! dB/dT-weighted LAYER PADE APPROXIMANTS
     3394      REAL(KIND=8) PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
     3395      REAL(KIND=8) PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
     3396      REAL(KIND=8) PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
     3397      REAL(KIND=8) PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
    35783398C
    35793399C-------------------------------------------------------------------------
    35803400C*  LOCAL VARIABLES:
    35813401      INTEGER INDB(KDLON),INDS(KDLON)
    3582       REAL*8 ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
    3583       REAL*8 ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)
     3402      REAL(KIND=8) ZBLAY(KDLON,KFLEV),ZBLEV(KDLON,KFLEV+1)
     3403      REAL(KIND=8) ZRES(KDLON),ZRES2(KDLON),ZTI(KDLON),ZTI2(KDLON)
    35843404c
    35853405      INTEGER jk, jl, ic, jnu, jf, jg
    … …  
    35873407      INTEGER k, j, ixtox, indto, ixtx, indt
    35883408      INTEGER indsu, indtp
    3589       REAL*8 zdsto1, zdstox, zdst1, zdstx
     3409      REAL(KIND=8) zdsto1, zdstox, zdst1, zdstx
    35903410c
    35913411C* Quelques parametres:
    3592       REAL*8 TSTAND
     3412      REAL(KIND=8) TSTAND
    35933413      PARAMETER (TSTAND=250.0)
    3594       REAL*8 TSTP
     3414      REAL(KIND=8) TSTP
    35953415      PARAMETER (TSTP=12.5)
    35963416      INTEGER MXIXT
    … …  
    35983418C
    35993419C* Used Data Block:
    3600       REAL*8 TINTP(11)
    3601       SAVE TINTP
    3602 c$OMP THREADPRIVATE(TINTP)
    3603       REAL*8 GA(11,16,3), GB(11,16,3)
    3604       SAVE GA, GB
    3605 c$OMP THREADPRIVATE(GA, GB)
    3606       REAL*8 XP(6,6)
    3607       SAVE XP
    3608 c$OMP THREADPRIVATE(XP)
     3420c     REAL*8 TINTP(11)
     3421c     SAVE TINTP
     3422cc$OMP THREADPRIVATE(TINTP)
     3423c     REAL*8 GA(11,16,3), GB(11,16,3)
     3424c     SAVE GA, GB
     3425cc$OMP THREADPRIVATE(GA, GB)
     3426c     REAL*8 XP(6,6)
     3427c     SAVE XP
     3428cc$OMP THREADPRIVATE(XP)
    36093429c
    3610       DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
    3611      S             262.5, 275., 287.5, 300., 312.5 /
     3430c     DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
     3431c    S             262.5, 275., 287.5, 300., 312.5 /
    36123432C-----------------------------------------------------------------------
    36133433C-- WATER VAPOR -- INT.1 -- 0- 500 CM-1 -- FROM ABS225 ----------------
    … …  
    36223442C
    36233443C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3624       DATA (GA( 1, 1,IC),IC=1,3) /
    3625      S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/
    3626       DATA (GB( 1, 1,IC),IC=1,3) /
    3627      S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/
    3628       DATA (GA( 1, 2,IC),IC=1,3) /
    3629      S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/
    3630       DATA (GB( 1, 2,IC),IC=1,3) /
    3631      S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/
     3444C     DATA (GA( 1, 1,IC),IC=1,3) /
     3445C    S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/
     3446C     DATA (GB( 1, 1,IC),IC=1,3) /
     3447C    S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/
     3448C     DATA (GA( 1, 2,IC),IC=1,3) /
     3449C    S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/
     3450C     DATA (GB( 1, 2,IC),IC=1,3) /
     3451C    S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/
    36323452C
    36333453C----- INTERVAL = 1 ----- T =  200.0
    36343454C
    36353455C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3636       DATA (GA( 2, 1,IC),IC=1,3) /
    3637      S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/
    3638       DATA (GB( 2, 1,IC),IC=1,3) /
    3639      S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/
    3640       DATA (GA( 2, 2,IC),IC=1,3) /
    3641      S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/
    3642       DATA (GB( 2, 2,IC),IC=1,3) /
    3643      S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/
     3456C     DATA (GA( 2, 1,IC),IC=1,3) /
     3457C    S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/
     3458C     DATA (GB( 2, 1,IC),IC=1,3) /
     3459C    S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/
     3460C     DATA (GA( 2, 2,IC),IC=1,3) /
     3461C    S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/
     3462C     DATA (GB( 2, 2,IC),IC=1,3) /
     3463C    S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/
    36443464C
    36453465C----- INTERVAL = 1 ----- T =  212.5
    36463466C
    36473467C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3648       DATA (GA( 3, 1,IC),IC=1,3) /
    3649      S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/
    3650       DATA (GB( 3, 1,IC),IC=1,3) /
    3651      S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/
    3652       DATA (GA( 3, 2,IC),IC=1,3) /
    3653      S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/
    3654       DATA (GB( 3, 2,IC),IC=1,3) /
    3655      S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/
     3468C     DATA (GA( 3, 1,IC),IC=1,3) /
     3469C    S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/
     3470C     DATA (GB( 3, 1,IC),IC=1,3) /
     3471C    S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/
     3472C     DATA (GA( 3, 2,IC),IC=1,3) /
     3473C    S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/
     3474C     DATA (GB( 3, 2,IC),IC=1,3) /
     3475C    S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/
    36563476C
    36573477C----- INTERVAL = 1 ----- T =  225.0
    36583478C
    36593479C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3660       DATA (GA( 4, 1,IC),IC=1,3) /
    3661      S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/
    3662       DATA (GB( 4, 1,IC),IC=1,3) /
    3663      S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/
    3664       DATA (GA( 4, 2,IC),IC=1,3) /
    3665      S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/
    3666       DATA (GB( 4, 2,IC),IC=1,3) /
    3667      S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/
     3480C     DATA (GA( 4, 1,IC),IC=1,3) /
     3481C    S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/
     3482C     DATA (GB( 4, 1,IC),IC=1,3) /
     3483C    S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/
     3484C     DATA (GA( 4, 2,IC),IC=1,3) /
     3485C    S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/
     3486C     DATA (GB( 4, 2,IC),IC=1,3) /
     3487C    S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/
    36683488C
    36693489C----- INTERVAL = 1 ----- T =  237.5
    36703490C
    36713491C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3672       DATA (GA( 5, 1,IC),IC=1,3) /
    3673      S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/
    3674       DATA (GB( 5, 1,IC),IC=1,3) /
    3675      S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/
    3676       DATA (GA( 5, 2,IC),IC=1,3) /
    3677      S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/
    3678       DATA (GB( 5, 2,IC),IC=1,3) /
    3679      S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/
     3492C     DATA (GA( 5, 1,IC),IC=1,3) /
     3493C    S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/
     3494C     DATA (GB( 5, 1,IC),IC=1,3) /
     3495C    S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/
     3496C     DATA (GA( 5, 2,IC),IC=1,3) /
     3497C    S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/
     3498C     DATA (GB( 5, 2,IC),IC=1,3) /
     3499C    S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/
    36803500C
    36813501C----- INTERVAL = 1 ----- T =  250.0
    36823502C
    36833503C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3684       DATA (GA( 6, 1,IC),IC=1,3) /
    3685      S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/
    3686       DATA (GB( 6, 1,IC),IC=1,3) /
    3687      S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/
    3688       DATA (GA( 6, 2,IC),IC=1,3) /
    3689      S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/
    3690       DATA (GB( 6, 2,IC),IC=1,3) /
    3691      S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/
     3504C     DATA (GA( 6, 1,IC),IC=1,3) /
     3505C    S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/
     3506C     DATA (GB( 6, 1,IC),IC=1,3) /
     3507C    S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/
     3508C     DATA (GA( 6, 2,IC),IC=1,3) /
     3509C    S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/
     3510C     DATA (GB( 6, 2,IC),IC=1,3) /
     3511C    S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/
    36923512C
    36933513C----- INTERVAL = 1 ----- T =  262.5
    36943514C
    36953515C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3696       DATA (GA( 7, 1,IC),IC=1,3) /
    3697      S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/
    3698       DATA (GB( 7, 1,IC),IC=1,3) /
    3699      S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/
    3700       DATA (GA( 7, 2,IC),IC=1,3) /
    3701      S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/
    3702       DATA (GB( 7, 2,IC),IC=1,3) /
    3703      S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/
     3516C     DATA (GA( 7, 1,IC),IC=1,3) /
     3517C    S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/
     3518C     DATA (GB( 7, 1,IC),IC=1,3) /
     3519C    S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/
     3520C     DATA (GA( 7, 2,IC),IC=1,3) /
     3521C    S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/
     3522C     DATA (GB( 7, 2,IC),IC=1,3) /
     3523C    S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/
    37043524C
    37053525C----- INTERVAL = 1 ----- T =  275.0
    37063526C
    37073527C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3708       DATA (GA( 8, 1,IC),IC=1,3) /
    3709      S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/
    3710       DATA (GB( 8, 1,IC),IC=1,3) /
    3711      S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/
    3712       DATA (GA( 8, 2,IC),IC=1,3) /
    3713      S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/
    3714       DATA (GB( 8, 2,IC),IC=1,3) /
    3715      S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/
     3528C     DATA (GA( 8, 1,IC),IC=1,3) /
     3529C    S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/
     3530C     DATA (GB( 8, 1,IC),IC=1,3) /
     3531C    S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/
     3532C     DATA (GA( 8, 2,IC),IC=1,3) /
     3533C    S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/
     3534C     DATA (GB( 8, 2,IC),IC=1,3) /
     3535C    S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/
    37163536C
    37173537C----- INTERVAL = 1 ----- T =  287.5
    37183538C
    37193539C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3720       DATA (GA( 9, 1,IC),IC=1,3) /
    3721      S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/
    3722       DATA (GB( 9, 1,IC),IC=1,3) /
    3723      S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/
    3724       DATA (GA( 9, 2,IC),IC=1,3) /
    3725      S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/
    3726       DATA (GB( 9, 2,IC),IC=1,3) /
    3727      S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/
     3540C     DATA (GA( 9, 1,IC),IC=1,3) /
     3541C    S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/
     3542C     DATA (GB( 9, 1,IC),IC=1,3) /
     3543C    S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/
     3544C     DATA (GA( 9, 2,IC),IC=1,3) /
     3545C    S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/
     3546C     DATA (GB( 9, 2,IC),IC=1,3) /
     3547C    S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/
    37283548C
    37293549C----- INTERVAL = 1 ----- T =  300.0
    37303550C
    37313551C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3732       DATA (GA(10, 1,IC),IC=1,3) /
    3733      S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/
    3734       DATA (GB(10, 1,IC),IC=1,3) /
    3735      S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/
    3736       DATA (GA(10, 2,IC),IC=1,3) /
    3737      S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/
    3738       DATA (GB(10, 2,IC),IC=1,3) /
    3739      S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/
     3552C     DATA (GA(10, 1,IC),IC=1,3) /
     3553C    S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/
     3554C     DATA (GB(10, 1,IC),IC=1,3) /
     3555C    S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/
     3556C     DATA (GA(10, 2,IC),IC=1,3) /
     3557C    S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/
     3558C     DATA (GB(10, 2,IC),IC=1,3) /
     3559C    S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/
    37403560C
    37413561C----- INTERVAL = 1 ----- T =  312.5
    37423562C
    37433563C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3744       DATA (GA(11, 1,IC),IC=1,3) /
    3745      S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/
    3746       DATA (GB(11, 1,IC),IC=1,3) /
    3747      S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/
    3748       DATA (GA(11, 2,IC),IC=1,3) /
    3749      S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/
    3750       DATA (GB(11, 2,IC),IC=1,3) /
    3751      S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/
     3564C     DATA (GA(11, 1,IC),IC=1,3) /
     3565C    S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/
     3566C     DATA (GB(11, 1,IC),IC=1,3) /
     3567C    S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/
     3568C     DATA (GA(11, 2,IC),IC=1,3) /
     3569C    S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/
     3570C     DATA (GB(11, 2,IC),IC=1,3) /
     3571C    S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/
    37523572C
    37533573C
    … …  
    37643584C
    37653585C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3766       DATA (GA( 1, 3,IC),IC=1,3) /
    3767      S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/
    3768       DATA (GB( 1, 3,IC),IC=1,3) /
    3769      S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/
    3770       DATA (GA( 1, 4,IC),IC=1,3) /
    3771      S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/
    3772       DATA (GB( 1, 4,IC),IC=1,3) /
    3773      S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/
     3586C     DATA (GA( 1, 3,IC),IC=1,3) /
     3587C    S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/
     3588C     DATA (GB( 1, 3,IC),IC=1,3) /
     3589C    S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/
     3590C     DATA (GA( 1, 4,IC),IC=1,3) /
     3591C    S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/
     3592C     DATA (GB( 1, 4,IC),IC=1,3) /
     3593C    S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/
    37743594C
    37753595C----- INTERVAL = 2 ----- T =  200.0
    37763596C
    37773597C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3778       DATA (GA( 2, 3,IC),IC=1,3) /
    3779      S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/
    3780       DATA (GB( 2, 3,IC),IC=1,3) /
    3781      S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/
    3782       DATA (GA( 2, 4,IC),IC=1,3) /
    3783      S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/
    3784       DATA (GB( 2, 4,IC),IC=1,3) /
    3785      S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/
     3598C     DATA (GA( 2, 3,IC),IC=1,3) /
     3599C    S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/
     3600C     DATA (GB( 2, 3,IC),IC=1,3) /
     3601C    S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/
     3602C     DATA (GA( 2, 4,IC),IC=1,3) /
     3603C    S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/
     3604C     DATA (GB( 2, 4,IC),IC=1,3) /
     3605C    S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/
    37863606C
    37873607C----- INTERVAL = 2 ----- T =  212.5
    37883608C
    37893609C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3790       DATA (GA( 3, 3,IC),IC=1,3) /
    3791      S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/
    3792       DATA (GB( 3, 3,IC),IC=1,3) /
    3793      S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/
    3794       DATA (GA( 3, 4,IC),IC=1,3) /
    3795      S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/
    3796       DATA (GB( 3, 4,IC),IC=1,3) /
    3797      S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/
     3610C     DATA (GA( 3, 3,IC),IC=1,3) /
     3611C    S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/
     3612C     DATA (GB( 3, 3,IC),IC=1,3) /
     3613C    S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/
     3614C     DATA (GA( 3, 4,IC),IC=1,3) /
     3615C    S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/
     3616C     DATA (GB( 3, 4,IC),IC=1,3) /
     3617C    S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/
    37983618C
    37993619C----- INTERVAL = 2 ----- T =  225.0
    38003620C
    38013621C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3802       DATA (GA( 4, 3,IC),IC=1,3) /
    3803      S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/
    3804       DATA (GB( 4, 3,IC),IC=1,3) /
    3805      S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/
    3806       DATA (GA( 4, 4,IC),IC=1,3) /
    3807      S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/
    3808       DATA (GB( 4, 4,IC),IC=1,3) /
    3809      S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/
     3622C     DATA (GA( 4, 3,IC),IC=1,3) /
     3623C    S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/
     3624C     DATA (GB( 4, 3,IC),IC=1,3) /
     3625C    S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/
     3626C     DATA (GA( 4, 4,IC),IC=1,3) /
     3627C    S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/
     3628C     DATA (GB( 4, 4,IC),IC=1,3) /
     3629C    S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/
    38103630C
    38113631C----- INTERVAL = 2 ----- T =  237.5
    38123632C
    38133633C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3814       DATA (GA( 5, 3,IC),IC=1,3) /
    3815      S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/
    3816       DATA (GB( 5, 3,IC),IC=1,3) /
    3817      S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/
    3818       DATA (GA( 5, 4,IC),IC=1,3) /
    3819      S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/
    3820       DATA (GB( 5, 4,IC),IC=1,3) /
    3821      S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/
     3634C     DATA (GA( 5, 3,IC),IC=1,3) /
     3635C    S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/
     3636C     DATA (GB( 5, 3,IC),IC=1,3) /
     3637C    S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/
     3638C     DATA (GA( 5, 4,IC),IC=1,3) /
     3639C    S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/
     3640C     DATA (GB( 5, 4,IC),IC=1,3) /
     3641C    S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/
    38223642C
    38233643C----- INTERVAL = 2 ----- T =  250.0
    38243644C
    38253645C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3826       DATA (GA( 6, 3,IC),IC=1,3) /
    3827      S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/
    3828       DATA (GB( 6, 3,IC),IC=1,3) /
    3829      S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/
    3830       DATA (GA( 6, 4,IC),IC=1,3) /
    3831      S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/
    3832       DATA (GB( 6, 4,IC),IC=1,3) /
    3833      S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/
     3646C     DATA (GA( 6, 3,IC),IC=1,3) /
     3647C    S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/
     3648C     DATA (GB( 6, 3,IC),IC=1,3) /
     3649C    S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/
     3650C     DATA (GA( 6, 4,IC),IC=1,3) /
     3651C    S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/
     3652C     DATA (GB( 6, 4,IC),IC=1,3) /
     3653C    S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/
    38343654C
    38353655C----- INTERVAL = 2 ----- T =  262.5
    38363656C
    38373657C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3838       DATA (GA( 7, 3,IC),IC=1,3) /
    3839      S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/
    3840       DATA (GB( 7, 3,IC),IC=1,3) /
    3841      S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/
    3842       DATA (GA( 7, 4,IC),IC=1,3) /
    3843      S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/
    3844       DATA (GB( 7, 4,IC),IC=1,3) /
    3845      S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/
     3658C     DATA (GA( 7, 3,IC),IC=1,3) /
     3659C    S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/
     3660C     DATA (GB( 7, 3,IC),IC=1,3) /
     3661C    S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/
     3662C     DATA (GA( 7, 4,IC),IC=1,3) /
     3663C    S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/
     3664C     DATA (GB( 7, 4,IC),IC=1,3) /
     3665C    S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/
    38463666C
    38473667C----- INTERVAL = 2 ----- T =  275.0
    38483668C
    38493669C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3850       DATA (GA( 8, 3,IC),IC=1,3) /
    3851      S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/
    3852       DATA (GB( 8, 3,IC),IC=1,3) /
    3853      S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/
    3854       DATA (GA( 8, 4,IC),IC=1,3) /
    3855      S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/
    3856       DATA (GB( 8, 4,IC),IC=1,3) /
    3857      S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/
     3670C     DATA (GA( 8, 3,IC),IC=1,3) /
     3671C    S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/
     3672C     DATA (GB( 8, 3,IC),IC=1,3) /
     3673C    S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/
     3674C     DATA (GA( 8, 4,IC),IC=1,3) /
     3675C    S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/
     3676C     DATA (GB( 8, 4,IC),IC=1,3) /
     3677C    S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/
    38583678C
    38593679C----- INTERVAL = 2 ----- T =  287.5
    38603680C
    38613681C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3862       DATA (GA( 9, 3,IC),IC=1,3) /
    3863      S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/
    3864       DATA (GB( 9, 3,IC),IC=1,3) /
    3865      S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/
    3866       DATA (GA( 9, 4,IC),IC=1,3) /
    3867      S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/
    3868       DATA (GB( 9, 4,IC),IC=1,3) /
    3869      S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/
     3682C     DATA (GA( 9, 3,IC),IC=1,3) /
     3683C    S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/
     3684C     DATA (GB( 9, 3,IC),IC=1,3) /
     3685C    S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/
     3686C     DATA (GA( 9, 4,IC),IC=1,3) /
     3687C    S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/
     3688C     DATA (GB( 9, 4,IC),IC=1,3) /
     3689C    S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/
    38703690C
    38713691C----- INTERVAL = 2 ----- T =  300.0
    38723692C
    38733693C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3874       DATA (GA(10, 3,IC),IC=1,3) /
    3875      S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/
    3876       DATA (GB(10, 3,IC),IC=1,3) /
    3877      S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/
    3878       DATA (GA(10, 4,IC),IC=1,3) /
    3879      S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/
    3880       DATA (GB(10, 4,IC),IC=1,3) /
    3881      S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/
     3694C     DATA (GA(10, 3,IC),IC=1,3) /
     3695C    S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/
     3696C     DATA (GB(10, 3,IC),IC=1,3) /
     3697C    S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/
     3698C     DATA (GA(10, 4,IC),IC=1,3) /
     3699C    S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/
     3700C     DATA (GB(10, 4,IC),IC=1,3) /
     3701C    S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/
    38823702C
    38833703C----- INTERVAL = 2 ----- T =  312.5
    38843704C
    38853705C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3886       DATA (GA(11, 3,IC),IC=1,3) /
    3887      S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/
    3888       DATA (GB(11, 3,IC),IC=1,3) /
    3889      S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/
    3890       DATA (GA(11, 4,IC),IC=1,3) /
    3891      S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/
    3892       DATA (GB(11, 4,IC),IC=1,3) /
    3893      S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/
     3706C     DATA (GA(11, 3,IC),IC=1,3) /
     3707C    S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/
     3708C     DATA (GB(11, 3,IC),IC=1,3) /
     3709C    S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/
     3710C     DATA (GA(11, 4,IC),IC=1,3) /
     3711C    S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/
     3712C     DATA (GB(11, 4,IC),IC=1,3) /
     3713C    S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/
    38943714C
    38953715C
    … …  
    39103730C
    39113731C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3912       DATA (GA( 1, 7,IC),IC=1,3) /
    3913      S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/
    3914       DATA (GB( 1, 7,IC),IC=1,3) /
    3915      S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/
    3916       DATA (GA( 1, 8,IC),IC=1,3) /
    3917      S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/
    3918       DATA (GB( 1, 8,IC),IC=1,3) /
    3919      S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/
     3732C     DATA (GA( 1, 7,IC),IC=1,3) /
     3733C    S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/
     3734C     DATA (GB( 1, 7,IC),IC=1,3) /
     3735C    S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/
     3736C     DATA (GA( 1, 8,IC),IC=1,3) /
     3737C    S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/
     3738C     DATA (GB( 1, 8,IC),IC=1,3) /
     3739C    S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/
    39203740C
    39213741C----- INTERVAL = 3 ----- T =  200.0
    39223742C
    39233743C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3924       DATA (GA( 2, 7,IC),IC=1,3) /
    3925      S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/
    3926       DATA (GB( 2, 7,IC),IC=1,3) /
    3927      S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/
    3928       DATA (GA( 2, 8,IC),IC=1,3) /
    3929      S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/
    3930       DATA (GB( 2, 8,IC),IC=1,3) /
    3931      S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/
     3744C     DATA (GA( 2, 7,IC),IC=1,3) /
     3745C    S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/
     3746C     DATA (GB( 2, 7,IC),IC=1,3) /
     3747C    S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/
     3748C     DATA (GA( 2, 8,IC),IC=1,3) /
     3749C    S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/
     3750C     DATA (GB( 2, 8,IC),IC=1,3) /
     3751C    S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/
    39323752C
    39333753C----- INTERVAL = 3 ----- T =  212.5
    39343754C
    39353755C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3936       DATA (GA( 3, 7,IC),IC=1,3) /
    3937      S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/
    3938       DATA (GB( 3, 7,IC),IC=1,3) /
    3939      S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/
    3940       DATA (GA( 3, 8,IC),IC=1,3) /
    3941      S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/
    3942       DATA (GB( 3, 8,IC),IC=1,3) /
    3943      S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/
     3756C     DATA (GA( 3, 7,IC),IC=1,3) /
     3757C    S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/
     3758C     DATA (GB( 3, 7,IC),IC=1,3) /
     3759C    S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/
     3760C     DATA (GA( 3, 8,IC),IC=1,3) /
     3761C    S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/
     3762C     DATA (GB( 3, 8,IC),IC=1,3) /
     3763C    S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/
    39443764C
    39453765C----- INTERVAL = 3 ----- T =  225.0
    39463766C
    39473767C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3948       DATA (GA( 4, 7,IC),IC=1,3) /
    3949      S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/
    3950       DATA (GB( 4, 7,IC),IC=1,3) /
    3951      S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/
    3952       DATA (GA( 4, 8,IC),IC=1,3) /
    3953      S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/
    3954       DATA (GB( 4, 8,IC),IC=1,3) /
    3955      S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/
     3768C     DATA (GA( 4, 7,IC),IC=1,3) /
     3769C    S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/
     3770C     DATA (GB( 4, 7,IC),IC=1,3) /
     3771C    S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/
     3772C     DATA (GA( 4, 8,IC),IC=1,3) /
     3773C    S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/
     3774C     DATA (GB( 4, 8,IC),IC=1,3) /
     3775C    S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/
    39563776C
    39573777C----- INTERVAL = 3 ----- T =  237.5
    39583778C
    39593779C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3960       DATA (GA( 5, 7,IC),IC=1,3) /
    3961      S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/
    3962       DATA (GB( 5, 7,IC),IC=1,3) /
    3963      S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/
    3964       DATA (GA( 5, 8,IC),IC=1,3) /
    3965      S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/
    3966       DATA (GB( 5, 8,IC),IC=1,3) /
    3967      S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/
     3780C     DATA (GA( 5, 7,IC),IC=1,3) /
     3781C    S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/
     3782C     DATA (GB( 5, 7,IC),IC=1,3) /
     3783C    S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/
     3784C     DATA (GA( 5, 8,IC),IC=1,3) /
     3785C    S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/
     3786C     DATA (GB( 5, 8,IC),IC=1,3) /
     3787C    S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/
    39683788C
    39693789C----- INTERVAL = 3 ----- T =  250.0
    39703790C
    39713791C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3972       DATA (GA( 6, 7,IC),IC=1,3) /
    3973      S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/
    3974       DATA (GB( 6, 7,IC),IC=1,3) /
    3975      S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/
    3976       DATA (GA( 6, 8,IC),IC=1,3) /
    3977      S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/
    3978       DATA (GB( 6, 8,IC),IC=1,3) /
    3979      S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/
     3792C     DATA (GA( 6, 7,IC),IC=1,3) /
     3793C    S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/
     3794C     DATA (GB( 6, 7,IC),IC=1,3) /
     3795C    S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/
     3796C     DATA (GA( 6, 8,IC),IC=1,3) /
     3797C    S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/
     3798C     DATA (GB( 6, 8,IC),IC=1,3) /
     3799C    S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/
    39803800C
    39813801C----- INTERVAL = 3 ----- T =  262.5
    39823802C
    39833803C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3984       DATA (GA( 7, 7,IC),IC=1,3) /
    3985      S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/
    3986       DATA (GB( 7, 7,IC),IC=1,3) /
    3987      S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/
    3988       DATA (GA( 7, 8,IC),IC=1,3) /
    3989      S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/
    3990       DATA (GB( 7, 8,IC),IC=1,3) /
    3991      S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/
     3804C     DATA (GA( 7, 7,IC),IC=1,3) /
     3805C    S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/
     3806C     DATA (GB( 7, 7,IC),IC=1,3) /
     3807C    S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/
     3808C     DATA (GA( 7, 8,IC),IC=1,3) /
     3809C    S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/
     3810C     DATA (GB( 7, 8,IC),IC=1,3) /
     3811C    S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/
    39923812C
    39933813C----- INTERVAL = 3 ----- T =  275.0
    39943814C
    39953815C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3996       DATA (GA( 8, 7,IC),IC=1,3) /
    3997      S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/
    3998       DATA (GB( 8, 7,IC),IC=1,3) /
    3999      S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/
    4000       DATA (GA( 8, 8,IC),IC=1,3) /
    4001      S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/
    4002       DATA (GB( 8, 8,IC),IC=1,3) /
    4003      S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/
     3816C     DATA (GA( 8, 7,IC),IC=1,3) /
     3817C    S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/
     3818C     DATA (GB( 8, 7,IC),IC=1,3) /
     3819C    S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/
     3820C     DATA (GA( 8, 8,IC),IC=1,3) /
     3821C    S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/
     3822C     DATA (GB( 8, 8,IC),IC=1,3) /
     3823C    S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/
    40043824C
    40053825C----- INTERVAL = 3 ----- T =  287.5
    40063826C
    40073827C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4008       DATA (GA( 9, 7,IC),IC=1,3) /
    4009      S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/
    4010       DATA (GB( 9, 7,IC),IC=1,3) /
    4011      S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/
    4012       DATA (GA( 9, 8,IC),IC=1,3) /
    4013      S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/
    4014       DATA (GB( 9, 8,IC),IC=1,3) /
    4015      S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/
     3828C     DATA (GA( 9, 7,IC),IC=1,3) /
     3829C    S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/
     3830C     DATA (GB( 9, 7,IC),IC=1,3) /
     3831C    S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/
     3832C     DATA (GA( 9, 8,IC),IC=1,3) /
     3833C    S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/
     3834C     DATA (GB( 9, 8,IC),IC=1,3) /
     3835C    S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/
    40163836C
    40173837C----- INTERVAL = 3 ----- T =  300.0
    40183838C
    40193839C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4020       DATA (GA(10, 7,IC),IC=1,3) /
    4021      S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/
    4022       DATA (GB(10, 7,IC),IC=1,3) /
    4023      S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/
    4024       DATA (GA(10, 8,IC),IC=1,3) /
    4025      S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/
    4026       DATA (GB(10, 8,IC),IC=1,3) /
    4027      S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/
     3840C     DATA (GA(10, 7,IC),IC=1,3) /
     3841C    S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/
     3842C     DATA (GB(10, 7,IC),IC=1,3) /
     3843C    S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/
     3844C     DATA (GA(10, 8,IC),IC=1,3) /
     3845C    S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/
     3846C     DATA (GB(10, 8,IC),IC=1,3) /
     3847C    S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/
    40283848C
    40293849C----- INTERVAL = 3 ----- T =  312.5
    40303850C
    40313851C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4032       DATA (GA(11, 7,IC),IC=1,3) /
    4033      S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/
    4034       DATA (GB(11, 7,IC),IC=1,3) /
    4035      S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/
    4036       DATA (GA(11, 8,IC),IC=1,3) /
    4037      S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/
    4038       DATA (GB(11, 8,IC),IC=1,3) /
    4039      S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/
     3852C     DATA (GA(11, 7,IC),IC=1,3) /
     3853C    S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/
     3854C     DATA (GB(11, 7,IC),IC=1,3) /
     3855C    S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/
     3856C     DATA (GA(11, 8,IC),IC=1,3) /
     3857C    S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/
     3858C     DATA (GB(11, 8,IC),IC=1,3) /
     3859C    S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/
    40403860C
    40413861C
    … …  
    40473867C
    40483868C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4049       DATA (GA( 1, 9,IC),IC=1,3) /
    4050      S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/
    4051       DATA (GB( 1, 9,IC),IC=1,3) /
    4052      S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/
    4053       DATA (GA( 1,10,IC),IC=1,3) /
    4054      S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/
    4055       DATA (GB( 1,10,IC),IC=1,3) /
    4056      S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/
     3869C     DATA (GA( 1, 9,IC),IC=1,3) /
     3870C    S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/
     3871C     DATA (GB( 1, 9,IC),IC=1,3) /
     3872C    S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/
     3873C     DATA (GA( 1,10,IC),IC=1,3) /
     3874C    S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/
     3875C     DATA (GB( 1,10,IC),IC=1,3) /
     3876C    S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/
    40573877C
    40583878C----- INTERVAL = 4 ----- T =  200.0
    40593879C
    40603880C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4061       DATA (GA( 2, 9,IC),IC=1,3) /
    4062      S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/
    4063       DATA (GB( 2, 9,IC),IC=1,3) /
    4064      S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/
    4065       DATA (GA( 2,10,IC),IC=1,3) /
    4066      S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/
    4067       DATA (GB( 2,10,IC),IC=1,3) /
    4068      S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/
     3881C     DATA (GA( 2, 9,IC),IC=1,3) /
     3882C    S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/
     3883C     DATA (GB( 2, 9,IC),IC=1,3) /
     3884C    S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/
     3885C     DATA (GA( 2,10,IC),IC=1,3) /
     3886C    S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/
     3887C     DATA (GB( 2,10,IC),IC=1,3) /
     3888C    S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/
    40693889C
    40703890C----- INTERVAL = 4 ----- T =  212.5
    40713891C
    40723892C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4073       DATA (GA( 3, 9,IC),IC=1,3) /
    4074      S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/
    4075       DATA (GB( 3, 9,IC),IC=1,3) /
    4076      S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/
    4077       DATA (GA( 3,10,IC),IC=1,3) /
    4078      S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/
    4079       DATA (GB( 3,10,IC),IC=1,3) /
    4080      S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/
     3893C     DATA (GA( 3, 9,IC),IC=1,3) /
     3894C    S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/
     3895C     DATA (GB( 3, 9,IC),IC=1,3) /
     3896C    S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/
     3897C     DATA (GA( 3,10,IC),IC=1,3) /
     3898C    S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/
     3899C     DATA (GB( 3,10,IC),IC=1,3) /
     3900C    S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/
    40813901C
    40823902C----- INTERVAL = 4 ----- T =  225.0
    40833903C
    40843904C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4085       DATA (GA( 4, 9,IC),IC=1,3) /
    4086      S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/
    4087       DATA (GB( 4, 9,IC),IC=1,3) /
    4088      S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/
    4089       DATA (GA( 4,10,IC),IC=1,3) /
    4090      S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/
    4091       DATA (GB( 4,10,IC),IC=1,3) /
    4092      S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/
     3905C     DATA (GA( 4, 9,IC),IC=1,3) /
     3906C    S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/
     3907C     DATA (GB( 4, 9,IC),IC=1,3) /
     3908C    S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/
     3909C     DATA (GA( 4,10,IC),IC=1,3) /
     3910C    S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/
     3911C     DATA (GB( 4,10,IC),IC=1,3) /
     3912C    S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/
    40933913C
    40943914C----- INTERVAL = 4 ----- T =  237.5
    40953915C
    40963916C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4097       DATA (GA( 5, 9,IC),IC=1,3) /
    4098      S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/
    4099       DATA (GB( 5, 9,IC),IC=1,3) /
    4100      S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/
    4101       DATA (GA( 5,10,IC),IC=1,3) /
    4102      S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/
    4103       DATA (GB( 5,10,IC),IC=1,3) /
    4104      S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/
     3917C     DATA (GA( 5, 9,IC),IC=1,3) /
     3918C    S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/
     3919C     DATA (GB( 5, 9,IC),IC=1,3) /
     3920C    S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/
     3921C     DATA (GA( 5,10,IC),IC=1,3) /
     3922C    S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/
     3923C     DATA (GB( 5,10,IC),IC=1,3) /
     3924C    S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/
    41053925C
    41063926C----- INTERVAL = 4 ----- T =  250.0
    41073927C
    41083928C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4109       DATA (GA( 6, 9,IC),IC=1,3) /
    4110      S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/
    4111       DATA (GB( 6, 9,IC),IC=1,3) /
    4112      S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/
    4113       DATA (GA( 6,10,IC),IC=1,3) /
    4114      S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/
    4115       DATA (GB( 6,10,IC),IC=1,3) /
    4116      S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/
     3929C     DATA (GA( 6, 9,IC),IC=1,3) /
     3930C    S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/
     3931C     DATA (GB( 6, 9,IC),IC=1,3) /
     3932C    S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/
     3933C     DATA (GA( 6,10,IC),IC=1,3) /
     3934C    S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/
     3935C     DATA (GB( 6,10,IC),IC=1,3) /
     3936C    S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/
    41173937C
    41183938C----- INTERVAL = 4 ----- T =  262.5
    41193939C
    41203940C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4121       DATA (GA( 7, 9,IC),IC=1,3) /
    4122      S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/
    4123       DATA (GB( 7, 9,IC),IC=1,3) /
    4124      S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/
    4125       DATA (GA( 7,10,IC),IC=1,3) /
    4126      S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/
    4127       DATA (GB( 7,10,IC),IC=1,3) /
    4128      S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/
     3941C     DATA (GA( 7, 9,IC),IC=1,3) /
     3942C    S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/
     3943C     DATA (GB( 7, 9,IC),IC=1,3) /
     3944C    S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/
     3945C     DATA (GA( 7,10,IC),IC=1,3) /
     3946C    S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/
     3947C     DATA (GB( 7,10,IC),IC=1,3) /
     3948C    S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/
    41293949C
    41303950C----- INTERVAL = 4 ----- T =  275.0
    41313951C
    41323952C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4133       DATA (GA( 8, 9,IC),IC=1,3) /
    4134      S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/
    4135       DATA (GB( 8, 9,IC),IC=1,3) /
    4136      S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/
    4137       DATA (GA( 8,10,IC),IC=1,3) /
    4138      S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/
    4139       DATA (GB( 8,10,IC),IC=1,3) /
    4140      S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/
     3953C     DATA (GA( 8, 9,IC),IC=1,3) /
     3954C    S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/
     3955C     DATA (GB( 8, 9,IC),IC=1,3) /
     3956C    S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/
     3957C     DATA (GA( 8,10,IC),IC=1,3) /
     3958C    S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/
     3959C     DATA (GB( 8,10,IC),IC=1,3) /
     3960C    S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/
    41413961C
    41423962C----- INTERVAL = 4 ----- T =  287.5
    41433963C
    41443964C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4145       DATA (GA( 9, 9,IC),IC=1,3) /
    4146      S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/
    4147       DATA (GB( 9, 9,IC),IC=1,3) /
    4148      S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/
    4149       DATA (GA( 9,10,IC),IC=1,3) /
    4150      S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/
    4151       DATA (GB( 9,10,IC),IC=1,3) /
    4152      S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/
     3965C     DATA (GA( 9, 9,IC),IC=1,3) /
     3966C    S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/
     3967C     DATA (GB( 9, 9,IC),IC=1,3) /
     3968C    S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/
     3969C     DATA (GA( 9,10,IC),IC=1,3) /
     3970C    S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/
     3971C     DATA (GB( 9,10,IC),IC=1,3) /
     3972C    S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/
    41533973C
    41543974C----- INTERVAL = 4 ----- T =  300.0
    41553975C
    41563976C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4157       DATA (GA(10, 9,IC),IC=1,3) /
    4158      S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/
    4159       DATA (GB(10, 9,IC),IC=1,3) /
    4160      S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/
    4161       DATA (GA(10,10,IC),IC=1,3) /
    4162      S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/
    4163       DATA (GB(10,10,IC),IC=1,3) /
    4164      S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/
     3977C     DATA (GA(10, 9,IC),IC=1,3) /
     3978C    S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/
     3979C     DATA (GB(10, 9,IC),IC=1,3) /
     3980C    S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/
     3981C     DATA (GA(10,10,IC),IC=1,3) /
     3982C    S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/
     3983C     DATA (GB(10,10,IC),IC=1,3) /
     3984C    S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/
    41653985C
    41663986C----- INTERVAL = 4 ----- T =  312.5
    41673987C
    41683988C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4169       DATA (GA(11, 9,IC),IC=1,3) /
    4170      S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/
    4171       DATA (GB(11, 9,IC),IC=1,3) /
    4172      S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/
    4173       DATA (GA(11,10,IC),IC=1,3) /
    4174      S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/
    4175       DATA (GB(11,10,IC),IC=1,3) /
    4176      S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/
     3989C     DATA (GA(11, 9,IC),IC=1,3) /
     3990C    S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/
     3991C     DATA (GB(11, 9,IC),IC=1,3) /
     3992C    S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/
     3993C     DATA (GA(11,10,IC),IC=1,3) /
     3994C    S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/
     3995C     DATA (GB(11,10,IC),IC=1,3) /
     3996C    S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/
    41773997C
    41783998C
    … …  
    41874007C
    41884008C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4189       DATA (GA( 1, 5,IC),IC=1,3) /
    4190      S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/
    4191       DATA (GB( 1, 5,IC),IC=1,3) /
    4192      S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/
    4193       DATA (GA( 1, 6,IC),IC=1,3) /
    4194      S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/
    4195       DATA (GB( 1, 6,IC),IC=1,3) /
    4196      S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/
     4009C     DATA (GA( 1, 5,IC),IC=1,3) /
     4010C    S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/
     4011C     DATA (GB( 1, 5,IC),IC=1,3) /
     4012C    S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/
     4013C     DATA (GA( 1, 6,IC),IC=1,3) /
     4014C    S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/
     4015C     DATA (GB( 1, 6,IC),IC=1,3) /
     4016C    S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/
    41974017C
    41984018C----- INTERVAL = 5 ----- T =  200.0
    41994019C
    42004020C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4201       DATA (GA( 2, 5,IC),IC=1,3) /
    4202      S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/
    4203       DATA (GB( 2, 5,IC),IC=1,3) /
    4204      S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/
    4205       DATA (GA( 2, 6,IC),IC=1,3) /
    4206      S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/
    4207       DATA (GB( 2, 6,IC),IC=1,3) /
    4208      S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/
     4021C     DATA (GA( 2, 5,IC),IC=1,3) /
     4022C    S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/
     4023C     DATA (GB( 2, 5,IC),IC=1,3) /
     4024C    S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/
     4025C     DATA (GA( 2, 6,IC),IC=1,3) /
     4026C    S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/
     4027C     DATA (GB( 2, 6,IC),IC=1,3) /
     4028C    S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/
    42094029C
    42104030C----- INTERVAL = 5 ----- T =  212.5
    42114031C
    42124032C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4213       DATA (GA( 3, 5,IC),IC=1,3) /
    4214      S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/
    4215       DATA (GB( 3, 5,IC),IC=1,3) /
    4216      S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/
    4217       DATA (GA( 3, 6,IC),IC=1,3) /
    4218      S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/
    4219       DATA (GB( 3, 6,IC),IC=1,3) /
    4220      S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/
     4033C     DATA (GA( 3, 5,IC),IC=1,3) /
     4034C    S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/
     4035C     DATA (GB( 3, 5,IC),IC=1,3) /
     4036C    S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/
     4037C     DATA (GA( 3, 6,IC),IC=1,3) /
     4038C    S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/
     4039C     DATA (GB( 3, 6,IC),IC=1,3) /
     4040C    S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/
    42214041C
    42224042C----- INTERVAL = 5 ----- T =  225.0
    42234043C
    42244044C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4225       DATA (GA( 4, 5,IC),IC=1,3) /
    4226      S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/
    4227       DATA (GB( 4, 5,IC),IC=1,3) /
    4228      S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/
    4229       DATA (GA( 4, 6,IC),IC=1,3) /
    4230      S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/
    4231       DATA (GB( 4, 6,IC),IC=1,3) /
    4232      S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/
     4045C     DATA (GA( 4, 5,IC),IC=1,3) /
     4046C    S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/
     4047C     DATA (GB( 4, 5,IC),IC=1,3) /
     4048C    S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/
     4049C     DATA (GA( 4, 6,IC),IC=1,3) /
     4050C    S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/
     4051C     DATA (GB( 4, 6,IC),IC=1,3) /
     4052C    S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/
    42334053C
    42344054C----- INTERVAL = 5 ----- T =  237.5
    42354055C
    42364056C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4237       DATA (GA( 5, 5,IC),IC=1,3) /
    4238      S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/
    4239       DATA (GB( 5, 5,IC),IC=1,3) /
    4240      S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/
    4241       DATA (GA( 5, 6,IC),IC=1,3) /
    4242      S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/
    4243       DATA (GB( 5, 6,IC),IC=1,3) /
    4244      S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/
     4057C     DATA (GA( 5, 5,IC),IC=1,3) /
     4058C    S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/
     4059C     DATA (GB( 5, 5,IC),IC=1,3) /
     4060C    S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/
     4061C     DATA (GA( 5, 6,IC),IC=1,3) /
     4062C    S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/
     4063C     DATA (GB( 5, 6,IC),IC=1,3) /
     4064C    S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/
    42454065C
    42464066C----- INTERVAL = 5 ----- T =  250.0
    42474067C
    42484068C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4249       DATA (GA( 6, 5,IC),IC=1,3) /
    4250      S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/
    4251       DATA (GB( 6, 5,IC),IC=1,3) /
    4252      S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/
    4253       DATA (GA( 6, 6,IC),IC=1,3) /
    4254      S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/
    4255       DATA (GB( 6, 6,IC),IC=1,3) /
    4256      S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/
     4069C     DATA (GA( 6, 5,IC),IC=1,3) /
     4070C    S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/
     4071C     DATA (GB( 6, 5,IC),IC=1,3) /
     4072C    S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/
     4073C     DATA (GA( 6, 6,IC),IC=1,3) /
     4074C    S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/
     4075C     DATA (GB( 6, 6,IC),IC=1,3) /
     4076C    S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/
    42574077C
    42584078C----- INTERVAL = 5 ----- T =  262.5
    42594079C
    42604080C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4261       DATA (GA( 7, 5,IC),IC=1,3) /
    4262      S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/
    4263       DATA (GB( 7, 5,IC),IC=1,3) /
    4264      S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/
    4265       DATA (GA( 7, 6,IC),IC=1,3) /
    4266      S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/
    4267       DATA (GB( 7, 6,IC),IC=1,3) /
    4268      S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/
     4081C     DATA (GA( 7, 5,IC),IC=1,3) /
     4082C    S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/
     4083C     DATA (GB( 7, 5,IC),IC=1,3) /
     4084C    S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/
     4085C     DATA (GA( 7, 6,IC),IC=1,3) /
     4086C    S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/
     4087C     DATA (GB( 7, 6,IC),IC=1,3) /
     4088C    S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/
    42694089C
    42704090C----- INTERVAL = 5 ----- T =  275.0
    42714091C
    42724092C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4273       DATA (GA( 8, 5,IC),IC=1,3) /
    4274      S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/
    4275       DATA (GB( 8, 5,IC),IC=1,3) /
    4276      S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/
    4277       DATA (GA( 8, 6,IC),IC=1,3) /
    4278      S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/
    4279       DATA (GB( 8, 6,IC),IC=1,3) /
    4280      S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/
     4093C     DATA (GA( 8, 5,IC),IC=1,3) /
     4094C    S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/
     4095C     DATA (GB( 8, 5,IC),IC=1,3) /
     4096C    S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/
     4097C     DATA (GA( 8, 6,IC),IC=1,3) /
     4098C    S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/
     4099C     DATA (GB( 8, 6,IC),IC=1,3) /
     4100C    S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/
    42814101C
    42824102C----- INTERVAL = 5 ----- T =  287.5
    42834103C
    42844104C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4285       DATA (GA( 9, 5,IC),IC=1,3) /
    4286      S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/
    4287       DATA (GB( 9, 5,IC),IC=1,3) /
    4288      S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/
    4289       DATA (GA( 9, 6,IC),IC=1,3) /
    4290      S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/
    4291       DATA (GB( 9, 6,IC),IC=1,3) /
    4292      S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/
     4105C     DATA (GA( 9, 5,IC),IC=1,3) /
     4106C    S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/
     4107C     DATA (GB( 9, 5,IC),IC=1,3) /
     4108C    S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/
     4109C     DATA (GA( 9, 6,IC),IC=1,3) /
     4110C    S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/
     4111C     DATA (GB( 9, 6,IC),IC=1,3) /
     4112C    S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/
    42934113C
    42944114C----- INTERVAL = 5 ----- T =  300.0
    42954115C
    42964116C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4297       DATA (GA(10, 5,IC),IC=1,3) /
    4298      S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/
    4299       DATA (GB(10, 5,IC),IC=1,3) /
    4300      S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/
    4301       DATA (GA(10, 6,IC),IC=1,3) /
    4302      S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/
    4303       DATA (GB(10, 6,IC),IC=1,3) /
    4304      S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/
     4117C     DATA (GA(10, 5,IC),IC=1,3) /
     4118C    S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/
     4119C     DATA (GB(10, 5,IC),IC=1,3) /
     4120C    S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/
     4121C     DATA (GA(10, 6,IC),IC=1,3) /
     4122C    S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/
     4123C     DATA (GB(10, 6,IC),IC=1,3) /
     4124C    S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/
    43054125C
    43064126C----- INTERVAL = 5 ----- T =  312.5
    43074127C
    43084128C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4309       DATA (GA(11, 5,IC),IC=1,3) /
    4310      S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/
    4311       DATA (GB(11, 5,IC),IC=1,3) /
    4312      S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/
    4313       DATA (GA(11, 6,IC),IC=1,3) /
    4314      S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/
    4315       DATA (GB(11, 6,IC),IC=1,3) /
    4316      S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/
     4129C     DATA (GA(11, 5,IC),IC=1,3) /
     4130C    S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/
     4131C     DATA (GB(11, 5,IC),IC=1,3) /
     4132C    S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/
     4133C     DATA (GA(11, 6,IC),IC=1,3) /
     4134C    S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/
     4135C     DATA (GB(11, 6,IC),IC=1,3) /
     4136C    S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/
    43174137C
    43184138C
    … …  
    43264146C
    43274147C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4328       DATA (GA( 1,11,IC),IC=1,3) /
    4329      S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/
    4330       DATA (GB( 1,11,IC),IC=1,3) /
    4331      S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/
    4332       DATA (GA( 1,12,IC),IC=1,3) /
    4333      S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/
    4334       DATA (GB( 1,12,IC),IC=1,3) /
    4335      S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/
     4148C     DATA (GA( 1,11,IC),IC=1,3) /
     4149C    S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/
     4150C     DATA (GB( 1,11,IC),IC=1,3) /
     4151C    S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/
     4152C     DATA (GA( 1,12,IC),IC=1,3) /
     4153C    S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/
     4154C     DATA (GB( 1,12,IC),IC=1,3) /
     4155C    S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/
    43364156C
    43374157C----- INTERVAL = 6 ----- T =  200.0
    43384158C
    43394159C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4340       DATA (GA( 2,11,IC),IC=1,3) /
    4341      S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/
    4342       DATA (GB( 2,11,IC),IC=1,3) /
    4343      S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/
    4344       DATA (GA( 2,12,IC),IC=1,3) /
    4345      S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/
    4346       DATA (GB( 2,12,IC),IC=1,3) /
    4347      S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/
     4160C     DATA (GA( 2,11,IC),IC=1,3) /
     4161C    S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/
     4162C     DATA (GB( 2,11,IC),IC=1,3) /
     4163C    S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/
     4164C     DATA (GA( 2,12,IC),IC=1,3) /
     4165C    S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/
     4166C     DATA (GB( 2,12,IC),IC=1,3) /
     4167C    S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/
    43484168C
    43494169C----- INTERVAL = 6 ----- T =  212.5
    43504170C
    43514171C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4352       DATA (GA( 3,11,IC),IC=1,3) /
    4353      S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/
    4354       DATA (GB( 3,11,IC),IC=1,3) /
    4355      S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/
    4356       DATA (GA( 3,12,IC),IC=1,3) /
    4357      S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/
    4358       DATA (GB( 3,12,IC),IC=1,3) /
    4359      S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/
     4172C     DATA (GA( 3,11,IC),IC=1,3) /
     4173C    S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/
     4174C     DATA (GB( 3,11,IC),IC=1,3) /
     4175C    S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/
     4176C     DATA (GA( 3,12,IC),IC=1,3) /
     4177C    S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/
     4178C     DATA (GB( 3,12,IC),IC=1,3) /
     4179C    S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/
    43604180C
    43614181C----- INTERVAL = 6 ----- T =  225.0
    43624182C
    43634183C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4364       DATA (GA( 4,11,IC),IC=1,3) /
    4365      S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/
    4366       DATA (GB( 4,11,IC),IC=1,3) /
    4367      S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/
    4368       DATA (GA( 4,12,IC),IC=1,3) /
    4369      S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/
    4370       DATA (GB( 4,12,IC),IC=1,3) /
    4371      S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/
     4184C     DATA (GA( 4,11,IC),IC=1,3) /
     4185C    S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/
     4186C     DATA (GB( 4,11,IC),IC=1,3) /
     4187C    S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/
     4188C     DATA (GA( 4,12,IC),IC=1,3) /
     4189C    S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/
     4190C     DATA (GB( 4,12,IC),IC=1,3) /
     4191C    S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/
    43724192C
    43734193C----- INTERVAL = 6 ----- T =  237.5
    43744194C
    43754195C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4376       DATA (GA( 5,11,IC),IC=1,3) /
    4377      S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/
    4378       DATA (GB( 5,11,IC),IC=1,3) /
    4379      S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/
    4380       DATA (GA( 5,12,IC),IC=1,3) /
    4381      S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/
    4382       DATA (GB( 5,12,IC),IC=1,3) /
    4383      S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/
     4196C     DATA (GA( 5,11,IC),IC=1,3) /
     4197C    S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/
     4198C     DATA (GB( 5,11,IC),IC=1,3) /
     4199C    S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/
     4200C     DATA (GA( 5,12,IC),IC=1,3) /
     4201C    S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/
     4202C     DATA (GB( 5,12,IC),IC=1,3) /
     4203C    S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/
    43844204C
    43854205C----- INTERVAL = 6 ----- T =  250.0
    43864206C
    43874207C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4388       DATA (GA( 6,11,IC),IC=1,3) /
    4389      S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/
    4390       DATA (GB( 6,11,IC),IC=1,3) /
    4391      S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/
    4392       DATA (GA( 6,12,IC),IC=1,3) /
    4393      S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/
    4394       DATA (GB( 6,12,IC),IC=1,3) /
    4395      S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/
     4208C     DATA (GA( 6,11,IC),IC=1,3) /
     4209C    S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/
     4210C     DATA (GB( 6,11,IC),IC=1,3) /
     4211C    S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/
     4212C     DATA (GA( 6,12,IC),IC=1,3) /
     4213C    S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/
     4214C     DATA (GB( 6,12,IC),IC=1,3) /
     4215C    S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/
    43964216C
    43974217C----- INTERVAL = 6 ----- T =  262.5
    43984218C
    43994219C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4400       DATA (GA( 7,11,IC),IC=1,3) /
    4401      S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/
    4402       DATA (GB( 7,11,IC),IC=1,3) /
    4403      S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/
    4404       DATA (GA( 7,12,IC),IC=1,3) /
    4405      S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/
    4406       DATA (GB( 7,12,IC),IC=1,3) /
    4407      S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/
     4220C     DATA (GA( 7,11,IC),IC=1,3) /
     4221C    S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/
     4222C     DATA (GB( 7,11,IC),IC=1,3) /
     4223C    S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/
     4224C     DATA (GA( 7,12,IC),IC=1,3) /
     4225C    S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/
     4226C     DATA (GB( 7,12,IC),IC=1,3) /
     4227C    S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/
    44084228C
    44094229C----- INTERVAL = 6 ----- T =  275.0
    44104230C
    44114231C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4412       DATA (GA( 8,11,IC),IC=1,3) /
    4413      S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/
    4414       DATA (GB( 8,11,IC),IC=1,3) /
    4415      S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/
    4416       DATA (GA( 8,12,IC),IC=1,3) /
    4417      S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/
    4418       DATA (GB( 8,12,IC),IC=1,3) /
    4419      S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/
     4232C     DATA (GA( 8,11,IC),IC=1,3) /
     4233C    S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/
     4234C     DATA (GB( 8,11,IC),IC=1,3) /
     4235C    S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/
     4236C     DATA (GA( 8,12,IC),IC=1,3) /
     4237C    S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/
     4238C     DATA (GB( 8,12,IC),IC=1,3) /
     4239C    S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/
    44204240C
    44214241C----- INTERVAL = 6 ----- T =  287.5
    44224242C
    44234243C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4424       DATA (GA( 9,11,IC),IC=1,3) /
    4425      S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/
    4426       DATA (GB( 9,11,IC),IC=1,3) /
    4427      S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/
    4428       DATA (GA( 9,12,IC),IC=1,3) /
    4429      S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/
    4430       DATA (GB( 9,12,IC),IC=1,3) /
    4431      S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/
     4244C     DATA (GA( 9,11,IC),IC=1,3) /
     4245C    S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/
     4246C     DATA (GB( 9,11,IC),IC=1,3) /
     4247C    S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/
     4248C     DATA (GA( 9,12,IC),IC=1,3) /
     4249C    S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/
     4250C     DATA (GB( 9,12,IC),IC=1,3) /
     4251C    S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/
    44324252C
    44334253C----- INTERVAL = 6 ----- T =  300.0
    44344254C
    44354255C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4436       DATA (GA(10,11,IC),IC=1,3) /
    4437      S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/
    4438       DATA (GB(10,11,IC),IC=1,3) /
    4439      S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/
    4440       DATA (GA(10,12,IC),IC=1,3) /
    4441      S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/
    4442       DATA (GB(10,12,IC),IC=1,3) /
    4443      S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/
     4256C     DATA (GA(10,11,IC),IC=1,3) /
     4257C    S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/
     4258C     DATA (GB(10,11,IC),IC=1,3) /
     4259C    S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/
     4260C     DATA (GA(10,12,IC),IC=1,3) /
     4261C    S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/
     4262C     DATA (GB(10,12,IC),IC=1,3) /
     4263C    S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/
    44444264C
    44454265C----- INTERVAL = 6 ----- T =  312.5
    44464266C
    44474267C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4448       DATA (GA(11,11,IC),IC=1,3) /
    4449      S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/
    4450       DATA (GB(11,11,IC),IC=1,3) /
    4451      S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/
    4452       DATA (GA(11,12,IC),IC=1,3) /
    4453      S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/
    4454       DATA (GB(11,12,IC),IC=1,3) /
    4455      S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/
     4268C     DATA (GA(11,11,IC),IC=1,3) /
     4269C    S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/
     4270C     DATA (GB(11,11,IC),IC=1,3) /
     4271C    S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/
     4272C     DATA (GA(11,12,IC),IC=1,3) /
     4273C    S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/
     4274C     DATA (GB(11,12,IC),IC=1,3) /
     4275C    S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/
    44564276C
    44574277C
    … …  
    44714291C
    44724292C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4473       DATA (GA( 1,13,IC),IC=1,3) /
    4474      S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/
    4475       DATA (GB( 1,13,IC),IC=1,3) /
    4476      S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/
    4477       DATA (GA( 1,14,IC),IC=1,3) /
    4478      S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/
    4479       DATA (GB( 1,14,IC),IC=1,3) /
    4480      S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/
     4293C     DATA (GA( 1,13,IC),IC=1,3) /
     4294C    S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/
     4295C     DATA (GB( 1,13,IC),IC=1,3) /
     4296C    S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/
     4297C     DATA (GA( 1,14,IC),IC=1,3) /
     4298C    S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/
     4299C     DATA (GB( 1,14,IC),IC=1,3) /
     4300C    S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/
    44814301C
    44824302C----- INTERVAL = 2 ----- T =  200.0
    44834303C
    44844304C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4485       DATA (GA( 2,13,IC),IC=1,3) /
    4486      S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/
    4487       DATA (GB( 2,13,IC),IC=1,3) /
    4488      S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/
    4489       DATA (GA( 2,14,IC),IC=1,3) /
    4490      S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/
    4491       DATA (GB( 2,14,IC),IC=1,3) /
    4492      S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/
     4305C     DATA (GA( 2,13,IC),IC=1,3) /
     4306C    S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/
     4307C     DATA (GB( 2,13,IC),IC=1,3) /
     4308C    S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/
     4309C     DATA (GA( 2,14,IC),IC=1,3) /
     4310C    S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/
     4311C     DATA (GB( 2,14,IC),IC=1,3) /
     4312C    S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/
    44934313C
    44944314C----- INTERVAL = 2 ----- T =  212.5
    44954315C
    44964316C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4497       DATA (GA( 3,13,IC),IC=1,3) /
    4498      S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/
    4499       DATA (GB( 3,13,IC),IC=1,3) /
    4500      S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/
    4501       DATA (GA( 3,14,IC),IC=1,3) /
    4502      S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/
    4503       DATA (GB( 3,14,IC),IC=1,3) /
    4504      S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/
     4317C     DATA (GA( 3,13,IC),IC=1,3) /
     4318C    S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/
     4319C     DATA (GB( 3,13,IC),IC=1,3) /
     4320C    S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/
     4321C     DATA (GA( 3,14,IC),IC=1,3) /
     4322C    S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/
     4323C     DATA (GB( 3,14,IC),IC=1,3) /
     4324C    S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/
    45054325C
    45064326C----- INTERVAL = 2 ----- T =  225.0
    45074327C
    45084328C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4509       DATA (GA( 4,13,IC),IC=1,3) /
    4510      S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/
    4511       DATA (GB( 4,13,IC),IC=1,3) /
    4512      S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/
    4513       DATA (GA( 4,14,IC),IC=1,3) /
    4514      S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/
    4515       DATA (GB( 4,14,IC),IC=1,3) /
    4516      S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/
     4329C     DATA (GA( 4,13,IC),IC=1,3) /
     4330C    S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/
     4331C     DATA (GB( 4,13,IC),IC=1,3) /
     4332C    S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/
     4333C     DATA (GA( 4,14,IC),IC=1,3) /
     4334C    S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/
     4335C     DATA (GB( 4,14,IC),IC=1,3) /
     4336C    S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/
    45174337C
    45184338C----- INTERVAL = 2 ----- T =  237.5
    45194339C
    45204340C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4521       DATA (GA( 5,13,IC),IC=1,3) /
    4522      S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/
    4523       DATA (GB( 5,13,IC),IC=1,3) /
    4524      S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/
    4525       DATA (GA( 5,14,IC),IC=1,3) /
    4526      S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/
    4527       DATA (GB( 5,14,IC),IC=1,3) /
    4528      S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/
     4341C     DATA (GA( 5,13,IC),IC=1,3) /
     4342C    S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/
     4343C     DATA (GB( 5,13,IC),IC=1,3) /
     4344C    S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/
     4345C     DATA (GA( 5,14,IC),IC=1,3) /
     4346C    S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/
     4347C     DATA (GB( 5,14,IC),IC=1,3) /
     4348C    S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/
    45294349C
    45304350C----- INTERVAL = 2 ----- T =  250.0
    45314351C
    45324352C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4533       DATA (GA( 6,13,IC),IC=1,3) /
    4534      S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/
    4535       DATA (GB( 6,13,IC),IC=1,3) /
    4536      S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/
    4537       DATA (GA( 6,14,IC),IC=1,3) /
    4538      S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/
    4539       DATA (GB( 6,14,IC),IC=1,3) /
    4540      S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/
     4353C     DATA (GA( 6,13,IC),IC=1,3) /
     4354C    S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/
     4355C     DATA (GB( 6,13,IC),IC=1,3) /
     4356C    S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/
     4357C     DATA (GA( 6,14,IC),IC=1,3) /
     4358C    S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/
     4359C     DATA (GB( 6,14,IC),IC=1,3) /
     4360C    S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/
    45414361C
    45424362C----- INTERVAL = 2 ----- T =  262.5
    45434363C
    45444364C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4545       DATA (GA( 7,13,IC),IC=1,3) /
    4546      S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/
    4547       DATA (GB( 7,13,IC),IC=1,3) /
    4548      S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/
    4549       DATA (GA( 7,14,IC),IC=1,3) /
    4550      S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/
    4551       DATA (GB( 7,14,IC),IC=1,3) /
    4552      S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/
     4365C     DATA (GA( 7,13,IC),IC=1,3) /
     4366C    S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/
     4367C     DATA (GB( 7,13,IC),IC=1,3) /
     4368C    S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/
     4369C     DATA (GA( 7,14,IC),IC=1,3) /
     4370C    S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/
     4371C     DATA (GB( 7,14,IC),IC=1,3) /
     4372C    S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/
    45534373C
    45544374C----- INTERVAL = 2 ----- T =  275.0
    45554375C
    45564376C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4557       DATA (GA( 8,13,IC),IC=1,3) /
    4558      S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/
    4559       DATA (GB( 8,13,IC),IC=1,3) /
    4560      S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/
    4561       DATA (GA( 8,14,IC),IC=1,3) /
    4562      S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/
    4563       DATA (GB( 8,14,IC),IC=1,3) /
    4564      S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/
     4377C     DATA (GA( 8,13,IC),IC=1,3) /
     4378C    S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/
     4379C     DATA (GB( 8,13,IC),IC=1,3) /
     4380C    S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/
     4381C     DATA (GA( 8,14,IC),IC=1,3) /
     4382C    S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/
     4383C     DATA (GB( 8,14,IC),IC=1,3) /
     4384C    S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/
    45654385C
    45664386C----- INTERVAL = 2 ----- T =  287.5
    45674387C
    45684388C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4569       DATA (GA( 9,13,IC),IC=1,3) /
    4570      S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/
    4571       DATA (GB( 9,13,IC),IC=1,3) /
    4572      S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/
    4573       DATA (GA( 9,14,IC),IC=1,3) /
    4574      S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/
    4575       DATA (GB( 9,14,IC),IC=1,3) /
    4576      S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/
     4389C     DATA (GA( 9,13,IC),IC=1,3) /
     4390C    S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/
     4391C     DATA (GB( 9,13,IC),IC=1,3) /
     4392C    S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/
     4393C     DATA (GA( 9,14,IC),IC=1,3) /
     4394C    S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/
     4395C     DATA (GB( 9,14,IC),IC=1,3) /
     4396C    S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/
    45774397C
    45784398C----- INTERVAL = 2 ----- T =  300.0
    45794399C
    45804400C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4581       DATA (GA(10,13,IC),IC=1,3) /
    4582      S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/
    4583       DATA (GB(10,13,IC),IC=1,3) /
    4584      S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/
    4585       DATA (GA(10,14,IC),IC=1,3) /
    4586      S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/
    4587       DATA (GB(10,14,IC),IC=1,3) /
    4588      S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/
     4401C     DATA (GA(10,13,IC),IC=1,3) /
     4402C    S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/
     4403C     DATA (GB(10,13,IC),IC=1,3) /
     4404C    S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/
     4405C     DATA (GA(10,14,IC),IC=1,3) /
     4406C    S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/
     4407C     DATA (GB(10,14,IC),IC=1,3) /
     4408C    S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/
    45894409C
    45904410C----- INTERVAL = 2 ----- T =  312.5
    45914411C
    45924412C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4593       DATA (GA(11,13,IC),IC=1,3) /
    4594      S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/
    4595       DATA (GB(11,13,IC),IC=1,3) /
    4596      S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/
    4597       DATA (GA(11,14,IC),IC=1,3) /
    4598      S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/
    4599       DATA (GB(11,14,IC),IC=1,3) /
    4600      S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/
     4413C     DATA (GA(11,13,IC),IC=1,3) /
     4414C    S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/
     4415C     DATA (GB(11,13,IC),IC=1,3) /
     4416C    S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/
     4417C     DATA (GA(11,14,IC),IC=1,3) /
     4418C    S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/
     4419C     DATA (GB(11,14,IC),IC=1,3) /
     4420C    S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/
    46014421C
    46024422C
    … …  
    46184438C
    46194439C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4620       DATA (GA( 1,15,IC),IC=1,3) /
    4621      S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/
    4622       DATA (GB( 1,15,IC),IC=1,3) /
    4623      S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/
    4624       DATA (GA( 1,16,IC),IC=1,3) /
    4625      S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/
    4626       DATA (GB( 1,16,IC),IC=1,3) /
    4627      S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/
     4440C     DATA (GA( 1,15,IC),IC=1,3) /
     4441C    S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/
     4442C     DATA (GB( 1,15,IC),IC=1,3) /
     4443C    S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/
     4444C     DATA (GA( 1,16,IC),IC=1,3) /
     4445C    S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/
     4446C     DATA (GB( 1,16,IC),IC=1,3) /
     4447C    S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/
    46284448C
    46294449C----- INTERVAL = 4 ----- T =  200.0
    46304450C
    46314451C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4632       DATA (GA( 2,15,IC),IC=1,3) /
    4633      S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/
    4634       DATA (GB( 2,15,IC),IC=1,3) /
    4635      S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/
    4636       DATA (GA( 2,16,IC),IC=1,3) /
    4637      S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/
    4638       DATA (GB( 2,16,IC),IC=1,3) /
    4639      S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/
     4452C     DATA (GA( 2,15,IC),IC=1,3) /
     4453C    S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/
     4454C     DATA (GB( 2,15,IC),IC=1,3) /
     4455C    S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/
     4456C     DATA (GA( 2,16,IC),IC=1,3) /
     4457C    S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/
     4458C     DATA (GB( 2,16,IC),IC=1,3) /
     4459C    S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/
    46404460C
    46414461C----- INTERVAL = 4 ----- T =  212.5
    46424462C
    46434463C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4644       DATA (GA( 3,15,IC),IC=1,3) /
    4645      S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/
    4646       DATA (GB( 3,15,IC),IC=1,3) /
    4647      S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/
    4648       DATA (GA( 3,16,IC),IC=1,3) /
    4649      S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/
    4650       DATA (GB( 3,16,IC),IC=1,3) /
    4651      S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/
     4464C     DATA (GA( 3,15,IC),IC=1,3) /
     4465C    S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/
     4466C     DATA (GB( 3,15,IC),IC=1,3) /
     4467C    S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/
     4468C     DATA (GA( 3,16,IC),IC=1,3) /
     4469C    S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/
     4470C     DATA (GB( 3,16,IC),IC=1,3) /
     4471C    S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/
    46524472C
    46534473C----- INTERVAL = 4 ----- T =  225.0
    46544474C
    46554475C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4656       DATA (GA( 4,15,IC),IC=1,3) /
    4657      S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/
    4658       DATA (GB( 4,15,IC),IC=1,3) /
    4659      S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/
    4660       DATA (GA( 4,16,IC),IC=1,3) /
    4661      S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/
    4662       DATA (GB( 4,16,IC),IC=1,3) /
    4663      S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/
     4476C     DATA (GA( 4,15,IC),IC=1,3) /
     4477C    S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/
     4478C     DATA (GB( 4,15,IC),IC=1,3) /
     4479C    S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/
     4480C     DATA (GA( 4,16,IC),IC=1,3) /
     4481C    S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/
     4482C     DATA (GB( 4,16,IC),IC=1,3) /
     4483C    S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/
    46644484C
    46654485C----- INTERVAL = 4 ----- T =  237.5
    46664486C
    46674487C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4668       DATA (GA( 5,15,IC),IC=1,3) /
    4669      S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/
    4670       DATA (GB( 5,15,IC),IC=1,3) /
    4671      S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/
    4672       DATA (GA( 5,16,IC),IC=1,3) /
    4673      S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/
    4674       DATA (GB( 5,16,IC),IC=1,3) /
    4675      S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/
     4488C     DATA (GA( 5,15,IC),IC=1,3) /
     4489C    S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/
     4490C     DATA (GB( 5,15,IC),IC=1,3) /
     4491C    S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/
     4492C     DATA (GA( 5,16,IC),IC=1,3) /
     4493C    S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/
     4494C     DATA (GB( 5,16,IC),IC=1,3) /
     4495C    S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/
    46764496C
    46774497C----- INTERVAL = 4 ----- T =  250.0
    46784498C
    46794499C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4680       DATA (GA( 6,15,IC),IC=1,3) /
    4681      S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/
    4682       DATA (GB( 6,15,IC),IC=1,3) /
    4683      S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/
    4684       DATA (GA( 6,16,IC),IC=1,3) /
    4685      S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/
    4686       DATA (GB( 6,16,IC),IC=1,3) /
    4687      S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/
     4500C     DATA (GA( 6,15,IC),IC=1,3) /
     4501C    S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/
     4502C     DATA (GB( 6,15,IC),IC=1,3) /
     4503C    S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/
     4504C     DATA (GA( 6,16,IC),IC=1,3) /
     4505C    S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/
     4506C     DATA (GB( 6,16,IC),IC=1,3) /
     4507C    S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/
    46884508C
    46894509C----- INTERVAL = 4 ----- T =  262.5
    46904510C
    46914511C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4692       DATA (GA( 7,15,IC),IC=1,3) /
    4693      S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/
    4694       DATA (GB( 7,15,IC),IC=1,3) /
    4695      S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/
    4696       DATA (GA( 7,16,IC),IC=1,3) /
    4697      S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/
    4698       DATA (GB( 7,16,IC),IC=1,3) /
    4699      S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/
     4512C     DATA (GA( 7,15,IC),IC=1,3) /
     4513C    S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/
     4514C     DATA (GB( 7,15,IC),IC=1,3) /
     4515C    S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/
     4516C     DATA (GA( 7,16,IC),IC=1,3) /
     4517C    S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/
     4518C     DATA (GB( 7,16,IC),IC=1,3) /
     4519C    S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/
    47004520C
    47014521C----- INTERVAL = 4 ----- T =  275.0
    47024522C
    47034523C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4704       DATA (GA( 8,15,IC),IC=1,3) /
    4705      S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/
    4706       DATA (GB( 8,15,IC),IC=1,3) /
    4707      S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/
    4708       DATA (GA( 8,16,IC),IC=1,3) /
    4709      S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/
    4710       DATA (GB( 8,16,IC),IC=1,3) /
    4711      S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/
     4524C     DATA (GA( 8,15,IC),IC=1,3) /
     4525C    S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/
     4526C     DATA (GB( 8,15,IC),IC=1,3) /
     4527C    S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/
     4528C     DATA (GA( 8,16,IC),IC=1,3) /
     4529C    S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/
     4530C     DATA (GB( 8,16,IC),IC=1,3) /
     4531C    S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/
    47124532C
    47134533C----- INTERVAL = 4 ----- T =  287.5
    47144534C
    47154535C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4716       DATA (GA( 9,15,IC),IC=1,3) /
    4717      S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/
    4718       DATA (GB( 9,15,IC),IC=1,3) /
    4719      S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/
    4720       DATA (GA( 9,16,IC),IC=1,3) /
    4721      S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/
    4722       DATA (GB( 9,16,IC),IC=1,3) /
    4723      S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/
     4536C     DATA (GA( 9,15,IC),IC=1,3) /
     4537C    S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/
     4538C     DATA (GB( 9,15,IC),IC=1,3) /
     4539C    S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/
     4540C     DATA (GA( 9,16,IC),IC=1,3) /
     4541C    S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/
     4542C     DATA (GB( 9,16,IC),IC=1,3) /
     4543C    S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/
    47244544C
    47254545C----- INTERVAL = 4 ----- T =  300.0
    47264546C
    47274547C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4728       DATA (GA(10,15,IC),IC=1,3) /
    4729      S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/
    4730       DATA (GB(10,15,IC),IC=1,3) /
    4731      S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/
    4732       DATA (GA(10,16,IC),IC=1,3) /
    4733      S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/
    4734       DATA (GB(10,16,IC),IC=1,3) /
    4735      S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/
     4548C     DATA (GA(10,15,IC),IC=1,3) /
     4549C    S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/
     4550C     DATA (GB(10,15,IC),IC=1,3) /
     4551C    S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/
     4552C     DATA (GA(10,16,IC),IC=1,3) /
     4553C    S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/
     4554C     DATA (GB(10,16,IC),IC=1,3) /
     4555C    S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/
    47364556C
    47374557C----- INTERVAL = 4 ----- T =  312.5
    47384558C
    47394559C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4740       DATA (GA(11,15,IC),IC=1,3) /
    4741      S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/
    4742       DATA (GB(11,15,IC),IC=1,3) /
    4743      S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/
    4744       DATA (GA(11,16,IC),IC=1,3) /
    4745      S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/
    4746       DATA (GB(11,16,IC),IC=1,3) /
    4747      S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/
     4560C     DATA (GA(11,15,IC),IC=1,3) /
     4561C    S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/
     4562C     DATA (GB(11,15,IC),IC=1,3) /
     4563C    S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/
     4564C     DATA (GA(11,16,IC),IC=1,3) /
     4565C    S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/
     4566C     DATA (GB(11,16,IC),IC=1,3) /
     4567C    S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/
     4568C
     4569C     ------------------------------------------------------------------
     4570C     DATA (( XP(  J,K),J=1,6),       K=1,6) /
     4571C    S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,
     4572C    S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,
     4573C    S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,
     4574C    S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,
     4575C    S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,
     4576C    S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,
     4577C    S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,
     4578C    S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,
     4579C    S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,
     4580C    S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,
     4581C    S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,
     4582C    S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /
    47484583
    4749 C     ------------------------------------------------------------------
    4750       DATA (( XP(  J,K),J=1,6),       K=1,6) /
    4751      S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,
    4752      S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,
    4753      S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,
    4754      S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,
    4755      S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,
    4756      S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,
    4757      S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,
    4758      S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,
    4759      S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,
    4760      S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,
    4761      S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,
    4762      S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /
    47634584C
    47644585C
    … …  
    47684589 100  CONTINUE
    47694590C
     4591!cdir collapse
    47704592      DO 102 JK = 1 , KFLEV+1
    47714593      DO 101 JL = 1, KDLON
    … …  
    49594781      INTEGER KUAER,KTRAER, KLIM
    49604782C
    4961       REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
    4962       REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
    4963       REAL*8 PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
    4964       REAL*8 PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
    4965       REAL*8 PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
    4966       REAL*8 PBTOP(KDLON,Ninter) ! T.O.A. SPECTRAL PLANCK FUNCTION
    4967       REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
    4968       REAL*8 PEMIS(KDLON) ! SURFACE EMISSIVITY
    4969       REAL*8 PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
    4970       REAL*8 PTAVE(KDLON,KFLEV) ! TEMPERATURE
    4971       REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    4972       REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    4973       REAL*8 PGASUR(KDLON,8,2) ! PADE APPROXIMANTS
    4974       REAL*8 PGBSUR(KDLON,8,2) ! PADE APPROXIMANTS
    4975       REAL*8 PGATOP(KDLON,8,2) ! PADE APPROXIMANTS
    4976       REAL*8 PGBTOP(KDLON,8,2) ! PADE APPROXIMANTS
    4977 C
    4978       REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
    4979       REAL*8 PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
    4980       REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
     4783      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! EFFECTIVE ABSORBER AMOUNTS
     4784      REAL(KIND=8) PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
     4785      REAL(KIND=8) PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
     4786      REAL(KIND=8) PBSUR(KDLON,Ninter) ! SURFACE SPECTRAL PLANCK FUNCTION
     4787      REAL(KIND=8) PBSUIN(KDLON) ! SURFACE PLANCK FUNCTION
     4788      REAL(KIND=8) PBTOP(KDLON,Ninter) ! T.O.A. SPECTRAL PLANCK FUNCTION
     4789      REAL(KIND=8) PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
     4790      REAL(KIND=8) PEMIS(KDLON) ! SURFACE EMISSIVITY
     4791      REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! HALF-LEVEL PRESSURE (MB)
     4792      REAL(KIND=8) PTAVE(KDLON,KFLEV) ! TEMPERATURE
     4793      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     4794      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     4795      REAL(KIND=8) PGASUR(KDLON,8,2) ! PADE APPROXIMANTS
     4796      REAL(KIND=8) PGBSUR(KDLON,8,2) ! PADE APPROXIMANTS
     4797      REAL(KIND=8) PGATOP(KDLON,8,2) ! PADE APPROXIMANTS
     4798      REAL(KIND=8) PGBTOP(KDLON,8,2) ! PADE APPROXIMANTS
     4799C
     4800      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
     4801      REAL(KIND=8) PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
     4802      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
    49814803C-----------------------------------------------------------------------
    49824804C LOCAL VARIABLES:
    4983       REAL*8 ZADJD(KDLON,KFLEV+1)
    4984       REAL*8 ZADJU(KDLON,KFLEV+1)
    4985       REAL*8 ZDBDT(KDLON,Ninter,KFLEV)
    4986       REAL*8 ZDISD(KDLON,KFLEV+1)
    4987       REAL*8 ZDISU(KDLON,KFLEV+1)
     4805      REAL(KIND=8) ZADJD(KDLON,KFLEV+1)
     4806      REAL(KIND=8) ZADJU(KDLON,KFLEV+1)
     4807      REAL(KIND=8) ZDBDT(KDLON,Ninter,KFLEV)
     4808      REAL(KIND=8) ZDISD(KDLON,KFLEV+1)
     4809      REAL(KIND=8) ZDISU(KDLON,KFLEV+1)
    49884810C
    49894811      INTEGER jk, jl
    … …  
    50764898      INTEGER KUAER,KTRAER, KLIM
    50774899C
    5078       REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
    5079       REAL*8 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
    5080       REAL*8 PADJU(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
    5081       REAL*8 PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
    5082       REAL*8 PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
    5083       REAL*8 PBSUR(KDLON,Ninter) ! SPECTRAL SURFACE PLANCK FUNCTION
    5084       REAL*8 PBSUI(KDLON) ! SURFACE PLANCK FUNCTION
    5085       REAL*8 PBTOP(KDLON,Ninter) ! SPECTRAL T.O.A. PLANCK FUNCTION
    5086       REAL*8 PDISD(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
    5087       REAL*8 PDISU(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
    5088       REAL*8 PEMIS(KDLON) ! SURFACE EMISSIVITY
    5089       REAL*8 PPMB(KDLON,KFLEV+1) ! PRESSURE MB
    5090       REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    5091       REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    5092       REAL*8 PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
    5093       REAL*8 PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
    5094       REAL*8 PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
    5095       REAL*8 PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
    5096 C
    5097       REAL*8 PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
    5098       REAL*8 PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
     4900      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
     4901      REAL(KIND=8) PADJD(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
     4902      REAL(KIND=8) PADJU(KDLON,KFLEV+1) ! CONTRIBUTION BY ADJACENT LAYERS
     4903      REAL(KIND=8) PB(KDLON,Ninter,KFLEV+1) ! SPECTRAL HALF-LEVEL PLANCK FUNCTIONS
     4904      REAL(KIND=8) PBINT(KDLON,KFLEV+1) ! HALF-LEVEL PLANCK FUNCTIONS
     4905      REAL(KIND=8) PBSUR(KDLON,Ninter) ! SPECTRAL SURFACE PLANCK FUNCTION
     4906      REAL(KIND=8) PBSUI(KDLON) ! SURFACE PLANCK FUNCTION
     4907      REAL(KIND=8) PBTOP(KDLON,Ninter) ! SPECTRAL T.O.A. PLANCK FUNCTION
     4908      REAL(KIND=8) PDISD(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
     4909      REAL(KIND=8) PDISU(KDLON,KFLEV+1) ! CONTRIBUTION BY DISTANT LAYERS
     4910      REAL(KIND=8) PEMIS(KDLON) ! SURFACE EMISSIVITY
     4911      REAL(KIND=8) PPMB(KDLON,KFLEV+1) ! PRESSURE MB
     4912      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     4913      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     4914      REAL(KIND=8) PGASUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
     4915      REAL(KIND=8) PGBSUR(KDLON,8,2) ! SURFACE PADE APPROXIMANTS
     4916      REAL(KIND=8) PGATOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
     4917      REAL(KIND=8) PGBTOP(KDLON,8,2) ! T.O.A. PADE APPROXIMANTS
     4918C
     4919      REAL(KIND=8) PFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES
     4920      REAL(KIND=8) PCTS(KDLON,KFLEV) ! COOLING-TO-SPACE TERM
    50994921C
    51004922C* LOCAL VARIABLES:
    51014923C
    5102       REAL*8 ZBGND(KDLON)
    5103       REAL*8 ZFD(KDLON)
    5104       REAL*8  ZFN10(KDLON)
    5105       REAL*8 ZFU(KDLON)
    5106       REAL*8  ZTT(KDLON,NTRA)
    5107       REAL*8 ZTT1(KDLON,NTRA)
    5108       REAL*8 ZTT2(KDLON,NTRA)
    5109       REAL*8  ZUU(KDLON,NUA)
    5110       REAL*8 ZCNSOL(KDLON)
    5111       REAL*8 ZCNTOP(KDLON)
     4924      REAL(KIND=8) ZBGND(KDLON)
     4925      REAL(KIND=8) ZFD(KDLON)
     4926      REAL(KIND=8)  ZFN10(KDLON)
     4927      REAL(KIND=8) ZFU(KDLON)
     4928      REAL(KIND=8)  ZTT(KDLON,NTRA)
     4929      REAL(KIND=8) ZTT1(KDLON,NTRA)
     4930      REAL(KIND=8) ZTT2(KDLON,NTRA)
     4931      REAL(KIND=8)  ZUU(KDLON,NUA)
     4932      REAL(KIND=8) ZCNSOL(KDLON)
     4933      REAL(KIND=8) ZCNTOP(KDLON)
    51124934C
    51134935      INTEGER jk, jl, ja
    51144936      INTEGER jstra, jstru
    51154937      INTEGER ind1, ind2, ind3, ind4, in, jlim
    5116       REAL*8 zctstr
     4938      REAL(KIND=8) zctstr
    51174939C-----------------------------------------------------------------------
    51184940C
    … …  
    54015223      INTEGER KUAER,KTRAER
    54025224C
    5403       REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
    5404       REAL*8 PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT
    5405       REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    5406       REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    5407 C
    5408       REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! ENERGY EXCHANGE MATRIX
    5409       REAL*8 PDISD(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
    5410       REAL*8 PDISU(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
     5225      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
     5226      REAL(KIND=8) PDBDT(KDLON,Ninter,KFLEV) ! LAYER PLANCK FUNCTION GRADIENT
     5227      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     5228      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     5229C
     5230      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! ENERGY EXCHANGE MATRIX
     5231      REAL(KIND=8) PDISD(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
     5232      REAL(KIND=8) PDISU(KDLON,KFLEV+1) !  CONTRIBUTION BY DISTANT LAYERS
    54115233C
    54125234C* LOCAL VARIABLES:
    54135235C
    5414       REAL*8 ZGLAYD(KDLON)
    5415       REAL*8 ZGLAYU(KDLON)
    5416       REAL*8 ZTT(KDLON,NTRA)
    5417       REAL*8 ZTT1(KDLON,NTRA)
    5418       REAL*8 ZTT2(KDLON,NTRA)
     5236      REAL(KIND=8) ZGLAYD(KDLON)
     5237      REAL(KIND=8) ZGLAYU(KDLON)
     5238      REAL(KIND=8) ZTT(KDLON,NTRA)
     5239      REAL(KIND=8) ZTT1(KDLON,NTRA)
     5240      REAL(KIND=8) ZTT2(KDLON,NTRA)
    54195241C
    54205242      INTEGER jl, jk, ja, ikp1, ikn, ikd1, jkj, ikd2
    54215243      INTEGER ikjp1, ikm1, ikj, jlk, iku1, ijkl, iku2
    54225244      INTEGER ind1, ind2, ind3, ind4, itt
    5423       REAL*8 zww, zdzxdg, zdzxmg
     5245      REAL(KIND=8) zww, zdzxdg, zdzxmg
    54245246C
    54255247C*         1.    INITIALIZATION
    … …  
    56195441     S  , PADJD,PADJU,PCNTRB,PDBDT)
    56205442       USE dimphy
     5443      USE radiation_AR4_param, only : WG1
    56215444      IMPLICIT none
    56225445cym#include "dimensions.h"
    … …  
    56565479      INTEGER KUAER,KTRAER
    56575480C
    5658       REAL*8 PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
    5659       REAL*8 PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
    5660       REAL*8 PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    5661       REAL*8 PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
    5662 C
    5663       REAL*8 PADJD(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
    5664       REAL*8 PADJU(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
    5665       REAL*8 PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
    5666       REAL*8 PDBDT(KDLON,Ninter,KFLEV) !  LAYER PLANCK FUNCTION GRADIENT
     5481      REAL(KIND=8) PABCU(KDLON,NUA,3*KFLEV+1) ! ABSORBER AMOUNTS
     5482      REAL(KIND=8) PDBSL(KDLON,Ninter,KFLEV*2) ! SUB-LAYER PLANCK FUNCTION GRADIENT
     5483      REAL(KIND=8) PGA(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     5484      REAL(KIND=8) PGB(KDLON,8,2,KFLEV) ! PADE APPROXIMANTS
     5485C
     5486      REAL(KIND=8) PADJD(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
     5487      REAL(KIND=8) PADJU(KDLON,KFLEV+1) ! CONTRIBUTION OF ADJACENT LAYERS
     5488      REAL(KIND=8) PCNTRB(KDLON,KFLEV+1,KFLEV+1) ! CLEAR-SKY ENERGY EXCHANGE MATRIX
     5489      REAL(KIND=8) PDBDT(KDLON,Ninter,KFLEV) !  LAYER PLANCK FUNCTION GRADIENT
    56675490C
    56685491C* LOCAL ARRAYS:
    56695492C
    5670       REAL*8 ZGLAYD(KDLON)
    5671       REAL*8 ZGLAYU(KDLON)
    5672       REAL*8 ZTT(KDLON,NTRA)
    5673       REAL*8 ZTT1(KDLON,NTRA)
    5674       REAL*8 ZTT2(KDLON,NTRA)
    5675       REAL*8 ZUU(KDLON,NUA)
     5493      REAL(KIND=8) ZGLAYD(KDLON)
     5494      REAL(KIND=8) ZGLAYU(KDLON)
     5495      REAL(KIND=8) ZTT(KDLON,NTRA)
     5496      REAL(KIND=8) ZTT1(KDLON,NTRA)
     5497      REAL(KIND=8) ZTT2(KDLON,NTRA)
     5498      REAL(KIND=8) ZUU(KDLON,NUA)
    56765499C
    56775500      INTEGER jk, jl, ja, im12, ind, inu, ixu, jg
    56785501      INTEGER ixd, ibs, idd, imu, jk1, jk2, jnu
    5679       REAL*8 zwtr
     5502      REAL(KIND=8) zwtr
    56805503c
    5681 C* Data Block:
    5682 c
    5683       REAL*8 WG1(2)
    5684       SAVE WG1
    5685 c$OMP THREADPRIVATE(WG1)
    5686       DATA (WG1(jk),jk=1,2) /1.0, 1.0/
     5504
    56875505C-----------------------------------------------------------------------
    56885506C
    … …  
    58655683C
    58665684C-----------------------------------------------------------------------
    5867       REAL*8 O1H, O2H
     5685      REAL(KIND=8) O1H, O2H
    58685686      PARAMETER (O1H=2230.)
    58695687      PARAMETER (O2H=100.)
    5870       REAL*8 RPIALF0
     5688      REAL(KIND=8) RPIALF0
    58715689      PARAMETER (RPIALF0=2.0)
    58725690C
    58735691C* ARGUMENTS:
    58745692C
    5875       REAL*8 PUU(KDLON,NUA)
    5876       REAL*8 PTT(KDLON,NTRA)
    5877       REAL*8 PGA(KDLON,8,2)
    5878       REAL*8 PGB(KDLON,8,2)
     5693      REAL(KIND=8) PUU(KDLON,NUA)
     5694      REAL(KIND=8) PTT(KDLON,NTRA)
     5695      REAL(KIND=8) PGA(KDLON,8,2)
     5696      REAL(KIND=8) PGB(KDLON,8,2)
    58795697C
    58805698C* LOCAL VARIABLES:
    58815699C
    5882       REAL*8 zz, zxd, zxn
    5883       REAL*8 zpu, zpu10, zpu11, zpu12, zpu13
    5884       REAL*8 zeu, zeu10, zeu11, zeu12, zeu13
    5885       REAL*8 zx, zy, zsq1, zsq2, zvxy, zuxy
    5886       REAL*8 zaercn, zto1, zto2, zxch4, zych4, zxn2o, zyn2o
    5887       REAL*8 zsqn21, zodn21, zsqh42, zodh42
    5888       REAL*8 zsqh41, zodh41, zsqn22, zodn22, zttf11, zttf12
    5889       REAL*8 zuu11, zuu12, za11, za12
     5700      REAL(KIND=8) zz, zxd, zxn
     5701      REAL(KIND=8) zpu, zpu10, zpu11, zpu12, zpu13
     5702      REAL(KIND=8) zeu, zeu10, zeu11, zeu12, zeu13
     5703      REAL(KIND=8) zx, zy, zsq1, zsq2, zvxy, zuxy
     5704      REAL(KIND=8) zaercn, zto1, zto2, zxch4, zych4, zxn2o, zyn2o
     5705      REAL(KIND=8) zsqn21, zodn21, zsqh42, zodh42
     5706      REAL(KIND=8) zsqh41, zodh41, zsqn22, zodn22, zttf11, zttf12
     5707      REAL(KIND=8) zuu11, zuu12, za11, za12
    58905708      INTEGER jl, ja
    58915709C     ------------------------------------------------------------------
    … …  
    58975715C
    58985716C
     5717!cdir collapse
    58995718      DO 130 JA = 1 , 8
    59005719      DO 120 JL = 1, KDLON
    … …  
    60415860C
    60425861C-----------------------------------------------------------------------
    6043       REAL*8 O1H, O2H
     5862      REAL(KIND=8) O1H, O2H
    60445863      PARAMETER (O1H=2230.)
    60455864      PARAMETER (O2H=100.)
    6046       REAL*8 RPIALF0
     5865      REAL(KIND=8) RPIALF0
    60475866      PARAMETER (RPIALF0=2.0)
    60485867C
    60495868C* ARGUMENTS:
    60505869C
    6051       REAL*8 PGA(KDLON,8,2) ! PADE APPROXIMANTS
    6052       REAL*8 PGB(KDLON,8,2) ! PADE APPROXIMANTS
    6053       REAL*8 PUU1(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 1
    6054       REAL*8 PUU2(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 2
    6055       REAL*8 PTT(KDLON,NTRA) ! TRANSMISSION FUNCTIONS
     5870      REAL(KIND=8) PGA(KDLON,8,2) ! PADE APPROXIMANTS
     5871      REAL(KIND=8) PGB(KDLON,8,2) ! PADE APPROXIMANTS
     5872      REAL(KIND=8) PUU1(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 1
     5873      REAL(KIND=8) PUU2(KDLON,NUA) ! ABSORBER AMOUNTS FROM TOP TO LEVEL 2
     5874      REAL(KIND=8) PTT(KDLON,NTRA) ! TRANSMISSION FUNCTIONS
    60565875C
    60575876C* LOCAL VARIABLES:
    60585877C
    60595878      INTEGER ja, jl
    6060       REAL*8 zz, zxd, zxn
    6061       REAL*8 zpu, zpu10, zpu11, zpu12, zpu13
    6062       REAL*8 zeu, zeu10, zeu11, zeu12, zeu13
    6063       REAL*8 zx, zy, zuxy, zsq1, zsq2, zvxy, zaercn, zto1, zto2
    6064       REAL*8 zxch4, zych4, zsqh41, zodh41
    6065       REAL*8 zxn2o, zyn2o, zsqn21, zodn21, zsqh42, zodh42
    6066       REAL*8 zsqn22, zodn22, za11, zttf11, za12, zttf12
    6067       REAL*8 zuu11, zuu12
     5879      REAL(KIND=8) zz, zxd, zxn
     5880      REAL(KIND=8) zpu, zpu10, zpu11, zpu12, zpu13
     5881      REAL(KIND=8) zeu, zeu10, zeu11, zeu12, zeu13
     5882      REAL(KIND=8) zx, zy, zuxy, zsq1, zsq2, zvxy, zaercn, zto1, zto2
     5883      REAL(KIND=8) zxch4, zych4, zsqh41, zodh41
     5884      REAL(KIND=8) zxn2o, zyn2o, zsqn21, zodn21, zsqh42, zodh42
     5885      REAL(KIND=8) zsqn22, zodn22, za11, zttf11, za12, zttf12
     5886      REAL(KIND=8) zuu11, zuu12
    60685887C     ------------------------------------------------------------------
    60695888C
    … …  
    60745893C
    60755894C
     5895
     5896!CDIR ON_ADB(PUU1)
     5897!CDIR ON_ADB(PUU2)
     5898!CDIR COLLAPSE
    60765899      DO 130 JA = 1 , 8
    60775900      DO 120 JL = 1, KDLON
Note: See TracChangeset for help on using the changeset viewer.