*=*=*=*= SULW.html =*=*=*=*
SUBROUTINE SULW C #include "dimensions.h" #include "dimphys.h" #include "dimradmars.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 c cinterval 1 DATA (ga(ie, 1, 1, 1),ie=1,4) s / 0.278355E-04, 0.164959E-01,-0.329261E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 1),ie=1,4) s / 0.278355E-04, 0.140436E-01, 0.525423E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 2),ie=1,4) s / 0.280173E-04, 0.166060E-01,-0.331191E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 2),ie=1,4) s / 0.280173E-04, 0.141373E-01, 0.529005E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 3),ie=1,4) s / 0.281289E-04, 0.166740E-01,-0.332372E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 3),ie=1,4) s / 0.281289E-04, 0.141950E-01, 0.531254E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 4),ie=1,4) s / 0.282457E-04, 0.167437E-01,-0.333606E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 4),ie=1,4) s / 0.282457E-04, 0.142545E-01, 0.533494E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 5),ie=1,4) s / 0.283684E-04, 0.168159E-01,-0.334902E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 5),ie=1,4) s / 0.283684E-04, 0.143163E-01, 0.535762E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 6),ie=1,4) s / 0.284905E-04, 0.168872E-01,-0.336193E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 6),ie=1,4) s / 0.284905E-04, 0.143774E-01, 0.537967E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 7),ie=1,4) s / 0.286079E-04, 0.169553E-01,-0.337435E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 7),ie=1,4) s / 0.286079E-04, 0.144359E-01, 0.540052E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 8),ie=1,4) s / 0.287189E-04, 0.170194E-01,-0.338610E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 8),ie=1,4) s / 0.287189E-04, 0.144910E-01, 0.542000E+00, 0.100000E+01/ DATA (ga(ie, 1, 1, 9),ie=1,4) s / 0.288231E-04, 0.170794E-01,-0.339714E-01, 0.000000E+00/ DATA (gb(ie, 1, 1, 9),ie=1,4) s / 0.288231E-04, 0.145426E-01, 0.543812E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,10),ie=1,4) s / 0.289207E-04, 0.171354E-01,-0.340749E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,10),ie=1,4) s / 0.289207E-04, 0.145910E-01, 0.545497E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,11),ie=1,4) s / 0.290123E-04, 0.171879E-01,-0.341720E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,11),ie=1,4) s / 0.290123E-04, 0.146362E-01, 0.547067E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,12),ie=1,4) s / 0.290984E-04, 0.172372E-01,-0.342634E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,12),ie=1,4) s / 0.290984E-04, 0.146787E-01, 0.548535E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,13),ie=1,4) s / 0.291796E-04, 0.172835E-01,-0.343495E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,13),ie=1,4) s / 0.291796E-04, 0.147187E-01, 0.549910E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,14),ie=1,4) s / 0.292561E-04, 0.173272E-01,-0.344307E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,14),ie=1,4) s / 0.292561E-04, 0.147564E-01, 0.551202E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,15),ie=1,4) s / 0.293283E-04, 0.173683E-01,-0.345075E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,15),ie=1,4) s / 0.293283E-04, 0.147920E-01, 0.552416E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,16),ie=1,4) s / 0.293965E-04, 0.174071E-01,-0.345799E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,16),ie=1,4) s / 0.293965E-04, 0.148255E-01, 0.553558E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,17),ie=1,4) s / 0.294606E-04, 0.174436E-01,-0.346481E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,17),ie=1,4) s / 0.294606E-04, 0.148570E-01, 0.554630E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,18),ie=1,4) s / 0.295209E-04, 0.174778E-01,-0.347121E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,18),ie=1,4) s / 0.295209E-04, 0.148866E-01, 0.555634E+00, 0.100000E+01/ DATA (ga(ie, 1, 1,19),ie=1,4) s / 0.295771E-04, 0.175098E-01,-0.347719E-01, 0.000000E+00/ DATA (gb(ie, 1, 1,19),ie=1,4) s / 0.295771E-04, 0.149142E-01, 0.556569E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 1),ie=1,4) s / 0.280566E-04, 0.166337E-01,-0.331615E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 1),ie=1,4) s / 0.280566E-04, 0.141601E-01, 0.530085E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 2),ie=1,4) s / 0.282097E-04, 0.167246E-01,-0.333227E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 2),ie=1,4) s / 0.282097E-04, 0.142376E-01, 0.532997E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 3),ie=1,4) s / 0.283806E-04, 0.168248E-01,-0.335030E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 3),ie=1,4) s / 0.283806E-04, 0.143235E-01, 0.536136E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 4),ie=1,4) s / 0.285475E-04, 0.169219E-01,-0.336794E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 4),ie=1,4) s / 0.285475E-04, 0.144067E-01, 0.539123E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 5),ie=1,4) s / 0.287052E-04, 0.170129E-01,-0.338462E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 5),ie=1,4) s / 0.287052E-04, 0.144851E-01, 0.541893E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 6),ie=1,4) s / 0.288529E-04, 0.170979E-01,-0.340028E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 6),ie=1,4) s / 0.288529E-04, 0.145583E-01, 0.544452E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 7),ie=1,4) s / 0.289912E-04, 0.171771E-01,-0.341494E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 7),ie=1,4) s / 0.289912E-04, 0.146266E-01, 0.546821E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 8),ie=1,4) s / 0.291206E-04, 0.172510E-01,-0.342866E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 8),ie=1,4) s / 0.291206E-04, 0.146904E-01, 0.549018E+00, 0.100000E+01/ DATA (ga(ie, 1, 2, 9),ie=1,4) s / 0.292415E-04, 0.173200E-01,-0.344150E-01, 0.000000E+00/ DATA (gb(ie, 1, 2, 9),ie=1,4) s / 0.292415E-04, 0.147499E-01, 0.551058E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,10),ie=1,4) s / 0.293546E-04, 0.173843E-01,-0.345351E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,10),ie=1,4) s / 0.293546E-04, 0.148055E-01, 0.552954E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,11),ie=1,4) s / 0.294601E-04, 0.174443E-01,-0.346472E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,11),ie=1,4) s / 0.294601E-04, 0.148574E-01, 0.554715E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,12),ie=1,4) s / 0.295585E-04, 0.175001E-01,-0.347518E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,12),ie=1,4) s / 0.295585E-04, 0.149057E-01, 0.556350E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,13),ie=1,4) s / 0.296500E-04, 0.175520E-01,-0.348491E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,13),ie=1,4) s / 0.296500E-04, 0.149505E-01, 0.557864E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,14),ie=1,4) s / 0.297350E-04, 0.176001E-01,-0.349395E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,14),ie=1,4) s / 0.297350E-04, 0.149922E-01, 0.559266E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,15),ie=1,4) s / 0.298138E-04, 0.176447E-01,-0.350234E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,15),ie=1,4) s / 0.298138E-04, 0.150308E-01, 0.560562E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,16),ie=1,4) s / 0.298869E-04, 0.176861E-01,-0.351013E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,16),ie=1,4) s / 0.298869E-04, 0.150666E-01, 0.561761E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,17),ie=1,4) s / 0.299549E-04, 0.177245E-01,-0.351736E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,17),ie=1,4) s / 0.299549E-04, 0.150999E-01, 0.562873E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,18),ie=1,4) s / 0.300184E-04, 0.177603E-01,-0.352412E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,18),ie=1,4) s / 0.300184E-04, 0.151309E-01, 0.563910E+00, 0.100000E+01/ DATA (ga(ie, 1, 2,19),ie=1,4) s / 0.300781E-04, 0.177941E-01,-0.353049E-01, 0.000000E+00/ DATA (gb(ie, 1, 2,19),ie=1,4) s / 0.300781E-04, 0.151602E-01, 0.564885E+00, 0.100000E+01/ DATA (ga(ie, 2, 1, 1),ie=1,4) s / 0.345270E-01, 0.224888E+01, 0.463477E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 1),ie=1,4) s / 0.345270E-01, 0.223719E+01, 0.682978E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 2),ie=1,4) s / 0.329361E-01, 0.214299E+01, 0.439620E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 2),ie=1,4) s / 0.329361E-01, 0.213139E+01, 0.657806E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 3),ie=1,4) s / 0.317652E-01, 0.206364E+01, 0.421747E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 3),ie=1,4) s / 0.317652E-01, 0.205213E+01, 0.638583E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 4),ie=1,4) s / 0.308749E-01, 0.200626E+01, 0.408637E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 4),ie=1,4) s / 0.308749E-01, 0.199480E+01, 0.624566E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 5),ie=1,4) s / 0.302088E-01, 0.196625E+01, 0.399303E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 5),ie=1,4) s / 0.302088E-01, 0.195482E+01, 0.614767E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 6),ie=1,4) s / 0.297183E-01, 0.193928E+01, 0.392826E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 6),ie=1,4) s / 0.297183E-01, 0.192786E+01, 0.608177E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 7),ie=1,4) s / 0.293620E-01, 0.192189E+01, 0.388469E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 7),ie=1,4) s / 0.293620E-01, 0.191045E+01, 0.603967E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 8),ie=1,4) s / 0.291074E-01, 0.191152E+01, 0.385680E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 8),ie=1,4) s / 0.291074E-01, 0.190006E+01, 0.601513E+01, 0.100000E+01/ DATA (ga(ie, 2, 1, 9),ie=1,4) s / 0.289299E-01, 0.190634E+01, 0.384061E+01, 0.000000E+00/ DATA (gb(ie, 2, 1, 9),ie=1,4) s / 0.289299E-01, 0.189485E+01, 0.600363E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,10),ie=1,4) s / 0.288113E-01, 0.190504E+01, 0.383322E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,10),ie=1,4) s / 0.288113E-01, 0.189351E+01, 0.600189E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,11),ie=1,4) s / 0.287379E-01, 0.190665E+01, 0.383252E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,11),ie=1,4) s / 0.287379E-01, 0.189508E+01, 0.600750E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,12),ie=1,4) s / 0.286996E-01, 0.191048E+01, 0.383694E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,12),ie=1,4) s / 0.286996E-01, 0.189886E+01, 0.601870E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,13),ie=1,4) s / 0.286884E-01, 0.191598E+01, 0.384530E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,13),ie=1,4) s / 0.286884E-01, 0.190432E+01, 0.603411E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,14),ie=1,4) s / 0.286983E-01, 0.192274E+01, 0.385665E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,14),ie=1,4) s / 0.286983E-01, 0.191103E+01, 0.605268E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,15),ie=1,4) s / 0.287244E-01, 0.193043E+01, 0.387026E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,15),ie=1,4) s / 0.287244E-01, 0.191867E+01, 0.607356E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,16),ie=1,4) s / 0.287629E-01, 0.193876E+01, 0.388551E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,16),ie=1,4) s / 0.287629E-01, 0.192696E+01, 0.609604E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,17),ie=1,4) s / 0.288104E-01, 0.194751E+01, 0.390186E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,17),ie=1,4) s / 0.288104E-01, 0.193566E+01, 0.611949E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,18),ie=1,4) s / 0.288641E-01, 0.195645E+01, 0.391884E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,18),ie=1,4) s / 0.288641E-01, 0.194456E+01, 0.614339E+01, 0.100000E+01/ DATA (ga(ie, 2, 1,19),ie=1,4) s / 0.289215E-01, 0.196539E+01, 0.393603E+01, 0.000000E+00/ DATA (gb(ie, 2, 1,19),ie=1,4) s / 0.289215E-01, 0.195346E+01, 0.616723E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 1),ie=1,4) s / 0.315461E-01, 0.204154E+01, 0.417134E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 1),ie=1,4) s / 0.315461E-01, 0.203010E+01, 0.632845E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 2),ie=1,4) s / 0.304413E-01, 0.197136E+01, 0.401057E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 2),ie=1,4) s / 0.304413E-01, 0.195999E+01, 0.615670E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 3),ie=1,4) s / 0.295810E-01, 0.192066E+01, 0.389180E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 3),ie=1,4) s / 0.295810E-01, 0.190932E+01, 0.603223E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 4),ie=1,4) s / 0.290035E-01, 0.189058E+01, 0.381829E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 4),ie=1,4) s / 0.290035E-01, 0.187923E+01, 0.595874E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 5),ie=1,4) s / 0.286421E-01, 0.187574E+01, 0.377854E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 5),ie=1,4) s / 0.286421E-01, 0.186436E+01, 0.592347E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 6),ie=1,4) s / 0.284332E-01, 0.187142E+01, 0.376225E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 6),ie=1,4) s / 0.284332E-01, 0.185999E+01, 0.591477E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 7),ie=1,4) s / 0.283304E-01, 0.187427E+01, 0.376212E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 7),ie=1,4) s / 0.283304E-01, 0.186277E+01, 0.592427E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 8),ie=1,4) s / 0.283015E-01, 0.188200E+01, 0.377314E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 8),ie=1,4) s / 0.283015E-01, 0.187043E+01, 0.594618E+01, 0.100000E+01/ DATA (ga(ie, 2, 2, 9),ie=1,4) s / 0.283247E-01, 0.189305E+01, 0.379189E+01, 0.000000E+00/ DATA (gb(ie, 2, 2, 9),ie=1,4) s / 0.283247E-01, 0.188142E+01, 0.597652E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,10),ie=1,4) s / 0.283846E-01, 0.190635E+01, 0.381599E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,10),ie=1,4) s / 0.283846E-01, 0.189464E+01, 0.601251E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,11),ie=1,4) s / 0.284703E-01, 0.192113E+01, 0.384373E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,11),ie=1,4) s / 0.284703E-01, 0.190934E+01, 0.605216E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,12),ie=1,4) s / 0.285741E-01, 0.193682E+01, 0.387387E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,12),ie=1,4) s / 0.285741E-01, 0.192496E+01, 0.609403E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,13),ie=1,4) s / 0.286900E-01, 0.195301E+01, 0.390549E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,13),ie=1,4) s / 0.286900E-01, 0.194108E+01, 0.613705E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,14),ie=1,4) s / 0.288136E-01, 0.196939E+01, 0.393785E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,14),ie=1,4) s / 0.288136E-01, 0.195739E+01, 0.618041E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,15),ie=1,4) s / 0.289414E-01, 0.198571E+01, 0.397040E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,15),ie=1,4) s / 0.289414E-01, 0.197363E+01, 0.622346E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,16),ie=1,4) s / 0.290705E-01, 0.200176E+01, 0.400267E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,16),ie=1,4) s / 0.290705E-01, 0.198962E+01, 0.626571E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,17),ie=1,4) s / 0.291985E-01, 0.201738E+01, 0.403426E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,17),ie=1,4) s / 0.291985E-01, 0.200519E+01, 0.630675E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,18),ie=1,4) s / 0.293230E-01, 0.203245E+01, 0.406486E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,18),ie=1,4) s / 0.293230E-01, 0.202020E+01, 0.634625E+01, 0.100000E+01/ DATA (ga(ie, 2, 2,19),ie=1,4) s / 0.294419E-01, 0.204684E+01, 0.409416E+01, 0.000000E+00/ DATA (gb(ie, 2, 2,19),ie=1,4) s / 0.294419E-01, 0.203453E+01, 0.638394E+01, 0.100000E+01/ DATA cst_voigt / 0.500E-02, 0.100E-01, 0.150E-01, 0.100E+00/ C 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. C RETURN END