Changeset 1250


Ignore:
Timestamp:
Oct 29, 2009, 2:55:23 PM (15 years ago)
Author:
yann meurdesoif
Message:

Optimisations SX9

YM

Location:
LMDZ4/branches/LMDZ4-dev
Files:
7 added
1 deleted
16 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/bld.cfg

    r1228 r1250  
    2727src::filtrez %SRC_PATH/filtrez
    2828src::bibio   %SRC_PATH/bibio
    29 src::ext_src %EXT_SRC
    3029
    3130bld::lib::dyn      %DYN
     
    3938bld::target              lib%{DYN}.a lib%{PHYS}.a libgrid.a libfiltrez.a libbibio.a
    4039bld::target              %EXEC%SUFF_NAME.e
    41 bld::exe_dep             %{DYN} %{PHYS} grid filtrez bibio ext_src
     40bld::exe_dep             %{DYN} %{PHYS} grid filtrez bibio
    4241
    4342
     
    6059bld::tool::fppkeys   %CPP_KEY %FPP_DEF
    6160
     61
     62#bld::tool::fflags::phys::readaerosol         %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt -pi auto
     63#bld::tool::fflags::phys::aeropt_2bands       %BASE_FFLAGS %PROD_FFLAGS  %INCDIR
     64#bld::tool::fflags::phys::radiation_AR4       %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -Wf,-O,extendreorder
     65#bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -f3
     66#bld::tool::fflags::phys::fisrtilp            %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
     67#bld::tool::fflags::phys::cv30_routines       %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -Wf,-O,extendreorder
     68#bld::tool::fflags::phys::cvltr               %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
     69#bld::tool::fflags::phys::clouds_gno          %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
     70#bld::tool::fflags::dyn::vlsplt_p             %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
     71#bld::tool::fflags::dyn::groupeun_p           %BASE_FFLAGS %PROD_FFLAGS  %INCDIR -C hopt
     72
     73
     74inc arch.opt
    6275
    6376# Pre-process code before analysing dependencies
  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F

    r1231 r1250  
    158158      REAL,ALLOCATABLE,SAVE :: pcvgt(:,:), pcvgq(:,:,:)
    159159c
    160 c      REAL,ALLOCATABLE,SAVE :: pvervel(:,:)
    161 c
    162160      REAL,ALLOCATABLE,SAVE :: zdufi(:,:),zdvfi(:,:)
    163161      REAL,ALLOCATABLE,SAVE :: zdtfi(:,:),zdqfi(:,:,:)
     
    175173      REAL,ALLOCATABLE,SAVE :: ztfi_omp(:,:)
    176174      REAL,ALLOCATABLE,SAVE :: zqfi_omp(:,:,:)
    177 c      REAL,ALLOCATABLE,SAVE :: pvervel_omp(:,:)
    178175      REAL,ALLOCATABLE,SAVE :: zdufi_omp(:,:)
    179176      REAL,ALLOCATABLE,SAVE :: zdvfi_omp(:,:)
     
    252249      ALLOCATE(pcvgu(klon,llm), pcvgv(klon,llm))
    253250      ALLOCATE(pcvgt(klon,llm), pcvgq(klon,llm,2))
    254 c      ALLOCATE(pvervel(klon,llm))
    255251      ALLOCATE(zdufi(klon,llm),zdvfi(klon,llm))
    256252      ALLOCATE(zdtfi(klon,llm),zdqfi(klon,llm,nqtot))
     
    278274
    279275c$OMP MASTER             
     276!CDIR ON_ADB(index_i)
     277!CDIR ON_ADB(index_j)
    280278      do ig0=1,klon
    281279        i=index_i(ig0)
     
    300298c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    301299      DO l = 1, llmp1
     300!CDIR ON_ADB(index_i)
     301!CDIR ON_ADB(index_j)
    302302        do ig0=1,klon
    303303          i=index_i(ig0)
     
    314314c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    315315      DO l=1,llm
    316 
     316!CDIR ON_ADB(index_i)
     317!CDIR ON_ADB(index_j)
    317318        do ig0=1,klon
    318319          i=index_i(ig0)
     
    321322          zplay(ig0,l)   = preff * pksurcp ** unskap
    322323          ztfi(ig0,l)    = pteta(i,j,l)  * pksurcp
    323 c          pcvgt(ig0,l)   = pdteta(i,j,l) * pksurcp / pmasse(i,j,l)
    324324        enddo
    325325
     
    335335c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    336336         DO l=1,llm
     337!CDIR ON_ADB(index_i)
     338!CDIR ON_ADB(index_j)
    337339           do ig0=1,klon
    338340             i=index_i(ig0)
     
    344346      ENDDO
    345347
    346 c   convergence dynamique pour les traceurs "EAU"
    347 ! Earth-specific treatment of first 2 tracers (water)
    348       if (planet_type=="earth") then
    349        DO iq=1,2
    350 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
    351          DO l=1,llm
    352            do ig0=1,klon
    353              i=index_i(ig0)
    354              j=index_j(ig0)
    355 c             pcvgq(ig0,l,iq) = pdq(i,j,l,iq) / pmasse(i,j,l)
    356            enddo
    357          ENDDO
    358 c$OMP END DO NOWAIT     
    359        ENDDO
    360       endif ! of if (planet_type=="earth")
    361 
    362348
    363349c   Geopotentiel calcule par rapport a la surface locale:
     
    378364c$OMP END DO NOWAIT
    379365     
    380 c   ....  Calcul de la vitesse  verticale  ( en Pa*m*s  ou Kg/s )  ....
    381 c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux
    382 c      de masse est calclue dans advtrac_p.F 
    383 c
    384 cc$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
    385 c      DO l=1,llm
    386 c        do ig0=1,klon
    387 c           i=index_i(ig0)
    388 c           j=index_j(ig0)
    389 c           pvervel(ig0,l) = pw(i,j,l)*g* unsaire(i,j)
    390 c        enddo
    391 c       if (is_north_pole) pvervel(1,l)=pw(1,1,l)*g /apoln
    392 c       if (is_south_pole) pvervel(klon,l)=pw(1,jjp1,l)*g/apols
    393 c      ENDDO
    394 cc$OMP END DO NOWAIT
    395366
    396367c
     
    406377c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
    407378      DO l=1,llm
     379!CDIR ON_ADB(index_i)
     380!CDIR ON_ADB(index_j)
     381!CDIR SPARSE
    408382        do ig0=kstart,kend
    409383          i=index_i(ig0)
     
    412386            zufi(ig0,l)= 0.5 *(  pucov(iim,j,l)/cu(iim,j)
    413387     $                         + pucov(1,j,l)/cu(1,j) )
    414 c            pcvgu(ig0,l)= 0.5*(  pducov(iim,j,l)/cu(iim,j)
    415 c     $                         + pducov(1,j,l)/cu(1,j) )
    416388          else
    417389            zufi(ig0,l)= 0.5*(  pucov(i-1,j,l)/cu(i-1,j)
    418390     $                       + pucov(i,j,l)/cu(i,j) )
    419 c            pcvgu(ig0,l)= 0.5*(  pducov(i-1,j,l)/cu(i-1,j)
    420 c     $                        + pducov(i,j,l)/cu(i,j) )
    421391          endif
    422392        enddo
    423393      ENDDO
    424394c$OMP END DO NOWAIT
     395
    425396c   46.champ v:
    426397c   -----------
     398
    427399c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    428400      DO l=1,llm
     401!CDIR ON_ADB(index_i)
     402!CDIR ON_ADB(index_j)
    429403        DO ig0=kstart,kend
    430404          i=index_i(ig0)
     
    433407     $                       + pvcov(i,j,l)/cv(i,j) )
    434408   
    435 c          pcvgv(ig0+i,l)= 0.5 * (  pdvcov(i,j-1,l)/cv(i,j-1)
    436 c     $                           + pdvcov(i,j,l)/cv(i,j) )
    437409         ENDDO
    438410      ENDDO
     
    449421
    450422           z1(1)   =(rlonu(1)-rlonu(iim)+2.*pi)*pvcov(1,1,l)/cv(1,1)
    451 c           z1bis(1)=(rlonu(1)-rlonu(iim)+2.*pi)*pdvcov(1,1,l)/cv(1,1)
    452423           DO i=2,iim
    453424              z1(i)   =(rlonu(i)-rlonu(i-1))*pvcov(i,1,l)/cv(i,1)
    454 c              z1bis(i)=(rlonu(i)-rlonu(i-1))*pdvcov(i,1,l)/cv(i,1)
    455425           ENDDO
    456426 
    457427           DO i=1,iim
    458428              zcos(i)   = COS(rlonv(i))*z1(i)
    459 c              zcosbis(i)= COS(rlonv(i))*z1bis(i)
    460429              zsin(i)   = SIN(rlonv(i))*z1(i)
    461 c              zsinbis(i)= SIN(rlonv(i))*z1bis(i)
    462430           ENDDO
    463431 
    464432           zufi(1,l)  = SSUM(iim,zcos,1)/pi
    465 c           pcvgu(1,l) = SSUM(iim,zcosbis,1)/pi
    466433           zvfi(1,l)  = SSUM(iim,zsin,1)/pi
    467 c           pcvgv(1,l) = SSUM(iim,zsinbis,1)/pi
    468434 
    469435        ENDDO
     
    482448 
    483449         z1(1)   =(rlonu(1)-rlonu(iim)+2.*pi)*pvcov(1,jjm,l)/cv(1,jjm)
    484 c         z1bis(1)=(rlonu(1)-rlonu(iim)+2.*pi)*pdvcov(1,jjm,l)/cv(1,jjm)
    485450           DO i=2,iim
    486            z1(i)   =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm)
    487 c           z1bis(i)=(rlonu(i)-rlonu(i-1))*pdvcov(i,jjm,l)/cv(i,jjm)
     451             z1(i)   =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm)
    488452           ENDDO
    489453 
    490454           DO i=1,iim
    491455              zcos(i)    = COS(rlonv(i))*z1(i)
    492 c              zcosbis(i) = COS(rlonv(i))*z1bis(i)
    493456              zsin(i)    = SIN(rlonv(i))*z1(i)
    494 c              zsinbis(i) = SIN(rlonv(i))*z1bis(i)
    495457           ENDDO
    496458 
    497459           zufi(klon,l)  = SSUM(iim,zcos,1)/pi
    498 c           pcvgu(klon,l) = SSUM(iim,zcosbis,1)/pi
    499460           zvfi(klon,l)  = SSUM(iim,zsin,1)/pi
    500 c           pcvgv(klon,l) = SSUM(iim,zsinbis,1)/pi
    501 
    502461        ENDDO
    503462c$OMP END DO NOWAIT       
     
    506465
    507466      IF (is_sequential) THEN
    508         if (planet_type=="earth") then
    509 #ifdef CPP_EARTH
    510467c
    511468cIM calcul PV a teta=350, 380, 405K
     
    514471     $           ntetaSTD,rtetaSTD,PVteta)
    515472c
    516 #endif
    517         endif
    518473      ENDIF
    519474
     
    525480c   ---------------------
    526481
    527 cc$OMP  PARALLEL DEFAULT(NONE)
    528 cc$OMP+ PRIVATE(i,l,offset,iq)
    529 cc$OMP+ SHARED(klon_omp_nb,nqtot,klon_omp_begin,
    530 cc$OMP+        debut,lafin,rdayvrai,heure,dtphys,zplev,zplay,
    531 cc$OMP+        zphi,zphis,presnivs,clesphy0,zufi,zvfi,ztfi,
    532 cc$OMP+        zqfi,pvervel,zdufi,zdvfi,zdtfi,zdqfi,zdpsrf)
    533 
    534 c PRIVATE(zplev_omp,zplay_omp,zphi_omp,zphis_omp,
    535 c c$OMP+                 presnivs_omp,zufi_omp,zvfi_omp,ztfi_omp,
    536 c c$OMP+                 zqfi_omp,pvervel_omp,zdufi_omp,zdvfi_omp,
    537 c c$OMP+                 zdtfi_omp,zdqfi_omp,zdpsrf_omp)
    538482
    539483c$OMP BARRIER
     
    550494        allocate(ztfi_omp(klon,llm))
    551495        allocate(zqfi_omp(klon,llm,nqtot))
    552 c        allocate(pvervel_omp(klon,llm))
    553496        allocate(zdufi_omp(klon,llm))
    554497        allocate(zdvfi_omp(klon,llm))
     
    617560      enddo
    618561       
    619 c      do l=1,llm
    620 c        do i=1,klon
    621 c         pvervel_omp(i,l)=pvervel(offset+i,l)
    622 c       enddo
    623 c      enddo
    624        
    625562      do l=1,llm
    626563        do i=1,klon
     
    660597     
    661598c$OMP BARRIER
    662 cym      call WriteField_phy_p('zdtfi_omp',zdtfi_omp(:,:),llm)
    663599     
    664600      if (planet_type=="earth") then
     
    681617     .             ztfi_omp,
    682618     .             zqfi_omp,
    683 c     .             pvervel_omp,
    684619c#ifdef INCA
    685620     .             flxwfi_omp,
     
    695630#endif
    696631      endif !of if (planet_type=="earth")
    697 
    698 cym      call WriteField_phy_p('zdtfi_omp',zdtfi_omp(:,:),llm)
    699 
    700632c$OMP BARRIER
    701633
     
    753685        enddo
    754686      enddo
    755        
    756 c      do l=1,llm
    757 c        do i=1,klon
    758 c         pvervel(offset+i,l)=pvervel_omp(i,l)
    759 c       enddo
    760 c      enddo
    761687       
    762688      do l=1,llm
     
    791717     
    792718
    793 cc$OMP END PARALLEL
    794719      klon=klon_mpi
    795720500   CONTINUE
     
    797722
    798723c$OMP MASTER
    799 cym      call WriteField_phy('zdtfi',zdtfi(:,:),llm)
    800724      call stop_timer(timer_physic)
    801725c$OMP END MASTER
     
    913837      DO l=1,llm
    914838
    915 !!cdir NODEP
     839!CDIR ON_ADB(index_i)
     840!CDIR ON_ADB(index_j)
     841!cdir NODEP
    916842        do ig0=kstart,kend
    917843          i=index_i(ig0)
     
    975901
    976902C
    977 
     903!cdir NODEP
    978904      DO iq=1,nqtot
    979905         iiq=niadv(iq)
    980906c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    981907         DO l=1,llm
    982 
    983 !!cdir NODEP           
     908!CDIR ON_ADB(index_i)
     909!CDIR ON_ADB(index_j)
     910!cdir NODEP           
    984911             DO ig0=kstart,kend
    985912              i=index_i(ig0)
     
    1009936c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    1010937      DO l=1,llm
    1011 !!cdir NODEP
     938!CDIR ON_ADB(index_i)
     939!CDIR ON_ADB(index_j)
     940!cdir NODEP
    1012941         do ig0=kstart,kend
    1013942           i=index_i(ig0)
     
    1052981c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)     
    1053982      DO l=1,llm
    1054 !!cdir NODEP
     983!CDIR ON_ADB(index_i)
     984!CDIR ON_ADB(index_j)
     985!cdir NODEP
    1055986        do ig0=kstart,kend
    1056987           i=index_i(ig0)
  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/cray.F

    r774 r1250  
    1313      real sx((n-1)*incx+1),sy((n-1)*incy+1)
    1414c
     15      if (incx.eq.1.and.incy.eq.1) then
     16      do 10 i=1,n
     17         sy(i)=sx(i)
     1810    continue
     19      else
    1520      iy=1
    1621      ix=1
    17       do 10 i=1,n
     22      do 11 i=1,n
    1823         sy(iy)=sx(ix)
    1924         ix=ix+incx
    2025         iy=iy+incy
    21 10    continue
     2611    continue
     27      endif
    2228c
    2329      return
     
    3238c
    3339      ssum=0.
     40      if (incx.eq.1) then
     41      do 10 i=1,n
     42         ssum=ssum+sx(i)
     4310    continue
     44      else
    3445      ix=1
    35       do 10 i=1,n
     46      do 11 i=1,n
    3647         ssum=ssum+sx(ix)
    3748         ix=ix+incx
    38 10    continue
     4911    continue
     50      endif
    3951c
    4052      return
  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/groupeun_p.F

    r1087 r1250  
    11      SUBROUTINE groupeun_p(jjmax,llmax,jjb,jje,q)
    22      USE parallel
     3      USE Write_Field_p
    34      IMPLICIT NONE
    45
     
    1718      REAL airecs,qs
    1819
    19       INTEGER i,j,l,ig,j1,j2,i0,jd
     20      INTEGER i,j,l,ig,ig2,j1,j2,i0,jd
    2021
    2122c--------------------------------------------------------------------c
     
    3738      LOGICAL, SAVE :: first = .TRUE.
    3839!$OMP THREADPRIVATE(first)
     40      INTEGER,SAVE :: i_index(iim,ngroup)
     41      INTEGER      :: offset
     42      REAL         :: qsum(iim/ngroup)
    3943
    4044      IF (first) THEN
     
    5458            j_start  = MAX(jjb, j1-jd)
    5559            j_finish = MIN(jje, j2-jd)
    56             DO j=j_start, j_finish
    57                DO i0=1,iim,2**(ngroup-ig+1)
    58                   qn=0.
    59                   DO i=i0,i0+2**(ngroup-ig+1)-1
    60                      qn=qn+q(i,j,l)
    61                   ENDDO
    62                   DO i=i0,i0+2**(ngroup-ig+1)-1
    63                      q(i,j,l)=qn*airen_tab(i,j,jd)
    64                   ENDDO
     60            DO ig2=1,ngroup-ig+1
     61              offset=2**(ig2-1)
     62              DO j=j_start, j_finish
     63!CDIR NODEP
     64!CDIR ON_ADB(q)
     65                 DO i0=1,iim,2**ig2
     66                   q(i0,j,l)=q(i0,j,l)+q(i0+offset,j,l)
     67                 ENDDO
     68              ENDDO
     69            ENDDO
     70           
     71            DO j=j_start, j_finish
     72!CDIR NODEP
     73!CDIR ON_ADB(q)
     74               DO i=1,iim
     75                 q(i,j,l)=q(i-MOD(i-1,2**(ngroup-ig+1)),j,l)
     76               ENDDO
     77            ENDDO
     78
     79            DO j=j_start, j_finish
     80!CDIR ON_ADB(airen_tab)
     81!CDIR ON_ADB(q)
     82               DO i=1,iim
     83                 q(i,j,l)=q(i,j,l)*airen_tab(i,j,jd)
    6584               ENDDO
    6685               q(iip1,j,l)=q(1,j,l)
    6786            ENDDO
    68         
     87       
    6988!c     Concerne le pole sud
    7089            j_start  = MAX(1+jjp1-jje-jd, j1-jd)
    7190            j_finish = MIN(1+jjp1-jjb-jd, j2-jd)
    72             DO j=j_start, j_finish
    73                DO i0=1,iim,2**(ngroup-ig+1)
    74                   qs=0.
    75                   DO i=i0,i0+2**(ngroup-ig+1)-1
    76                      qs=qs+q(i,jjp1-j+1-jd,l)
    77                   ENDDO
    78                   DO i=i0,i0+2**(ngroup-ig+1)-1
    79                      q(i,jjp1-j+1-jd,l)=qs*aires_tab(i,jjp1-j+1,jd)
    80                   ENDDO
     91            DO ig2=1,ngroup-ig+1
     92              offset=2**(ig2-1)
     93              DO j=j_start, j_finish
     94!CDIR NODEP
     95!CDIR ON_ADB(q)
     96                 DO i0=1,iim,2**ig2
     97                   q(i0,jjp1-j+1-jd,l)= q(i0,jjp1-j+1-jd,l)
     98     &                                 +q(i0+offset,jjp1-j+1-jd,l)
     99                 ENDDO
     100              ENDDO
     101            ENDDO
     102
     103
     104            DO j=j_start, j_finish
     105!CDIR NODEP
     106!CDIR ON_ADB(q)
     107               DO i=1,iim
     108                 q(i,jjp1-j+1-jd,l)=q(i-MOD(i-1,2**(ngroup-ig+1)),
     109     &                                jjp1-j+1-jd,l)
     110               ENDDO
     111            ENDDO
     112
     113            DO j=j_start, j_finish
     114!CDIR ON_ADB(aires_tab)
     115!CDIR ON_ADB(q)
     116               DO i=1,iim
     117                 q(i,jjp1-j+1-jd,l)=q(i,jjp1-j+1-jd,l)* 
     118     &                              aires_tab(i,jjp1-j+1,jd)
    81119               ENDDO
    82120               q(iip1,jjp1-j+1-jd,l)=q(1,jjp1-j+1-jd,l)
    83121            ENDDO
     122
    84123       
    85124            j1=j2+1
  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/leapfrog_p.F

    r1247 r1250  
    327327c$OMP MASTER
    328328      ItCount=ItCount+1
    329       if (MOD(ItCount,1)==1) then
     329      if (MOD(ItCount,12)==0) then
    330330        debug=.true.
    331331      else
  • LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_fft_wrapper.F90

    r986 r1250  
    1919    INTEGER,INTENT(IN)  :: nb_vect
    2020    REAL,INTENT(IN)     :: vect(vsize+inc,nb_vect)
    21     COMPLEX,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect)
     21    COMPLEX*16,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect)
    2222   
    2323    STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
     
    2929    INTEGER,INTENT(IN)  :: nb_vect
    3030    REAL,INTENT(INOUT)    :: vect(vsize+inc,nb_vect)
    31     COMPLEX,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)
     31    COMPLEX*16,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)
    3232 
    3333    STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique"
  • LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_filtre_fft.F90

    r1249 r1250  
    199199
    200200    REAL               :: vect(iim+inc,jj_end-jj_begin+1,nbniv)
    201     COMPLEX            :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv)
     201    COMPLEX*16         :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv)
    202202    INTEGER            :: nb_vect
    203203    INTEGER :: i,j,l
     
    260260    REAL,INTENT(INOUT) :: vect_inout(iim+1,nlat,nbniv)
    261261
    262      REAL               :: vect(iim+inc,jj_end-jj_begin+1,nbniv)
    263     COMPLEX            :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv)
     262    REAL               :: vect(iim+inc,jj_end-jj_begin+1,nbniv)
     263    COMPLEX*16         :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv)
    264264    INTEGER            :: nb_vect
    265265    INTEGER :: i,j,l
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90

    r1249 r1250  
    5151  REAL, DIMENSION(klon,klev,naero_tot,nbands) ::  cg_ae
    5252  LOGICAL ::  soluble
    53   INTEGER :: i, k, ierr, inu, m, mrfspecies
     53  INTEGER :: i, k,n, ierr, inu, m, mrfspecies
    5454  INTEGER :: spsol, spinsol, spss
    55   INTEGER :: RH_num
     55  INTEGER :: RH_num(klon,klev)
    5656  INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA
    5757
     
    122122  REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./)
    123123  REAL, PARAMETER :: RH_MAX=95.
    124   REAL:: DELTA, rh, H
    125   REAL:: tau_ae2b_int(KLON,KLEV,nbands)   ! Intermediate computation of epaisseur optique aerosol
    126   REAL:: piz_ae2b_int(KLON,KLEV,nbands)   ! Intermediate computation of Single scattering albedo
    127   REAL:: cg_ae2b_int(KLON,KLEV,nbands)    ! Intermediate computation of Assymetry parameter
     124  REAL:: DELTA(klon,klev), rh(klon,klev), H
     125  REAL:: tau_ae2b_int   ! Intermediate computation of epaisseur optique aerosol
     126  REAL:: piz_ae2b_int   ! Intermediate computation of Single scattering albedo
     127  REAL:: cg_ae2b_int    ! Intermediate computation of Assymetry parameter
     128  REAL :: Fact_RH(nbre_RH)
    128129  REAL :: zrho
    129130  REAL :: fac
     
    146147  REAL:: piz_aeri_2bands(nbands,naero_insoluble)           !-- unit
    147148
     149  INTEGER :: id
     150  LOGICAL :: used_aer(naero_tot)
     151  REAL :: tmp_var, tmp_var_pi
    148152
    149153  DATA presnivs_19/&
     
    586590
    587591  DO k=1, klev
    588      DO i=1, klon
    589 !        IF (t_seri(i,k).EQ.0.) THEN
    590 !           WRITE(lunout,*) 't_seri(i,k)=0 for i=',i,'k=',k
    591 !           CALL abort_gcm('aeropt_2bands','t_seri=0',1)
    592 !        END IF
    593 !        IF (pplay(i,k).EQ.0.) THEN
    594 !           WRITE(lunout,*) 'pplay(i,k)=0 for i=',i,'k=',k
    595 !           CALL abort_gcm('aeropt_2bands','pplay=0',1)
    596 !        END IF
    597         zrho=pplay(i,k)/t_seri(i,k)/RD                    ! kg/m3
    598         mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9
     592    DO i=1, klon
     593      zrho=pplay(i,k)/t_seri(i,k)/RD                  ! kg/m3
     594!CDIR UNROLL=naero_spc
     595      mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9
    599596!RAF zrho
    600         mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9
    601      ENDDO
     597!CDIR UNROLL=naero_spc
     598      mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9
     599      zdp1(i,k)=pdel(i,k)/(gravit*delt)      ! air mass auxiliary  variable --> zdp1 [kg/(m^2 *s)]
     600    ENDDO
    602601  ENDDO
    603602
    604603  IF (flag_aerosol .EQ. 1) THEN
    605      nb_aer = 1
     604     nb_aer = 2
    606605     ALLOCATE (aerosol_name(nb_aer))
    607606     aerosol_name(1) = id_ASSO4M
    608 
     607     aerosol_name(2) = id_CSSO4M
    609608  ELSEIF (flag_aerosol .EQ. 2) THEN
    610609     nb_aer = 2
     
    650649  !      compute optical_thickness_at_gridpoint_per_species
    651650
    652   tau_ae(:,:,:,:)=0.
    653 !RAF
    654   tau_ae_pi(:,:,:,:)=0.
    655   piz_ae(:,:,:,:)=0.
    656   cg_ae(:,:,:,:)=0.
    657   tau_allaer(:,:,:,:)=0.
    658   piz_allaer(:,:,:,:)=0.
    659   cg_allaer(:,:,:,:)=0.
    660 
    661   !
    662   ! Calculations that need to be done since we are not in the subroutines INCA
    663   !     
    664   ! air mass auxiliary  variable --> zdp1 [kg/(m^2 *s)]
    665   zdp1(:,:)=pdel(:,:)/(gravit*delt)
    666 
    667 
     651
     652
     653!!CDIR ON_ADB(RH_tab)
     654!CDIR ON_ADB(fact_RH)
     655!CDIR SHORTLOOP
     656  DO n=1,nbre_RH-1
     657    fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n))
     658  ENDDO
     659   
     660  DO k=1, KLEV
     661!!CDIR ON_ADB(RH_tab)
     662!CDIR ON_ADB(fact_RH)
     663    DO i=1, KLON
     664      rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
     665      RH_num(i,k) = INT( rh(i,k)/10. + 1.)
     666      IF (rh(i,k).GT.85.) RH_num(i,k)=10
     667      IF (rh(i,k).GT.90.) RH_num(i,k)=11
     668     
     669      DELTA(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k))
     670    ENDDO
     671  ENDDO
     672
     673  used_aer(:)=.FALSE.
     674   
    668675  DO m=1,nb_aer   ! tau is only computed for each mass
    669 
    670      fac=1.0
     676    fac=1.0
    671677     IF (aerosol_name(m).EQ.id_ASBCM) THEN
    672678         soluble=.TRUE.
     
    715721     ENDIF
    716722
    717 
    718      tau_ae2b_int(:,:,:)=0.
    719      piz_ae2b_int(:,:,:)=0.
    720      cg_ae2b_int(:,:,:)=0.
    721 
    722      DO inu=1,nbands
     723    id=aerosol_name(m)
     724    used_aer(id)=.TRUE.
     725
     726     
     727    IF (soluble) THEN
     728
     729      IF (spss.NE.0) THEN
     730
     731         IF (spss.EQ.1) THEN !accumulation mode
     732            DO k=1, KLEV
     733!CDIR ON_ADB(A1_ASSSM_b1)
     734!CDIR ON_ADB(A2_ASSSM_b1)
     735!CDIR ON_ADB(A3_ASSSM_b1)
     736!CDIR ON_ADB(B1_ASSSM_b1)
     737!CDIR ON_ADB(B2_ASSSM_b1)
     738!CDIR ON_ADB(C1_ASSSM_b1)
     739!CDIR ON_ADB(C2_ASSSM_b2)
     740!CDIR ON_ADB(A1_ASSSM_b2)
     741!CDIR ON_ADB(A2_ASSSM_b2)
     742!CDIR ON_ADB(A3_ASSSM_b2)
     743!CDIR ON_ADB(B1_ASSSM_b2)
     744!CDIR ON_ADB(B2_ASSSM_b2)
     745!CDIR ON_ADB(C1_ASSSM_b2)
     746!CDIR ON_ADB(C2_ASSSM_b2)
     747              DO i=1, KLON
     748                H=rh(i,k)/100
     749                tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     750                tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     751
     752                ! band 1
     753                tau_ae2b_int=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05)
     754                piz_ae2b_int=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H
     755                cg_ae2b_int=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H
     756
     757                tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int
     758                tau_ae_pi(i,k,id,1) =  tmp_var_pi* tau_ae2b_int
     759                piz_ae(i,k,id,1) = piz_ae2b_int
     760                cg_ae(i,k,id,1)= cg_ae2b_int
     761               
     762                !band 2
     763                tau_ae2b_int=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05)
     764                piz_ae2b_int=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H
     765                cg_ae2b_int=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H
     766
     767                tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int
     768                tau_ae_pi(i,k,id,2) =  tmp_var_pi* tau_ae2b_int
     769                piz_ae(i,k,id,2) = piz_ae2b_int
     770                cg_ae(i,k,id,2)= cg_ae2b_int
     771
     772              ENDDO
     773            ENDDO
     774          ENDIF
     775
     776          IF (spss.EQ.2) THEN !coarse mode
     777            DO k=1, KLEV
     778!CDIR ON_ADB(A1_CSSSM_b1)
     779!CDIR ON_ADB(A2_CSSSM_b1)
     780!CDIR ON_ADB(A3_CSSSM_b1)
     781!CDIR ON_ADB(B1_CSSSM_b1)
     782!CDIR ON_ADB(B2_CSSSM_b1)
     783!CDIR ON_ADB(C1_CSSSM_b1)
     784!CDIR ON_ADB(C2_CSSSM_b2)
     785!CDIR ON_ADB(A1_CSSSM_b2)
     786!CDIR ON_ADB(A2_CSSSM_b2)
     787!CDIR ON_ADB(A3_CSSSM_b2)
     788!CDIR ON_ADB(B1_CSSSM_b2)
     789!CDIR ON_ADB(B2_CSSSM_b2)
     790!CDIR ON_ADB(C1_CSSSM_b2)
     791!CDIR ON_ADB(C2_CSSSM_b2)
     792              DO i=1, KLON
     793                H=rh(i,k)/100
     794                tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     795                tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     796                ! band 1
     797                tau_ae2b_int=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05)
     798                piz_ae2b_int=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H
     799                cg_ae2b_int=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H
     800
     801                tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int
     802                tau_ae_pi(i,k,id,1) =  tmp_var_pi* tau_ae2b_int
     803                piz_ae(i,k,id,1) = piz_ae2b_int
     804                cg_ae(i,k,id,1)= cg_ae2b_int
     805
     806                ! band 2
     807                tau_ae2b_int=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05)
     808                piz_ae2b_int=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H
     809                cg_ae2b_int=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H
     810
     811                tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int
     812                tau_ae_pi(i,k,id,2) =  tmp_var_pi* tau_ae2b_int
     813                piz_ae(i,k,id,2) = piz_ae2b_int
     814                cg_ae(i,k,id,2)= cg_ae2b_int
     815
     816             ENDDO
     817           ENDDO
     818         ENDIF
     819
     820         IF (spss.EQ.3) THEN !super coarse mode
     821            DO k=1, KLEV
     822!CDIR ON_ADB(A1_SSSSM_b1)
     823!CDIR ON_ADB(A2_SSSSM_b1)
     824!CDIR ON_ADB(A3_SSSSM_b1)
     825!CDIR ON_ADB(B1_SSSSM_b1)
     826!CDIR ON_ADB(B2_SSSSM_b1)
     827!CDIR ON_ADB(C1_SSSSM_b1)
     828!CDIR ON_ADB(C2_SSSSM_b2)
     829!CDIR ON_ADB(A1_SSSSM_b2)
     830!CDIR ON_ADB(A2_SSSSM_b2)
     831!CDIR ON_ADB(A3_SSSSM_b2)
     832!CDIR ON_ADB(B1_SSSSM_b2)
     833!CDIR ON_ADB(B2_SSSSM_b2)
     834!CDIR ON_ADB(C1_SSSSM_b2)
     835!CDIR ON_ADB(C2_SSSSM_b2)
     836              DO i=1, KLON
     837                H=rh(i,k)/100
     838                tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     839                tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     840
     841                ! band 1
     842                tau_ae2b_int=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05)
     843                piz_ae2b_int=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H
     844                cg_ae2b_int=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H
     845
     846                tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int
     847                tau_ae_pi(i,k,id,1) =  tmp_var_pi* tau_ae2b_int
     848                piz_ae(i,k,id,1) = piz_ae2b_int
     849                cg_ae(i,k,id,1)= cg_ae2b_int
     850
     851                ! band 2
     852                tau_ae2b_int=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05)
     853                piz_ae2b_int=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H
     854                cg_ae2b_int=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H
     855
     856                tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int
     857                tau_ae_pi(i,k,id,2) =  tmp_var_pi* tau_ae2b_int
     858                piz_ae(i,k,id,2) = piz_ae2b_int
     859                cg_ae(i,k,id,2)= cg_ae2b_int
     860
     861              ENDDO
     862            ENDDO
     863          ENDIF
     864
     865        ELSE
     866                       
     867!CDIR ON_ADB(alpha_aers_2bands)
     868!CDIR ON_ADB(piz_aers_2bands)
     869!CDIR ON_ADB(cg_aers_2bands)
     870          DO k=1, KLEV
     871            DO i=1, KLON
     872              tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     873              tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac
     874!CDIR UNROLL=nbands
     875              DO inu=1,nbands
     876
     877                tau_ae2b_int= alpha_aers_2bands(RH_num(i,k),inu,spsol)+ &
     878                              DELTA(i,k)* (alpha_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
     879                              alpha_aers_2bands(RH_num(i,k),inu,spsol))
     880                     
     881                piz_ae2b_int = piz_aers_2bands(RH_num(i,k),inu,spsol) + &
     882                               DELTA(i,k)* (piz_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
     883                               piz_aers_2bands(RH_num(i,k),inu,spsol))
     884                     
     885                cg_ae2b_int = cg_aers_2bands(RH_num(i,k),inu,spsol) + &
     886                              DELTA(i,k)* (cg_aers_2bands(RH_num(i,k)+1,inu,spsol) - &
     887                              cg_aers_2bands(RH_num(i,k),inu,spsol))
     888
     889                tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
     890                tau_ae_pi(i,k,id,inu) =  tmp_var_pi* tau_ae2b_int
     891                piz_ae(i,k,id,inu) = piz_ae2b_int
     892                cg_ae(i,k,id,inu)= cg_ae2b_int
     893                         
     894              ENDDO
     895            ENDDO
     896          ENDDO
     897       
     898        ENDIF                     
     899
     900      ELSE                                                    ! For all aerosol insoluble components
     901
     902!CDIR ON_ADB(alpha_aers_2bands)
     903!CDIR ON_ADB(piz_aers_2bands)
     904!CDIR ON_ADB(cg_aers_2bands)
    723905        DO k=1, KLEV
    724            DO i=1, KLON
    725 
    726               rh=MIN(RHcl(i,k)*100.,RH_MAX)
    727               RH_num = INT( rh/10. + 1.)
    728 
    729               IF (rh.GT.85.) RH_num=10
    730               IF (rh.GT.90.) RH_num=11
    731               DELTA=(rh-RH_tab(RH_num))/(RH_tab(RH_num+1)-RH_tab(RH_num))
    732 
    733 !           DO inu=1,nbands
    734               IF (soluble) THEN
    735 
    736               ! First optical parameters are computed for seasalt
    737                   IF (spss.NE.0) THEN
    738                       H=rh/100
    739                       IF (spss.EQ.1) THEN !accumulation mode
    740                           ! band 1
    741                           tau_ae2b_int(i,k,1)=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05)
    742                           piz_ae2b_int(i,k,1)=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H
    743                           cg_ae2b_int(i,k,1)=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H
    744                           !band 2
    745                           tau_ae2b_int(i,k,2)=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05)
    746                           piz_ae2b_int(i,k,2)=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H
    747                           cg_ae2b_int(i,k,2)=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H
    748                       ENDIF
    749                       IF (spss.EQ.2) THEN !coarse mode
    750                           ! band 1
    751                           tau_ae2b_int(i,k,1)=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05)
    752                           piz_ae2b_int(i,k,1)=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H
    753                           cg_ae2b_int(i,k,1)=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H
    754                           ! band 2
    755                           tau_ae2b_int(i,k,2)=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05)
    756                           piz_ae2b_int(i,k,2)=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H
    757                           cg_ae2b_int(i,k,2)=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H
    758                       ENDIF
    759                       IF (spss.EQ.3) THEN !super coarse mode
    760                           ! band 1
    761                           tau_ae2b_int(i,k,1)=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05)
    762                           piz_ae2b_int(i,k,1)=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H
    763                           cg_ae2b_int(i,k,1)=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H
    764                           ! band 2
    765                           tau_ae2b_int(i,k,2)=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05)
    766                           piz_ae2b_int(i,k,2)=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H
    767                           cg_ae2b_int(i,k,2)=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H
    768                       ENDIF
    769                   ELSE
    770                       tau_ae2b_int(i,k,inu)= &
    771                          alpha_aers_2bands(RH_num,inu,spsol)+ &
    772                          DELTA* (alpha_aers_2bands(RH_num+1,inu,spsol) - &
    773                          alpha_aers_2bands(RH_num,inu,spsol))
    774                      
    775                       piz_ae2b_int(i,k,inu)= &
    776                          piz_aers_2bands(RH_num,inu,spsol) + &
    777                          DELTA* (piz_aers_2bands(RH_num+1,inu,spsol) - &
    778                          piz_aers_2bands(RH_num,inu,spsol))
    779                      
    780                       cg_ae2b_int(i,k,inu)= &
    781                          cg_aers_2bands(RH_num,inu,spsol) + &
    782                          DELTA* (cg_aers_2bands(RH_num+1,inu,spsol) - &
    783                          cg_aers_2bands(RH_num,inu,spsol))
    784                   ENDIF
    785 
    786                  tau_ae(i,k,aerosol_name(m),inu) = &
    787                       mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt* &
    788                                          tau_ae2b_int(i,k,inu)*fac
    789 !RAF fac
    790                   tau_ae_pi(i,k,aerosol_name(m),inu) = &
    791                      mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt* &
    792                                            tau_ae2b_int(i,k,inu)*fac
    793 
    794               ELSE                                                    ! For all aerosol insoluble components
    795                  tau_ae2b_int(i,k,inu) = alpha_aeri_2bands(inu,spinsol)
    796                  piz_ae2b_int(i,k,inu) = piz_aeri_2bands(inu,spinsol)
    797                  cg_ae2b_int(i,k,inu) = cg_aeri_2bands(inu,spinsol)
    798 
    799                  tau_ae(i,k,aerosol_name(m),inu) = &
    800                       mass_temp(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)* &
    801                                                delt*tau_ae2b_int(i,k,inu)*fac
    802 !RAF fac
    803                  tau_ae_pi(i,k,aerosol_name(m),inu) = &
    804                     mass_temp_pi(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)* &
    805                                                  delt*tau_ae2b_int(i,k,inu)*fac
    806               ENDIF
    807 
    808               piz_ae(i,k,aerosol_name(m),inu) = piz_ae2b_int(i,k,inu)
    809 
    810               cg_ae(i,k,aerosol_name(m),inu)= cg_ae2b_int(i,k,inu)
    811 
    812            ENDDO    ! nbands : boucle sur les bandes spectrale
    813         ENDDO     ! klon : Boucle sur les points geographiques (grille horizontale)
    814      ENDDO     ! klev : Boucle sur les niveaux verticaux
    815   ENDDO     ! nb_aer : Boucle  sur les masses de traceurs
    816 
     906          DO i=1, KLON
     907            tmp_var=mass_temp(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)*delt*fac
     908            tmp_var_pi=mass_temp_pi(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)*delt*fac
     909!CDIR UNROLL=nbands
     910            DO inu=1,nbands
     911              tau_ae2b_int = alpha_aeri_2bands(inu,spinsol)
     912              piz_ae2b_int = piz_aeri_2bands(inu,spinsol)
     913              cg_ae2b_int = cg_aeri_2bands(inu,spinsol)
     914
     915              tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int
     916              tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int
     917              piz_ae(i,k,id,inu) = piz_ae2b_int
     918              cg_ae(i,k,id,inu)= cg_ae2b_int
     919            ENDDO
     920          ENDDO
     921        ENDDO
     922
     923      ENDIF ! soluble
     924
     925    ENDDO  ! nb_aer 
     926
     927  DO m=1,nb_aer   
     928    IF (.NOT. used_aer(m)) THEN
     929      tau_ae(:,:,:,:)=0.
     930      tau_ae_pi(:,:,:,:)=0.
     931      piz_ae(:,:,:,:)=0.
     932      cg_ae(:,:,:,:)=0.
     933    ENDIF
     934  ENDDO
    817935
    818936  DO inu=1, nbands
    819      DO mrfspecies=1,naero_grp
     937    DO mrfspecies=1,naero_grp
     938      IF (mrfspecies .EQ. 2) THEN             ! = total aerosol AER     
    820939        DO k=1, KLEV
    821            DO i=1, KLON
    822               IF (mrfspecies .EQ. 2) THEN             ! = total aerosol AER     
    823                  tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)+ &
    824                       tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+ &                                                   
    825                       tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ &   
    826                       tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)+ &
    827                       tau_ae(i,k,id_CIDUSTM,inu)
    828                  tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
     940          DO i=1, KLON
     941            tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)+ &
     942                                           tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+   &                                                   
     943                                           tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ &       
     944                                           tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+   &
     945                                           tau_ae(i,k,id_SSSSM,inu)+ tau_ae(i,k,id_CIDUSTM,inu)
     946             tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
    829947                 
    830                  piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
    831                       tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &
    832                       tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ &
    833                       tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ &
    834                       tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &
    835                       tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &   
    836                       tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
    837                       tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
    838                       tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ &
    839                       tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu))/tau_allaer(i,k,mrfspecies,inu)
    840                  piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
    841 
    842                  cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ &
     948             piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ &
     949                                             tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ &
     950                                             tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ &
     951                                             tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ &
     952                                             tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ &
     953                                             tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ &     
     954                                             tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ &
     955                                             tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ &
     956                                             tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ &
     957                                             tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) &
     958                                            /tau_allaer(i,k,mrfspecies,inu)
     959             piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
     960
     961             cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ &
    843962                      tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ &
    844963                      tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ &
     
    851970                      tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ &
    852971                      (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
    853 
    854               ELSEIF (mrfspecies .EQ. 3) THEN             ! = natural aerosol NAT
     972          ENDDO   
     973        ENDDO
     974
     975      ELSEIF (mrfspecies .EQ. 3) THEN             ! = natural aerosol NAT
     976
     977        DO k=1, KLEV
     978          DO i=1, KLON
    855979!RAF
    856980                 tau_allaer(i,k,mrfspecies,inu)=tau_ae_pi(i,k,id_ASSO4M,inu)+ &
     
    8921016                      cg_ae(i,k,id_CIDUSTM,inu))/ &
    8931017                      (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
    894 
    895 !END RAF
    896               ELSEIF (mrfspecies .EQ. 4) THEN             ! = BC
    897                  tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)
    898                  tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
    899                  piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) &
     1018          ENDDO
     1019        ENDDO
     1020                   
     1021      ELSEIF (mrfspecies .EQ. 4) THEN             ! = BC
     1022        DO k=1, KLEV
     1023          DO i=1, KLON
     1024            tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)
     1025            tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
     1026            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) &
    9001027                      +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ &
    9011028                      tau_allaer(i,k,mrfspecies,inu)
    902                  piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
    903                 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) *cg_ae(i,k,id_ASBCM,inu)&
     1029            piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
     1030            cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) *cg_ae(i,k,id_ASBCM,inu)&
    9041031                      +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu))/ &
    9051032                      (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
    906 
    907               ELSEIF (mrfspecies .EQ. 5) THEN             ! = SO4
    908                  tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)
    909                  tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
    910                  piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) &
     1033          ENDDO
     1034        ENDDO
     1035             
     1036      ELSEIF (mrfspecies .EQ. 5) THEN             ! = SO4
     1037
     1038        DO k=1, KLEV
     1039          DO i=1, KLON
     1040            tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)
     1041            tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
     1042            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) &
    9111043                      +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ &
    9121044                      tau_allaer(i,k,mrfspecies,inu)
    913                  piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
    914                 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) *cg_ae(i,k,id_CSSO4M,inu)&
     1045            piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
     1046            cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) *cg_ae(i,k,id_CSSO4M,inu)&
    9151047                      +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu))/ &
    9161048                      (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
    917 
    918               ELSEIF (mrfspecies .EQ. 6) THEN             ! = POM
    919                  tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)
    920                  tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
    921                  piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) &
     1049          ENDDO
     1050        ENDDO
     1051
     1052      ELSEIF (mrfspecies .EQ. 6) THEN             ! = POM
     1053
     1054        DO k=1, KLEV
     1055          DO i=1, KLON
     1056            tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)
     1057            tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
     1058            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) &
    9221059                      +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ &
    9231060                      tau_allaer(i,k,mrfspecies,inu)
    924                  piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
    925                 cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) *cg_ae(i,k,id_ASPOMM,inu)&
     1061            piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
     1062            cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) *cg_ae(i,k,id_ASPOMM,inu)&
    9261063                      +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu))/ &
    9271064                      (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
    928 
    929               ELSEIF (mrfspecies .EQ. 7) THEN             ! = DUST
    930                  tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu)
    931                  tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
    932                  piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu)
    933                  cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu)
    934 
    935               ELSEIF (mrfspecies .EQ. 8) THEN             ! = SS
    936                  tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)
    937                  tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
    938                  piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) &
    939                       +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) &
    940                       +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ &
    941                       tau_allaer(i,k,mrfspecies,inu)
    942                  piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
    943                  cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) *cg_ae(i,k,id_ASSSM,inu)&
    944                       +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) &
    945                       +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ &
    946                       (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
    947 
    948               ELSEIF (mrfspecies .EQ. 9) THEN             ! = NO3
    949                  tau_allaer(i,k,mrfspecies,inu)=0.   ! preliminary
    950                  piz_allaer(i,k,mrfspecies,inu)=0.
    951                  cg_allaer(i,k,mrfspecies,inu)=0.
    952               ENDIF
    953            ENDDO
     1065          ENDDO
    9541066        ENDDO
    955      ENDDO
     1067             
     1068      ELSEIF (mrfspecies .EQ. 7) THEN             ! = DUST
     1069
     1070        DO k=1, KLEV
     1071          DO i=1, KLON
     1072            tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu)
     1073            tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
     1074            piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu)
     1075            cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu)
     1076          ENDDO
     1077        ENDDO
     1078
     1079      ELSEIF (mrfspecies .EQ. 8) THEN             ! = SS
     1080
     1081        DO k=1, KLEV
     1082          DO i=1, KLON
     1083            tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)
     1084            tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20)
     1085            piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) &
     1086                    +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) &
     1087                    +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ &
     1088                    tau_allaer(i,k,mrfspecies,inu)
     1089            piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20)
     1090            cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) *cg_ae(i,k,id_ASSSM,inu)&
     1091                    +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) &
     1092                    +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ &
     1093                    (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu))
     1094          ENDDO
     1095        ENDDO
     1096     
     1097      ELSEIF (mrfspecies .EQ. 9) THEN             ! = NO3
     1098     
     1099        DO k=1, KLEV
     1100          DO i=1, KLON
     1101            tau_allaer(i,k,mrfspecies,inu)=0.   ! preliminary
     1102            piz_allaer(i,k,mrfspecies,inu)=0.
     1103            cg_allaer(i,k,mrfspecies,inu)=0.
     1104          ENDDO
     1105        ENDDO
     1106     
     1107      ELSE
     1108
     1109        DO k=1, KLEV
     1110          DO i=1, KLON
     1111            tau_allaer(i,k,mrfspecies,inu)=0. 
     1112            piz_allaer(i,k,mrfspecies,inu)=0.
     1113            cg_allaer(i,k,mrfspecies,inu)=0.
     1114          ENDDO
     1115        ENDDO
     1116           
     1117      ENDIF
     1118
     1119    ENDDO
    9561120  ENDDO
    9571121
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90

    r1249 r1250  
    126126
    127127  REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./)
    128   REAL :: DELTA, rh, H
    129   REAL :: tau_ae5wv_int(KLON,KLEV,las) ! Intermediate computation of epaisseur optique aerosol
    130   REAL :: piz_ae5wv_int(KLON,KLEV,las) ! Intermediate single scattering albedo aerosol
    131   REAL :: cg_ae5wv_int(KLON,KLEV,las)  ! Intermediate asymmetry parameter aerosol
     128  REAL :: DELTA(klon,klev), rh(klon,klev), H
     129  REAL :: tau_ae5wv_int ! Intermediate computation of epaisseur optique aerosol
     130  REAL :: piz_ae5wv_int ! Intermediate single scattering albedo aerosol
     131  REAL :: cg_ae5wv_int  ! Intermediate asymmetry parameter aerosol
    132132  REAL, PARAMETER :: RH_MAX=95.
    133133  REAL :: taue670(KLON)       ! epaisseur optique aerosol absorption 550 nm
     
    158158  ! Proprietes optiques
    159159  !
    160   REAL :: radry = 287.054                     ! dry air mass constant
    161 
     160  REAL :: radry = 287.054
     161  REAL :: tau_tmp                     ! dry air mass constant
     162  REAL :: fact_RH(nbre_RH)
     163  LOGICAL :: used_tau(naero_spc)
     164 
    162165  DATA presnivs_19/&
    163166       100426.5,  98327.6, 95346.5, 90966.8, 84776.9, &
     
    477480  ! Initialisations
    478481  ai(:) = 0.
    479   tau_ae5wv_int(:,:,:) = 0.
    480   piz_ae5wv_int(:,:,:) = 0.
    481   cg_ae5wv_int(:,:,:)  = 0.
    482482  tausum(:,:,:) = 0.
    483   tau(:,:,:,:)  = 0.
    484483
    485484
     
    489488!      IF (pplay(i,k).EQ.0) stop  'stop aeropt_5wv p '
    490489      zrho=pplay(i,k)/t_seri(i,k)/RD                  ! kg/m3
     490!CDIR UNROLL=naero_spc
    491491      mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9
     492      zdp1(i,k)=pdel(i,k)/(gravit*delt)     ! air mass auxiliary  variable --> zdp1 [kg/(m^2 *s)]
     493
    492494    ENDDO
    493495  ENDDO
     
    545547  ! Calculations that need to be done since we are not in the subroutines INCA
    546548  !     
    547   ! air mass auxiliary  variable --> zdp1 [kg/(m^2 *s)]
    548   zdp1=pdel(:,:)/(gravit*delt)
    549  
     549
     550!!CDIR ON_ADB(RH_tab)
     551!CDIR ON_ADB(fact_RH)
     552!CDIR NOVECTOR
     553  DO RH_Num=1,nbre_RH-1
     554    fact_RH(RH_num)=1./(RH_tab(RH_num+1)-RH_tab(RH_num))
     555  ENDDO
     556   
     557  DO k=1, KLEV
     558!!CDIR ON_ADB(RH_tab)
     559!CDIR ON_ADB(fact_RH)
     560    DO i=1, KLON
     561      rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX)
     562      RH_num = INT( rh(i,k)/10. + 1.)
     563      IF (rh(i,k).GT.85.) RH_num=10
     564      IF (rh(i,k).GT.90.) RH_num=11
     565      DELTA(i,k)=(rh(i,k)-RH_tab(RH_num))*fact_RH(RH_num)
     566    ENDDO
     567  ENDDO
     568
     569!CDIR SHORTLOOP 
     570  used_tau(:)=.FALSE.
     571   
    550572  DO m=1,nb_aer   ! tau is only computed for each mass   
    551573    fac=1.0
     
    595617        CYCLE
    596618    ENDIF
    597    
     619
     620    used_tau(spsol)=.TRUE.
    598621    DO la=1,las
    599       tau3d(:,:)=0.
    600       piz3d(:,:)=0.
    601       cg3d(:,:)=0.
    602       abs3d(:,:)=0.
    603      
    604       DO k=1, KLEV
    605         DO i=1, KLON
    606          
    607           rh=MIN(RHcl(i,k)*100.,RH_MAX)
    608           RH_num = INT( rh/10. + 1.)
    609          
    610           IF (rh.GT.85.) RH_num=10
    611           IF (rh.GT.90.) RH_num=11
    612           DELTA=(rh-RH_tab(RH_num))/(RH_tab(RH_num+1)-RH_tab(RH_num))
    613          
    614           IF (soluble) THEN
    615 
    616               IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm
    617                   H=rh/100
    618                   IF (spss.EQ.1) THEN !accumulation mode
    619                       tau_ae5wv_int(i,k,la)=A1_ASSSM(k)+A2_ASSSM(k)*H+A3_ASSSM(k)/(H-1.05)
    620                       piz_ae5wv_int(i,k,la)=1-B1_ASSSM(k)-B2_ASSSM(k)*H
    621                       cg_ae5wv_int(i,k,la)=C1_ASSSM(k)+C2_ASSSM(k)*H
    622                   ENDIF
    623                   IF (spss.EQ.2) THEN !coarse mode
    624                       tau_ae5wv_int(i,k,la)=A1_CSSSM(k)+A2_CSSSM(k)*H+A3_CSSSM(k)/(H-1.05)
    625                       piz_ae5wv_int(i,k,la)=1-B1_CSSSM(k)-B2_CSSSM(k)*H
    626                       cg_ae5wv_int(i,k,la)=C1_CSSSM(k)+C2_CSSSM(k)*H
    627                   ENDIF
    628                   IF (spss.EQ.3) THEN !super coarse mode
    629                       tau_ae5wv_int(i,k,la)=A1_SSSSM(k)+A2_SSSSM(k)*H+A3_SSSSM(k)/(H-1.05)
    630                       piz_ae5wv_int(i,k,la)=1-B1_SSSSM(k)-B2_SSSSM(k)*H
    631                       cg_ae5wv_int(i,k,la)=C1_SSSSM(k)+C2_SSSSM(k)*H
    632                   ENDIF
    633               ELSE
    634                  
    635                   tau_ae5wv_int(i,k,la) = &
    636                      alpha_aers_5wv(RH_num,la,spsol)+DELTA* &
    637                      (alpha_aers_5wv(RH_num+1,la,spsol) - &
    638                      alpha_aers_5wv(RH_num,la,spsol))
    639                  
    640                   piz_ae5wv_int(i,k,la) = &
    641                      piz_aers_5wv(RH_num,la,spsol)+DELTA* &
    642                      (piz_aers_5wv(RH_num+1,la,spsol) - &
    643                      piz_aers_5wv(RH_num,la,spsol))
    644                  
    645                   cg_ae5wv_int(i,k,la) = &
    646                      cg_aers_5wv(RH_num,la,spsol)+DELTA* &
    647                      (cg_aers_5wv(RH_num+1,la,spsol) - &
    648                      cg_aers_5wv(RH_num,la,spsol))
    649               ENDIF
    650 
    651               tau3d(i,k) = &
    652                  mass_temp(i,k,spsol)*1000.*zdp1(i,k)*tau_ae5wv_int(i,k,la)*delt*fac
    653 
    654           ELSE                                                  ! For insoluble aerosol
    655               tau_ae5wv_int(i,k,la) = alpha_aeri_5wv(la,spinsol)
    656               piz_ae5wv_int(i,k,la) = piz_aeri_5wv(la,spinsol)
    657               cg_ae5wv_int(i,k,la)  = cg_aeri_5wv(la,spinsol)
    658 
    659               tau3d(i,k) = &
    660                  mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* &
    661                       tau_ae5wv_int(i,k,la)*delt*fac
     622
     623      IF (soluble) THEN
     624
     625        IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm
     626          IF (spss.EQ.1) THEN !accumulation mode
     627            DO k=1, KLEV
     628!CDIR ON_ADB(A1_ASSSM)
     629!CDIR ON_ADB(A2_ASSSM)
     630!CDIR ON_ADB(A3_ASSSM)
     631              DO i=1, KLON
     632                H=rh(i,k)/100
     633                tau_ae5wv_int=A1_ASSSM(k)+A2_ASSSM(k)*H+A3_ASSSM(k)/(H-1.05)
     634                tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k)   &
     635                                   *tau_ae5wv_int*delt*fac
     636                tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol)
     637              ENDDO
     638            ENDDO
    662639          ENDIF
    663          
    664         ENDDO     ! Boucle sur les points géographiques (grille horizontale)
    665       ENDDO     ! Boucle sur les niveaux verticaux
    666 
    667       IF (soluble) THEN
    668 
    669           tau(:,:,la,spsol)=tau3d(:,:)
    670      
     640 
     641          IF (spss.EQ.2) THEN !coarse mode
     642            DO k=1, KLEV
     643!CDIR ON_ADB(A1_CSSSM)
     644!CDIR ON_ADB(A2_CSSSM)
     645!CDIR ON_ADB(A3_CSSSM)
     646              DO i=1, KLON
     647                H=rh(i,k)/100
     648                tau_ae5wv_int=A1_CSSSM(k)+A2_CSSSM(k)*H+A3_CSSSM(k)/(H-1.05)
     649                tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k)  &
     650                                   *tau_ae5wv_int*delt*fac
     651                tausum(i,la,spsol) = tausum(i,la,spsol)+tau(i,k,la,spsol)
     652              ENDDO
     653            ENDDO
     654          ENDIF
     655
     656          IF (spss.EQ.3) THEN !super coarse mode
     657            DO k=1, KLEV
     658!CDIR ON_ADB(A1_SSSSM)
     659!CDIR ON_ADB(A2_SSSSM)
     660!CDIR ON_ADB(A3_SSSSM)
     661              DO i=1, KLON
     662                H=rh(i,k)/100
     663                tau_ae5wv_int=A1_SSSSM(k)+A2_SSSSM(k)*H+A3_SSSSM(k)/(H-1.05)
     664                tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k)  &
     665                                   *tau_ae5wv_int*delt*fac
     666                tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol)
     667              ENDDO
     668            ENDDO
     669          ENDIF
     670
     671        ELSE
    671672          DO k=1, KLEV
    672             DO i=1,KLON
    673               tausum(i,la,spsol)=tausum(i,la,spsol)+tau3d(i,k)
     673!CDIR ON_ADB(alpha_aers_5wv)
     674            DO i=1, KLON
     675              tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* &
     676                             (alpha_aers_5wv(RH_num+1,la,spsol) - &
     677                              alpha_aers_5wv(RH_num,la,spsol))
     678
     679              tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k)   &
     680                                 *tau_ae5wv_int*delt*fac
     681              tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol)
    674682            ENDDO
    675683          ENDDO
    676       ELSE                                            ! For insoluble aerosol
    677           tau(:,:,la,naero_soluble+spinsol)=tau3d(:,:)
    678      
    679           DO k=1, KLEV
    680             DO i=1,KLON
    681               tausum(i,la,naero_soluble+spinsol)= &
    682                     tausum(i,la,naero_soluble+spinsol)+tau3d(i,k)
    683             ENDDO     ! Boucle sur les points géographiques (grille horizontale)
    684           ENDDO     ! Boucle sur les niveaux verticaux
     684        ENDIF
     685
     686      ELSE                                                  ! For insoluble aerosol
     687        DO k=1, KLEV
     688!CDIR ON_ADB(alpha_aeri_5wv)
     689          DO i=1, KLON
     690            tau_ae5wv_int = alpha_aeri_5wv(la,spinsol)
     691            tau(i,k,la,naero_soluble+spinsol) = mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* &
     692                                                tau_ae5wv_int*delt*fac
     693            tausum(i,la,naero_soluble+spinsol)= tausum(i,la,naero_soluble+spinsol)  &
     694                                               +tau(i,k,la,naero_soluble+spinsol)
     695          ENDDO
     696        ENDDO
    685697      ENDIF
    686 
    687698    ENDDO   ! boucle sur les longueurs d'onde
    688699  ENDDO     ! Boucle  sur les masses de traceurs
     700
     701  DO m=1,naero_spc
     702    IF (.NOT.used_tau(m)) tau(:,:,:,m)=0.
     703  ENDDO 
    689704!
    690705!
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/carbon_cycle_mod.F90

    r1249 r1250  
    99
    1010! Variables read from parmeter file physiq.def
    11   LOGICAL, SAVE, PUBLIC :: carbon_cycle_tr        ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys
     11  LOGICAL, PUBLIC :: carbon_cycle_tr        ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys
    1212!$OMP THREADPRIVATE(carbon_cycle_tr)
    13   LOGICAL, SAVE,PUBLIC :: carbon_cycle_cpl       ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES)
     13  LOGICAL, PUBLIC :: carbon_cycle_cpl       ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES)
    1414!$OMP THREADPRIVATE(carbon_cycle_cpl)
    1515  LOGICAL :: carbon_cycle_emis_comp=.FALSE. ! Calculation of emission compatible
     
    1717! Scalare values when no transport, from physiq.def
    1818  REAL :: fos_fuel_s  ! carbon_cycle_fos_fuel dans physiq.def
     19!$OMP THREADPRIVATE(fos_fuel_s)
    1920  REAL :: emis_land_s ! not yet implemented
     21!$OMP THREADPRIVATE(emis_land_s)
    2022
    2123  INTEGER :: ntr_co2                ! Number of tracers concerning the carbon cycle
     
    2527  INTEGER :: id_fco2_land_use       !  - " -
    2628  INTEGER :: id_fco2_fos_fuel       !  - " -
    27  
     29!$OMP THREADPRIVATE(ntr_co2, id_fco2_tot, id_fco2_ocn, id_fco2_land, id_fco2_land_use, id_fco2_fos_fuel) 
     30
    2831  REAL, DIMENSION(:), ALLOCATABLE :: fos_fuel        ! CO2 fossil fuel emission from file [gC/m2/d]
     32!$OMP THREADPRIVATE(fos_fuel)
    2933  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocn_day ! flux CO2 from ocean for 1 day (cumulated) [gC/m2/d]
     34!$OMP THREADPRIVATE(fco2_ocn_day)
    3035  REAL, DIMENSION(:), ALLOCATABLE :: fco2_land_day   ! flux CO2 from land for 1 day (cumulated)  [gC/m2/d]
     36!$OMP THREADPRIVATE(fco2_land_day)
    3137  REAL, DIMENSION(:), ALLOCATABLE :: fco2_lu_day     ! Emission from land use change for 1 day (cumulated) [gC/m2/d]
     38!$OMP THREADPRIVATE(fco2_lu_day)
    3239
    3340! Following 2 fields will be initialized in surf_land_orchidee at each time step
    3441  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_inst  ! flux CO2 from land at one time step
     42!$OMP THREADPRIVATE(fco2_land_inst)
    3543  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_lu_inst    ! Emission from land use change at one time step
     44!$OMP THREADPRIVATE(fco2_lu_inst)
    3645
    3746! Calculated co2 field to be send to the ocean via the coupler and to ORCHIDEE
    3847  REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: co2_send
    39 
     48!$OMP THREADPRIVATE(co2_send)
    4049
    4150CONTAINS
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/clouds_gno.F

    r1162 r1250  
    4747     
    4848      INTEGER i,K, n, m
    49       REAL mu(klon), qsat(klon), delta(klon), beta(klon)
    50       real zu2(klon),zv2(klon)
    51       REAL xx(klon), aux(klon), coeff(klon), block(klon)
    52       REAL  dist(klon), fprime(klon), det(klon)
    53       REAL pi, u(klon), v(klon), erfcu(klon), erfcv(klon)
    54       REAL  xx1(klon), xx2(klon)
    55       real erf,kkk
     49      REAL mu(klon), qsat, delta(klon), beta(klon)
     50      real zu2,zv2
     51      REAL xx(klon), aux(klon), coeff, block
     52      REAL  dist, fprime, det
     53      REAL pi, u, v, erfcu, erfcv
     54      REAL  xx1, xx2
     55      real erf,hsqrtlog_2,v2
    5656      real sqrtpi,sqrt2,zx1,zx2,exdel
    5757c lconv = true si le calcul a converge (entre autre si qsub < min_q)
    5858       LOGICAL lconv(klon)
    5959
    60 cym
    61       cldf(:,:)=0.0
     60!cdir arraycomb
     61      cldf  (1:klon,1:ND)=0.0        ! cym
     62      ratqsc(1:klon,1:ND)=0.0
     63      ptconv(1:klon,1:ND)=.false.
     64!cdir end arraycomb
    6265     
    6366      pi = ACOS(-1.)
    6467      sqrtpi=sqrt(pi)
    6568      sqrt2=sqrt(2.)
    66 
    67       ptconv=.false.
    68       ratqsc=0.
    69 
     69      hsqrtlog_2=0.5*SQRT(log(2.))
    7070
    7171      DO 500 K = 1, ND
     
    7474      mu(i) = R(i,K)
    7575      mu(i) = MAX(mu(i),min_mu)
    76       qsat(i) = RS(i,K)
    77       qsat(i) = MAX(qsat(i),min_mu)
    78       delta(i) = log(mu(i)/qsat(i))
    79                                     enddo ! vector
     76      qsat = RS(i,K)
     77      qsat = MAX(qsat,min_mu)
     78      delta(i) = log(mu(i)/qsat)
     79c                                   enddo ! vector
    8080
    8181C
     
    106106c  suffisamment d'eau nuageuse.
    107107
    108                                     do i=1,klon ! vector
     108c                                   do i=1,klon ! vector
    109109
    110110      IF ( QSUB(i,K) .lt. min_Q ) THEN
     
    124124c --  roots of equation v > vmax:
    125125
    126         det(i) = delta(i) + vmax(i)**2.
    127         if (det(i).LE.0.0) vmax(i) = vmax0 + 1.0
    128         det(i) = delta(i) + vmax(i)**2.
    129 
    130         if (det(i).LE.0.) then
     126        det = delta(i) + vmax(i)*vmax(i)
     127        if (det.LE.0.0) vmax(i) = vmax0 + 1.0
     128        det = delta(i) + vmax(i)*vmax(i)
     129
     130        if (det.LE.0.) then
    131131          xx(i) = -0.0001
    132132        else
    133133         zx1=-sqrt2*vmax(i)
    134          zx2=SQRT(1.0+delta(i)/(vmax(i)**2.))
    135          xx1(i)=zx1*(1.0-zx2)
    136          xx2(i)=zx1*(1.0+zx2)
    137          xx(i) = 1.01 * xx1(i)
    138          if ( xx1(i) .GE. 0.0 ) xx(i) = 0.5*xx2(i)
     134         zx2=SQRT(1.0+delta(i)/(vmax(i)*vmax(i)))
     135         xx1=zx1*(1.0-zx2)
     136         xx2=zx1*(1.0+zx2)
     137         xx(i) = 1.01 * xx1
     138         if ( xx1 .GE. 0.0 ) xx(i) = 0.5*xx2
    139139        endif
    140         if (delta(i).LT.0.) xx(i) = -0.5*SQRT(log(2.))
     140        if (delta(i).LT.0.) xx(i) = -hsqrtlog_2
    141141
    142142      ENDIF
     
    153153        if (.not.lconv(i)) then
    154154
    155           u(i) = delta(i)/(xx(i)*sqrt2) + xx(i)/(2.*sqrt2)
    156           v(i) = delta(i)/(xx(i)*sqrt2) - xx(i)/(2.*sqrt2)
    157 
    158           IF ( v(i) .GT. vmax(i) ) THEN
    159 
    160             IF (     ABS(u(i))  .GT. vmax(i)
     155          u = delta(i)/(xx(i)*sqrt2) + xx(i)/(2.*sqrt2)
     156          v = delta(i)/(xx(i)*sqrt2) - xx(i)/(2.*sqrt2)
     157          v2 = v*v
     158
     159          IF ( v .GT. vmax(i) ) THEN
     160
     161            IF (     ABS(u)  .GT. vmax(i)
    161162     :          .AND.  delta(i) .LT. 0. ) THEN
    162163
     
    171172             endif
    172173             xx(i) = -SQRT(aux(i))
    173              block(i) = EXP(-v(i)*v(i)) / v(i) / sqrtpi
    174              dist(i) = 0.0
    175              fprime(i) = 1.0
     174             block = EXP(-v*v) / v / sqrtpi
     175             dist = 0.0
     176             fprime = 1.0
    176177
    177178            ELSE
     
    179180c -- erfv -> 1.0, use an asymptotic expression of erfv for v large:
    180181
    181              erfcu(i) = 1.0-ERF(u(i))
     182             erfcu = 1.0-ERF(u)
    182183c  !!! ATTENTION : rajout d'un seuil pour l'exponentiel
    183              aux(i) = sqrtpi*erfcu(i)*EXP(min(v(i)*v(i),100.))
    184              coeff(i) = 1.0 - 1./2./(v(i)**2.) + 3./4./(v(i)**4.)
    185              block(i) = coeff(i) * EXP(-v(i)*v(i)) / v(i) / sqrtpi
    186              dist(i) = v(i) * aux(i) / coeff(i) - beta(i)
    187              fprime(i) = 2.0 / xx(i) * (v(i)**2.)
    188      :           * ( coeff(i)*EXP(-delta(i)) - u(i) * aux(i) )
    189      :           / coeff(i) / coeff(i)
     184             aux(i) = sqrtpi*erfcu*EXP(min(v2,100.))
     185             coeff = 1.0 - 0.5/(v2) + 0.75/(v2*v2)
     186             block = coeff * EXP(-v2) / v / sqrtpi
     187             dist = v * aux(i) / coeff - beta(i)
     188             fprime = 2.0 / xx(i) * (v2)
     189     :           * ( EXP(-delta(i)) - u * aux(i) / coeff )
     190     :           / coeff
    190191           
    191192            ENDIF ! ABS(u)
     
    195196c -- general case:
    196197
    197            erfcu(i) = 1.0-ERF(u(i))
    198            erfcv(i) = 1.0-ERF(v(i))
    199            block(i) = erfcv(i)
    200            dist(i) = erfcu(i) / erfcv(i) - beta(i)
    201            zu2(i)=u(i)*u(i)
    202            zv2(i)=v(i)*v(i)
    203            if(zu2(i).gt.20..or. zv2(i).gt.20.) then
     198           erfcu = 1.0-ERF(u)
     199           erfcv = 1.0-ERF(v)
     200           block = erfcv
     201           dist = erfcu / erfcv - beta(i)
     202           zu2=u*u
     203           zv2=v2
     204           if(zu2.gt.20..or. zv2.gt.20.) then
    204205c              print*,'ATTENTION !!! xx(',i,') =', xx(i)
    205206c           print*,'ATTENTION !!! klon,ND,R,RS,QSUB,PTCONV,RATQSC,CLDF',
     
    207208c     .CLDF(i,k)
    208209c              print*,'ATTENTION !!! zu2 zv2 =',zu2(i),zv2(i)
    209               zu2(i)=20.
    210               zv2(i)=20.
    211              fprime(i) = 0.
     210              zu2=20.
     211              zv2=20.
     212             fprime = 0.
    212213           else
    213              fprime(i) = 2. /sqrtpi /xx(i) /erfcv(i)**2.
    214      :           * (   erfcv(i)*v(i)*EXP(-zu2(i))
    215      :               - erfcu(i)*u(i)*EXP(-zv2(i)) )
     214             fprime = 2. /sqrtpi /xx(i) /(erfcv*erfcv)
     215     :           * (   erfcv*v*EXP(-zu2)
     216     :               - erfcu*u*EXP(-zv2) )
    216217           endif
    217218          ENDIF ! x
     
    223224!              stop
    224225!          endif
    225           if (abs(fprime(i)).lt.1.e-11) then
     226          if (abs(fprime).lt.1.e-11) then
    226227!              print*,'avant test fprime<.e-11 '
    227228!     s        ,i,k,lconv(i),u(i),v(i),beta(i),fprime(i)
    228229!              print*,'klon,ND,R,RS,QSUB',
    229230!     s        klon,ND,R(i,k),rs(i,k),qsub(i,k)
    230               fprime(i)=sign(1.e-11,fprime(i))
     231              fprime=sign(1.e-11,fprime)
    231232          endif
    232233
    233234
    234           if ( ABS(dist(i)/beta(i)) .LT. epsilon ) then
     235          if ( ABS(dist/beta(i)) .LT. epsilon ) then
    235236c           print*,'v-u **2',(v(i)-u(i))**2
    236237c           print*,'exp v-u **2',exp((v(i)-u(i))**2)
     
    238239            lconv(i)=.true.
    239240c  borne pour l'exponentielle
    240             ratqsc(i,k)=min(2.*(v(i)-u(i))**2,20.)
     241            ratqsc(i,k)=min(2.*(v-u)*(v-u),20.)
    241242            ratqsc(i,k)=sqrt(exp(ratqsc(i,k))-1.)
    242             CLDF(i,K) = 0.5 * block(i)
     243            CLDF(i,K) = 0.5 * block
    243244          else
    244             xx(i) = xx(i) - dist(i)/fprime(i)
     245            xx(i) = xx(i) - dist/fprime
    245246          endif
    246247c         print*,'apres test ',i,k,lconv(i)
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/cvltr.F90

    r1191 r1250  
    2828
    2929! Variables locales     
    30   REAL,DIMENSION(klon,klev)       :: zed
     30! REAL,DIMENSION(klon,klev)       :: zed
    3131  REAL,DIMENSION(klon,klev,klev)  :: zmd
    3232  REAL,DIMENSION(klon,klev,klev)  :: za
     
    3434  REAL,DIMENSION(klon,klev)       :: zmfp,zmfu
    3535  INTEGER                         :: i,k,j
     36  REAL                            :: pdtimeRG
    3637
    3738! =========================================
    3839! calcul des tendances liees au downdraft
    3940! =========================================
    40   zed(:,:)=0.
    41   zmfd(:,:)=0.
    42   zmfa(:,:)=0.
    43   zmfu(:,:)=0.
    44   zmfp(:,:)=0.
    45   zmd(:,:,:)=0.
    46   za(:,:,:)=0.
     41!cdir collapse
     42  DO j=1,klev
     43  DO i=1,klon
     44!   zed(i,j)=0.
     45    zmfd(i,j)=0.
     46    zmfa(i,j)=0.
     47    zmfu(i,j)=0.
     48    zmfp(i,j)=0.
     49  END DO
     50  END DO
     51!cdir collapse
     52  DO k=1,klev
     53  DO j=1,klev
     54  DO i=1,klon
     55    zmd(i,j,k)=0.
     56    za (i,j,k)=0.
     57  END DO
     58  END DO
     59  END DO
    4760! entrainement
    48   DO k=1,klev-1
    49      DO i=1,klon
    50         zed(i,k)=max(0.,mp(i,k)-mp(i,k+1))
    51      END DO
    52   END DO
     61! DO k=1,klev-1
     62!    DO i=1,klon
     63!       zed(i,k)=max(0.,mp(i,k)-mp(i,k+1))
     64!    END DO
     65! END DO
    5366
    5467! calcul de la matrice d echange
    5568! matrice de distribution de la masse entrainee en k
    5669
    57   DO k=1,klev
     70  DO k=1,klev-1
    5871     DO i=1,klon
    59         zmd(i,k,k)=zed(i,k)
     72        zmd(i,k,k)=max(0.,mp(i,k)-mp(i,k+1))
    6073     END DO
    6174  END DO
     
    126139  DO k=1, klev
    127140     DO i=1, klon
     141        dx(i,k)=paprs(i,k)-paprs(i,k+1)
     142     ENDDO
     143  ENDDO
     144  pdtimeRG=pdtime*RG
     145!cdir collapse
     146  DO k=1, klev
     147     DO i=1, klon
    128148        dx(i,k)=(zmfd(i,k)+zmfu(i,k)       &
    129              +zmfa(i,k)+zmfp(i,k))*pdtime  &
    130              *RG/(paprs(i,k)-paprs(i,k+1))
     149                +zmfa(i,k)+zmfp(i,k))*pdtimeRG/dx(i,k)
    131150        !          print*,'dx',k,dx(i,k)
    132151     ENDDO
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/fisrtilp.F

    r1091 r1250  
    5454      REAL frac_impa(klon,klev)
    5555      REAL frac_nucl(klon,klev)
    56       real zct(klon),zcl(klon)
     56      real zct      ,zcl
    5757cAA
    5858c
     
    8787      REAL ztglace, zt(klon)
    8888      INTEGER nexpo ! exponentiel pour glace/eau
    89       REAL zdz(klon),zrho(klon),ztot(klon), zrhol(klon)
    90       REAL zchau(klon),zfroi(klon),zfice(klon),zneb(klon)
     89      REAL zdz(klon),zrho(klon),ztot      , zrhol(klon)
     90      REAL zchau      ,zfroi      ,zfice(klon),zneb(klon)
    9191c
    9292      LOGICAL appel1er
     
    150150cAA Initialisation a 1 des coefs des fractions lessivees
    151151c
     152!cdir collapse
    152153      DO k = 1, klev
    153154       DO i = 1, klon
     
    161162c
    162163cMAf Initialisation a 0 de zoliq
    163        DO i = 1, klon
    164           zoliq(i)=0.
    165        ENDDO
     164c      DO i = 1, klon
     165c         zoliq(i)=0.
     166c      ENDDO
    166167c Determiner les nuages froids par leur temperature
    167168c  nexpo regle la raideur de la transition eau liquide / eau glace.
     
    173174c Initialiser les sorties:
    174175c
     176!cdir collapse
    175177      DO k = 1, klev+1
    176178      DO i = 1, klon
     
    180182      ENDDO
    181183
     184!cdir collapse
    182185      DO k = 1, klev
    183186      DO i = 1, klon
     
    194197         rain(i) = 0.0
    195198         snow(i) = 0.0
    196       ENDDO
     199         zoliq(i)=0.
     200c     ENDDO
    197201c
    198202c Initialiser le flux de precipitation a zero
    199203c
    200       DO i = 1, klon
     204c     DO i = 1, klon
    201205         zrfl(i) = 0.0
    202206         zneb(i) = seuil_neb
     
    441445         zrhol(i) = zrho(i) * zoliq(i) / zneb(i)
    442446
    443          if (ptconv(i,k)) then
    444             zcl(i)=cld_lc_con
    445             zct(i)=1./cld_tau_con
    446          else
    447             zcl(i)=cld_lc_lsc
    448             zct(i)=1./cld_tau_lsc
    449          endif
    450 c  quantit�d'eau ��minier.
    451          zchau(i) = zct(i)*dtime/FLOAT(ninter) * zoliq(i)
    452      .         *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl(i))**2)) *(1.-zfice(i))
    453 c  meme chose pour la glace.
    454          if (ptconv(i,k)) then
    455             zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
     447         IF (zneb(i).EQ.seuil_neb) THEN
     448             ztot = 0.0
     449         ELSE
     450c  quantite d'eau a eliminer: zchau
     451c  meme chose pour la glace: zfroi
     452             if (ptconv(i,k)) then
     453                zcl   =cld_lc_con
     454                zct   =1./cld_tau_con
     455                zfroi    = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
    456456     .              *fallvc(zrhol(i)) * zfice(i)
    457          else
    458             zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
     457             else
     458                zcl   =cld_lc_lsc
     459                zct   =1./cld_tau_lsc
     460                zfroi    = dtime/FLOAT(ninter)/zdz(i)*zoliq(i)
    459461     .              *fallvs(zrhol(i)) * zfice(i)
    460          endif
    461          ztot(i) = zchau(i) + zfroi(i)
    462          IF (zneb(i).EQ.seuil_neb) ztot(i) = 0.0
    463          ztot(i) = MIN(MAX(ztot(i),0.0),zoliq(i))
    464          zoliq(i) = MAX(zoliq(i)-ztot(i), 0.0)
     462             endif
     463             zchau    = zct   *dtime/FLOAT(ninter) * zoliq(i)
     464     .         *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl   )**2)) *(1.-zfice(i))
     465             ztot    = zchau    + zfroi
     466             ztot    = MAX(ztot   ,0.0)
     467         ENDIF
     468         ztot    = MIN(ztot,zoliq(i))
     469         zoliq(i) = MAX(zoliq(i)-ztot   , 0.0)
    465470         radliq(i,k) = radliq(i,k) + zoliq(i)/FLOAT(ninter+1)
    466471      ENDIF
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/hbtm.F

    r776 r1250  
    762762            endif
    763763c
     764            qsatbef(i) = qqsat    ! bug dans la version orig ???
    764765          endif
    765           qsatbef(i) = qqsat
    766766camn ???? cette ligne a deja ete faite normalement ?
    767767        endif
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F

    r1231 r1250  
    155155c$OMP THREADPRIVATE(initialized)
    156156cjq-end
     157      INTEGER JK
     158      REAL tmp_
    157159      if(.not.initialized) then
    158160        flag_aer=0.
     
    162164        allocate(ZFSUPAI(KDLON,KFLEV+1))
    163165        allocate(ZFSDNAI(KDLON,KFLEV+1))
    164         ZFSUPAD(:,:)=0.
    165         ZFSDNAD(:,:)=0.
    166         ZFSUPAI(:,:)=0.
    167         ZFSDNAI(:,:)=0.
    168        
     166        DO JK = 1 , KDLON*(KFLEV+1)
     167          ZFSUPAD(JK,1) = 0.0     ! ZFSUPAD(:,:)=0.
     168          ZFSDNAD(JK,1) = 0.0     ! ZFSDNAD(:,:)=0.
     169          ZFSUPAI(JK,1) = 0.0     ! ZFSUPAI(:,:)=0.
     170          ZFSDNAI(JK,1) = 0.0     ! ZFSDNAI(:,:)=0.
     171        END DO
    169172      endif
    170173!rv
     
    179182      IF (MOD(itapsw,swpas).EQ.0) THEN
    180183c
     184      tmp_ = 1./( dobson_u * 1e3 * RG)
     185!cdir collapse
    181186      DO JK = 1 , KFLEV
    182       DO JL = 1, KDLON
    183          ZCLDSW0(JL,JK) = 0.0
    184          ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG * PDP(JL,JK)
    185       ENDDO
     187        DO JL = 1, KDLON
     188          ZCLDSW0(JL,JK) = 0.0
     189          ZOZ(JL,JK) = POZON(JL,JK)*tmp_*PDP(JL,JK)
     190        ENDDO
    186191      ENDDO
    187192C
     
    342347     S                PRMU,PSEC,PUD)
    343348      USE dimphy
     349      USE radiation_AR4_param, only :
     350     S     ZPDH2O,ZPDUMG,ZPRH2O,ZPRUMG,RTDH2O,RTDUMG,RTH2O,RTUMG
    344351      IMPLICIT none
    345352cym#include "dimensions.h"
     
    389396      INTEGER jl, jk, jkp1, jkl, jklp1, ja
    390397C
    391 C* Prescribed Data:
    392 c
    393       REAL(KIND=8) ZPDH2O,ZPDUMG
    394       SAVE ZPDH2O,ZPDUMG
    395 c$OMP THREADPRIVATE(ZPDH2O,ZPDUMG)
    396       REAL(KIND=8) ZPRH2O,ZPRUMG
    397       SAVE ZPRH2O,ZPRUMG
    398 c$OMP THREADPRIVATE(ZPRH2O,ZPRUMG)
    399       REAL(KIND=8) RTDH2O,RTDUMG
    400       SAVE RTDH2O,RTDUMG
    401 c$OMP THREADPRIVATE(RTDH2O,RTDUMG)
    402       REAL(KIND=8) RTH2O ,RTUMG
    403       SAVE RTH2O ,RTUMG
    404 c$OMP THREADPRIVATE(RTH2O ,RTUMG)
    405       DATA ZPDH2O,ZPDUMG / 0.8   , 0.75 /
    406       DATA ZPRH2O,ZPRUMG / 30000., 30000. /
    407       DATA RTDH2O,RTDUMG /  0.40  , 0.375 /
    408       DATA RTH2O ,RTUMG  /  240.  , 240.  /
    409398C     ------------------------------------------------------------------
    410399C
     
    537526     S  ,  PFD   , PFU)
    538527      USE dimphy
     528      USE radiation_AR4_param, only : RSUN, RRAY
    539529      IMPLICIT none
    540530cym#include "dimensions.h"
     
    623613C
    624614      INTEGER jl, jk, k, jaj, ikm1, ikl
    625 c
    626 c Prescribed Data:
    627 c
    628       REAL(KIND=8) RSUN(2)
    629       SAVE RSUN
    630 c$OMP THREADPRIVATE(RSUN)
    631       REAL(KIND=8) RRAY(2,6)
    632       SAVE RRAY
    633 c$OMP THREADPRIVATE(RRAY)
    634       DATA RSUN(1) / 0.441676 /
    635       DATA RSUN(2) / 0.558324 /
    636       DATA (RRAY(1,K),K=1,6) /
    637      S .428937E-01, .890743E+00,-.288555E+01,
    638      S .522744E+01,-.469173E+01, .161645E+01/
    639       DATA (RRAY(2,K),K=1,6) /
    640      S .697200E-02, .173297E-01,-.850903E-01,
    641      S .248261E+00,-.302031E+00, .129662E+00/
     615
    642616C     ------------------------------------------------------------------
    643617C
     
    780754     S  ,  PFDOWN,PFUP                                            )
    781755      USE dimphy
     756      USE radiation_AR4_param, only : RSUN, RRAY
    782757      IMPLICIT none
    783758cym#include "dimensions.h"
     
    900875      REAL(KIND=8) ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11
    901876C
    902 C* Prescribed Data:
    903 C
    904       REAL(KIND=8) RSUN(2)
    905       SAVE RSUN
    906 c$OMP THREADPRIVATE(RSUN)
    907       REAL(KIND=8) RRAY(2,6)
    908       SAVE RRAY
    909 c$OMP THREADPRIVATE(RRAY)
    910       DATA RSUN(1) / 0.441676 /
    911       DATA RSUN(2) / 0.558324 /
    912       DATA (RRAY(1,K),K=1,6) /
    913      S .428937E-01, .890743E+00,-.288555E+01,
    914      S .522744E+01,-.469173E+01, .161645E+01/
    915       DATA (RRAY(2,K),K=1,6) /
    916      S .697200E-02, .173297E-01,-.850903E-01,
    917      S .248261E+00,-.302031E+00, .129662E+00/
     877
    918878C
    919879C     ------------------------------------------------------------------
     
    13291289     S  , PRK   , PRMU0 , PTAUAZ, PTRA1 , PTRA2                   )
    13301290      USE dimphy
     1291      USE radiation_AR4_param, only : TAUA, RPIZA, RCGA
    13311292      IMPLICIT none
    13321293cym#include "dimensions.h"
     
    13981359      REAL(KIND=8) ZBMU0, ZBMU1, ZRE11
    13991360C
    1400 C* Prescribed Data for Aerosols:
    1401 C
    1402       REAL(KIND=8) TAUA(2,5), RPIZA(2,5), RCGA(2,5)
    1403       SAVE TAUA, RPIZA, RCGA
    1404 c$OMP THREADPRIVATE(TAUA, RPIZA, RCGA)
    1405       DATA ((TAUA(IN,JA),JA=1,5),IN=1,2) /
    1406      S .730719, .912819, .725059, .745405, .682188 ,
    1407      S .730719, .912819, .725059, .745405, .682188 /
    1408       DATA ((RPIZA(IN,JA),JA=1,5),IN=1,2) /
    1409      S .872212, .982545, .623143, .944887, .997975 ,
    1410      S .872212, .982545, .623143, .944887, .997975 /
    1411       DATA ((RCGA (IN,JA),JA=1,5),IN=1,2) /
    1412      S .647596, .739002, .580845, .662657, .624246 ,
    1413      S .647596, .739002, .580845, .662657, .624246 /
     1361
    14141362C     ------------------------------------------------------------------
    14151363C
     
    14191367 100  CONTINUE
    14201368C
     1369!cdir collapse
    14211370      DO 103 JK = 1 , KFLEV+1
    14221371      DO 102 JA = 1 , 6
     
    21792128      SUBROUTINE SWTT_LMDAR4 (KNU,KA,PU,PTR)
    21802129      USE dimphy
     2130      USE radiation_AR4_param, only : APAD, BPAD, D
    21812131      IMPLICIT none
    21822132cym#include "dimensions.h"
     
    22252175      INTEGER jl, i,j
    22262176C
    2227 C* Prescribed Data:
    2228 C
    2229       REAL(KIND=8) APAD(2,3,7), BPAD(2,3,7), D(2,3)
    2230       SAVE APAD, BPAD, D
    2231 c$OMP THREADPRIVATE(APAD, BPAD, D)
    2232       DATA ((APAD(1,I,J),I=1,3),J=1,7) /
    2233      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2234      S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01,
    2235      S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00,
    2236      S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02,
    2237      S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02,
    2238      S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02,
    2239      S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 /
    2240       DATA ((APAD(2,I,J),I=1,3),J=1,7) /
    2241      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2242      S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02,
    2243      S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00,
    2244      S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00,
    2245      S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00,
    2246      S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00,
    2247      S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /
    2248 C
    2249       DATA ((BPAD(1,I,J),I=1,3),J=1,7) /
    2250      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2251      S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01,
    2252      S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00,
    2253      S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02,
    2254      S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02,
    2255      S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02,
    2256      S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 /
    2257       DATA ((BPAD(2,I,J),I=1,3),J=1,7) /
    2258      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2259      S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02,
    2260      S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01,
    2261      S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00,
    2262      S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00,
    2263      S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00,
    2264      S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /
    2265 c
    2266       DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 /
    2267       DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 /
     2177
    22682178C
    22692179C-----------------------------------------------------------------------
     
    22972207      SUBROUTINE SWTT1_LMDAR4(KNU,KABS,KIND, PU, PTR)
    22982208      USE dimphy
     2209      USE radiation_AR4_param, only : APAD, BPAD, D
    22992210      IMPLICIT none
    23002211cym#include "dimensions.h"
     
    23452256      INTEGER jl, ja, i, j, ia
    23462257C
    2347 C* Prescribed Data:
    2348 C
    2349       REAL(KIND=8) APAD(2,3,7), BPAD(2,3,7), D(2,3)
    2350       SAVE APAD, BPAD, D
    2351 c$OMP THREADPRIVATE(APAD, BPAD, D)
    2352       DATA ((APAD(1,I,J),I=1,3),J=1,7) /
    2353      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2354      S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01,
    2355      S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00,
    2356      S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02,
    2357      S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02,
    2358      S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02,
    2359      S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 /
    2360       DATA ((APAD(2,I,J),I=1,3),J=1,7) /
    2361      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2362      S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02,
    2363      S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00,
    2364      S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00,
    2365      S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00,
    2366      S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00,
    2367      S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /
    2368 C
    2369       DATA ((BPAD(1,I,J),I=1,3),J=1,7) /
    2370      S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04,
    2371      S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01,
    2372      S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00,
    2373      S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02,
    2374      S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02,
    2375      S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02,
    2376      S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 /
    2377       DATA ((BPAD(2,I,J),I=1,3),J=1,7) /
    2378      S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03,
    2379      S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02,
    2380      S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01,
    2381      S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00,
    2382      S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00,
    2383      S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00,
    2384      S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /
    2385 c
    2386       DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 /
    2387       DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 /
     2258
    23882259C-----------------------------------------------------------------------
    23892260C
     
    26142485     S               PABCU)
    26152486      USE dimphy
     2487      USE radiation_AR4_param, only : TREF, RT1, RAER, AT, BT, OCT
    26162488      IMPLICIT none
    26172489cym#include "dimensions.h"
     
    27072579      REAL(KIND=8) R10E  ! DECIMAL/NATURAL LOG.FACTOR
    27082580      PARAMETER (R10E=0.4342945)
    2709 c
    2710 c Used Data Block:
    2711 c
    2712       REAL(KIND=8) TREF
    2713       SAVE TREF
    2714 c$OMP THREADPRIVATE(TREF)
    2715       REAL(KIND=8) RT1(2)
    2716       SAVE RT1
    2717 c$OMP THREADPRIVATE(RT1)
    2718       REAL(KIND=8) RAER(5,5)
    2719       SAVE RAER
    2720 c$OMP THREADPRIVATE(RAER)
    2721       REAL(KIND=8) AT(8,3), BT(8,3)
    2722       SAVE AT, BT
    2723 c$OMP THREADPRIVATE(AT, BT)
    2724       REAL(KIND=8) OCT(4)
    2725       SAVE OCT
    2726 c$OMP THREADPRIVATE(OCT)
    2727       DATA TREF /250.0/
    2728       DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 /
    2729       DATA RAER / .038520, .037196, .040532, .054934, .038520
    2730      1          , .12613 , .18313 , .10357 , .064106, .126130
    2731      2          , .012579, .013649, .018652, .025181, .012579
    2732      3          , .011890, .016142, .021105, .028908, .011890
    2733      4          , .013792, .026810, .052203, .066338, .013792 /
    2734       DATA (AT(1,IR),IR=1,3) /
    2735      S 0.298199E-02,-.394023E-03,0.319566E-04 /
    2736       DATA (BT(1,IR),IR=1,3) /
    2737      S-0.106432E-04,0.660324E-06,0.174356E-06 /
    2738       DATA (AT(2,IR),IR=1,3) /
    2739      S 0.143676E-01,0.366501E-02,-.160822E-02 /
    2740       DATA (BT(2,IR),IR=1,3) /
    2741      S-0.553979E-04,-.101701E-04,0.920868E-05 /
    2742       DATA (AT(3,IR),IR=1,3) /
    2743      S 0.197861E-01,0.315541E-02,-.174547E-02 /
    2744       DATA (BT(3,IR),IR=1,3) /
    2745      S-0.877012E-04,0.513302E-04,0.523138E-06 /
    2746       DATA (AT(4,IR),IR=1,3) /
    2747      S 0.289560E-01,-.208807E-02,-.121943E-02 /
    2748       DATA (BT(4,IR),IR=1,3) /
    2749      S-0.165960E-03,0.157704E-03,-.146427E-04 /
    2750       DATA (AT(5,IR),IR=1,3) /
    2751      S 0.103800E-01,0.436296E-02,-.161431E-02 /
    2752       DATA (BT(5,IR),IR=1,3) /
    2753      S -.276744E-04,-.327381E-04,0.127646E-04 /
    2754       DATA (AT(6,IR),IR=1,3) /
    2755      S 0.868859E-02,-.972752E-03,0.000000E-00 /
    2756       DATA (BT(6,IR),IR=1,3) /
    2757      S -.278412E-04,-.713940E-06,0.117469E-05 /
    2758       DATA (AT(7,IR),IR=1,3) /
    2759      S 0.250073E-03,0.455875E-03,0.109242E-03 /
    2760       DATA (BT(7,IR),IR=1,3) /
    2761      S 0.199846E-05,-.216313E-05,0.175991E-06 /
    2762       DATA (AT(8,IR),IR=1,3) /
    2763      S 0.307423E-01,0.110879E-02,-.322172E-03 /
    2764       DATA (BT(8,IR),IR=1,3) /
    2765      S-0.108482E-03,0.258096E-05,-.814575E-06 /
    2766 c
    2767       DATA OCT /-.326E-03, -.102E-05, .137E-02, -.535E-05/
     2581
    27682582C-----------------------------------------------------------------------
    27692583c
     
    35103324     S  , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP)
    35113325      USE dimphy
     3326      USE radiation_AR4_param, only : TINTP, XP, GA, GB
    35123327      IMPLICIT none
    35133328cym#include "dimensions.h"
     
    36043419C
    36053420C* Used Data Block:
    3606       REAL(KIND=8) TINTP(11)
    3607       SAVE TINTP
     3421c     REAL*8 TINTP(11)
     3422c     SAVE TINTP
    36083423c$OMP THREADPRIVATE(TINTP)
    3609       REAL(KIND=8) GA(11,16,3), GB(11,16,3)
    3610       SAVE GA, GB
     3424c     REAL*8 GA(11,16,3), GB(11,16,3)
     3425c     SAVE GA, GB
    36113426c$OMP THREADPRIVATE(GA, GB)
    3612       REAL(KIND=8) XP(6,6)
    3613       SAVE XP
     3427c     REAL*8 XP(6,6)
     3428c     SAVE XP
    36143429c$OMP THREADPRIVATE(XP)
    36153430c
    3616       DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
    3617      S             262.5, 275., 287.5, 300., 312.5 /
     3431c     DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250.,
     3432c    S             262.5, 275., 287.5, 300., 312.5 /
    36183433C-----------------------------------------------------------------------
    36193434C-- WATER VAPOR -- INT.1 -- 0- 500 CM-1 -- FROM ABS225 ----------------
     
    36283443C
    36293444C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3630       DATA (GA( 1, 1,IC),IC=1,3) /
    3631      S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/
    3632       DATA (GB( 1, 1,IC),IC=1,3) /
    3633      S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/
    3634       DATA (GA( 1, 2,IC),IC=1,3) /
    3635      S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/
    3636       DATA (GB( 1, 2,IC),IC=1,3) /
    3637      S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/
     3445C     DATA (GA( 1, 1,IC),IC=1,3) /
     3446C    S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/
     3447C     DATA (GB( 1, 1,IC),IC=1,3) /
     3448C    S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/
     3449C     DATA (GA( 1, 2,IC),IC=1,3) /
     3450C    S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/
     3451C     DATA (GB( 1, 2,IC),IC=1,3) /
     3452C    S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/
    36383453C
    36393454C----- INTERVAL = 1 ----- T =  200.0
    36403455C
    36413456C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3642       DATA (GA( 2, 1,IC),IC=1,3) /
    3643      S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/
    3644       DATA (GB( 2, 1,IC),IC=1,3) /
    3645      S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/
    3646       DATA (GA( 2, 2,IC),IC=1,3) /
    3647      S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/
    3648       DATA (GB( 2, 2,IC),IC=1,3) /
    3649      S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/
     3457C     DATA (GA( 2, 1,IC),IC=1,3) /
     3458C    S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/
     3459C     DATA (GB( 2, 1,IC),IC=1,3) /
     3460C    S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/
     3461C     DATA (GA( 2, 2,IC),IC=1,3) /
     3462C    S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/
     3463C     DATA (GB( 2, 2,IC),IC=1,3) /
     3464C    S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/
    36503465C
    36513466C----- INTERVAL = 1 ----- T =  212.5
    36523467C
    36533468C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3654       DATA (GA( 3, 1,IC),IC=1,3) /
    3655      S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/
    3656       DATA (GB( 3, 1,IC),IC=1,3) /
    3657      S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/
    3658       DATA (GA( 3, 2,IC),IC=1,3) /
    3659      S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/
    3660       DATA (GB( 3, 2,IC),IC=1,3) /
    3661      S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/
     3469C     DATA (GA( 3, 1,IC),IC=1,3) /
     3470C    S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/
     3471C     DATA (GB( 3, 1,IC),IC=1,3) /
     3472C    S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/
     3473C     DATA (GA( 3, 2,IC),IC=1,3) /
     3474C    S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/
     3475C     DATA (GB( 3, 2,IC),IC=1,3) /
     3476C    S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/
    36623477C
    36633478C----- INTERVAL = 1 ----- T =  225.0
    36643479C
    36653480C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3666       DATA (GA( 4, 1,IC),IC=1,3) /
    3667      S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/
    3668       DATA (GB( 4, 1,IC),IC=1,3) /
    3669      S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/
    3670       DATA (GA( 4, 2,IC),IC=1,3) /
    3671      S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/
    3672       DATA (GB( 4, 2,IC),IC=1,3) /
    3673      S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/
     3481C     DATA (GA( 4, 1,IC),IC=1,3) /
     3482C    S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/
     3483C     DATA (GB( 4, 1,IC),IC=1,3) /
     3484C    S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/
     3485C     DATA (GA( 4, 2,IC),IC=1,3) /
     3486C    S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/
     3487C     DATA (GB( 4, 2,IC),IC=1,3) /
     3488C    S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/
    36743489C
    36753490C----- INTERVAL = 1 ----- T =  237.5
    36763491C
    36773492C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3678       DATA (GA( 5, 1,IC),IC=1,3) /
    3679      S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/
    3680       DATA (GB( 5, 1,IC),IC=1,3) /
    3681      S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/
    3682       DATA (GA( 5, 2,IC),IC=1,3) /
    3683      S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/
    3684       DATA (GB( 5, 2,IC),IC=1,3) /
    3685      S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/
     3493C     DATA (GA( 5, 1,IC),IC=1,3) /
     3494C    S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/
     3495C     DATA (GB( 5, 1,IC),IC=1,3) /
     3496C    S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/
     3497C     DATA (GA( 5, 2,IC),IC=1,3) /
     3498C    S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/
     3499C     DATA (GB( 5, 2,IC),IC=1,3) /
     3500C    S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/
    36863501C
    36873502C----- INTERVAL = 1 ----- T =  250.0
    36883503C
    36893504C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3690       DATA (GA( 6, 1,IC),IC=1,3) /
    3691      S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/
    3692       DATA (GB( 6, 1,IC),IC=1,3) /
    3693      S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/
    3694       DATA (GA( 6, 2,IC),IC=1,3) /
    3695      S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/
    3696       DATA (GB( 6, 2,IC),IC=1,3) /
    3697      S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/
     3505C     DATA (GA( 6, 1,IC),IC=1,3) /
     3506C    S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/
     3507C     DATA (GB( 6, 1,IC),IC=1,3) /
     3508C    S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/
     3509C     DATA (GA( 6, 2,IC),IC=1,3) /
     3510C    S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/
     3511C     DATA (GB( 6, 2,IC),IC=1,3) /
     3512C    S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/
    36983513C
    36993514C----- INTERVAL = 1 ----- T =  262.5
    37003515C
    37013516C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3702       DATA (GA( 7, 1,IC),IC=1,3) /
    3703      S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/
    3704       DATA (GB( 7, 1,IC),IC=1,3) /
    3705      S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/
    3706       DATA (GA( 7, 2,IC),IC=1,3) /
    3707      S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/
    3708       DATA (GB( 7, 2,IC),IC=1,3) /
    3709      S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/
     3517C     DATA (GA( 7, 1,IC),IC=1,3) /
     3518C    S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/
     3519C     DATA (GB( 7, 1,IC),IC=1,3) /
     3520C    S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/
     3521C     DATA (GA( 7, 2,IC),IC=1,3) /
     3522C    S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/
     3523C     DATA (GB( 7, 2,IC),IC=1,3) /
     3524C    S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/
    37103525C
    37113526C----- INTERVAL = 1 ----- T =  275.0
    37123527C
    37133528C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3714       DATA (GA( 8, 1,IC),IC=1,3) /
    3715      S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/
    3716       DATA (GB( 8, 1,IC),IC=1,3) /
    3717      S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/
    3718       DATA (GA( 8, 2,IC),IC=1,3) /
    3719      S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/
    3720       DATA (GB( 8, 2,IC),IC=1,3) /
    3721      S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/
     3529C     DATA (GA( 8, 1,IC),IC=1,3) /
     3530C    S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/
     3531C     DATA (GB( 8, 1,IC),IC=1,3) /
     3532C    S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/
     3533C     DATA (GA( 8, 2,IC),IC=1,3) /
     3534C    S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/
     3535C     DATA (GB( 8, 2,IC),IC=1,3) /
     3536C    S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/
    37223537C
    37233538C----- INTERVAL = 1 ----- T =  287.5
    37243539C
    37253540C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3726       DATA (GA( 9, 1,IC),IC=1,3) /
    3727      S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/
    3728       DATA (GB( 9, 1,IC),IC=1,3) /
    3729      S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/
    3730       DATA (GA( 9, 2,IC),IC=1,3) /
    3731      S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/
    3732       DATA (GB( 9, 2,IC),IC=1,3) /
    3733      S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/
     3541C     DATA (GA( 9, 1,IC),IC=1,3) /
     3542C    S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/
     3543C     DATA (GB( 9, 1,IC),IC=1,3) /
     3544C    S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/
     3545C     DATA (GA( 9, 2,IC),IC=1,3) /
     3546C    S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/
     3547C     DATA (GB( 9, 2,IC),IC=1,3) /
     3548C    S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/
    37343549C
    37353550C----- INTERVAL = 1 ----- T =  300.0
    37363551C
    37373552C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3738       DATA (GA(10, 1,IC),IC=1,3) /
    3739      S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/
    3740       DATA (GB(10, 1,IC),IC=1,3) /
    3741      S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/
    3742       DATA (GA(10, 2,IC),IC=1,3) /
    3743      S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/
    3744       DATA (GB(10, 2,IC),IC=1,3) /
    3745      S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/
     3553C     DATA (GA(10, 1,IC),IC=1,3) /
     3554C    S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/
     3555C     DATA (GB(10, 1,IC),IC=1,3) /
     3556C    S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/
     3557C     DATA (GA(10, 2,IC),IC=1,3) /
     3558C    S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/
     3559C     DATA (GB(10, 2,IC),IC=1,3) /
     3560C    S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/
    37463561C
    37473562C----- INTERVAL = 1 ----- T =  312.5
    37483563C
    37493564C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    3750       DATA (GA(11, 1,IC),IC=1,3) /
    3751      S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/
    3752       DATA (GB(11, 1,IC),IC=1,3) /
    3753      S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/
    3754       DATA (GA(11, 2,IC),IC=1,3) /
    3755      S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/
    3756       DATA (GB(11, 2,IC),IC=1,3) /
    3757      S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/
     3565C     DATA (GA(11, 1,IC),IC=1,3) /
     3566C    S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/
     3567C     DATA (GB(11, 1,IC),IC=1,3) /
     3568C    S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/
     3569C     DATA (GA(11, 2,IC),IC=1,3) /
     3570C    S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/
     3571C     DATA (GB(11, 2,IC),IC=1,3) /
     3572C    S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/
    37583573C
    37593574C
     
    37703585C
    37713586C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3772       DATA (GA( 1, 3,IC),IC=1,3) /
    3773      S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/
    3774       DATA (GB( 1, 3,IC),IC=1,3) /
    3775      S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/
    3776       DATA (GA( 1, 4,IC),IC=1,3) /
    3777      S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/
    3778       DATA (GB( 1, 4,IC),IC=1,3) /
    3779      S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/
     3587C     DATA (GA( 1, 3,IC),IC=1,3) /
     3588C    S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/
     3589C     DATA (GB( 1, 3,IC),IC=1,3) /
     3590C    S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/
     3591C     DATA (GA( 1, 4,IC),IC=1,3) /
     3592C    S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/
     3593C     DATA (GB( 1, 4,IC),IC=1,3) /
     3594C    S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/
    37803595C
    37813596C----- INTERVAL = 2 ----- T =  200.0
    37823597C
    37833598C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3784       DATA (GA( 2, 3,IC),IC=1,3) /
    3785      S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/
    3786       DATA (GB( 2, 3,IC),IC=1,3) /
    3787      S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/
    3788       DATA (GA( 2, 4,IC),IC=1,3) /
    3789      S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/
    3790       DATA (GB( 2, 4,IC),IC=1,3) /
    3791      S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/
     3599C     DATA (GA( 2, 3,IC),IC=1,3) /
     3600C    S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/
     3601C     DATA (GB( 2, 3,IC),IC=1,3) /
     3602C    S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/
     3603C     DATA (GA( 2, 4,IC),IC=1,3) /
     3604C    S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/
     3605C     DATA (GB( 2, 4,IC),IC=1,3) /
     3606C    S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/
    37923607C
    37933608C----- INTERVAL = 2 ----- T =  212.5
    37943609C
    37953610C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3796       DATA (GA( 3, 3,IC),IC=1,3) /
    3797      S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/
    3798       DATA (GB( 3, 3,IC),IC=1,3) /
    3799      S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/
    3800       DATA (GA( 3, 4,IC),IC=1,3) /
    3801      S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/
    3802       DATA (GB( 3, 4,IC),IC=1,3) /
    3803      S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/
     3611C     DATA (GA( 3, 3,IC),IC=1,3) /
     3612C    S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/
     3613C     DATA (GB( 3, 3,IC),IC=1,3) /
     3614C    S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/
     3615C     DATA (GA( 3, 4,IC),IC=1,3) /
     3616C    S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/
     3617C     DATA (GB( 3, 4,IC),IC=1,3) /
     3618C    S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/
    38043619C
    38053620C----- INTERVAL = 2 ----- T =  225.0
    38063621C
    38073622C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3808       DATA (GA( 4, 3,IC),IC=1,3) /
    3809      S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/
    3810       DATA (GB( 4, 3,IC),IC=1,3) /
    3811      S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/
    3812       DATA (GA( 4, 4,IC),IC=1,3) /
    3813      S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/
    3814       DATA (GB( 4, 4,IC),IC=1,3) /
    3815      S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/
     3623C     DATA (GA( 4, 3,IC),IC=1,3) /
     3624C    S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/
     3625C     DATA (GB( 4, 3,IC),IC=1,3) /
     3626C    S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/
     3627C     DATA (GA( 4, 4,IC),IC=1,3) /
     3628C    S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/
     3629C     DATA (GB( 4, 4,IC),IC=1,3) /
     3630C    S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/
    38163631C
    38173632C----- INTERVAL = 2 ----- T =  237.5
    38183633C
    38193634C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3820       DATA (GA( 5, 3,IC),IC=1,3) /
    3821      S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/
    3822       DATA (GB( 5, 3,IC),IC=1,3) /
    3823      S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/
    3824       DATA (GA( 5, 4,IC),IC=1,3) /
    3825      S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/
    3826       DATA (GB( 5, 4,IC),IC=1,3) /
    3827      S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/
     3635C     DATA (GA( 5, 3,IC),IC=1,3) /
     3636C    S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/
     3637C     DATA (GB( 5, 3,IC),IC=1,3) /
     3638C    S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/
     3639C     DATA (GA( 5, 4,IC),IC=1,3) /
     3640C    S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/
     3641C     DATA (GB( 5, 4,IC),IC=1,3) /
     3642C    S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/
    38283643C
    38293644C----- INTERVAL = 2 ----- T =  250.0
    38303645C
    38313646C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3832       DATA (GA( 6, 3,IC),IC=1,3) /
    3833      S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/
    3834       DATA (GB( 6, 3,IC),IC=1,3) /
    3835      S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/
    3836       DATA (GA( 6, 4,IC),IC=1,3) /
    3837      S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/
    3838       DATA (GB( 6, 4,IC),IC=1,3) /
    3839      S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/
     3647C     DATA (GA( 6, 3,IC),IC=1,3) /
     3648C    S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/
     3649C     DATA (GB( 6, 3,IC),IC=1,3) /
     3650C    S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/
     3651C     DATA (GA( 6, 4,IC),IC=1,3) /
     3652C    S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/
     3653C     DATA (GB( 6, 4,IC),IC=1,3) /
     3654C    S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/
    38403655C
    38413656C----- INTERVAL = 2 ----- T =  262.5
    38423657C
    38433658C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3844       DATA (GA( 7, 3,IC),IC=1,3) /
    3845      S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/
    3846       DATA (GB( 7, 3,IC),IC=1,3) /
    3847      S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/
    3848       DATA (GA( 7, 4,IC),IC=1,3) /
    3849      S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/
    3850       DATA (GB( 7, 4,IC),IC=1,3) /
    3851      S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/
     3659C     DATA (GA( 7, 3,IC),IC=1,3) /
     3660C    S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/
     3661C     DATA (GB( 7, 3,IC),IC=1,3) /
     3662C    S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/
     3663C     DATA (GA( 7, 4,IC),IC=1,3) /
     3664C    S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/
     3665C     DATA (GB( 7, 4,IC),IC=1,3) /
     3666C    S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/
    38523667C
    38533668C----- INTERVAL = 2 ----- T =  275.0
    38543669C
    38553670C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3856       DATA (GA( 8, 3,IC),IC=1,3) /
    3857      S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/
    3858       DATA (GB( 8, 3,IC),IC=1,3) /
    3859      S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/
    3860       DATA (GA( 8, 4,IC),IC=1,3) /
    3861      S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/
    3862       DATA (GB( 8, 4,IC),IC=1,3) /
    3863      S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/
     3671C     DATA (GA( 8, 3,IC),IC=1,3) /
     3672C    S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/
     3673C     DATA (GB( 8, 3,IC),IC=1,3) /
     3674C    S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/
     3675C     DATA (GA( 8, 4,IC),IC=1,3) /
     3676C    S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/
     3677C     DATA (GB( 8, 4,IC),IC=1,3) /
     3678C    S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/
    38643679C
    38653680C----- INTERVAL = 2 ----- T =  287.5
    38663681C
    38673682C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3868       DATA (GA( 9, 3,IC),IC=1,3) /
    3869      S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/
    3870       DATA (GB( 9, 3,IC),IC=1,3) /
    3871      S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/
    3872       DATA (GA( 9, 4,IC),IC=1,3) /
    3873      S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/
    3874       DATA (GB( 9, 4,IC),IC=1,3) /
    3875      S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/
     3683C     DATA (GA( 9, 3,IC),IC=1,3) /
     3684C    S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/
     3685C     DATA (GB( 9, 3,IC),IC=1,3) /
     3686C    S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/
     3687C     DATA (GA( 9, 4,IC),IC=1,3) /
     3688C    S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/
     3689C     DATA (GB( 9, 4,IC),IC=1,3) /
     3690C    S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/
    38763691C
    38773692C----- INTERVAL = 2 ----- T =  300.0
    38783693C
    38793694C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3880       DATA (GA(10, 3,IC),IC=1,3) /
    3881      S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/
    3882       DATA (GB(10, 3,IC),IC=1,3) /
    3883      S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/
    3884       DATA (GA(10, 4,IC),IC=1,3) /
    3885      S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/
    3886       DATA (GB(10, 4,IC),IC=1,3) /
    3887      S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/
     3695C     DATA (GA(10, 3,IC),IC=1,3) /
     3696C    S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/
     3697C     DATA (GB(10, 3,IC),IC=1,3) /
     3698C    S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/
     3699C     DATA (GA(10, 4,IC),IC=1,3) /
     3700C    S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/
     3701C     DATA (GB(10, 4,IC),IC=1,3) /
     3702C    S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/
    38883703C
    38893704C----- INTERVAL = 2 ----- T =  312.5
    38903705C
    38913706C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3892       DATA (GA(11, 3,IC),IC=1,3) /
    3893      S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/
    3894       DATA (GB(11, 3,IC),IC=1,3) /
    3895      S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/
    3896       DATA (GA(11, 4,IC),IC=1,3) /
    3897      S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/
    3898       DATA (GB(11, 4,IC),IC=1,3) /
    3899      S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/
     3707C     DATA (GA(11, 3,IC),IC=1,3) /
     3708C    S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/
     3709C     DATA (GB(11, 3,IC),IC=1,3) /
     3710C    S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/
     3711C     DATA (GA(11, 4,IC),IC=1,3) /
     3712C    S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/
     3713C     DATA (GB(11, 4,IC),IC=1,3) /
     3714C    S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/
    39003715C
    39013716C
     
    39163731C
    39173732C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3918       DATA (GA( 1, 7,IC),IC=1,3) /
    3919      S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/
    3920       DATA (GB( 1, 7,IC),IC=1,3) /
    3921      S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/
    3922       DATA (GA( 1, 8,IC),IC=1,3) /
    3923      S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/
    3924       DATA (GB( 1, 8,IC),IC=1,3) /
    3925      S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/
     3733C     DATA (GA( 1, 7,IC),IC=1,3) /
     3734C    S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/
     3735C     DATA (GB( 1, 7,IC),IC=1,3) /
     3736C    S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/
     3737C     DATA (GA( 1, 8,IC),IC=1,3) /
     3738C    S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/
     3739C     DATA (GB( 1, 8,IC),IC=1,3) /
     3740C    S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/
    39263741C
    39273742C----- INTERVAL = 3 ----- T =  200.0
    39283743C
    39293744C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3930       DATA (GA( 2, 7,IC),IC=1,3) /
    3931      S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/
    3932       DATA (GB( 2, 7,IC),IC=1,3) /
    3933      S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/
    3934       DATA (GA( 2, 8,IC),IC=1,3) /
    3935      S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/
    3936       DATA (GB( 2, 8,IC),IC=1,3) /
    3937      S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/
     3745C     DATA (GA( 2, 7,IC),IC=1,3) /
     3746C    S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/
     3747C     DATA (GB( 2, 7,IC),IC=1,3) /
     3748C    S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/
     3749C     DATA (GA( 2, 8,IC),IC=1,3) /
     3750C    S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/
     3751C     DATA (GB( 2, 8,IC),IC=1,3) /
     3752C    S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/
    39383753C
    39393754C----- INTERVAL = 3 ----- T =  212.5
    39403755C
    39413756C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3942       DATA (GA( 3, 7,IC),IC=1,3) /
    3943      S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/
    3944       DATA (GB( 3, 7,IC),IC=1,3) /
    3945      S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/
    3946       DATA (GA( 3, 8,IC),IC=1,3) /
    3947      S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/
    3948       DATA (GB( 3, 8,IC),IC=1,3) /
    3949      S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/
     3757C     DATA (GA( 3, 7,IC),IC=1,3) /
     3758C    S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/
     3759C     DATA (GB( 3, 7,IC),IC=1,3) /
     3760C    S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/
     3761C     DATA (GA( 3, 8,IC),IC=1,3) /
     3762C    S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/
     3763C     DATA (GB( 3, 8,IC),IC=1,3) /
     3764C    S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/
    39503765C
    39513766C----- INTERVAL = 3 ----- T =  225.0
    39523767C
    39533768C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3954       DATA (GA( 4, 7,IC),IC=1,3) /
    3955      S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/
    3956       DATA (GB( 4, 7,IC),IC=1,3) /
    3957      S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/
    3958       DATA (GA( 4, 8,IC),IC=1,3) /
    3959      S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/
    3960       DATA (GB( 4, 8,IC),IC=1,3) /
    3961      S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/
     3769C     DATA (GA( 4, 7,IC),IC=1,3) /
     3770C    S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/
     3771C     DATA (GB( 4, 7,IC),IC=1,3) /
     3772C    S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/
     3773C     DATA (GA( 4, 8,IC),IC=1,3) /
     3774C    S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/
     3775C     DATA (GB( 4, 8,IC),IC=1,3) /
     3776C    S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/
    39623777C
    39633778C----- INTERVAL = 3 ----- T =  237.5
    39643779C
    39653780C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3966       DATA (GA( 5, 7,IC),IC=1,3) /
    3967      S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/
    3968       DATA (GB( 5, 7,IC),IC=1,3) /
    3969      S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/
    3970       DATA (GA( 5, 8,IC),IC=1,3) /
    3971      S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/
    3972       DATA (GB( 5, 8,IC),IC=1,3) /
    3973      S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/
     3781C     DATA (GA( 5, 7,IC),IC=1,3) /
     3782C    S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/
     3783C     DATA (GB( 5, 7,IC),IC=1,3) /
     3784C    S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/
     3785C     DATA (GA( 5, 8,IC),IC=1,3) /
     3786C    S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/
     3787C     DATA (GB( 5, 8,IC),IC=1,3) /
     3788C    S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/
    39743789C
    39753790C----- INTERVAL = 3 ----- T =  250.0
    39763791C
    39773792C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3978       DATA (GA( 6, 7,IC),IC=1,3) /
    3979      S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/
    3980       DATA (GB( 6, 7,IC),IC=1,3) /
    3981      S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/
    3982       DATA (GA( 6, 8,IC),IC=1,3) /
    3983      S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/
    3984       DATA (GB( 6, 8,IC),IC=1,3) /
    3985      S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/
     3793C     DATA (GA( 6, 7,IC),IC=1,3) /
     3794C    S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/
     3795C     DATA (GB( 6, 7,IC),IC=1,3) /
     3796C    S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/
     3797C     DATA (GA( 6, 8,IC),IC=1,3) /
     3798C    S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/
     3799C     DATA (GB( 6, 8,IC),IC=1,3) /
     3800C    S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/
    39863801C
    39873802C----- INTERVAL = 3 ----- T =  262.5
    39883803C
    39893804C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    3990       DATA (GA( 7, 7,IC),IC=1,3) /
    3991      S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/
    3992       DATA (GB( 7, 7,IC),IC=1,3) /
    3993      S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/
    3994       DATA (GA( 7, 8,IC),IC=1,3) /
    3995      S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/
    3996       DATA (GB( 7, 8,IC),IC=1,3) /
    3997      S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/
     3805C     DATA (GA( 7, 7,IC),IC=1,3) /
     3806C    S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/
     3807C     DATA (GB( 7, 7,IC),IC=1,3) /
     3808C    S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/
     3809C     DATA (GA( 7, 8,IC),IC=1,3) /
     3810C    S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/
     3811C     DATA (GB( 7, 8,IC),IC=1,3) /
     3812C    S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/
    39983813C
    39993814C----- INTERVAL = 3 ----- T =  275.0
    40003815C
    40013816C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4002       DATA (GA( 8, 7,IC),IC=1,3) /
    4003      S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/
    4004       DATA (GB( 8, 7,IC),IC=1,3) /
    4005      S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/
    4006       DATA (GA( 8, 8,IC),IC=1,3) /
    4007      S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/
    4008       DATA (GB( 8, 8,IC),IC=1,3) /
    4009      S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/
     3817C     DATA (GA( 8, 7,IC),IC=1,3) /
     3818C    S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/
     3819C     DATA (GB( 8, 7,IC),IC=1,3) /
     3820C    S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/
     3821C     DATA (GA( 8, 8,IC),IC=1,3) /
     3822C    S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/
     3823C     DATA (GB( 8, 8,IC),IC=1,3) /
     3824C    S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/
    40103825C
    40113826C----- INTERVAL = 3 ----- T =  287.5
    40123827C
    40133828C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4014       DATA (GA( 9, 7,IC),IC=1,3) /
    4015      S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/
    4016       DATA (GB( 9, 7,IC),IC=1,3) /
    4017      S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/
    4018       DATA (GA( 9, 8,IC),IC=1,3) /
    4019      S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/
    4020       DATA (GB( 9, 8,IC),IC=1,3) /
    4021      S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/
     3829C     DATA (GA( 9, 7,IC),IC=1,3) /
     3830C    S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/
     3831C     DATA (GB( 9, 7,IC),IC=1,3) /
     3832C    S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/
     3833C     DATA (GA( 9, 8,IC),IC=1,3) /
     3834C    S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/
     3835C     DATA (GB( 9, 8,IC),IC=1,3) /
     3836C    S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/
    40223837C
    40233838C----- INTERVAL = 3 ----- T =  300.0
    40243839C
    40253840C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4026       DATA (GA(10, 7,IC),IC=1,3) /
    4027      S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/
    4028       DATA (GB(10, 7,IC),IC=1,3) /
    4029      S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/
    4030       DATA (GA(10, 8,IC),IC=1,3) /
    4031      S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/
    4032       DATA (GB(10, 8,IC),IC=1,3) /
    4033      S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/
     3841C     DATA (GA(10, 7,IC),IC=1,3) /
     3842C    S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/
     3843C     DATA (GB(10, 7,IC),IC=1,3) /
     3844C    S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/
     3845C     DATA (GA(10, 8,IC),IC=1,3) /
     3846C    S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/
     3847C     DATA (GB(10, 8,IC),IC=1,3) /
     3848C    S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/
    40343849C
    40353850C----- INTERVAL = 3 ----- T =  312.5
    40363851C
    40373852C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4038       DATA (GA(11, 7,IC),IC=1,3) /
    4039      S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/
    4040       DATA (GB(11, 7,IC),IC=1,3) /
    4041      S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/
    4042       DATA (GA(11, 8,IC),IC=1,3) /
    4043      S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/
    4044       DATA (GB(11, 8,IC),IC=1,3) /
    4045      S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/
     3853C     DATA (GA(11, 7,IC),IC=1,3) /
     3854C    S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/
     3855C     DATA (GB(11, 7,IC),IC=1,3) /
     3856C    S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/
     3857C     DATA (GA(11, 8,IC),IC=1,3) /
     3858C    S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/
     3859C     DATA (GB(11, 8,IC),IC=1,3) /
     3860C    S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/
    40463861C
    40473862C
     
    40533868C
    40543869C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4055       DATA (GA( 1, 9,IC),IC=1,3) /
    4056      S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/
    4057       DATA (GB( 1, 9,IC),IC=1,3) /
    4058      S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/
    4059       DATA (GA( 1,10,IC),IC=1,3) /
    4060      S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/
    4061       DATA (GB( 1,10,IC),IC=1,3) /
    4062      S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/
     3870C     DATA (GA( 1, 9,IC),IC=1,3) /
     3871C    S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/
     3872C     DATA (GB( 1, 9,IC),IC=1,3) /
     3873C    S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/
     3874C     DATA (GA( 1,10,IC),IC=1,3) /
     3875C    S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/
     3876C     DATA (GB( 1,10,IC),IC=1,3) /
     3877C    S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/
    40633878C
    40643879C----- INTERVAL = 4 ----- T =  200.0
    40653880C
    40663881C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4067       DATA (GA( 2, 9,IC),IC=1,3) /
    4068      S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/
    4069       DATA (GB( 2, 9,IC),IC=1,3) /
    4070      S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/
    4071       DATA (GA( 2,10,IC),IC=1,3) /
    4072      S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/
    4073       DATA (GB( 2,10,IC),IC=1,3) /
    4074      S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/
     3882C     DATA (GA( 2, 9,IC),IC=1,3) /
     3883C    S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/
     3884C     DATA (GB( 2, 9,IC),IC=1,3) /
     3885C    S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/
     3886C     DATA (GA( 2,10,IC),IC=1,3) /
     3887C    S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/
     3888C     DATA (GB( 2,10,IC),IC=1,3) /
     3889C    S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/
    40753890C
    40763891C----- INTERVAL = 4 ----- T =  212.5
    40773892C
    40783893C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4079       DATA (GA( 3, 9,IC),IC=1,3) /
    4080      S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/
    4081       DATA (GB( 3, 9,IC),IC=1,3) /
    4082      S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/
    4083       DATA (GA( 3,10,IC),IC=1,3) /
    4084      S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/
    4085       DATA (GB( 3,10,IC),IC=1,3) /
    4086      S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/
     3894C     DATA (GA( 3, 9,IC),IC=1,3) /
     3895C    S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/
     3896C     DATA (GB( 3, 9,IC),IC=1,3) /
     3897C    S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/
     3898C     DATA (GA( 3,10,IC),IC=1,3) /
     3899C    S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/
     3900C     DATA (GB( 3,10,IC),IC=1,3) /
     3901C    S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/
    40873902C
    40883903C----- INTERVAL = 4 ----- T =  225.0
    40893904C
    40903905C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4091       DATA (GA( 4, 9,IC),IC=1,3) /
    4092      S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/
    4093       DATA (GB( 4, 9,IC),IC=1,3) /
    4094      S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/
    4095       DATA (GA( 4,10,IC),IC=1,3) /
    4096      S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/
    4097       DATA (GB( 4,10,IC),IC=1,3) /
    4098      S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/
     3906C     DATA (GA( 4, 9,IC),IC=1,3) /
     3907C    S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/
     3908C     DATA (GB( 4, 9,IC),IC=1,3) /
     3909C    S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/
     3910C     DATA (GA( 4,10,IC),IC=1,3) /
     3911C    S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/
     3912C     DATA (GB( 4,10,IC),IC=1,3) /
     3913C    S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/
    40993914C
    41003915C----- INTERVAL = 4 ----- T =  237.5
    41013916C
    41023917C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4103       DATA (GA( 5, 9,IC),IC=1,3) /
    4104      S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/
    4105       DATA (GB( 5, 9,IC),IC=1,3) /
    4106      S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/
    4107       DATA (GA( 5,10,IC),IC=1,3) /
    4108      S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/
    4109       DATA (GB( 5,10,IC),IC=1,3) /
    4110      S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/
     3918C     DATA (GA( 5, 9,IC),IC=1,3) /
     3919C    S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/
     3920C     DATA (GB( 5, 9,IC),IC=1,3) /
     3921C    S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/
     3922C     DATA (GA( 5,10,IC),IC=1,3) /
     3923C    S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/
     3924C     DATA (GB( 5,10,IC),IC=1,3) /
     3925C    S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/
    41113926C
    41123927C----- INTERVAL = 4 ----- T =  250.0
    41133928C
    41143929C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4115       DATA (GA( 6, 9,IC),IC=1,3) /
    4116      S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/
    4117       DATA (GB( 6, 9,IC),IC=1,3) /
    4118      S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/
    4119       DATA (GA( 6,10,IC),IC=1,3) /
    4120      S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/
    4121       DATA (GB( 6,10,IC),IC=1,3) /
    4122      S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/
     3930C     DATA (GA( 6, 9,IC),IC=1,3) /
     3931C    S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/
     3932C     DATA (GB( 6, 9,IC),IC=1,3) /
     3933C    S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/
     3934C     DATA (GA( 6,10,IC),IC=1,3) /
     3935C    S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/
     3936C     DATA (GB( 6,10,IC),IC=1,3) /
     3937C    S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/
    41233938C
    41243939C----- INTERVAL = 4 ----- T =  262.5
    41253940C
    41263941C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4127       DATA (GA( 7, 9,IC),IC=1,3) /
    4128      S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/
    4129       DATA (GB( 7, 9,IC),IC=1,3) /
    4130      S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/
    4131       DATA (GA( 7,10,IC),IC=1,3) /
    4132      S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/
    4133       DATA (GB( 7,10,IC),IC=1,3) /
    4134      S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/
     3942C     DATA (GA( 7, 9,IC),IC=1,3) /
     3943C    S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/
     3944C     DATA (GB( 7, 9,IC),IC=1,3) /
     3945C    S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/
     3946C     DATA (GA( 7,10,IC),IC=1,3) /
     3947C    S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/
     3948C     DATA (GB( 7,10,IC),IC=1,3) /
     3949C    S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/
    41353950C
    41363951C----- INTERVAL = 4 ----- T =  275.0
    41373952C
    41383953C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4139       DATA (GA( 8, 9,IC),IC=1,3) /
    4140      S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/
    4141       DATA (GB( 8, 9,IC),IC=1,3) /
    4142      S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/
    4143       DATA (GA( 8,10,IC),IC=1,3) /
    4144      S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/
    4145       DATA (GB( 8,10,IC),IC=1,3) /
    4146      S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/
     3954C     DATA (GA( 8, 9,IC),IC=1,3) /
     3955C    S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/
     3956C     DATA (GB( 8, 9,IC),IC=1,3) /
     3957C    S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/
     3958C     DATA (GA( 8,10,IC),IC=1,3) /
     3959C    S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/
     3960C     DATA (GB( 8,10,IC),IC=1,3) /
     3961C    S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/
    41473962C
    41483963C----- INTERVAL = 4 ----- T =  287.5
    41493964C
    41503965C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4151       DATA (GA( 9, 9,IC),IC=1,3) /
    4152      S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/
    4153       DATA (GB( 9, 9,IC),IC=1,3) /
    4154      S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/
    4155       DATA (GA( 9,10,IC),IC=1,3) /
    4156      S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/
    4157       DATA (GB( 9,10,IC),IC=1,3) /
    4158      S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/
     3966C     DATA (GA( 9, 9,IC),IC=1,3) /
     3967C    S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/
     3968C     DATA (GB( 9, 9,IC),IC=1,3) /
     3969C    S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/
     3970C     DATA (GA( 9,10,IC),IC=1,3) /
     3971C    S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/
     3972C     DATA (GB( 9,10,IC),IC=1,3) /
     3973C    S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/
    41593974C
    41603975C----- INTERVAL = 4 ----- T =  300.0
    41613976C
    41623977C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4163       DATA (GA(10, 9,IC),IC=1,3) /
    4164      S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/
    4165       DATA (GB(10, 9,IC),IC=1,3) /
    4166      S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/
    4167       DATA (GA(10,10,IC),IC=1,3) /
    4168      S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/
    4169       DATA (GB(10,10,IC),IC=1,3) /
    4170      S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/
     3978C     DATA (GA(10, 9,IC),IC=1,3) /
     3979C    S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/
     3980C     DATA (GB(10, 9,IC),IC=1,3) /
     3981C    S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/
     3982C     DATA (GA(10,10,IC),IC=1,3) /
     3983C    S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/
     3984C     DATA (GB(10,10,IC),IC=1,3) /
     3985C    S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/
    41713986C
    41723987C----- INTERVAL = 4 ----- T =  312.5
    41733988C
    41743989C-- INDICES FOR PADE APPROXIMATION     1   28   37   45
    4175       DATA (GA(11, 9,IC),IC=1,3) /
    4176      S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/
    4177       DATA (GB(11, 9,IC),IC=1,3) /
    4178      S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/
    4179       DATA (GA(11,10,IC),IC=1,3) /
    4180      S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/
    4181       DATA (GB(11,10,IC),IC=1,3) /
    4182      S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/
     3990C     DATA (GA(11, 9,IC),IC=1,3) /
     3991C    S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/
     3992C     DATA (GB(11, 9,IC),IC=1,3) /
     3993C    S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/
     3994C     DATA (GA(11,10,IC),IC=1,3) /
     3995C    S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/
     3996C     DATA (GB(11,10,IC),IC=1,3) /
     3997C    S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/
    41833998C
    41843999C
     
    41934008C
    41944009C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4195       DATA (GA( 1, 5,IC),IC=1,3) /
    4196      S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/
    4197       DATA (GB( 1, 5,IC),IC=1,3) /
    4198      S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/
    4199       DATA (GA( 1, 6,IC),IC=1,3) /
    4200      S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/
    4201       DATA (GB( 1, 6,IC),IC=1,3) /
    4202      S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/
     4010C     DATA (GA( 1, 5,IC),IC=1,3) /
     4011C    S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/
     4012C     DATA (GB( 1, 5,IC),IC=1,3) /
     4013C    S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/
     4014C     DATA (GA( 1, 6,IC),IC=1,3) /
     4015C    S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/
     4016C     DATA (GB( 1, 6,IC),IC=1,3) /
     4017C    S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/
    42034018C
    42044019C----- INTERVAL = 5 ----- T =  200.0
    42054020C
    42064021C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4207       DATA (GA( 2, 5,IC),IC=1,3) /
    4208      S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/
    4209       DATA (GB( 2, 5,IC),IC=1,3) /
    4210      S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/
    4211       DATA (GA( 2, 6,IC),IC=1,3) /
    4212      S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/
    4213       DATA (GB( 2, 6,IC),IC=1,3) /
    4214      S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/
     4022C     DATA (GA( 2, 5,IC),IC=1,3) /
     4023C    S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/
     4024C     DATA (GB( 2, 5,IC),IC=1,3) /
     4025C    S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/
     4026C     DATA (GA( 2, 6,IC),IC=1,3) /
     4027C    S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/
     4028C     DATA (GB( 2, 6,IC),IC=1,3) /
     4029C    S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/
    42154030C
    42164031C----- INTERVAL = 5 ----- T =  212.5
    42174032C
    42184033C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4219       DATA (GA( 3, 5,IC),IC=1,3) /
    4220      S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/
    4221       DATA (GB( 3, 5,IC),IC=1,3) /
    4222      S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/
    4223       DATA (GA( 3, 6,IC),IC=1,3) /
    4224      S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/
    4225       DATA (GB( 3, 6,IC),IC=1,3) /
    4226      S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/
     4034C     DATA (GA( 3, 5,IC),IC=1,3) /
     4035C    S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/
     4036C     DATA (GB( 3, 5,IC),IC=1,3) /
     4037C    S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/
     4038C     DATA (GA( 3, 6,IC),IC=1,3) /
     4039C    S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/
     4040C     DATA (GB( 3, 6,IC),IC=1,3) /
     4041C    S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/
    42274042C
    42284043C----- INTERVAL = 5 ----- T =  225.0
    42294044C
    42304045C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4231       DATA (GA( 4, 5,IC),IC=1,3) /
    4232      S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/
    4233       DATA (GB( 4, 5,IC),IC=1,3) /
    4234      S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/
    4235       DATA (GA( 4, 6,IC),IC=1,3) /
    4236      S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/
    4237       DATA (GB( 4, 6,IC),IC=1,3) /
    4238      S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/
     4046C     DATA (GA( 4, 5,IC),IC=1,3) /
     4047C    S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/
     4048C     DATA (GB( 4, 5,IC),IC=1,3) /
     4049C    S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/
     4050C     DATA (GA( 4, 6,IC),IC=1,3) /
     4051C    S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/
     4052C     DATA (GB( 4, 6,IC),IC=1,3) /
     4053C    S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/
    42394054C
    42404055C----- INTERVAL = 5 ----- T =  237.5
    42414056C
    42424057C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4243       DATA (GA( 5, 5,IC),IC=1,3) /
    4244      S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/
    4245       DATA (GB( 5, 5,IC),IC=1,3) /
    4246      S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/
    4247       DATA (GA( 5, 6,IC),IC=1,3) /
    4248      S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/
    4249       DATA (GB( 5, 6,IC),IC=1,3) /
    4250      S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/
     4058C     DATA (GA( 5, 5,IC),IC=1,3) /
     4059C    S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/
     4060C     DATA (GB( 5, 5,IC),IC=1,3) /
     4061C    S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/
     4062C     DATA (GA( 5, 6,IC),IC=1,3) /
     4063C    S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/
     4064C     DATA (GB( 5, 6,IC),IC=1,3) /
     4065C    S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/
    42514066C
    42524067C----- INTERVAL = 5 ----- T =  250.0
    42534068C
    42544069C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4255       DATA (GA( 6, 5,IC),IC=1,3) /
    4256      S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/
    4257       DATA (GB( 6, 5,IC),IC=1,3) /
    4258      S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/
    4259       DATA (GA( 6, 6,IC),IC=1,3) /
    4260      S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/
    4261       DATA (GB( 6, 6,IC),IC=1,3) /
    4262      S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/
     4070C     DATA (GA( 6, 5,IC),IC=1,3) /
     4071C    S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/
     4072C     DATA (GB( 6, 5,IC),IC=1,3) /
     4073C    S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/
     4074C     DATA (GA( 6, 6,IC),IC=1,3) /
     4075C    S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/
     4076C     DATA (GB( 6, 6,IC),IC=1,3) /
     4077C    S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/
    42634078C
    42644079C----- INTERVAL = 5 ----- T =  262.5
    42654080C
    42664081C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4267       DATA (GA( 7, 5,IC),IC=1,3) /
    4268      S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/
    4269       DATA (GB( 7, 5,IC),IC=1,3) /
    4270      S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/
    4271       DATA (GA( 7, 6,IC),IC=1,3) /
    4272      S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/
    4273       DATA (GB( 7, 6,IC),IC=1,3) /
    4274      S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/
     4082C     DATA (GA( 7, 5,IC),IC=1,3) /
     4083C    S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/
     4084C     DATA (GB( 7, 5,IC),IC=1,3) /
     4085C    S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/
     4086C     DATA (GA( 7, 6,IC),IC=1,3) /
     4087C    S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/
     4088C     DATA (GB( 7, 6,IC),IC=1,3) /
     4089C    S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/
    42754090C
    42764091C----- INTERVAL = 5 ----- T =  275.0
    42774092C
    42784093C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4279       DATA (GA( 8, 5,IC),IC=1,3) /
    4280      S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/
    4281       DATA (GB( 8, 5,IC),IC=1,3) /
    4282      S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/
    4283       DATA (GA( 8, 6,IC),IC=1,3) /
    4284      S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/
    4285       DATA (GB( 8, 6,IC),IC=1,3) /
    4286      S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/
     4094C     DATA (GA( 8, 5,IC),IC=1,3) /
     4095C    S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/
     4096C     DATA (GB( 8, 5,IC),IC=1,3) /
     4097C    S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/
     4098C     DATA (GA( 8, 6,IC),IC=1,3) /
     4099C    S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/
     4100C     DATA (GB( 8, 6,IC),IC=1,3) /
     4101C    S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/
    42874102C
    42884103C----- INTERVAL = 5 ----- T =  287.5
    42894104C
    42904105C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4291       DATA (GA( 9, 5,IC),IC=1,3) /
    4292      S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/
    4293       DATA (GB( 9, 5,IC),IC=1,3) /
    4294      S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/
    4295       DATA (GA( 9, 6,IC),IC=1,3) /
    4296      S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/
    4297       DATA (GB( 9, 6,IC),IC=1,3) /
    4298      S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/
     4106C     DATA (GA( 9, 5,IC),IC=1,3) /
     4107C    S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/
     4108C     DATA (GB( 9, 5,IC),IC=1,3) /
     4109C    S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/
     4110C     DATA (GA( 9, 6,IC),IC=1,3) /
     4111C    S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/
     4112C     DATA (GB( 9, 6,IC),IC=1,3) /
     4113C    S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/
    42994114C
    43004115C----- INTERVAL = 5 ----- T =  300.0
    43014116C
    43024117C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4303       DATA (GA(10, 5,IC),IC=1,3) /
    4304      S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/
    4305       DATA (GB(10, 5,IC),IC=1,3) /
    4306      S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/
    4307       DATA (GA(10, 6,IC),IC=1,3) /
    4308      S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/
    4309       DATA (GB(10, 6,IC),IC=1,3) /
    4310      S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/
     4118C     DATA (GA(10, 5,IC),IC=1,3) /
     4119C    S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/
     4120C     DATA (GB(10, 5,IC),IC=1,3) /
     4121C    S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/
     4122C     DATA (GA(10, 6,IC),IC=1,3) /
     4123C    S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/
     4124C     DATA (GB(10, 6,IC),IC=1,3) /
     4125C    S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/
    43114126C
    43124127C----- INTERVAL = 5 ----- T =  312.5
    43134128C
    43144129C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4315       DATA (GA(11, 5,IC),IC=1,3) /
    4316      S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/
    4317       DATA (GB(11, 5,IC),IC=1,3) /
    4318      S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/
    4319       DATA (GA(11, 6,IC),IC=1,3) /
    4320      S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/
    4321       DATA (GB(11, 6,IC),IC=1,3) /
    4322      S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/
     4130C     DATA (GA(11, 5,IC),IC=1,3) /
     4131C    S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/
     4132C     DATA (GB(11, 5,IC),IC=1,3) /
     4133C    S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/
     4134C     DATA (GA(11, 6,IC),IC=1,3) /
     4135C    S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/
     4136C     DATA (GB(11, 6,IC),IC=1,3) /
     4137C    S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/
    43234138C
    43244139C
     
    43324147C
    43334148C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4334       DATA (GA( 1,11,IC),IC=1,3) /
    4335      S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/
    4336       DATA (GB( 1,11,IC),IC=1,3) /
    4337      S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/
    4338       DATA (GA( 1,12,IC),IC=1,3) /
    4339      S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/
    4340       DATA (GB( 1,12,IC),IC=1,3) /
    4341      S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/
     4149C     DATA (GA( 1,11,IC),IC=1,3) /
     4150C    S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/
     4151C     DATA (GB( 1,11,IC),IC=1,3) /
     4152C    S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/
     4153C     DATA (GA( 1,12,IC),IC=1,3) /
     4154C    S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/
     4155C     DATA (GB( 1,12,IC),IC=1,3) /
     4156C    S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/
    43424157C
    43434158C----- INTERVAL = 6 ----- T =  200.0
    43444159C
    43454160C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4346       DATA (GA( 2,11,IC),IC=1,3) /
    4347      S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/
    4348       DATA (GB( 2,11,IC),IC=1,3) /
    4349      S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/
    4350       DATA (GA( 2,12,IC),IC=1,3) /
    4351      S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/
    4352       DATA (GB( 2,12,IC),IC=1,3) /
    4353      S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/
     4161C     DATA (GA( 2,11,IC),IC=1,3) /
     4162C    S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/
     4163C     DATA (GB( 2,11,IC),IC=1,3) /
     4164C    S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/
     4165C     DATA (GA( 2,12,IC),IC=1,3) /
     4166C    S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/
     4167C     DATA (GB( 2,12,IC),IC=1,3) /
     4168C    S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/
    43544169C
    43554170C----- INTERVAL = 6 ----- T =  212.5
    43564171C
    43574172C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4358       DATA (GA( 3,11,IC),IC=1,3) /
    4359      S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/
    4360       DATA (GB( 3,11,IC),IC=1,3) /
    4361      S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/
    4362       DATA (GA( 3,12,IC),IC=1,3) /
    4363      S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/
    4364       DATA (GB( 3,12,IC),IC=1,3) /
    4365      S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/
     4173C     DATA (GA( 3,11,IC),IC=1,3) /
     4174C    S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/
     4175C     DATA (GB( 3,11,IC),IC=1,3) /
     4176C    S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/
     4177C     DATA (GA( 3,12,IC),IC=1,3) /
     4178C    S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/
     4179C     DATA (GB( 3,12,IC),IC=1,3) /
     4180C    S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/
    43664181C
    43674182C----- INTERVAL = 6 ----- T =  225.0
    43684183C
    43694184C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4370       DATA (GA( 4,11,IC),IC=1,3) /
    4371      S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/
    4372       DATA (GB( 4,11,IC),IC=1,3) /
    4373      S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/
    4374       DATA (GA( 4,12,IC),IC=1,3) /
    4375      S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/
    4376       DATA (GB( 4,12,IC),IC=1,3) /
    4377      S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/
     4185C     DATA (GA( 4,11,IC),IC=1,3) /
     4186C    S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/
     4187C     DATA (GB( 4,11,IC),IC=1,3) /
     4188C    S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/
     4189C     DATA (GA( 4,12,IC),IC=1,3) /
     4190C    S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/
     4191C     DATA (GB( 4,12,IC),IC=1,3) /
     4192C    S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/
    43784193C
    43794194C----- INTERVAL = 6 ----- T =  237.5
    43804195C
    43814196C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4382       DATA (GA( 5,11,IC),IC=1,3) /
    4383      S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/
    4384       DATA (GB( 5,11,IC),IC=1,3) /
    4385      S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/
    4386       DATA (GA( 5,12,IC),IC=1,3) /
    4387      S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/
    4388       DATA (GB( 5,12,IC),IC=1,3) /
    4389      S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/
     4197C     DATA (GA( 5,11,IC),IC=1,3) /
     4198C    S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/
     4199C     DATA (GB( 5,11,IC),IC=1,3) /
     4200C    S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/
     4201C     DATA (GA( 5,12,IC),IC=1,3) /
     4202C    S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/
     4203C     DATA (GB( 5,12,IC),IC=1,3) /
     4204C    S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/
    43904205C
    43914206C----- INTERVAL = 6 ----- T =  250.0
    43924207C
    43934208C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4394       DATA (GA( 6,11,IC),IC=1,3) /
    4395      S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/
    4396       DATA (GB( 6,11,IC),IC=1,3) /
    4397      S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/
    4398       DATA (GA( 6,12,IC),IC=1,3) /
    4399      S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/
    4400       DATA (GB( 6,12,IC),IC=1,3) /
    4401      S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/
     4209C     DATA (GA( 6,11,IC),IC=1,3) /
     4210C    S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/
     4211C     DATA (GB( 6,11,IC),IC=1,3) /
     4212C    S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/
     4213C     DATA (GA( 6,12,IC),IC=1,3) /
     4214C    S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/
     4215C     DATA (GB( 6,12,IC),IC=1,3) /
     4216C    S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/
    44024217C
    44034218C----- INTERVAL = 6 ----- T =  262.5
    44044219C
    44054220C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4406       DATA (GA( 7,11,IC),IC=1,3) /
    4407      S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/
    4408       DATA (GB( 7,11,IC),IC=1,3) /
    4409      S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/
    4410       DATA (GA( 7,12,IC),IC=1,3) /
    4411      S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/
    4412       DATA (GB( 7,12,IC),IC=1,3) /
    4413      S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/
     4221C     DATA (GA( 7,11,IC),IC=1,3) /
     4222C    S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/
     4223C     DATA (GB( 7,11,IC),IC=1,3) /
     4224C    S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/
     4225C     DATA (GA( 7,12,IC),IC=1,3) /
     4226C    S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/
     4227C     DATA (GB( 7,12,IC),IC=1,3) /
     4228C    S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/
    44144229C
    44154230C----- INTERVAL = 6 ----- T =  275.0
    44164231C
    44174232C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4418       DATA (GA( 8,11,IC),IC=1,3) /
    4419      S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/
    4420       DATA (GB( 8,11,IC),IC=1,3) /
    4421      S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/
    4422       DATA (GA( 8,12,IC),IC=1,3) /
    4423      S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/
    4424       DATA (GB( 8,12,IC),IC=1,3) /
    4425      S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/
     4233C     DATA (GA( 8,11,IC),IC=1,3) /
     4234C    S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/
     4235C     DATA (GB( 8,11,IC),IC=1,3) /
     4236C    S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/
     4237C     DATA (GA( 8,12,IC),IC=1,3) /
     4238C    S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/
     4239C     DATA (GB( 8,12,IC),IC=1,3) /
     4240C    S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/
    44264241C
    44274242C----- INTERVAL = 6 ----- T =  287.5
    44284243C
    44294244C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4430       DATA (GA( 9,11,IC),IC=1,3) /
    4431      S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/
    4432       DATA (GB( 9,11,IC),IC=1,3) /
    4433      S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/
    4434       DATA (GA( 9,12,IC),IC=1,3) /
    4435      S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/
    4436       DATA (GB( 9,12,IC),IC=1,3) /
    4437      S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/
     4245C     DATA (GA( 9,11,IC),IC=1,3) /
     4246C    S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/
     4247C     DATA (GB( 9,11,IC),IC=1,3) /
     4248C    S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/
     4249C     DATA (GA( 9,12,IC),IC=1,3) /
     4250C    S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/
     4251C     DATA (GB( 9,12,IC),IC=1,3) /
     4252C    S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/
    44384253C
    44394254C----- INTERVAL = 6 ----- T =  300.0
    44404255C
    44414256C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4442       DATA (GA(10,11,IC),IC=1,3) /
    4443      S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/
    4444       DATA (GB(10,11,IC),IC=1,3) /
    4445      S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/
    4446       DATA (GA(10,12,IC),IC=1,3) /
    4447      S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/
    4448       DATA (GB(10,12,IC),IC=1,3) /
    4449      S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/
     4257C     DATA (GA(10,11,IC),IC=1,3) /
     4258C    S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/
     4259C     DATA (GB(10,11,IC),IC=1,3) /
     4260C    S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/
     4261C     DATA (GA(10,12,IC),IC=1,3) /
     4262C    S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/
     4263C     DATA (GB(10,12,IC),IC=1,3) /
     4264C    S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/
    44504265C
    44514266C----- INTERVAL = 6 ----- T =  312.5
    44524267C
    44534268C-- INDICES FOR PADE APPROXIMATION   1 35 40 45
    4454       DATA (GA(11,11,IC),IC=1,3) /
    4455      S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/
    4456       DATA (GB(11,11,IC),IC=1,3) /
    4457      S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/
    4458       DATA (GA(11,12,IC),IC=1,3) /
    4459      S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/
    4460       DATA (GB(11,12,IC),IC=1,3) /
    4461      S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/
     4269C     DATA (GA(11,11,IC),IC=1,3) /
     4270C    S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/
     4271C     DATA (GB(11,11,IC),IC=1,3) /
     4272C    S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/
     4273C     DATA (GA(11,12,IC),IC=1,3) /
     4274C    S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/
     4275C     DATA (GB(11,12,IC),IC=1,3) /
     4276C    S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/
    44624277C
    44634278C
     
    44774292C
    44784293C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4479       DATA (GA( 1,13,IC),IC=1,3) /
    4480      S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/
    4481       DATA (GB( 1,13,IC),IC=1,3) /
    4482      S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/
    4483       DATA (GA( 1,14,IC),IC=1,3) /
    4484      S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/
    4485       DATA (GB( 1,14,IC),IC=1,3) /
    4486      S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/
     4294C     DATA (GA( 1,13,IC),IC=1,3) /
     4295C    S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/
     4296C     DATA (GB( 1,13,IC),IC=1,3) /
     4297C    S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/
     4298C     DATA (GA( 1,14,IC),IC=1,3) /
     4299C    S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/
     4300C     DATA (GB( 1,14,IC),IC=1,3) /
     4301C    S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/
    44874302C
    44884303C----- INTERVAL = 2 ----- T =  200.0
    44894304C
    44904305C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4491       DATA (GA( 2,13,IC),IC=1,3) /
    4492      S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/
    4493       DATA (GB( 2,13,IC),IC=1,3) /
    4494      S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/
    4495       DATA (GA( 2,14,IC),IC=1,3) /
    4496      S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/
    4497       DATA (GB( 2,14,IC),IC=1,3) /
    4498      S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/
     4306C     DATA (GA( 2,13,IC),IC=1,3) /
     4307C    S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/
     4308C     DATA (GB( 2,13,IC),IC=1,3) /
     4309C    S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/
     4310C     DATA (GA( 2,14,IC),IC=1,3) /
     4311C    S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/
     4312C     DATA (GB( 2,14,IC),IC=1,3) /
     4313C    S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/
    44994314C
    45004315C----- INTERVAL = 2 ----- T =  212.5
    45014316C
    45024317C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4503       DATA (GA( 3,13,IC),IC=1,3) /
    4504      S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/
    4505       DATA (GB( 3,13,IC),IC=1,3) /
    4506      S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/
    4507       DATA (GA( 3,14,IC),IC=1,3) /
    4508      S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/
    4509       DATA (GB( 3,14,IC),IC=1,3) /
    4510      S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/
     4318C     DATA (GA( 3,13,IC),IC=1,3) /
     4319C    S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/
     4320C     DATA (GB( 3,13,IC),IC=1,3) /
     4321C    S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/
     4322C     DATA (GA( 3,14,IC),IC=1,3) /
     4323C    S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/
     4324C     DATA (GB( 3,14,IC),IC=1,3) /
     4325C    S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/
    45114326C
    45124327C----- INTERVAL = 2 ----- T =  225.0
    45134328C
    45144329C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4515       DATA (GA( 4,13,IC),IC=1,3) /
    4516      S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/
    4517       DATA (GB( 4,13,IC),IC=1,3) /
    4518      S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/
    4519       DATA (GA( 4,14,IC),IC=1,3) /
    4520      S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/
    4521       DATA (GB( 4,14,IC),IC=1,3) /
    4522      S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/
     4330C     DATA (GA( 4,13,IC),IC=1,3) /
     4331C    S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/
     4332C     DATA (GB( 4,13,IC),IC=1,3) /
     4333C    S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/
     4334C     DATA (GA( 4,14,IC),IC=1,3) /
     4335C    S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/
     4336C     DATA (GB( 4,14,IC),IC=1,3) /
     4337C    S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/
    45234338C
    45244339C----- INTERVAL = 2 ----- T =  237.5
    45254340C
    45264341C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4527       DATA (GA( 5,13,IC),IC=1,3) /
    4528      S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/
    4529       DATA (GB( 5,13,IC),IC=1,3) /
    4530      S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/
    4531       DATA (GA( 5,14,IC),IC=1,3) /
    4532      S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/
    4533       DATA (GB( 5,14,IC),IC=1,3) /
    4534      S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/
     4342C     DATA (GA( 5,13,IC),IC=1,3) /
     4343C    S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/
     4344C     DATA (GB( 5,13,IC),IC=1,3) /
     4345C    S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/
     4346C     DATA (GA( 5,14,IC),IC=1,3) /
     4347C    S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/
     4348C     DATA (GB( 5,14,IC),IC=1,3) /
     4349C    S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/
    45354350C
    45364351C----- INTERVAL = 2 ----- T =  250.0
    45374352C
    45384353C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4539       DATA (GA( 6,13,IC),IC=1,3) /
    4540      S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/
    4541       DATA (GB( 6,13,IC),IC=1,3) /
    4542      S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/
    4543       DATA (GA( 6,14,IC),IC=1,3) /
    4544      S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/
    4545       DATA (GB( 6,14,IC),IC=1,3) /
    4546      S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/
     4354C     DATA (GA( 6,13,IC),IC=1,3) /
     4355C    S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/
     4356C     DATA (GB( 6,13,IC),IC=1,3) /
     4357C    S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/
     4358C     DATA (GA( 6,14,IC),IC=1,3) /
     4359C    S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/
     4360C     DATA (GB( 6,14,IC),IC=1,3) /
     4361C    S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/
    45474362C
    45484363C----- INTERVAL = 2 ----- T =  262.5
    45494364C
    45504365C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4551       DATA (GA( 7,13,IC),IC=1,3) /
    4552      S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/
    4553       DATA (GB( 7,13,IC),IC=1,3) /
    4554      S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/
    4555       DATA (GA( 7,14,IC),IC=1,3) /
    4556      S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/
    4557       DATA (GB( 7,14,IC),IC=1,3) /
    4558      S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/
     4366C     DATA (GA( 7,13,IC),IC=1,3) /
     4367C    S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/
     4368C     DATA (GB( 7,13,IC),IC=1,3) /
     4369C    S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/
     4370C     DATA (GA( 7,14,IC),IC=1,3) /
     4371C    S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/
     4372C     DATA (GB( 7,14,IC),IC=1,3) /
     4373C    S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/
    45594374C
    45604375C----- INTERVAL = 2 ----- T =  275.0
    45614376C
    45624377C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4563       DATA (GA( 8,13,IC),IC=1,3) /
    4564      S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/
    4565       DATA (GB( 8,13,IC),IC=1,3) /
    4566      S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/
    4567       DATA (GA( 8,14,IC),IC=1,3) /
    4568      S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/
    4569       DATA (GB( 8,14,IC),IC=1,3) /
    4570      S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/
     4378C     DATA (GA( 8,13,IC),IC=1,3) /
     4379C    S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/
     4380C     DATA (GB( 8,13,IC),IC=1,3) /
     4381C    S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/
     4382C     DATA (GA( 8,14,IC),IC=1,3) /
     4383C    S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/
     4384C     DATA (GB( 8,14,IC),IC=1,3) /
     4385C    S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/
    45714386C
    45724387C----- INTERVAL = 2 ----- T =  287.5
    45734388C
    45744389C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4575       DATA (GA( 9,13,IC),IC=1,3) /
    4576      S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/
    4577       DATA (GB( 9,13,IC),IC=1,3) /
    4578      S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/
    4579       DATA (GA( 9,14,IC),IC=1,3) /
    4580      S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/
    4581       DATA (GB( 9,14,IC),IC=1,3) /
    4582      S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/
     4390C     DATA (GA( 9,13,IC),IC=1,3) /
     4391C    S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/
     4392C     DATA (GB( 9,13,IC),IC=1,3) /
     4393C    S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/
     4394C     DATA (GA( 9,14,IC),IC=1,3) /
     4395C    S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/
     4396C     DATA (GB( 9,14,IC),IC=1,3) /
     4397C    S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/
    45834398C
    45844399C----- INTERVAL = 2 ----- T =  300.0
    45854400C
    45864401C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4587       DATA (GA(10,13,IC),IC=1,3) /
    4588      S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/
    4589       DATA (GB(10,13,IC),IC=1,3) /
    4590      S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/
    4591       DATA (GA(10,14,IC),IC=1,3) /
    4592      S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/
    4593       DATA (GB(10,14,IC),IC=1,3) /
    4594      S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/
     4402C     DATA (GA(10,13,IC),IC=1,3) /
     4403C    S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/
     4404C     DATA (GB(10,13,IC),IC=1,3) /
     4405C    S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/
     4406C     DATA (GA(10,14,IC),IC=1,3) /
     4407C    S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/
     4408C     DATA (GB(10,14,IC),IC=1,3) /
     4409C    S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/
    45954410C
    45964411C----- INTERVAL = 2 ----- T =  312.5
    45974412C
    45984413C-- INDICES FOR PADE APPROXIMATION   1 30 38 45
    4599       DATA (GA(11,13,IC),IC=1,3) /
    4600      S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/
    4601       DATA (GB(11,13,IC),IC=1,3) /
    4602      S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/
    4603       DATA (GA(11,14,IC),IC=1,3) /
    4604      S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/
    4605       DATA (GB(11,14,IC),IC=1,3) /
    4606      S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/
     4414C     DATA (GA(11,13,IC),IC=1,3) /
     4415C    S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/
     4416C     DATA (GB(11,13,IC),IC=1,3) /
     4417C    S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/
     4418C     DATA (GA(11,14,IC),IC=1,3) /
     4419C    S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/
     4420C     DATA (GB(11,14,IC),IC=1,3) /
     4421C    S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/
    46074422C
    46084423C
     
    46244439C
    46254440C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4626       DATA (GA( 1,15,IC),IC=1,3) /
    4627      S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/
    4628       DATA (GB( 1,15,IC),IC=1,3) /
    4629      S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/
    4630       DATA (GA( 1,16,IC),IC=1,3) /
    4631      S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/
    4632       DATA (GB( 1,16,IC),IC=1,3) /
    4633      S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/
     4441C     DATA (GA( 1,15,IC),IC=1,3) /
     4442C    S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/
     4443C     DATA (GB( 1,15,IC),IC=1,3) /
     4444C    S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/
     4445C     DATA (GA( 1,16,IC),IC=1,3) /
     4446C    S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/
     4447C     DATA (GB( 1,16,IC),IC=1,3) /
     4448C    S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/
    46344449C
    46354450C----- INTERVAL = 4 ----- T =  200.0
    46364451C
    46374452C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4638       DATA (GA( 2,15,IC),IC=1,3) /
    4639      S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/
    4640       DATA (GB( 2,15,IC),IC=1,3) /
    4641      S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/
    4642       DATA (GA( 2,16,IC),IC=1,3) /
    4643      S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/
    4644       DATA (GB( 2,16,IC),IC=1,3) /
    4645      S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/
     4453C     DATA (GA( 2,15,IC),IC=1,3) /
     4454C    S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/
     4455C     DATA (GB( 2,15,IC),IC=1,3) /
     4456C    S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/
     4457C     DATA (GA( 2,16,IC),IC=1,3) /
     4458C    S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/
     4459C     DATA (GB( 2,16,IC),IC=1,3) /
     4460C    S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/
    46464461C
    46474462C----- INTERVAL = 4 ----- T =  212.5
    46484463C
    46494464C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4650       DATA (GA( 3,15,IC),IC=1,3) /
    4651      S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/
    4652       DATA (GB( 3,15,IC),IC=1,3) /
    4653      S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/
    4654       DATA (GA( 3,16,IC),IC=1,3) /
    4655      S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/
    4656       DATA (GB( 3,16,IC),IC=1,3) /
    4657      S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/
     4465C     DATA (GA( 3,15,IC),IC=1,3) /
     4466C    S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/
     4467C     DATA (GB( 3,15,IC),IC=1,3) /
     4468C    S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/
     4469C     DATA (GA( 3,16,IC),IC=1,3) /
     4470C    S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/
     4471C     DATA (GB( 3,16,IC),IC=1,3) /
     4472C    S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/
    46584473C
    46594474C----- INTERVAL = 4 ----- T =  225.0
    46604475C
    46614476C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4662       DATA (GA( 4,15,IC),IC=1,3) /
    4663      S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/
    4664       DATA (GB( 4,15,IC),IC=1,3) /
    4665      S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/
    4666       DATA (GA( 4,16,IC),IC=1,3) /
    4667      S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/
    4668       DATA (GB( 4,16,IC),IC=1,3) /
    4669      S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/
     4477C     DATA (GA( 4,15,IC),IC=1,3) /
     4478C    S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/
     4479C     DATA (GB( 4,15,IC),IC=1,3) /
     4480C    S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/
     4481C     DATA (GA( 4,16,IC),IC=1,3) /
     4482C    S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/
     4483C     DATA (GB( 4,16,IC),IC=1,3) /
     4484C    S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/
    46704485C
    46714486C----- INTERVAL = 4 ----- T =  237.5
    46724487C
    46734488C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4674       DATA (GA( 5,15,IC),IC=1,3) /
    4675      S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/
    4676       DATA (GB( 5,15,IC),IC=1,3) /
    4677      S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/
    4678       DATA (GA( 5,16,IC),IC=1,3) /
    4679      S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/
    4680       DATA (GB( 5,16,IC),IC=1,3) /
    4681      S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/
     4489C     DATA (GA( 5,15,IC),IC=1,3) /
     4490C    S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/
     4491C     DATA (GB( 5,15,IC),IC=1,3) /
     4492C    S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/
     4493C     DATA (GA( 5,16,IC),IC=1,3) /
     4494C    S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/
     4495C     DATA (GB( 5,16,IC),IC=1,3) /
     4496C    S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/
    46824497C
    46834498C----- INTERVAL = 4 ----- T =  250.0
    46844499C
    46854500C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4686       DATA (GA( 6,15,IC),IC=1,3) /
    4687      S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/
    4688       DATA (GB( 6,15,IC),IC=1,3) /
    4689      S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/
    4690       DATA (GA( 6,16,IC),IC=1,3) /
    4691      S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/
    4692       DATA (GB( 6,16,IC),IC=1,3) /
    4693      S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/
     4501C     DATA (GA( 6,15,IC),IC=1,3) /
     4502C    S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/
     4503C     DATA (GB( 6,15,IC),IC=1,3) /
     4504C    S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/
     4505C     DATA (GA( 6,16,IC),IC=1,3) /
     4506C    S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/
     4507C     DATA (GB( 6,16,IC),IC=1,3) /
     4508C    S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/
    46944509C
    46954510C----- INTERVAL = 4 ----- T =  262.5
    46964511C
    46974512C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4698       DATA (GA( 7,15,IC),IC=1,3) /
    4699      S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/
    4700       DATA (GB( 7,15,IC),IC=1,3) /
    4701      S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/
    4702       DATA (GA( 7,16,IC),IC=1,3) /
    4703      S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/
    4704       DATA (GB( 7,16,IC),IC=1,3) /
    4705      S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/
     4513C     DATA (GA( 7,15,IC),IC=1,3) /
     4514C    S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/
     4515C     DATA (GB( 7,15,IC),IC=1,3) /
     4516C    S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/
     4517C     DATA (GA( 7,16,IC),IC=1,3) /
     4518C    S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/
     4519C     DATA (GB( 7,16,IC),IC=1,3) /
     4520C    S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/
    47064521C
    47074522C----- INTERVAL = 4 ----- T =  275.0
    47084523C
    47094524C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4710       DATA (GA( 8,15,IC),IC=1,3) /
    4711      S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/
    4712       DATA (GB( 8,15,IC),IC=1,3) /
    4713      S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/
    4714       DATA (GA( 8,16,IC),IC=1,3) /
    4715      S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/
    4716       DATA (GB( 8,16,IC),IC=1,3) /
    4717      S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/
     4525C     DATA (GA( 8,15,IC),IC=1,3) /
     4526C    S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/
     4527C     DATA (GB( 8,15,IC),IC=1,3) /
     4528C    S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/
     4529C     DATA (GA( 8,16,IC),IC=1,3) /
     4530C    S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/
     4531C     DATA (GB( 8,16,IC),IC=1,3) /
     4532C    S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/
    47184533C
    47194534C----- INTERVAL = 4 ----- T =  287.5
    47204535C
    47214536C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4722       DATA (GA( 9,15,IC),IC=1,3) /
    4723      S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/
    4724       DATA (GB( 9,15,IC),IC=1,3) /
    4725      S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/
    4726       DATA (GA( 9,16,IC),IC=1,3) /
    4727      S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/
    4728       DATA (GB( 9,16,IC),IC=1,3) /
    4729      S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/
     4537C     DATA (GA( 9,15,IC),IC=1,3) /
     4538C    S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/
     4539C     DATA (GB( 9,15,IC),IC=1,3) /
     4540C    S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/
     4541C     DATA (GA( 9,16,IC),IC=1,3) /
     4542C    S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/
     4543C     DATA (GB( 9,16,IC),IC=1,3) /
     4544C    S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/
    47304545C
    47314546C----- INTERVAL = 4 ----- T =  300.0
    47324547C
    47334548C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4734       DATA (GA(10,15,IC),IC=1,3) /
    4735      S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/
    4736       DATA (GB(10,15,IC),IC=1,3) /
    4737      S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/
    4738       DATA (GA(10,16,IC),IC=1,3) /
    4739      S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/
    4740       DATA (GB(10,16,IC),IC=1,3) /
    4741      S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/
     4549C     DATA (GA(10,15,IC),IC=1,3) /
     4550C    S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/
     4551C     DATA (GB(10,15,IC),IC=1,3) /
     4552C    S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/
     4553C     DATA (GA(10,16,IC),IC=1,3) /
     4554C    S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/
     4555C     DATA (GB(10,16,IC),IC=1,3) /
     4556C    S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/
    47424557C
    47434558C----- INTERVAL = 4 ----- T =  312.5
    47444559C
    47454560C-- INDICES FOR PADE APPROXIMATION     1   15   29   45
    4746       DATA (GA(11,15,IC),IC=1,3) /
    4747      S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/
    4748       DATA (GB(11,15,IC),IC=1,3) /
    4749      S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/
    4750       DATA (GA(11,16,IC),IC=1,3) /
    4751      S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/
    4752       DATA (GB(11,16,IC),IC=1,3) /
    4753      S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/
     4561C     DATA (GA(11,15,IC),IC=1,3) /
     4562C    S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/
     4563C     DATA (GB(11,15,IC),IC=1,3) /
     4564C    S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/
     4565C     DATA (GA(11,16,IC),IC=1,3) /
     4566C    S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/
     4567C     DATA (GB(11,16,IC),IC=1,3) /
     4568C    S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/
     4569C
     4570C     ------------------------------------------------------------------
     4571C     DATA (( XP(  J,K),J=1,6),       K=1,6) /
     4572C    S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,
     4573C    S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,
     4574C    S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,
     4575C    S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,
     4576C    S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,
     4577C    S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,
     4578C    S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,
     4579C    S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,
     4580C    S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,
     4581C    S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,
     4582C    S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,
     4583C    S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /
    47544584
    4755 C     ------------------------------------------------------------------
    4756       DATA (( XP(  J,K),J=1,6),       K=1,6) /
    4757      S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,
    4758      S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,
    4759      S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,
    4760      S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,
    4761      S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,
    4762      S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,
    4763      S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,
    4764      S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,
    4765      S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,
    4766      S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,
    4767      S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,
    4768      S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /
    47694585C
    47704586C
     
    47744590 100  CONTINUE
    47754591C
     4592!cdir collapse
    47764593      DO 102 JK = 1 , KFLEV+1
    47774594      DO 101 JL = 1, KDLON
     
    56255442     S  , PADJD,PADJU,PCNTRB,PDBDT)
    56265443       USE dimphy
     5444      USE radiation_AR4_param, only : WG1
    56275445      IMPLICIT none
    56285446cym#include "dimensions.h"
     
    56855503      REAL(KIND=8) zwtr
    56865504c
    5687 C* Data Block:
    5688 c
    5689       REAL(KIND=8) WG1(2)
    5690       SAVE WG1
    5691 c$OMP THREADPRIVATE(WG1)
    5692       DATA (WG1(jk),jk=1,2) /1.0, 1.0/
     5505
    56935506C-----------------------------------------------------------------------
    56945507C
     
    59035716C
    59045717C
     5718!cdir collapse
    59055719      DO 130 JA = 1 , 8
    59065720      DO 120 JL = 1, KDLON
     
    60805894C
    60815895C
     5896
     5897!CDIR ON_ADB(PUU1)
     5898!CDIR ON_ADB(PUU2)
     5899!CDIR COLLAPSE
    60825900      DO 130 JA = 1 , 8
    60835901      DO 120 JL = 1, KDLON
  • LMDZ4/branches/LMDZ4-dev/makelmdz_fcm

    r1249 r1250  
    162162# lecture des chemins propres à l'architecture de la machine #
    163163###############################################################
    164 
    165 rm -f ./arch.path
    166 ln -s ./arch/arch-${arch}.path ./arch.path
     164rm -f .void_file
     165echo > .void_file
     166rm -f arch.path
     167ln -s arch/arch-${arch}.path ./arch.path
    167168source arch.path
    168169
     
    360361rm -f bin/${code}${SUFF_NAME}.e
    361362rm -f arch.fcm
     363rm -f arch.opt
    362364
    363365echo "%ARCH          $arch"          >> $config_fcm
     
    380382
    381383ln -s arch/arch-${arch}.fcm arch.fcm
     384if test -f arch/arch-${arch}.opt &&  [ $compil_mod = "prod" ]
     385  then
     386  ln -s arch/arch-${arch}.opt arch.opt
     387else
     388  ln -s .void_file arch.opt
     389fi
     390
     391
    382392rm -f $LIBOGCM/${arch}${SUFF_NAME}/.config/fcm.bld.lock
    383393./build_gcm
     
    387397ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config config
    388398ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config/tmp tmp_src
    389 #\rm -f $libf/grid/dimensions.h
Note: See TracChangeset for help on using the changeset viewer.