| 1 | SUBROUTINE SULW |
|---|
| 2 | C |
|---|
| 3 | #include "dimensions.h" |
|---|
| 4 | #include "dimphys.h" |
|---|
| 5 | #include "dimradmars.h" |
|---|
| 6 | #include "callkeys.h" |
|---|
| 7 | #include "yomlw.h" |
|---|
| 8 | C |
|---|
| 9 | C PURPOSE. |
|---|
| 10 | C -------- |
|---|
| 11 | C INITIALIZE YOMLW, THE COMMON THAT CONTAINS COEFFICIENTS |
|---|
| 12 | C NEEDED TO RUN THE LONGWAVE RADIATION ROUTINES |
|---|
| 13 | C |
|---|
| 14 | C AUTHOR. |
|---|
| 15 | C ------- |
|---|
| 16 | C JEAN-JACQUES MORCRETTE *ECMWF* |
|---|
| 17 | c Simplified F. Forget 01/2000 |
|---|
| 18 | C |
|---|
| 19 | C ------------------------------------------------------------------ |
|---|
| 20 | C ---------------------------------------------------------------- |
|---|
| 21 | C |
|---|
| 22 | C * 1. SET VALUES. |
|---|
| 23 | C ----------- |
|---|
| 24 | C |
|---|
| 25 | C |
|---|
| 26 | C-- REFERENCE TEMPERATURES FOR THE PLANCK FUNCTIONS ----- TSTAND |
|---|
| 27 | C TEMPERATURE DEPENDENCE OF THE ABSORPTION ----- TREF |
|---|
| 28 | DATA TSTAND,TREF / 200.0, 200.0 / |
|---|
| 29 | C |
|---|
| 30 | C-- ROOTS AND WEIGHTS FOR THE 2-POINT GAUSSIAN QUADRATURE |
|---|
| 31 | c DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 / |
|---|
| 32 | c DATA (WG1(IG1),IG1=1,2) / 1.0 , 1.0 / |
|---|
| 33 | |
|---|
| 34 | c----------------------------------------------------------------------- |
|---|
| 35 | C-- COEFFICIENTS OF THE POLYNOMIALS GIVING THE PLANCK FUNCTIONS |
|---|
| 36 | C |
|---|
| 37 | DATA xp/ |
|---|
| 38 | |
|---|
| 39 | s 0.63849788E+01, 0.30969419E+02, 0.44790835E+02, |
|---|
| 40 | s 0.52651048E+01,-0.18799237E+02, 0.92836181E+01, |
|---|
| 41 | |
|---|
| 42 | s 0.26166790E+02, 0.12348011E+03, 0.17868306E+03, |
|---|
| 43 | s 0.33657659E+02,-0.66869343E+02, 0.21017507E+02, |
|---|
| 44 | c |
|---|
| 45 | c s 0.47045285E+02, 0.12234377E+03, 0.61873275E+02, |
|---|
| 46 | c s -0.31971883E+02, 0.59168472E+01, 0.91927407E+01, |
|---|
| 47 | c |
|---|
| 48 | c s 0.11101254E+02, 0.86037325E+02, 0.25892695E+03, |
|---|
| 49 | c s 0.35582991E+03, 0.16958020E+03,-0.41311413E+02/ |
|---|
| 50 | c |
|---|
| 51 | s 0.11101254E+02, 0.86037325E+02, 0.25892695E+03, |
|---|
| 52 | s 0.35582991E+03, 0.16958020E+03,-0.41311413E+02, |
|---|
| 53 | |
|---|
| 54 | s 0.47045285E+02, 0.12234377E+03, 0.61873275E+02, |
|---|
| 55 | s -0.31971883E+02, 0.59168472E+01, 0.91927407E+01/ |
|---|
| 56 | |
|---|
| 57 | c |
|---|
| 58 | c----------------------------------------------------------------------- |
|---|
| 59 | c temperature dependency of absorber amounts: |
|---|
| 60 | c ------------------------------------------- |
|---|
| 61 | |
|---|
| 62 | DATA at |
|---|
| 63 | s / 0.694E-03, 0.272E-02, 0.275E-02, 0.178E-01/ |
|---|
| 64 | DATA bt |
|---|
| 65 | s / 0.328E-05, 0.298E-05,-0.705E-04,-0.163E-04/ |
|---|
| 66 | |
|---|
| 67 | c----------------------------------------------------------------------- |
|---|
| 68 | DATA (ga(ie, 1),ie=1,4) |
|---|
| 69 | s / 0.288231E-04, 0.170794E-01,-0.339714E-01, 0.000000E+00/ |
|---|
| 70 | DATA (gb(ie, 1),ie=1,4) |
|---|
| 71 | s / 0.288231E-04, 0.145426E-01, 0.543812E+00, 0.100000E+01/ |
|---|
| 72 | |
|---|
| 73 | DATA (ga(ie, 2),ie=1,4) |
|---|
| 74 | s / 0.289299E-01, 0.190634E+01, 0.384061E+01, 0.000000E+00/ |
|---|
| 75 | DATA (gb(ie, 2),ie=1,4) |
|---|
| 76 | s / 0.289299E-01, 0.189485E+01, 0.600363E+01, 0.100000E+01/ |
|---|
| 77 | |
|---|
| 78 | DATA cst_voigt / 0.500E-02, 0.100E-01, 0.150E-01, 0.100E+00/ |
|---|
| 79 | |
|---|
| 80 | return |
|---|
| 81 | end |
|---|
