!    SUBROUTINE QSAT_WAT-----------------------------------------------
!  
!    SUITABLE FOR SINGLE COLUMN MODEL USE
!  
!    ORIGINAL VERSION WRITTEN FOR CYBER/ETA-10 BY M.FISHER DURING
!    SPRING AND SUMMER 1989
!    CODE REWORKED FOR CYBER/ETA10 BY D.GREGORY SPRING/SUMMER 1989
!  
!    CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90
!  
!    CODE REVIEWED BY C.WILSON   07/09/90
!  
!    VERSION NO. 1 DATED 31/1/90
!  
!    COSMOS DATA SET 'M20.UMCON(QSAT)'
!  
!    PROGRAMMING STANDARDS : UNIFIED MODEL DOCUMENTATION PAPER NO. 4
!    VERSION NO. 1
!  
!    SYSTEM TASK : P27
!  
!    PURPOSE : RETURNS A SATURATION MIXING RATIO GIVEN
!              A TEMPERATURE AND PRESSURE USING SATURATION
!              VAPOUR PRESSURES CALCULATED USING THE
!              GOFF-GRATCH FORMULAE, ADOPTED BY THE WMO AS
!              TAKEN FROM LANDOLT-BORNSTEIN, 1987 NUMERICAL
!              DATA AND FUNCTIONAL RELATIONSHIPS IN SCIENCE
!              AND TECHNOLOGY. GROUP V/VOL 4B METEOROLOGY.
!              PHYSICAL AND CHEMICAL PROPERTIES OF AIR, P35
!  
!              VALUES IN THE LOOKUP TABLE ARE OVER WATER ABOVE AND
!              BELOW 0 DEG C.
!  
!              NOTE : FOR VAPOUR PRESSURE OVER WATER THIS FORMULA
!                     IS VALID FOR TEMPERATURES BETWEEN 373K AND
!                     223K.  THE VALUES FOR SATURATED VAPOUR
!                     OVER WATER IN THE LOOKUP TABLE BELOW ARE
!                     OUT OF THE LOWER END OF THIS RANGE.  HOWEVER
!                     IT IS STANDARD WMO PRACTICE TO USE THE FORMULA
!                     BELOW ITS ACCEPTED RANGE FOR USE WITH THE
!                     CALCULATION OF DEW POINTS IN THE UPPER
!                     ATMOSPHERE
!  
!    DOCUMENTATION : UNIFIED MODEL DOCUMENTATION PAPER
!  
!  END-----------------------------------------------------------------
!
!    ARGUMENTS---------------------------------------------------------
!
      SUBROUTINE QSAT_WAT (QS,T,P,NPNTS)
!
      IMPLICIT NONE
!
!----------------------------------------------------------------------
! MODEL CONSTANTS
!----------------------------------------------------------------------
!
!  ------------------COMDECK C_EPSLON-----------------------------------
! EPSILON IS RATIO OF MOLECULAR WEIGHTS OF WATER AND DRY AIR
      REAL*8 EPSILON, C_VIRTUAL

      PARAMETER(EPSILON=0.62198D+00,
     &          C_VIRTUAL=1.0D+00/EPSILON-1.0D+00)
!*----------------------------------------------------------------------

!
!----------------------------------------------------------------------
! LOCAL CONSTANTS
!----------------------------------------------------------------------
!
      REAL*8 T_LOW           ! LOWEST TEMPERATURE FOR WHICH LOOK-UP
                           ! TABLE OF SATURATION WATER VAPOUR
                           ! PRESSURE IS VALID (K)
!
      REAL*8 T_HIGH          ! HIGHEST TEMPERATURE FOR WHICH LOOK-UP
                           ! TABLE OF SATURATION WATER VAPOUR
                           ! PRESSURES IS VALID (K)
!
      REAL*8 DELTA_T         ! TEMPERATURE INCREMENT OF THE LOOK-UP
                           ! TABLE OF SATURATION VAPOUR PRESSURES
!
      INTEGER N            ! SIZE OF LOOK-UP TABLE OF SATURATION
                           ! WATER VAPOUR PRESSURES
!
      PARAMETER ( T_LOW = 183.150D+00,
     *            T_HIGH = 338.150D+00,
     *            DELTA_T = 0.10D+00,
     *            N = ((T_HIGH - T_LOW + (DELTA_T*0.5))/DELTA_T) +
     *            1.00D+00
     *          )
!
!
!----------------------------------------------------------------------
! VECTOR LENGTHS AND LOOP COUNTERS
!----------------------------------------------------------------------
!
      INTEGER NPNTS        ! VECTOR LENGTH
!
      INTEGER I            ! LOOP COUNTER
!
      INTEGER IES          ! LOOP COUNTER FOR DATA STATEMENT
                           ! LOOK-UP TABLE
!
!
!----------------------------------------------------------------------
! VARIABLES WHICH ARE INPUT
!----------------------------------------------------------------------
!
      REAL*8 T(NPNTS)        ! IN TEMPERATURE (K)
!
      REAL*8 P(NPNTS)        ! IN PRESSURE (PA)
!
!
!----------------------------------------------------------------------
! VARIABLES WHICH ARE OUTPUT
!----------------------------------------------------------------------
!
      REAL*8 QS(NPNTS)       ! OUT SATURATION MIXING RATIO AT TEMPERATURE
                           !     T AND PRESSURE P (KG/KG)
!
!
!----------------------------------------------------------------------
! VARIABLES WHICH ARE DEFINED LOCALLY
!----------------------------------------------------------------------
!
      REAL*8 ES(N)           ! TABLE OF SATURATION WATER VAPOUR
                           ! PRESSURE (PA) - SET BY DATA STATEMENT
                           ! CALCULATED FROM THE GOFF-GRATCH FORMULAE
                           ! AS TAKEN FROM LANDOLT-BORNSTEIN, 1987
                           ! NUMERICAL DATA AND FUNCTIONAL RELATIONSHIPS
                           ! IN SCIENCE AND TECHNOLOGY. GROUP V/ VOL 4B
                           ! METEOROLOGY. PHYSICAL AND CHEMICAL
                           ! PROPERTIES OF AIR, P35
!
      REAL*8 ATABLE          ! WORK VARIABLES
!
      INTEGER ITABLE       ! WORK VARIABLES
!
!*---------------------------------------------------------------------
!L
!L---------------------------------------------------------------------
!L NO SIGNIFICANT STRUCTURE
!L---------------------------------------------------------------------
!L
!
!---------------------------------------------------------------------
! SATURATION WATER VAPOUR PRESSURE
!
! VALUES ABOVE AND BELOW 0 DEG C ARE OVER WATER
!
! VALUES BELOW -50 DEC C ARE OUTSIDE FORMAL RANGE OF
! WMO ADOPTED FORMULA (SEE ABOVE).  HOWEVER STANDARD
! PRACTICE TO USE THESE VALUES
!---------------------------------------------------------------------
!
      DATA (ES(IES),IES=    1, 95) /
     *0.186905D-01,0.190449D-01,0.194059D-01,0.197727D-01,0.201462D-01,
     *0.205261D-01,0.209122D-01,0.213052D-01,0.217050D-01,0.221116D-01,
     *0.225252D-01,0.229463D-01,0.233740D-01,0.238090D-01,0.242518D-01,
     *0.247017D-01,0.251595D-01,0.256252D-01,0.260981D-01,0.265795D-01,
     *0.270691D-01,0.275667D-01,0.280733D-01,0.285876D-01,0.291105D-01,
     *0.296429D-01,0.301835D-01,0.307336D-01,0.312927D-01,0.318611D-01,
     *0.324390D-01,0.330262D-01,0.336232D-01,0.342306D-01,0.348472D-01,
     *0.354748D-01,0.361117D-01,0.367599D-01,0.374185D-01,0.380879D-01,
     *0.387689D-01,0.394602D-01,0.401626D-01,0.408771D-01,0.416033D-01,
     *0.423411D-01,0.430908D-01,0.438524D-01,0.446263D-01,0.454124D-01,
     *0.462122D-01,0.470247D-01,0.478491D-01,0.486874D-01,0.495393D-01,
     *0.504057D-01,0.512847D-01,0.521784D-01,0.530853D-01,0.540076D-01,
     *0.549444D-01,0.558959D-01,0.568633D-01,0.578448D-01,0.588428D-01,
     *0.598566D-01,0.608858D-01,0.619313D-01,0.629926D-01,0.640706D-01,
     *0.651665D-01,0.662795D-01,0.674095D-01,0.685570D-01,0.697219D-01,
     *0.709063D-01,0.721076D-01,0.733284D-01,0.745679D-01,0.758265D-01,
     *0.771039D-01,0.784026D-01,0.797212D-01,0.810577D-01,0.824164D-01,
     *0.837971D-01,0.851970D-01,0.866198D-01,0.880620D-01,0.895281D-01,
     *0.910178D-01,0.925278D-01,0.940622D-01,0.956177D-01,0.971984D-01/
      DATA (ES(IES),IES= 96,190) /
     *0.988051D-01,0.100433D+00,0.102085D+00,0.103764D+00,0.105467D+00,
     *0.107196D+00,0.108953D+00,0.110732D+00,0.112541D+00,0.114376D+00,
     *0.116238D+00,0.118130D+00,0.120046D+00,0.121993D+00,0.123969D+00,
     *0.125973D+00,0.128009D+00,0.130075D+00,0.132167D+00,0.134296D+00,
     *0.136452D+00,0.138642D+00,0.140861D+00,0.143115D+00,0.145404D+00,
     *0.147723D+00,0.150078D+00,0.152466D+00,0.154889D+00,0.157346D+00,
     *0.159841D+00,0.162372D+00,0.164939D+00,0.167545D+00,0.170185D+00,
     *0.172866D+00,0.175584D+00,0.178340D+00,0.181139D+00,0.183977D+00,
     *0.186855D+00,0.189773D+00,0.192737D+00,0.195736D+00,0.198783D+00,
     *0.201875D+00,0.205007D+00,0.208186D+00,0.211409D+00,0.214676D+00,
     *0.217993D+00,0.221355D+00,0.224764D+00,0.228220D+00,0.231728D+00,
     *0.235284D+00,0.238888D+00,0.242542D+00,0.246251D+00,0.250010D+00,
     *0.253821D+00,0.257688D+00,0.261602D+00,0.265575D+00,0.269607D+00,
     *0.273689D+00,0.277830D+00,0.282027D+00,0.286287D+00,0.290598D+00,
     *0.294972D+00,0.299405D+00,0.303904D+00,0.308462D+00,0.313082D+00,
     *0.317763D+00,0.322512D+00,0.327324D+00,0.332201D+00,0.337141D+00,
     *0.342154D+00,0.347234D+00,0.352387D+00,0.357601D+00,0.362889D+00,
     *0.368257D+00,0.373685D+00,0.379194D+00,0.384773D+00,0.390433D+00,
     *0.396159D+00,0.401968D+00,0.407861D+00,0.413820D+00,0.419866D+00/
      DATA (ES(IES),IES=191,285) /
     *0.425999D+00,0.432203D+00,0.438494D+00,0.444867D+00,0.451332D+00,
     *0.457879D+00,0.464510D+00,0.471226D+00,0.478037D+00,0.484935D+00,
     *0.491920D+00,0.499005D+00,0.506181D+00,0.513447D+00,0.520816D+00,
     *0.528279D+00,0.535835D+00,0.543497D+00,0.551256D+00,0.559113D+00,
     *0.567081D+00,0.575147D+00,0.583315D+00,0.591585D+00,0.599970D+00,
     *0.608472D+00,0.617069D+00,0.625785D+00,0.634609D+00,0.643556D+00,
     *0.652611D+00,0.661782D+00,0.671077D+00,0.680487D+00,0.690015D+00,
     *0.699656D+00,0.709433D+00,0.719344D+00,0.729363D+00,0.739518D+00,
     *0.749795D+00,0.760217D+00,0.770763D+00,0.781454D+00,0.792258D+00,
     *0.803208D+00,0.814309D+00,0.825528D+00,0.836914D+00,0.848422D+00,
     *0.860086D+00,0.871891D+00,0.883837D+00,0.895944D+00,0.908214D+00,
     *0.920611D+00,0.933175D+00,0.945890D+00,0.958776D+00,0.971812D+00,
     *0.985027D+00,0.998379D+00,0.101193D+01,0.102561D+01,0.103949D+01,
     *0.105352D+01,0.106774D+01,0.108213D+01,0.109669D+01,0.111144D+01,
     *0.112636D+01,0.114148D+01,0.115676D+01,0.117226D+01,0.118791D+01,
     *0.120377D+01,0.121984D+01,0.123608D+01,0.125252D+01,0.126919D+01,
     *0.128604D+01,0.130309D+01,0.132036D+01,0.133782D+01,0.135549D+01,
     *0.137339D+01,0.139150D+01,0.140984D+01,0.142839D+01,0.144715D+01,
     *0.146616D+01,0.148538D+01,0.150482D+01,0.152450D+01,0.154445D+01/
      DATA (ES(IES),IES=286,380) /
     *0.156459D+01,0.158502D+01,0.160564D+01,0.162654D+01,0.164766D+01,
     *0.166906D+01,0.169070D+01,0.171257D+01,0.173473D+01,0.175718D+01,
     *0.177984D+01,0.180282D+01,0.182602D+01,0.184951D+01,0.187327D+01,
     *0.189733D+01,0.192165D+01,0.194629D+01,0.197118D+01,0.199636D+01,
     *0.202185D+01,0.204762D+01,0.207372D+01,0.210010D+01,0.212678D+01,
     *0.215379D+01,0.218109D+01,0.220873D+01,0.223668D+01,0.226497D+01,
     *0.229357D+01,0.232249D+01,0.235176D+01,0.238134D+01,0.241129D+01,
     *0.244157D+01,0.247217D+01,0.250316D+01,0.253447D+01,0.256617D+01,
     *0.259821D+01,0.263064D+01,0.266341D+01,0.269661D+01,0.273009D+01,
     *0.276403D+01,0.279834D+01,0.283302D+01,0.286811D+01,0.290358D+01,
     *0.293943D+01,0.297571D+01,0.301236D+01,0.304946D+01,0.308702D+01,
     *0.312491D+01,0.316326D+01,0.320208D+01,0.324130D+01,0.328092D+01,
     *0.332102D+01,0.336162D+01,0.340264D+01,0.344407D+01,0.348601D+01,
     *0.352838D+01,0.357118D+01,0.361449D+01,0.365834D+01,0.370264D+01,
     *0.374737D+01,0.379265D+01,0.383839D+01,0.388469D+01,0.393144D+01,
     *0.397876D+01,0.402656D+01,0.407492D+01,0.412378D+01,0.417313D+01,
     *0.422306D+01,0.427359D+01,0.432454D+01,0.437617D+01,0.442834D+01,
     *0.448102D+01,0.453433D+01,0.458816D+01,0.464253D+01,0.469764D+01,
     *0.475321D+01,0.480942D+01,0.486629D+01,0.492372D+01,0.498173D+01/
      DATA (ES(IES),IES=381,475) /
     *0.504041D+01,0.509967D+01,0.515962D+01,0.522029D+01,0.528142D+01,
     *0.534337D+01,0.540595D+01,0.546912D+01,0.553292D+01,0.559757D+01,
     *0.566273D+01,0.572864D+01,0.579532D+01,0.586266D+01,0.593075D+01,
     *0.599952D+01,0.606895D+01,0.613918D+01,0.621021D+01,0.628191D+01,
     *0.635433D+01,0.642755D+01,0.650162D+01,0.657639D+01,0.665188D+01,
     *0.672823D+01,0.680532D+01,0.688329D+01,0.696198D+01,0.704157D+01,
     *0.712206D+01,0.720319D+01,0.728534D+01,0.736829D+01,0.745204D+01,
     *0.753671D+01,0.762218D+01,0.770860D+01,0.779588D+01,0.788408D+01,
     *0.797314D+01,0.806318D+01,0.815408D+01,0.824599D+01,0.833874D+01,
     *0.843254D+01,0.852721D+01,0.862293D+01,0.871954D+01,0.881724D+01,
     *0.891579D+01,0.901547D+01,0.911624D+01,0.921778D+01,0.932061D+01,
     *0.942438D+01,0.952910D+01,0.963497D+01,0.974181D+01,0.984982D+01,
     *0.995887D+01,0.100690D+02,0.101804D+02,0.102926D+02,0.104063D+02,
     *0.105210D+02,0.106367D+02,0.107536D+02,0.108719D+02,0.109912D+02,
     *0.111116D+02,0.112333D+02,0.113563D+02,0.114804D+02,0.116056D+02,
     *0.117325D+02,0.118602D+02,0.119892D+02,0.121197D+02,0.122513D+02,
     *0.123844D+02,0.125186D+02,0.126543D+02,0.127912D+02,0.129295D+02,
     *0.130691D+02,0.132101D+02,0.133527D+02,0.134965D+02,0.136415D+02,
     *0.137882D+02,0.139361D+02,0.140855D+02,0.142366D+02,0.143889D+02/
      DATA (ES(IES),IES=476,570) /
     *0.145429D+02,0.146982D+02,0.148552D+02,0.150135D+02,0.151735D+02,
     *0.153349D+02,0.154979D+02,0.156624D+02,0.158286D+02,0.159965D+02,
     *0.161659D+02,0.163367D+02,0.165094D+02,0.166838D+02,0.168597D+02,
     *0.170375D+02,0.172168D+02,0.173979D+02,0.175806D+02,0.177651D+02,
     *0.179513D+02,0.181394D+02,0.183293D+02,0.185210D+02,0.187146D+02,
     *0.189098D+02,0.191066D+02,0.193059D+02,0.195065D+02,0.197095D+02,
     *0.199142D+02,0.201206D+02,0.203291D+02,0.205397D+02,0.207522D+02,
     *0.209664D+02,0.211831D+02,0.214013D+02,0.216221D+02,0.218448D+02,
     *0.220692D+02,0.222959D+02,0.225250D+02,0.227559D+02,0.229887D+02,
     *0.232239D+02,0.234614D+02,0.237014D+02,0.239428D+02,0.241872D+02,
     *0.244335D+02,0.246824D+02,0.249332D+02,0.251860D+02,0.254419D+02,
     *0.256993D+02,0.259600D+02,0.262225D+02,0.264873D+02,0.267552D+02,
     *0.270248D+02,0.272970D+02,0.275719D+02,0.278497D+02,0.281295D+02,
     *0.284117D+02,0.286965D+02,0.289843D+02,0.292743D+02,0.295671D+02,
     *0.298624D+02,0.301605D+02,0.304616D+02,0.307650D+02,0.310708D+02,
     *0.313803D+02,0.316915D+02,0.320064D+02,0.323238D+02,0.326437D+02,
     *0.329666D+02,0.332928D+02,0.336215D+02,0.339534D+02,0.342885D+02,
     *0.346263D+02,0.349666D+02,0.353109D+02,0.356572D+02,0.360076D+02,
     *0.363606D+02,0.367164D+02,0.370757D+02,0.374383D+02,0.378038D+02/
      DATA (ES(IES),IES=571,665) /
     *0.381727D+02,0.385453D+02,0.389206D+02,0.392989D+02,0.396807D+02,
     *0.400663D+02,0.404555D+02,0.408478D+02,0.412428D+02,0.416417D+02,
     *0.420445D+02,0.424502D+02,0.428600D+02,0.432733D+02,0.436900D+02,
     *0.441106D+02,0.445343D+02,0.449620D+02,0.453930D+02,0.458280D+02,
     *0.462672D+02,0.467096D+02,0.471561D+02,0.476070D+02,0.480610D+02,
     *0.485186D+02,0.489813D+02,0.494474D+02,0.499170D+02,0.503909D+02,
     *0.508693D+02,0.513511D+02,0.518376D+02,0.523277D+02,0.528232D+02,
     *0.533213D+02,0.538240D+02,0.543315D+02,0.548437D+02,0.553596D+02,
     *0.558802D+02,0.564046D+02,0.569340D+02,0.574672D+02,0.580061D+02,
     *0.585481D+02,0.590963D+02,0.596482D+02,0.602041D+02,0.607649D+02,
     *0.613311D+02,0.619025D+02,0.624779D+02,0.630574D+02,0.636422D+02,
     *0.642324D+02,0.648280D+02,0.654278D+02,0.660332D+02,0.666426D+02,
     *0.672577D+02,0.678771D+02,0.685034D+02,0.691328D+02,0.697694D+02,
     *0.704103D+02,0.710556D+02,0.717081D+02,0.723639D+02,0.730269D+02,
     *0.736945D+02,0.743681D+02,0.750463D+02,0.757309D+02,0.764214D+02,
     *0.771167D+02,0.778182D+02,0.785246D+02,0.792373D+02,0.799564D+02,
     *0.806804D+02,0.814109D+02,0.821479D+02,0.828898D+02,0.836384D+02,
     *0.843922D+02,0.851525D+02,0.859198D+02,0.866920D+02,0.874712D+02,
     *0.882574D+02,0.890486D+02,0.898470D+02,0.906525D+02,0.914634D+02/
      DATA (ES(IES),IES=666,760) /
     *0.922814D+02,0.931048D+02,0.939356D+02,0.947736D+02,0.956171D+02,
     *0.964681D+02,0.973246D+02,0.981907D+02,0.990605D+02,0.999399D+02,
     *0.100825D+03,0.101718D+03,0.102617D+03,0.103523D+03,0.104438D+03,
     *0.105358D+03,0.106287D+03,0.107221D+03,0.108166D+03,0.109115D+03,
     *0.110074D+03,0.111039D+03,0.112012D+03,0.112992D+03,0.113981D+03,
     *0.114978D+03,0.115981D+03,0.116993D+03,0.118013D+03,0.119041D+03,
     *0.120077D+03,0.121122D+03,0.122173D+03,0.123234D+03,0.124301D+03,
     *0.125377D+03,0.126463D+03,0.127556D+03,0.128657D+03,0.129769D+03,
     *0.130889D+03,0.132017D+03,0.133152D+03,0.134299D+03,0.135453D+03,
     *0.136614D+03,0.137786D+03,0.138967D+03,0.140158D+03,0.141356D+03,
     *0.142565D+03,0.143781D+03,0.145010D+03,0.146247D+03,0.147491D+03,
     *0.148746D+03,0.150011D+03,0.151284D+03,0.152571D+03,0.153862D+03,
     *0.155168D+03,0.156481D+03,0.157805D+03,0.159137D+03,0.160478D+03,
     *0.161832D+03,0.163198D+03,0.164569D+03,0.165958D+03,0.167348D+03,
     *0.168757D+03,0.170174D+03,0.171599D+03,0.173037D+03,0.174483D+03,
     *0.175944D+03,0.177414D+03,0.178892D+03,0.180387D+03,0.181886D+03,
     *0.183402D+03,0.184930D+03,0.186463D+03,0.188012D+03,0.189571D+03,
     *0.191146D+03,0.192730D+03,0.194320D+03,0.195930D+03,0.197546D+03,
     *0.199175D+03,0.200821D+03,0.202473D+03,0.204142D+03,0.205817D+03/
      DATA (ES(IES),IES=761,855) /
     *0.207510D+03,0.209216D+03,0.210928D+03,0.212658D+03,0.214398D+03,
     *0.216152D+03,0.217920D+03,0.219698D+03,0.221495D+03,0.223297D+03,
     *0.225119D+03,0.226951D+03,0.228793D+03,0.230654D+03,0.232522D+03,
     *0.234413D+03,0.236311D+03,0.238223D+03,0.240151D+03,0.242090D+03,
     *0.244049D+03,0.246019D+03,0.248000D+03,0.249996D+03,0.252009D+03,
     *0.254037D+03,0.256077D+03,0.258128D+03,0.260200D+03,0.262284D+03,
     *0.264384D+03,0.266500D+03,0.268629D+03,0.270779D+03,0.272936D+03,
     *0.275110D+03,0.277306D+03,0.279509D+03,0.281734D+03,0.283966D+03,
     *0.286227D+03,0.288494D+03,0.290780D+03,0.293083D+03,0.295398D+03,
     *0.297737D+03,0.300089D+03,0.302453D+03,0.304841D+03,0.307237D+03,
     *0.309656D+03,0.312095D+03,0.314541D+03,0.317012D+03,0.319496D+03,
     *0.322005D+03,0.324527D+03,0.327063D+03,0.329618D+03,0.332193D+03,
     *0.334788D+03,0.337396D+03,0.340025D+03,0.342673D+03,0.345329D+03,
     *0.348019D+03,0.350722D+03,0.353440D+03,0.356178D+03,0.358938D+03,
     *0.361718D+03,0.364513D+03,0.367322D+03,0.370160D+03,0.373012D+03,
     *0.375885D+03,0.378788D+03,0.381691D+03,0.384631D+03,0.387579D+03,
     *0.390556D+03,0.393556D+03,0.396563D+03,0.399601D+03,0.402646D+03,
     *0.405730D+03,0.408829D+03,0.411944D+03,0.415083D+03,0.418236D+03,
     *0.421422D+03,0.424632D+03,0.427849D+03,0.431099D+03,0.434365D+03/
      DATA (ES(IES),IES=856,950) /
     *0.437655D+03,0.440970D+03,0.444301D+03,0.447666D+03,0.451038D+03,
     *0.454445D+03,0.457876D+03,0.461316D+03,0.464790D+03,0.468281D+03,
     *0.471798D+03,0.475342D+03,0.478902D+03,0.482497D+03,0.486101D+03,
     *0.489741D+03,0.493408D+03,0.497083D+03,0.500804D+03,0.504524D+03,
     *0.508290D+03,0.512074D+03,0.515877D+03,0.519717D+03,0.523566D+03,
     *0.527462D+03,0.531367D+03,0.535301D+03,0.539264D+03,0.543245D+03,
     *0.547265D+03,0.551305D+03,0.555363D+03,0.559462D+03,0.563579D+03,
     *0.567727D+03,0.571905D+03,0.576102D+03,0.580329D+03,0.584576D+03,
     *0.588865D+03,0.593185D+03,0.597514D+03,0.601885D+03,0.606276D+03,
     *0.610699D+03,0.615151D+03,0.619625D+03,0.624140D+03,0.628671D+03,
     *0.633243D+03,0.637845D+03,0.642465D+03,0.647126D+03,0.651806D+03,
     *0.656527D+03,0.661279D+03,0.666049D+03,0.670861D+03,0.675692D+03,
     *0.680566D+03,0.685471D+03,0.690396D+03,0.695363D+03,0.700350D+03,
     *0.705381D+03,0.710444D+03,0.715527D+03,0.720654D+03,0.725801D+03,
     *0.730994D+03,0.736219D+03,0.741465D+03,0.746756D+03,0.752068D+03,
     *0.757426D+03,0.762819D+03,0.768231D+03,0.773692D+03,0.779172D+03,
     *0.784701D+03,0.790265D+03,0.795849D+03,0.801483D+03,0.807137D+03,
     *0.812842D+03,0.818582D+03,0.824343D+03,0.830153D+03,0.835987D+03,
     *0.841871D+03,0.847791D+03,0.853733D+03,0.859727D+03,0.865743D+03/
      DATA (ES(IES),IES=951,1045) /
     *0.871812D+03,0.877918D+03,0.884046D+03,0.890228D+03,0.896433D+03,
     *0.902690D+03,0.908987D+03,0.915307D+03,0.921681D+03,0.928078D+03,
     *0.934531D+03,0.941023D+03,0.947539D+03,0.954112D+03,0.960708D+03,
     *0.967361D+03,0.974053D+03,0.980771D+03,0.987545D+03,0.994345D+03,
     *0.100120D+04,0.100810D+04,0.101502D+04,0.102201D+04,0.102902D+04,
     *0.103608D+04,0.104320D+04,0.105033D+04,0.105753D+04,0.106475D+04,
     *0.107204D+04,0.107936D+04,0.108672D+04,0.109414D+04,0.110158D+04,
     *0.110908D+04,0.111663D+04,0.112421D+04,0.113185D+04,0.113952D+04,
     *0.114725D+04,0.115503D+04,0.116284D+04,0.117071D+04,0.117861D+04,
     *0.118658D+04,0.119459D+04,0.120264D+04,0.121074D+04,0.121888D+04,
     *0.122709D+04,0.123534D+04,0.124362D+04,0.125198D+04,0.126036D+04,
     *0.126881D+04,0.127731D+04,0.128584D+04,0.129444D+04,0.130307D+04,
     *0.131177D+04,0.132053D+04,0.132931D+04,0.133817D+04,0.134705D+04,
     *0.135602D+04,0.136503D+04,0.137407D+04,0.138319D+04,0.139234D+04,
     *0.140156D+04,0.141084D+04,0.142015D+04,0.142954D+04,0.143896D+04,
     *0.144845D+04,0.145800D+04,0.146759D+04,0.147725D+04,0.148694D+04,
     *0.149672D+04,0.150655D+04,0.151641D+04,0.152635D+04,0.153633D+04,
     *0.154639D+04,0.155650D+04,0.156665D+04,0.157688D+04,0.158715D+04,
     *0.159750D+04,0.160791D+04,0.161836D+04,0.162888D+04,0.163945D+04/
      DATA (ES(IES),IES=1046,1140) /
     *0.165010D+04,0.166081D+04,0.167155D+04,0.168238D+04,0.169325D+04,
     *0.170420D+04,0.171522D+04,0.172627D+04,0.173741D+04,0.174859D+04,
     *0.175986D+04,0.177119D+04,0.178256D+04,0.179402D+04,0.180552D+04,
     *0.181711D+04,0.182877D+04,0.184046D+04,0.185224D+04,0.186407D+04,
     *0.187599D+04,0.188797D+04,0.190000D+04,0.191212D+04,0.192428D+04,
     *0.193653D+04,0.194886D+04,0.196122D+04,0.197368D+04,0.198618D+04,
     *0.199878D+04,0.201145D+04,0.202416D+04,0.203698D+04,0.204983D+04,
     *0.206278D+04,0.207580D+04,0.208887D+04,0.210204D+04,0.211525D+04,
     *0.212856D+04,0.214195D+04,0.215538D+04,0.216892D+04,0.218249D+04,
     *0.219618D+04,0.220994D+04,0.222375D+04,0.223766D+04,0.225161D+04,
     *0.226567D+04,0.227981D+04,0.229399D+04,0.230829D+04,0.232263D+04,
     *0.233708D+04,0.235161D+04,0.236618D+04,0.238087D+04,0.239560D+04,
     *0.241044D+04,0.242538D+04,0.244035D+04,0.245544D+04,0.247057D+04,
     *0.248583D+04,0.250116D+04,0.251654D+04,0.253204D+04,0.254759D+04,
     *0.256325D+04,0.257901D+04,0.259480D+04,0.261073D+04,0.262670D+04,
     *0.264279D+04,0.265896D+04,0.267519D+04,0.269154D+04,0.270794D+04,
     *0.272447D+04,0.274108D+04,0.275774D+04,0.277453D+04,0.279137D+04,
     *0.280834D+04,0.282540D+04,0.284251D+04,0.285975D+04,0.287704D+04,
     *0.289446D+04,0.291198D+04,0.292954D+04,0.294725D+04,0.296499D+04/
      DATA (ES(IES),IES=1141,1235) /
     *0.298288D+04,0.300087D+04,0.301890D+04,0.303707D+04,0.305529D+04,
     *0.307365D+04,0.309211D+04,0.311062D+04,0.312927D+04,0.314798D+04,
     *0.316682D+04,0.318577D+04,0.320477D+04,0.322391D+04,0.324310D+04,
     *0.326245D+04,0.328189D+04,0.330138D+04,0.332103D+04,0.334073D+04,
     *0.336058D+04,0.338053D+04,0.340054D+04,0.342069D+04,0.344090D+04,
     *0.346127D+04,0.348174D+04,0.350227D+04,0.352295D+04,0.354369D+04,
     *0.356458D+04,0.358559D+04,0.360664D+04,0.362787D+04,0.364914D+04,
     *0.367058D+04,0.369212D+04,0.371373D+04,0.373548D+04,0.375731D+04,
     *0.377929D+04,0.380139D+04,0.382355D+04,0.384588D+04,0.386826D+04,
     *0.389081D+04,0.391348D+04,0.393620D+04,0.395910D+04,0.398205D+04,
     *0.400518D+04,0.402843D+04,0.405173D+04,0.407520D+04,0.409875D+04,
     *0.412246D+04,0.414630D+04,0.417019D+04,0.419427D+04,0.421840D+04,
     *0.424272D+04,0.426715D+04,0.429165D+04,0.431634D+04,0.434108D+04,
     *0.436602D+04,0.439107D+04,0.441618D+04,0.444149D+04,0.446685D+04,
     *0.449241D+04,0.451810D+04,0.454385D+04,0.456977D+04,0.459578D+04,
     *0.462197D+04,0.464830D+04,0.467468D+04,0.470127D+04,0.472792D+04,
     *0.475477D+04,0.478175D+04,0.480880D+04,0.483605D+04,0.486336D+04,
     *0.489087D+04,0.491853D+04,0.494623D+04,0.497415D+04,0.500215D+04,
     *0.503034D+04,0.505867D+04,0.508707D+04,0.511568D+04,0.514436D+04/
      DATA (ES(IES),IES=1236,1330) /
     *0.517325D+04,0.520227D+04,0.523137D+04,0.526068D+04,0.529005D+04,
     *0.531965D+04,0.534939D+04,0.537921D+04,0.540923D+04,0.543932D+04,
     *0.546965D+04,0.550011D+04,0.553064D+04,0.556139D+04,0.559223D+04,
     *0.562329D+04,0.565449D+04,0.568577D+04,0.571727D+04,0.574884D+04,
     *0.578064D+04,0.581261D+04,0.584464D+04,0.587692D+04,0.590924D+04,
     *0.594182D+04,0.597455D+04,0.600736D+04,0.604039D+04,0.607350D+04,
     *0.610685D+04,0.614036D+04,0.617394D+04,0.620777D+04,0.624169D+04,
     *0.627584D+04,0.631014D+04,0.634454D+04,0.637918D+04,0.641390D+04,
     *0.644887D+04,0.648400D+04,0.651919D+04,0.655467D+04,0.659021D+04,
     *0.662599D+04,0.666197D+04,0.669800D+04,0.673429D+04,0.677069D+04,
     *0.680735D+04,0.684415D+04,0.688104D+04,0.691819D+04,0.695543D+04,
     *0.699292D+04,0.703061D+04,0.706837D+04,0.710639D+04,0.714451D+04,
     *0.718289D+04,0.722143D+04,0.726009D+04,0.729903D+04,0.733802D+04,
     *0.737729D+04,0.741676D+04,0.745631D+04,0.749612D+04,0.753602D+04,
     *0.757622D+04,0.761659D+04,0.765705D+04,0.769780D+04,0.773863D+04,
     *0.777975D+04,0.782106D+04,0.786246D+04,0.790412D+04,0.794593D+04,
     *0.798802D+04,0.803028D+04,0.807259D+04,0.811525D+04,0.815798D+04,
     *0.820102D+04,0.824427D+04,0.828757D+04,0.833120D+04,0.837493D+04,
     *0.841895D+04,0.846313D+04,0.850744D+04,0.855208D+04,0.859678D+04/
      DATA (ES(IES),IES=1331,1425) /
     *0.864179D+04,0.868705D+04,0.873237D+04,0.877800D+04,0.882374D+04,
     *0.886979D+04,0.891603D+04,0.896237D+04,0.900904D+04,0.905579D+04,
     *0.910288D+04,0.915018D+04,0.919758D+04,0.924529D+04,0.929310D+04,
     *0.934122D+04,0.938959D+04,0.943804D+04,0.948687D+04,0.953575D+04,
     *0.958494D+04,0.963442D+04,0.968395D+04,0.973384D+04,0.978383D+04,
     *0.983412D+04,0.988468D+04,0.993534D+04,0.998630D+04,0.100374D+05,
     *0.100888D+05,0.101406D+05,0.101923D+05,0.102444D+05,0.102966D+05,
     *0.103492D+05,0.104020D+05,0.104550D+05,0.105082D+05,0.105616D+05,
     *0.106153D+05,0.106693D+05,0.107234D+05,0.107779D+05,0.108325D+05,
     *0.108874D+05,0.109425D+05,0.109978D+05,0.110535D+05,0.111092D+05,
     *0.111653D+05,0.112217D+05,0.112782D+05,0.113350D+05,0.113920D+05,
     *0.114493D+05,0.115070D+05,0.115646D+05,0.116228D+05,0.116809D+05,
     *0.117396D+05,0.117984D+05,0.118574D+05,0.119167D+05,0.119762D+05,
     *0.120360D+05,0.120962D+05,0.121564D+05,0.122170D+05,0.122778D+05,
     *0.123389D+05,0.124004D+05,0.124619D+05,0.125238D+05,0.125859D+05,
     *0.126484D+05,0.127111D+05,0.127739D+05,0.128372D+05,0.129006D+05,
     *0.129644D+05,0.130285D+05,0.130927D+05,0.131573D+05,0.132220D+05,
     *0.132872D+05,0.133526D+05,0.134182D+05,0.134842D+05,0.135503D+05,
     *0.136168D+05,0.136836D+05,0.137505D+05,0.138180D+05,0.138854D+05/
      DATA (ES(IES),IES=1426,1520) /
     *0.139534D+05,0.140216D+05,0.140900D+05,0.141588D+05,0.142277D+05,
     *0.142971D+05,0.143668D+05,0.144366D+05,0.145069D+05,0.145773D+05,
     *0.146481D+05,0.147192D+05,0.147905D+05,0.148622D+05,0.149341D+05,
     *0.150064D+05,0.150790D+05,0.151517D+05,0.152250D+05,0.152983D+05,
     *0.153721D+05,0.154462D+05,0.155205D+05,0.155952D+05,0.156701D+05,
     *0.157454D+05,0.158211D+05,0.158969D+05,0.159732D+05,0.160496D+05,
     *0.161265D+05,0.162037D+05,0.162811D+05,0.163589D+05,0.164369D+05,
     *0.165154D+05,0.165942D+05,0.166732D+05,0.167526D+05,0.168322D+05,
     *0.169123D+05,0.169927D+05,0.170733D+05,0.171543D+05,0.172356D+05,
     *0.173173D+05,0.173993D+05,0.174815D+05,0.175643D+05,0.176471D+05,
     *0.177305D+05,0.178143D+05,0.178981D+05,0.179826D+05,0.180671D+05,
     *0.181522D+05,0.182377D+05,0.183232D+05,0.184093D+05,0.184955D+05,
     *0.185823D+05,0.186695D+05,0.187568D+05,0.188447D+05,0.189326D+05,
     *0.190212D+05,0.191101D+05,0.191991D+05,0.192887D+05,0.193785D+05,
     *0.194688D+05,0.195595D+05,0.196503D+05,0.197417D+05,0.198332D+05,
     *0.199253D+05,0.200178D+05,0.201105D+05,0.202036D+05,0.202971D+05,
     *0.203910D+05,0.204853D+05,0.205798D+05,0.206749D+05,0.207701D+05,
     *0.208659D+05,0.209621D+05,0.210584D+05,0.211554D+05,0.212524D+05,
     *0.213501D+05,0.214482D+05,0.215465D+05,0.216452D+05,0.217442D+05/
      DATA (ES(IES),IES=1521,1551) /
     *0.218439D+05,0.219439D+05,0.220440D+05,0.221449D+05,0.222457D+05,
     *0.223473D+05,0.224494D+05,0.225514D+05,0.226542D+05,0.227571D+05,
     *0.228606D+05,0.229646D+05,0.230687D+05,0.231734D+05,0.232783D+05,
     *0.233839D+05,0.234898D+05,0.235960D+05,0.237027D+05,0.238097D+05,
     *0.239173D+05,0.240254D+05,0.241335D+05,0.242424D+05,0.243514D+05,
     *0.244611D+05,0.245712D+05,0.246814D+05,0.247923D+05,0.249034D+05,
     *0.250152D+05/
!
      DO 10 I=1,NPNTS
!
       IF (T(I) .LE. T_LOW) THEN
          QS(I) = EPSILON * ES(1) / P(I)
!
       ELSE IF (T(I) .GE. T_HIGH) THEN
          QS(I) = EPSILON * ES(N) / P(I)
!
       ELSE
!
          ATABLE = (T(I) - T_LOW + DELTA_T) / DELTA_T
          ITABLE = ATABLE
          ATABLE = ATABLE - ITABLE
!
          QS(I)  = EPSILON * ((1.0D+00-ATABLE) * ES(ITABLE)
     *            +     ATABLE  * ES(ITABLE+1)) / P(I)
!
       END IF
!
 10   CONTINUE
!
      RETURN
      END
