*=*=*=*= SULW.html =*=*=*=*
SUBROUTINE SULW

SUBROUTINE SULW


      SUBROUTINE SULW
C
#include "dimensions.h"
#include "dimphys.h"
#include "dimradmars.h"
#include "callkeys.h"
C
C**** *SULW*   - INITIALIZE COMMON YOMLW
C
C     PURPOSE.
C     --------
C           INITIALIZE YOMLW, THE COMMON THAT CONTAINS COEFFICIENTS
C           NEEDED TO RUN THE LONGWAVE RADIATION ROUTINES
C
C**   INTERFACE.
C     ----------
C        *CALL* *SULW
C
C        EXPLICIT ARGUMENTS :
C        --------------------
C        NONE
C
C        IMPLICIT ARGUMENTS :
C        --------------------
C        COMMON YOMLW
C
C     METHOD.
C     -------
C        SEE DOCUMENTATION
C
C     EXTERNALS.
C     ----------
C
C     REFERENCE.
C     ----------
C        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
C     "IN CORE MODEL"
C
C     AUTHOR.
C     -------
C        JEAN-JACQUES MORCRETTE  *ECMWF*
C
C     MODIFICATIONS.
C     --------------
C        ORIGINAL : 88-12-15
C     ------------------------------------------------------------------
#include "yomlw.h"
C      ----------------------------------------------------------------
C
C*       1.    SET VALUES.
C              -----------
C
C
C-- REFERENCE PRESSION FOR THE TRANSMISSION
C
C
C-- REFERENCE TEMPERATURES FOR THE PLANCK FUNCTIONS ----- TSTAND
C          TEMPERATURE DEPENDENCE OF THE ABSORPTION ----- TREF
      DATA TSTAND,TREF / 200.0, 200.0 /
C
C-- ROOTS AND WEIGHTS FOR THE 2-POINT GAUSSIAN QUADRATURE
      DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 /
      DATA (WG1(IG1),IG1=1,2) /  1.0        ,  1.0         /

c-----------------------------------------------------------------------
C-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE PLANCK FUNCTIONS
C
       DATA xp/

     s  0.63849788E+01, 0.30969419E+02, 0.44790835E+02,
     s  0.52651048E+01,-0.18799237E+02, 0.92836181E+01,

     s  0.26166790E+02, 0.12348011E+03, 0.17868306E+03,
     s  0.33657659E+02,-0.66869343E+02, 0.21017507E+02,
c
c     s  0.47045285E+02, 0.12234377E+03, 0.61873275E+02,
c     s -0.31971883E+02, 0.59168472E+01, 0.91927407E+01,
c
c     s  0.11101254E+02, 0.86037325E+02, 0.25892695E+03,
c     s  0.35582991E+03, 0.16958020E+03,-0.41311413E+02/
c
     s  0.11101254E+02, 0.86037325E+02, 0.25892695E+03,
     s  0.35582991E+03, 0.16958020E+03,-0.41311413E+02,

     s  0.47045285E+02, 0.12234377E+03, 0.61873275E+02,
     s -0.31971883E+02, 0.59168472E+01, 0.91927407E+01/
c
c-----------------------------------------------------------------------
c   temperature dependency of absorber amounts:
c   -------------------------------------------

      DATA at
     s / 0.694E-03, 0.272E-02, 0.275E-02, 0.178E-01/
      DATA bt
     s / 0.328E-05, 0.298E-05,-0.705E-04,-0.163E-04/

c-----------------------------------------------------------------------
C-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE TEMPERATURE DEPENDENCE OF
C     THE ABSORPTION

c     DATA (ga(ie, 1, 1, 9),ie=1,4)
c    s / 0.288231E-04, 0.170794E-01,-0.339714E-01, 0.000000E+00/
c     DATA (gb(ie, 1, 1, 9),ie=1,4)
c    s / 0.288231E-04, 0.145426E-01, 0.543812E+00, 0.100000E+01/

c     DATA (ga(ie, 2, 1, 9),ie=1,4)
c    s / 0.289299E-01, 0.190634E+01, 0.384061E+01, 0.000000E+00/
c     DATA (gb(ie, 2, 1, 9),ie=1,4)
c    s / 0.289299E-01, 0.189485E+01, 0.600363E+01, 0.100000E+01/

      DATA (ga(ie, 1),ie=1,4)
     s / 0.288231E-04, 0.170794E-01,-0.339714E-01, 0.000000E+00/
      DATA (gb(ie, 1),ie=1,4)
     s / 0.288231E-04, 0.145426E-01, 0.543812E+00, 0.100000E+01/

      DATA (ga(ie, 2),ie=1,4)
     s / 0.289299E-01, 0.190634E+01, 0.384061E+01, 0.000000E+00/
      DATA (gb(ie, 2),ie=1,4)
     s / 0.289299E-01, 0.189485E+01, 0.600363E+01, 0.100000E+01/

      DATA cst_voigt / 0.500E-02, 0.100E-01, 0.150E-01, 0.100E+00/

      NG1=2
      NG1P1=NG1+1
      Nnir=nir
      NIPD=2
      NIPD2=2*NIPD
      NTR=11
      NNTRA=NTRA
      NNUA=NUA
      npade=npademx
      nt_pade=nt_pademx
      nabs=nabsmx
      dt_pade=10.
      tmin_pade=125.

      return
      end