*IF DEF,C90_2B QSAT2B.2
C (c) CROWN COPYRIGHT 1997, METEOROLOGICAL OFFICE, All Rights Reserved. QSAT2B.3
C QSAT2B.4
C Use, duplication or disclosure of this code is subject to the QSAT2B.5
C restrictions as set forth in the contract. QSAT2B.6
C QSAT2B.7
C Meteorological Office QSAT2B.8
C London Road QSAT2B.9
C BRACKNELL QSAT2B.10
C Berkshire UK QSAT2B.11
C RG12 2SZ QSAT2B.12
C QSAT2B.13
C If no contract has been raised with this copy of the code, the use, QSAT2B.14
C duplication or disclosure of it is strictly prohibited. Permission QSAT2B.15
C to do so must first be obtained in writing from the Head of Numerical QSAT2B.16
C Modelling at the above address. QSAT2B.17
C ******************************COPYRIGHT****************************** QSAT2B.18
C QSAT2B.19
CLL SUBROUTINES QSAT_WAT AND QSAT ---------------------------------- QSAT2B.20
CLL QSAT2B.21
CLL PURPOSE: QSAT2B.22
CLL RETURNS A SATURATION MIXING RATIO GIVEN A TEMPERATURE QSAT2B.23
CLL AND PRESSURE USING SATURATION VAPOUR PRESSURES CALCULATED USING THE QSAT2B.24
CLL GOFF-GRATCH FORMULAE, ADOPTED BY THE WMO AS TAKEN FROM QSAT2B.25
CLL LANDOLT-BORNSTEIN, 1987 NUMERICAL DATA AND FUNCTIONAL RELATIONSHIPS QSAT2B.26
CLL IN SCIENCE AND TECHNOLOGY. GROUP V/VOL 4B METEOROLOGY. PHYSICAL AND QSAT2B.27
CLL CHEMICAL PROPERTIES OF AIR, P35 QSAT2B.28
CLL QSAT2B.29
CLL VALUES IN THE LOOKUP TABLE ARE OVER WATER ABOVE AND BELOW 0 DEG C. QSAT2B.30
CLL QSAT2B.31
CLL NOTE : FOR VAPOUR PRESSURE OVER WATER THIS FORMULA IS VALID FOR QSAT2B.32
CLL TEMPERATURES BETWEEN 373K AND 223K. THE VALUES FOR SATURATED VAPOUR QSAT2B.33
CLL OVER WATER IN THE LOOKUP TABLE BELOW ARE OUT OF THE LOWER END OF QSAT2B.34
CLL THIS RANGE. HOWEVER IT IS STANDARD WMO PRACTICE TO USE THE FORMULA QSAT2B.35
CLL BELOW ITS ACCEPTED RANGE FOR USE WITH THE CALCULATION OF DEW POINTS QSAT2B.36
CLL IN THE UPPER ATMOSPHERE QSAT2B.37
CLL QSAT2B.38
CLL SUITABLE FOR SINGLE COLUMN MODEL USE QSAT2B.39
CLL QSAT2B.40
CLL CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90 QSAT2B.41
CLL QSAT2B.42
CLL MODEL MODIFICATION HISTORY: QSAT2B.43
CLL VERSION DATE QSAT2B.44
CLL ###### NOTE: VERSION IS T3E OPTIMISED, IMPROVED FORMULATION (2B) GHM1F405.8
CLL QSAT2B.45
!LL 4.4 11/08/97 New version optimised for T3E. QSAT2B.46
!LL Not bit-reproducible with QSAT2A. QSAT2B.47
CLL 4.4 11/8/97 Replace IF statements by MAX or MIN to improve QSAT2B.48
CLL speed of code on T3E. QSAT2B.49
CLL D. Salmond/A. Dickinson QSAT2B.50
CLL 4.5 25/06/98 Correct potential failure introduced by optimising: ARR6F405.1
CLL Extend lookup array from (1:) to (0:) to cater for ARR6F405.2
CLL special case of extreme low temperatures ARR6F405.3
CLL (.LE.T_LOW) for which the array index is rounded ARR6F405.4
CLL down due to machine precision. R.Rawlins ARR6F405.5
CLL QSAT2B.51
CLL QSAT2B.52
CLL PROGRAMMING STANDARDS : QSAT2B.53
CLL QSAT2B.54
CLL LOGICAL COMPONENTS COVERED: P27 QSAT2B.55
CLL QSAT2B.56
CLL DOCUMENTATION : GHM1F405.9
CLL QSAT2B.58
CLLEND----------------------------------------------------------------- QSAT2B.59
C QSAT2B.60
C*L ARGUMENTS--------------------------------------------------------- QSAT2B.61
C QSAT2B.62
SUBROUTINE QSAT_WAT (QS,T,P,NPNTS) 19QSAT2B.63
C QSAT2B.64
IMPLICIT NONE QSAT2B.65
C QSAT2B.66
C---------------------------------------------------------------------- QSAT2B.67
C MODEL CONSTANTS QSAT2B.68
C---------------------------------------------------------------------- QSAT2B.69
C QSAT2B.70
*CALL C_EPSLON 
QSAT2B.71
*CALL C_0_DG_C QSAT2B.72
C QSAT2B.73
C---------------------------------------------------------------------- QSAT2B.74
C LOCAL CONSTANTS QSAT2B.75
C---------------------------------------------------------------------- QSAT2B.76
C QSAT2B.77
REAL T_LOW ! LOWEST TEMPERATURE FOR WHICH LOOK-UP QSAT2B.78
! TABLE OF SATURATION WATER VAPOUR QSAT2B.79
! PRESSURE IS VALID (K) QSAT2B.80
C QSAT2B.81
REAL T_HIGH ! HIGHEST TEMPERATURE FOR WHICH LOOK-UP QSAT2B.82
! TABLE OF SATURATION WATER VAPOUR QSAT2B.83
! PRESSURES IS VALID (K) QSAT2B.84
C QSAT2B.85
REAL DELTA_T ! TEMPERATURE INCREMENT OF THE LOOK-UP QSAT2B.86
! TABLE OF SATURATION VAPOUR PRESSURES QSAT2B.87
C QSAT2B.88
INTEGER N ! SIZE OF LOOK-UP TABLE OF SATURATION QSAT2B.89
! WATER VAPOUR PRESSURES QSAT2B.90
C QSAT2B.91
PARAMETER ( T_LOW = 183.15, QSAT2B.92
* T_HIGH = 338.15, QSAT2B.93
* DELTA_T = 0.1, QSAT2B.94
* N = ((T_HIGH - T_LOW + (DELTA_T*0.5))/DELTA_T) + 1.0 QSAT2B.95
* ) ! gives N=1551 ARR6F405.6
C QSAT2B.97
C QSAT2B.98
C---------------------------------------------------------------------- QSAT2B.99
C VECTOR LENGTHS AND LOOP COUNTERS QSAT2B.100
C---------------------------------------------------------------------- QSAT2B.101
C QSAT2B.102
INTEGER NPNTS ! VECTOR LENGTH QSAT2B.103
C QSAT2B.104
INTEGER I ! LOOP COUNTER QSAT2B.105
C QSAT2B.106
INTEGER IES ! LOOP COUNTER FOR DATA STATEMENT QSAT2B.107
! LOOK-UP TABLE QSAT2B.108
C QSAT2B.109
C QSAT2B.110
C---------------------------------------------------------------------- QSAT2B.111
C VARIABLES WHICH ARE INPUT QSAT2B.112
C---------------------------------------------------------------------- QSAT2B.113
C QSAT2B.114
REAL T(NPNTS) ! IN TEMPERATURE (K) QSAT2B.115
C QSAT2B.116
REAL P(NPNTS) ! IN PRESSURE (PA) QSAT2B.117
C QSAT2B.118
C QSAT2B.119
C---------------------------------------------------------------------- QSAT2B.120
C VARIABLES WHICH ARE OUTPUT QSAT2B.121
C---------------------------------------------------------------------- QSAT2B.122
C QSAT2B.123
REAL QS(NPNTS) ! OUT SATURATION MIXING RATIO AT TEMPERATURE QSAT2B.124
! T AND PRESSURE P (KG/KG) QSAT2B.125
C QSAT2B.126
C QSAT2B.127
C---------------------------------------------------------------------- QSAT2B.128
C VARIABLES WHICH ARE DEFINED LOCALLY QSAT2B.129
C---------------------------------------------------------------------- QSAT2B.130
C QSAT2B.131
REAL ES(0:N+1) ! TABLE OF SATURATION WATER VAPOUR ARR6F405.7
! PRESSURE (PA) - SET BY DATA STATEMENT QSAT2B.133
! CALCULATED FROM THE GOFF-GRATCH FORMULAE QSAT2B.134
! AS TAKEN FROM LANDOLT-BORNSTEIN, 1987 QSAT2B.135
! NUMERICAL DATA AND FUNCTIONAL RELATIONSHIPS QSAT2B.136
! IN SCIENCE AND TECHNOLOGY. GROUP V/ VOL 4B QSAT2B.137
! METEOROLOGY. PHYSICAL AND CHEMICAL QSAT2B.138
! PROPERTIES OF AIR, P35 QSAT2B.139
C QSAT2B.140
REAL ATABLE ! WORK VARIABLES QSAT2B.141
C QSAT2B.142
INTEGER ITABLE ! WORK VARIABLES QSAT2B.143
C QSAT2B.144
C VARIABLES INTRODUCED BY DLR. QSAT2B.145
C QSAT2B.146
REAL FSUBW ! FACTOR THAT CONVERTS FROM SAT VAPOUR QSAT2B.147
! PRESSURE IN A PURE WATER SYSTEM TO QSAT2B.148
! SAT VAPOUR PRESSURE IN AIR. QSAT2B.149
C QSAT2B.150
REAL ONE_MINUS_EPSILON ! ONE MINUS THE RATIO OF THE MOLECULAR QSAT2B.151
! WEIGHTS OF WATER AND DRY AIR QSAT2B.152
REAL TT QSAT2B.153
C QSAT2B.154
PARAMETER(ONE_MINUS_EPSILON = 1.0 - EPSILON) QSAT2B.155
C QSAT2B.156
C*--------------------------------------------------------------------- QSAT2B.157
CL QSAT2B.158
CL--------------------------------------------------------------------- QSAT2B.159
CL NO SIGNIFICANT STRUCTURE QSAT2B.160
CL--------------------------------------------------------------------- QSAT2B.161
CL QSAT2B.162
C QSAT2B.163
C--------------------------------------------------------------------- QSAT2B.164
C SATURATION WATER VAPOUR PRESSURE QSAT2B.165
C QSAT2B.166
C VALUES ABOVE AND BELOW 0 DEG C ARE OVER WATER QSAT2B.167
C QSAT2B.168
C VALUES BELOW -50 DEC C ARE OUTSIDE FORMAL RANGE OF QSAT2B.169
C WMO ADOPTED FORMULA (SEE ABOVE). HOWEVER STANDARD QSAT2B.170
C PRACTICE TO USE THESE VALUES QSAT2B.171
C--------------------------------------------------------------------- QSAT2B.172
C QSAT2B.173
C Note: 0 element is a repeat of 1st element to cater for special case ARR6F405.8
C of low temperatures (.LE.T_LOW) for which the array index is ARR6F405.9
C rounded down due to machine precision. ARR6F405.10
DATA (ES(IES),IES= 0, 95) / 0.186905E-01, ARR6F405.11
*0.186905E-01,0.190449E-01,0.194059E-01,0.197727E-01,0.201462E-01, QSAT2B.175
*0.205261E-01,0.209122E-01,0.213052E-01,0.217050E-01,0.221116E-01, QSAT2B.176
*0.225252E-01,0.229463E-01,0.233740E-01,0.238090E-01,0.242518E-01, QSAT2B.177
*0.247017E-01,0.251595E-01,0.256252E-01,0.260981E-01,0.265795E-01, QSAT2B.178
*0.270691E-01,0.275667E-01,0.280733E-01,0.285876E-01,0.291105E-01, QSAT2B.179
*0.296429E-01,0.301835E-01,0.307336E-01,0.312927E-01,0.318611E-01, QSAT2B.180
*0.324390E-01,0.330262E-01,0.336232E-01,0.342306E-01,0.348472E-01, QSAT2B.181
*0.354748E-01,0.361117E-01,0.367599E-01,0.374185E-01,0.380879E-01, QSAT2B.182
*0.387689E-01,0.394602E-01,0.401626E-01,0.408771E-01,0.416033E-01, QSAT2B.183
*0.423411E-01,0.430908E-01,0.438524E-01,0.446263E-01,0.454124E-01, QSAT2B.184
*0.462122E-01,0.470247E-01,0.478491E-01,0.486874E-01,0.495393E-01, QSAT2B.185
*0.504057E-01,0.512847E-01,0.521784E-01,0.530853E-01,0.540076E-01, QSAT2B.186
*0.549444E-01,0.558959E-01,0.568633E-01,0.578448E-01,0.588428E-01, QSAT2B.187
*0.598566E-01,0.608858E-01,0.619313E-01,0.629926E-01,0.640706E-01, QSAT2B.188
*0.651665E-01,0.662795E-01,0.674095E-01,0.685570E-01,0.697219E-01, QSAT2B.189
*0.709063E-01,0.721076E-01,0.733284E-01,0.745679E-01,0.758265E-01, QSAT2B.190
*0.771039E-01,0.784026E-01,0.797212E-01,0.810577E-01,0.824164E-01, QSAT2B.191
*0.837971E-01,0.851970E-01,0.866198E-01,0.880620E-01,0.895281E-01, QSAT2B.192
*0.910178E-01,0.925278E-01,0.940622E-01,0.956177E-01,0.971984E-01/ QSAT2B.193
DATA (ES(IES),IES= 96,190) / QSAT2B.194
*0.988051E-01,0.100433E+00,0.102085E+00,0.103764E+00,0.105467E+00, QSAT2B.195
*0.107196E+00,0.108953E+00,0.110732E+00,0.112541E+00,0.114376E+00, QSAT2B.196
*0.116238E+00,0.118130E+00,0.120046E+00,0.121993E+00,0.123969E+00, QSAT2B.197
*0.125973E+00,0.128009E+00,0.130075E+00,0.132167E+00,0.134296E+00, QSAT2B.198
*0.136452E+00,0.138642E+00,0.140861E+00,0.143115E+00,0.145404E+00, QSAT2B.199
*0.147723E+00,0.150078E+00,0.152466E+00,0.154889E+00,0.157346E+00, QSAT2B.200
*0.159841E+00,0.162372E+00,0.164939E+00,0.167545E+00,0.170185E+00, QSAT2B.201
*0.172866E+00,0.175584E+00,0.178340E+00,0.181139E+00,0.183977E+00, QSAT2B.202
*0.186855E+00,0.189773E+00,0.192737E+00,0.195736E+00,0.198783E+00, QSAT2B.203
*0.201875E+00,0.205007E+00,0.208186E+00,0.211409E+00,0.214676E+00, QSAT2B.204
*0.217993E+00,0.221355E+00,0.224764E+00,0.228220E+00,0.231728E+00, QSAT2B.205
*0.235284E+00,0.238888E+00,0.242542E+00,0.246251E+00,0.250010E+00, QSAT2B.206
*0.253821E+00,0.257688E+00,0.261602E+00,0.265575E+00,0.269607E+00, QSAT2B.207
*0.273689E+00,0.277830E+00,0.282027E+00,0.286287E+00,0.290598E+00, QSAT2B.208
*0.294972E+00,0.299405E+00,0.303904E+00,0.308462E+00,0.313082E+00, QSAT2B.209
*0.317763E+00,0.322512E+00,0.327324E+00,0.332201E+00,0.337141E+00, QSAT2B.210
*0.342154E+00,0.347234E+00,0.352387E+00,0.357601E+00,0.362889E+00, QSAT2B.211
*0.368257E+00,0.373685E+00,0.379194E+00,0.384773E+00,0.390433E+00, QSAT2B.212
*0.396159E+00,0.401968E+00,0.407861E+00,0.413820E+00,0.419866E+00/ QSAT2B.213
DATA (ES(IES),IES=191,285) / QSAT2B.214
*0.425999E+00,0.432203E+00,0.438494E+00,0.444867E+00,0.451332E+00, QSAT2B.215
*0.457879E+00,0.464510E+00,0.471226E+00,0.478037E+00,0.484935E+00, QSAT2B.216
*0.491920E+00,0.499005E+00,0.506181E+00,0.513447E+00,0.520816E+00, QSAT2B.217
*0.528279E+00,0.535835E+00,0.543497E+00,0.551256E+00,0.559113E+00, QSAT2B.218
*0.567081E+00,0.575147E+00,0.583315E+00,0.591585E+00,0.599970E+00, QSAT2B.219
*0.608472E+00,0.617069E+00,0.625785E+00,0.634609E+00,0.643556E+00, QSAT2B.220
*0.652611E+00,0.661782E+00,0.671077E+00,0.680487E+00,0.690015E+00, QSAT2B.221
*0.699656E+00,0.709433E+00,0.719344E+00,0.729363E+00,0.739518E+00, QSAT2B.222
*0.749795E+00,0.760217E+00,0.770763E+00,0.781454E+00,0.792258E+00, QSAT2B.223
*0.803208E+00,0.814309E+00,0.825528E+00,0.836914E+00,0.848422E+00, QSAT2B.224
*0.860086E+00,0.871891E+00,0.883837E+00,0.895944E+00,0.908214E+00, QSAT2B.225
*0.920611E+00,0.933175E+00,0.945890E+00,0.958776E+00,0.971812E+00, QSAT2B.226
*0.985027E+00,0.998379E+00,0.101193E+01,0.102561E+01,0.103949E+01, QSAT2B.227
*0.105352E+01,0.106774E+01,0.108213E+01,0.109669E+01,0.111144E+01, QSAT2B.228
*0.112636E+01,0.114148E+01,0.115676E+01,0.117226E+01,0.118791E+01, QSAT2B.229
*0.120377E+01,0.121984E+01,0.123608E+01,0.125252E+01,0.126919E+01, QSAT2B.230
*0.128604E+01,0.130309E+01,0.132036E+01,0.133782E+01,0.135549E+01, QSAT2B.231
*0.137339E+01,0.139150E+01,0.140984E+01,0.142839E+01,0.144715E+01, QSAT2B.232
*0.146616E+01,0.148538E+01,0.150482E+01,0.152450E+01,0.154445E+01/ QSAT2B.233
DATA (ES(IES),IES=286,380) / QSAT2B.234
*0.156459E+01,0.158502E+01,0.160564E+01,0.162654E+01,0.164766E+01, QSAT2B.235
*0.166906E+01,0.169070E+01,0.171257E+01,0.173473E+01,0.175718E+01, QSAT2B.236
*0.177984E+01,0.180282E+01,0.182602E+01,0.184951E+01,0.187327E+01, QSAT2B.237
*0.189733E+01,0.192165E+01,0.194629E+01,0.197118E+01,0.199636E+01, QSAT2B.238
*0.202185E+01,0.204762E+01,0.207372E+01,0.210010E+01,0.212678E+01, QSAT2B.239
*0.215379E+01,0.218109E+01,0.220873E+01,0.223668E+01,0.226497E+01, QSAT2B.240
*0.229357E+01,0.232249E+01,0.235176E+01,0.238134E+01,0.241129E+01, QSAT2B.241
*0.244157E+01,0.247217E+01,0.250316E+01,0.253447E+01,0.256617E+01, QSAT2B.242
*0.259821E+01,0.263064E+01,0.266341E+01,0.269661E+01,0.273009E+01, QSAT2B.243
*0.276403E+01,0.279834E+01,0.283302E+01,0.286811E+01,0.290358E+01, QSAT2B.244
*0.293943E+01,0.297571E+01,0.301236E+01,0.304946E+01,0.308702E+01, QSAT2B.245
*0.312491E+01,0.316326E+01,0.320208E+01,0.324130E+01,0.328092E+01, QSAT2B.246
*0.332102E+01,0.336162E+01,0.340264E+01,0.344407E+01,0.348601E+01, QSAT2B.247
*0.352838E+01,0.357118E+01,0.361449E+01,0.365834E+01,0.370264E+01, QSAT2B.248
*0.374737E+01,0.379265E+01,0.383839E+01,0.388469E+01,0.393144E+01, QSAT2B.249
*0.397876E+01,0.402656E+01,0.407492E+01,0.412378E+01,0.417313E+01, QSAT2B.250
*0.422306E+01,0.427359E+01,0.432454E+01,0.437617E+01,0.442834E+01, QSAT2B.251
*0.448102E+01,0.453433E+01,0.458816E+01,0.464253E+01,0.469764E+01, QSAT2B.252
*0.475321E+01,0.480942E+01,0.486629E+01,0.492372E+01,0.498173E+01/ QSAT2B.253
DATA (ES(IES),IES=381,475) / QSAT2B.254
*0.504041E+01,0.509967E+01,0.515962E+01,0.522029E+01,0.528142E+01, QSAT2B.255
*0.534337E+01,0.540595E+01,0.546912E+01,0.553292E+01,0.559757E+01, QSAT2B.256
*0.566273E+01,0.572864E+01,0.579532E+01,0.586266E+01,0.593075E+01, QSAT2B.257
*0.599952E+01,0.606895E+01,0.613918E+01,0.621021E+01,0.628191E+01, QSAT2B.258
*0.635433E+01,0.642755E+01,0.650162E+01,0.657639E+01,0.665188E+01, QSAT2B.259
*0.672823E+01,0.680532E+01,0.688329E+01,0.696198E+01,0.704157E+01, QSAT2B.260
*0.712206E+01,0.720319E+01,0.728534E+01,0.736829E+01,0.745204E+01, QSAT2B.261
*0.753671E+01,0.762218E+01,0.770860E+01,0.779588E+01,0.788408E+01, QSAT2B.262
*0.797314E+01,0.806318E+01,0.815408E+01,0.824599E+01,0.833874E+01, QSAT2B.263
*0.843254E+01,0.852721E+01,0.862293E+01,0.871954E+01,0.881724E+01, QSAT2B.264
*0.891579E+01,0.901547E+01,0.911624E+01,0.921778E+01,0.932061E+01, QSAT2B.265
*0.942438E+01,0.952910E+01,0.963497E+01,0.974181E+01,0.984982E+01, QSAT2B.266
*0.995887E+01,0.100690E+02,0.101804E+02,0.102926E+02,0.104063E+02, QSAT2B.267
*0.105210E+02,0.106367E+02,0.107536E+02,0.108719E+02,0.109912E+02, QSAT2B.268
*0.111116E+02,0.112333E+02,0.113563E+02,0.114804E+02,0.116056E+02, QSAT2B.269
*0.117325E+02,0.118602E+02,0.119892E+02,0.121197E+02,0.122513E+02, QSAT2B.270
*0.123844E+02,0.125186E+02,0.126543E+02,0.127912E+02,0.129295E+02, QSAT2B.271
*0.130691E+02,0.132101E+02,0.133527E+02,0.134965E+02,0.136415E+02, QSAT2B.272
*0.137882E+02,0.139361E+02,0.140855E+02,0.142366E+02,0.143889E+02/ QSAT2B.273
DATA (ES(IES),IES=476,570) / QSAT2B.274
*0.145429E+02,0.146982E+02,0.148552E+02,0.150135E+02,0.151735E+02, QSAT2B.275
*0.153349E+02,0.154979E+02,0.156624E+02,0.158286E+02,0.159965E+02, QSAT2B.276
*0.161659E+02,0.163367E+02,0.165094E+02,0.166838E+02,0.168597E+02, QSAT2B.277
*0.170375E+02,0.172168E+02,0.173979E+02,0.175806E+02,0.177651E+02, QSAT2B.278
*0.179513E+02,0.181394E+02,0.183293E+02,0.185210E+02,0.187146E+02, QSAT2B.279
*0.189098E+02,0.191066E+02,0.193059E+02,0.195065E+02,0.197095E+02, QSAT2B.280
*0.199142E+02,0.201206E+02,0.203291E+02,0.205397E+02,0.207522E+02, QSAT2B.281
*0.209664E+02,0.211831E+02,0.214013E+02,0.216221E+02,0.218448E+02, QSAT2B.282
*0.220692E+02,0.222959E+02,0.225250E+02,0.227559E+02,0.229887E+02, QSAT2B.283
*0.232239E+02,0.234614E+02,0.237014E+02,0.239428E+02,0.241872E+02, QSAT2B.284
*0.244335E+02,0.246824E+02,0.249332E+02,0.251860E+02,0.254419E+02, QSAT2B.285
*0.256993E+02,0.259600E+02,0.262225E+02,0.264873E+02,0.267552E+02, QSAT2B.286
*0.270248E+02,0.272970E+02,0.275719E+02,0.278497E+02,0.281295E+02, QSAT2B.287
*0.284117E+02,0.286965E+02,0.289843E+02,0.292743E+02,0.295671E+02, QSAT2B.288
*0.298624E+02,0.301605E+02,0.304616E+02,0.307650E+02,0.310708E+02, QSAT2B.289
*0.313803E+02,0.316915E+02,0.320064E+02,0.323238E+02,0.326437E+02, QSAT2B.290
*0.329666E+02,0.332928E+02,0.336215E+02,0.339534E+02,0.342885E+02, QSAT2B.291
*0.346263E+02,0.349666E+02,0.353109E+02,0.356572E+02,0.360076E+02, QSAT2B.292
*0.363606E+02,0.367164E+02,0.370757E+02,0.374383E+02,0.378038E+02/ QSAT2B.293
DATA (ES(IES),IES=571,665) / QSAT2B.294
*0.381727E+02,0.385453E+02,0.389206E+02,0.392989E+02,0.396807E+02, QSAT2B.295
*0.400663E+02,0.404555E+02,0.408478E+02,0.412428E+02,0.416417E+02, QSAT2B.296
*0.420445E+02,0.424502E+02,0.428600E+02,0.432733E+02,0.436900E+02, QSAT2B.297
*0.441106E+02,0.445343E+02,0.449620E+02,0.453930E+02,0.458280E+02, QSAT2B.298
*0.462672E+02,0.467096E+02,0.471561E+02,0.476070E+02,0.480610E+02, QSAT2B.299
*0.485186E+02,0.489813E+02,0.494474E+02,0.499170E+02,0.503909E+02, QSAT2B.300
*0.508693E+02,0.513511E+02,0.518376E+02,0.523277E+02,0.528232E+02, QSAT2B.301
*0.533213E+02,0.538240E+02,0.543315E+02,0.548437E+02,0.553596E+02, QSAT2B.302
*0.558802E+02,0.564046E+02,0.569340E+02,0.574672E+02,0.580061E+02, QSAT2B.303
*0.585481E+02,0.590963E+02,0.596482E+02,0.602041E+02,0.607649E+02, QSAT2B.304
*0.613311E+02,0.619025E+02,0.624779E+02,0.630574E+02,0.636422E+02, QSAT2B.305
*0.642324E+02,0.648280E+02,0.654278E+02,0.660332E+02,0.666426E+02, QSAT2B.306
*0.672577E+02,0.678771E+02,0.685034E+02,0.691328E+02,0.697694E+02, QSAT2B.307
*0.704103E+02,0.710556E+02,0.717081E+02,0.723639E+02,0.730269E+02, QSAT2B.308
*0.736945E+02,0.743681E+02,0.750463E+02,0.757309E+02,0.764214E+02, QSAT2B.309
*0.771167E+02,0.778182E+02,0.785246E+02,0.792373E+02,0.799564E+02, QSAT2B.310
*0.806804E+02,0.814109E+02,0.821479E+02,0.828898E+02,0.836384E+02, QSAT2B.311
*0.843922E+02,0.851525E+02,0.859198E+02,0.866920E+02,0.874712E+02, QSAT2B.312
*0.882574E+02,0.890486E+02,0.898470E+02,0.906525E+02,0.914634E+02/ QSAT2B.313
DATA (ES(IES),IES=666,760) / QSAT2B.314
*0.922814E+02,0.931048E+02,0.939356E+02,0.947736E+02,0.956171E+02, QSAT2B.315
*0.964681E+02,0.973246E+02,0.981907E+02,0.990605E+02,0.999399E+02, QSAT2B.316
*0.100825E+03,0.101718E+03,0.102617E+03,0.103523E+03,0.104438E+03, QSAT2B.317
*0.105358E+03,0.106287E+03,0.107221E+03,0.108166E+03,0.109115E+03, QSAT2B.318
*0.110074E+03,0.111039E+03,0.112012E+03,0.112992E+03,0.113981E+03, QSAT2B.319
*0.114978E+03,0.115981E+03,0.116993E+03,0.118013E+03,0.119041E+03, QSAT2B.320
*0.120077E+03,0.121122E+03,0.122173E+03,0.123234E+03,0.124301E+03, QSAT2B.321
*0.125377E+03,0.126463E+03,0.127556E+03,0.128657E+03,0.129769E+03, QSAT2B.322
*0.130889E+03,0.132017E+03,0.133152E+03,0.134299E+03,0.135453E+03, QSAT2B.323
*0.136614E+03,0.137786E+03,0.138967E+03,0.140158E+03,0.141356E+03, QSAT2B.324
*0.142565E+03,0.143781E+03,0.145010E+03,0.146247E+03,0.147491E+03, QSAT2B.325
*0.148746E+03,0.150011E+03,0.151284E+03,0.152571E+03,0.153862E+03, QSAT2B.326
*0.155168E+03,0.156481E+03,0.157805E+03,0.159137E+03,0.160478E+03, QSAT2B.327
*0.161832E+03,0.163198E+03,0.164569E+03,0.165958E+03,0.167348E+03, QSAT2B.328
*0.168757E+03,0.170174E+03,0.171599E+03,0.173037E+03,0.174483E+03, QSAT2B.329
*0.175944E+03,0.177414E+03,0.178892E+03,0.180387E+03,0.181886E+03, QSAT2B.330
*0.183402E+03,0.184930E+03,0.186463E+03,0.188012E+03,0.189571E+03, QSAT2B.331
*0.191146E+03,0.192730E+03,0.194320E+03,0.195930E+03,0.197546E+03, QSAT2B.332
*0.199175E+03,0.200821E+03,0.202473E+03,0.204142E+03,0.205817E+03/ QSAT2B.333
DATA (ES(IES),IES=761,855) / QSAT2B.334
*0.207510E+03,0.209216E+03,0.210928E+03,0.212658E+03,0.214398E+03, QSAT2B.335
*0.216152E+03,0.217920E+03,0.219698E+03,0.221495E+03,0.223297E+03, QSAT2B.336
*0.225119E+03,0.226951E+03,0.228793E+03,0.230654E+03,0.232522E+03, QSAT2B.337
*0.234413E+03,0.236311E+03,0.238223E+03,0.240151E+03,0.242090E+03, QSAT2B.338
*0.244049E+03,0.246019E+03,0.248000E+03,0.249996E+03,0.252009E+03, QSAT2B.339
*0.254037E+03,0.256077E+03,0.258128E+03,0.260200E+03,0.262284E+03, QSAT2B.340
*0.264384E+03,0.266500E+03,0.268629E+03,0.270779E+03,0.272936E+03, QSAT2B.341
*0.275110E+03,0.277306E+03,0.279509E+03,0.281734E+03,0.283966E+03, QSAT2B.342
*0.286227E+03,0.288494E+03,0.290780E+03,0.293083E+03,0.295398E+03, QSAT2B.343
*0.297737E+03,0.300089E+03,0.302453E+03,0.304841E+03,0.307237E+03, QSAT2B.344
*0.309656E+03,0.312095E+03,0.314541E+03,0.317012E+03,0.319496E+03, QSAT2B.345
*0.322005E+03,0.324527E+03,0.327063E+03,0.329618E+03,0.332193E+03, QSAT2B.346
*0.334788E+03,0.337396E+03,0.340025E+03,0.342673E+03,0.345329E+03, QSAT2B.347
*0.348019E+03,0.350722E+03,0.353440E+03,0.356178E+03,0.358938E+03, QSAT2B.348
*0.361718E+03,0.364513E+03,0.367322E+03,0.370160E+03,0.373012E+03, QSAT2B.349
*0.375885E+03,0.378788E+03,0.381691E+03,0.384631E+03,0.387579E+03, QSAT2B.350
*0.390556E+03,0.393556E+03,0.396563E+03,0.399601E+03,0.402646E+03, QSAT2B.351
*0.405730E+03,0.408829E+03,0.411944E+03,0.415083E+03,0.418236E+03, QSAT2B.352
*0.421422E+03,0.424632E+03,0.427849E+03,0.431099E+03,0.434365E+03/ QSAT2B.353
DATA (ES(IES),IES=856,950) / QSAT2B.354
*0.437655E+03,0.440970E+03,0.444301E+03,0.447666E+03,0.451038E+03, QSAT2B.355
*0.454445E+03,0.457876E+03,0.461316E+03,0.464790E+03,0.468281E+03, QSAT2B.356
*0.471798E+03,0.475342E+03,0.478902E+03,0.482497E+03,0.486101E+03, QSAT2B.357
*0.489741E+03,0.493408E+03,0.497083E+03,0.500804E+03,0.504524E+03, QSAT2B.358
*0.508290E+03,0.512074E+03,0.515877E+03,0.519717E+03,0.523566E+03, QSAT2B.359
*0.527462E+03,0.531367E+03,0.535301E+03,0.539264E+03,0.543245E+03, QSAT2B.360
*0.547265E+03,0.551305E+03,0.555363E+03,0.559462E+03,0.563579E+03, QSAT2B.361
*0.567727E+03,0.571905E+03,0.576102E+03,0.580329E+03,0.584576E+03, QSAT2B.362
*0.588865E+03,0.593185E+03,0.597514E+03,0.601885E+03,0.606276E+03, QSAT2B.363
*0.610699E+03,0.615151E+03,0.619625E+03,0.624140E+03,0.628671E+03, QSAT2B.364
*0.633243E+03,0.637845E+03,0.642465E+03,0.647126E+03,0.651806E+03, QSAT2B.365
*0.656527E+03,0.661279E+03,0.666049E+03,0.670861E+03,0.675692E+03, QSAT2B.366
*0.680566E+03,0.685471E+03,0.690396E+03,0.695363E+03,0.700350E+03, QSAT2B.367
*0.705381E+03,0.710444E+03,0.715527E+03,0.720654E+03,0.725801E+03, QSAT2B.368
*0.730994E+03,0.736219E+03,0.741465E+03,0.746756E+03,0.752068E+03, QSAT2B.369
*0.757426E+03,0.762819E+03,0.768231E+03,0.773692E+03,0.779172E+03, QSAT2B.370
*0.784701E+03,0.790265E+03,0.795849E+03,0.801483E+03,0.807137E+03, QSAT2B.371
*0.812842E+03,0.818582E+03,0.824343E+03,0.830153E+03,0.835987E+03, QSAT2B.372
*0.841871E+03,0.847791E+03,0.853733E+03,0.859727E+03,0.865743E+03/ QSAT2B.373
DATA (ES(IES),IES=951,1045) / QSAT2B.374
*0.871812E+03,0.877918E+03,0.884046E+03,0.890228E+03,0.896433E+03, QSAT2B.375
*0.902690E+03,0.908987E+03,0.915307E+03,0.921681E+03,0.928078E+03, QSAT2B.376
*0.934531E+03,0.941023E+03,0.947539E+03,0.954112E+03,0.960708E+03, QSAT2B.377
*0.967361E+03,0.974053E+03,0.980771E+03,0.987545E+03,0.994345E+03, QSAT2B.378
*0.100120E+04,0.100810E+04,0.101502E+04,0.102201E+04,0.102902E+04, QSAT2B.379
*0.103608E+04,0.104320E+04,0.105033E+04,0.105753E+04,0.106475E+04, QSAT2B.380
*0.107204E+04,0.107936E+04,0.108672E+04,0.109414E+04,0.110158E+04, QSAT2B.381
*0.110908E+04,0.111663E+04,0.112421E+04,0.113185E+04,0.113952E+04, QSAT2B.382
*0.114725E+04,0.115503E+04,0.116284E+04,0.117071E+04,0.117861E+04, QSAT2B.383
*0.118658E+04,0.119459E+04,0.120264E+04,0.121074E+04,0.121888E+04, QSAT2B.384
*0.122709E+04,0.123534E+04,0.124362E+04,0.125198E+04,0.126036E+04, QSAT2B.385
*0.126881E+04,0.127731E+04,0.128584E+04,0.129444E+04,0.130307E+04, QSAT2B.386
*0.131177E+04,0.132053E+04,0.132931E+04,0.133817E+04,0.134705E+04, QSAT2B.387
*0.135602E+04,0.136503E+04,0.137407E+04,0.138319E+04,0.139234E+04, QSAT2B.388
*0.140156E+04,0.141084E+04,0.142015E+04,0.142954E+04,0.143896E+04, QSAT2B.389
*0.144845E+04,0.145800E+04,0.146759E+04,0.147725E+04,0.148694E+04, QSAT2B.390
*0.149672E+04,0.150655E+04,0.151641E+04,0.152635E+04,0.153633E+04, QSAT2B.391
*0.154639E+04,0.155650E+04,0.156665E+04,0.157688E+04,0.158715E+04, QSAT2B.392
*0.159750E+04,0.160791E+04,0.161836E+04,0.162888E+04,0.163945E+04/ QSAT2B.393
DATA (ES(IES),IES=1046,1140) / QSAT2B.394
*0.165010E+04,0.166081E+04,0.167155E+04,0.168238E+04,0.169325E+04, QSAT2B.395
*0.170420E+04,0.171522E+04,0.172627E+04,0.173741E+04,0.174859E+04, QSAT2B.396
*0.175986E+04,0.177119E+04,0.178256E+04,0.179402E+04,0.180552E+04, QSAT2B.397
*0.181711E+04,0.182877E+04,0.184046E+04,0.185224E+04,0.186407E+04, QSAT2B.398
*0.187599E+04,0.188797E+04,0.190000E+04,0.191212E+04,0.192428E+04, QSAT2B.399
*0.193653E+04,0.194886E+04,0.196122E+04,0.197368E+04,0.198618E+04, QSAT2B.400
*0.199878E+04,0.201145E+04,0.202416E+04,0.203698E+04,0.204983E+04, QSAT2B.401
*0.206278E+04,0.207580E+04,0.208887E+04,0.210204E+04,0.211525E+04, QSAT2B.402
*0.212856E+04,0.214195E+04,0.215538E+04,0.216892E+04,0.218249E+04, QSAT2B.403
*0.219618E+04,0.220994E+04,0.222375E+04,0.223766E+04,0.225161E+04, QSAT2B.404
*0.226567E+04,0.227981E+04,0.229399E+04,0.230829E+04,0.232263E+04, QSAT2B.405
*0.233708E+04,0.235161E+04,0.236618E+04,0.238087E+04,0.239560E+04, QSAT2B.406
*0.241044E+04,0.242538E+04,0.244035E+04,0.245544E+04,0.247057E+04, QSAT2B.407
*0.248583E+04,0.250116E+04,0.251654E+04,0.253204E+04,0.254759E+04, QSAT2B.408
*0.256325E+04,0.257901E+04,0.259480E+04,0.261073E+04,0.262670E+04, QSAT2B.409
*0.264279E+04,0.265896E+04,0.267519E+04,0.269154E+04,0.270794E+04, QSAT2B.410
*0.272447E+04,0.274108E+04,0.275774E+04,0.277453E+04,0.279137E+04, QSAT2B.411
*0.280834E+04,0.282540E+04,0.284251E+04,0.285975E+04,0.287704E+04, QSAT2B.412
*0.289446E+04,0.291198E+04,0.292954E+04,0.294725E+04,0.296499E+04/ QSAT2B.413
DATA (ES(IES),IES=1141,1235) / QSAT2B.414
*0.298288E+04,0.300087E+04,0.301890E+04,0.303707E+04,0.305529E+04, QSAT2B.415
*0.307365E+04,0.309211E+04,0.311062E+04,0.312927E+04,0.314798E+04, QSAT2B.416
*0.316682E+04,0.318577E+04,0.320477E+04,0.322391E+04,0.324310E+04, QSAT2B.417
*0.326245E+04,0.328189E+04,0.330138E+04,0.332103E+04,0.334073E+04, QSAT2B.418
*0.336058E+04,0.338053E+04,0.340054E+04,0.342069E+04,0.344090E+04, QSAT2B.419
*0.346127E+04,0.348174E+04,0.350227E+04,0.352295E+04,0.354369E+04, QSAT2B.420
*0.356458E+04,0.358559E+04,0.360664E+04,0.362787E+04,0.364914E+04, QSAT2B.421
*0.367058E+04,0.369212E+04,0.371373E+04,0.373548E+04,0.375731E+04, QSAT2B.422
*0.377929E+04,0.380139E+04,0.382355E+04,0.384588E+04,0.386826E+04, QSAT2B.423
*0.389081E+04,0.391348E+04,0.393620E+04,0.395910E+04,0.398205E+04, QSAT2B.424
*0.400518E+04,0.402843E+04,0.405173E+04,0.407520E+04,0.409875E+04, QSAT2B.425
*0.412246E+04,0.414630E+04,0.417019E+04,0.419427E+04,0.421840E+04, QSAT2B.426
*0.424272E+04,0.426715E+04,0.429165E+04,0.431634E+04,0.434108E+04, QSAT2B.427
*0.436602E+04,0.439107E+04,0.441618E+04,0.444149E+04,0.446685E+04, QSAT2B.428
*0.449241E+04,0.451810E+04,0.454385E+04,0.456977E+04,0.459578E+04, QSAT2B.429
*0.462197E+04,0.464830E+04,0.467468E+04,0.470127E+04,0.472792E+04, QSAT2B.430
*0.475477E+04,0.478175E+04,0.480880E+04,0.483605E+04,0.486336E+04, QSAT2B.431
*0.489087E+04,0.491853E+04,0.494623E+04,0.497415E+04,0.500215E+04, QSAT2B.432
*0.503034E+04,0.505867E+04,0.508707E+04,0.511568E+04,0.514436E+04/ QSAT2B.433
DATA (ES(IES),IES=1236,1330) / QSAT2B.434
*0.517325E+04,0.520227E+04,0.523137E+04,0.526068E+04,0.529005E+04, QSAT2B.435
*0.531965E+04,0.534939E+04,0.537921E+04,0.540923E+04,0.543932E+04, QSAT2B.436
*0.546965E+04,0.550011E+04,0.553064E+04,0.556139E+04,0.559223E+04, QSAT2B.437
*0.562329E+04,0.565449E+04,0.568577E+04,0.571727E+04,0.574884E+04, QSAT2B.438
*0.578064E+04,0.581261E+04,0.584464E+04,0.587692E+04,0.590924E+04, QSAT2B.439
*0.594182E+04,0.597455E+04,0.600736E+04,0.604039E+04,0.607350E+04, QSAT2B.440
*0.610685E+04,0.614036E+04,0.617394E+04,0.620777E+04,0.624169E+04, QSAT2B.441
*0.627584E+04,0.631014E+04,0.634454E+04,0.637918E+04,0.641390E+04, QSAT2B.442
*0.644887E+04,0.648400E+04,0.651919E+04,0.655467E+04,0.659021E+04, QSAT2B.443
*0.662599E+04,0.666197E+04,0.669800E+04,0.673429E+04,0.677069E+04, QSAT2B.444
*0.680735E+04,0.684415E+04,0.688104E+04,0.691819E+04,0.695543E+04, QSAT2B.445
*0.699292E+04,0.703061E+04,0.706837E+04,0.710639E+04,0.714451E+04, QSAT2B.446
*0.718289E+04,0.722143E+04,0.726009E+04,0.729903E+04,0.733802E+04, QSAT2B.447
*0.737729E+04,0.741676E+04,0.745631E+04,0.749612E+04,0.753602E+04, QSAT2B.448
*0.757622E+04,0.761659E+04,0.765705E+04,0.769780E+04,0.773863E+04, QSAT2B.449
*0.777975E+04,0.782106E+04,0.786246E+04,0.790412E+04,0.794593E+04, QSAT2B.450
*0.798802E+04,0.803028E+04,0.807259E+04,0.811525E+04,0.815798E+04, QSAT2B.451
*0.820102E+04,0.824427E+04,0.828757E+04,0.833120E+04,0.837493E+04, QSAT2B.452
*0.841895E+04,0.846313E+04,0.850744E+04,0.855208E+04,0.859678E+04/ QSAT2B.453
DATA (ES(IES),IES=1331,1425) / QSAT2B.454
*0.864179E+04,0.868705E+04,0.873237E+04,0.877800E+04,0.882374E+04, QSAT2B.455
*0.886979E+04,0.891603E+04,0.896237E+04,0.900904E+04,0.905579E+04, QSAT2B.456
*0.910288E+04,0.915018E+04,0.919758E+04,0.924529E+04,0.929310E+04, QSAT2B.457
*0.934122E+04,0.938959E+04,0.943804E+04,0.948687E+04,0.953575E+04, QSAT2B.458
*0.958494E+04,0.963442E+04,0.968395E+04,0.973384E+04,0.978383E+04, QSAT2B.459
*0.983412E+04,0.988468E+04,0.993534E+04,0.998630E+04,0.100374E+05, QSAT2B.460
*0.100888E+05,0.101406E+05,0.101923E+05,0.102444E+05,0.102966E+05, QSAT2B.461
*0.103492E+05,0.104020E+05,0.104550E+05,0.105082E+05,0.105616E+05, QSAT2B.462
*0.106153E+05,0.106693E+05,0.107234E+05,0.107779E+05,0.108325E+05, QSAT2B.463
*0.108874E+05,0.109425E+05,0.109978E+05,0.110535E+05,0.111092E+05, QSAT2B.464
*0.111653E+05,0.112217E+05,0.112782E+05,0.113350E+05,0.113920E+05, QSAT2B.465
*0.114493E+05,0.115070E+05,0.115646E+05,0.116228E+05,0.116809E+05, QSAT2B.466
*0.117396E+05,0.117984E+05,0.118574E+05,0.119167E+05,0.119762E+05, QSAT2B.467
*0.120360E+05,0.120962E+05,0.121564E+05,0.122170E+05,0.122778E+05, QSAT2B.468
*0.123389E+05,0.124004E+05,0.124619E+05,0.125238E+05,0.125859E+05, QSAT2B.469
*0.126484E+05,0.127111E+05,0.127739E+05,0.128372E+05,0.129006E+05, QSAT2B.470
*0.129644E+05,0.130285E+05,0.130927E+05,0.131573E+05,0.132220E+05, QSAT2B.471
*0.132872E+05,0.133526E+05,0.134182E+05,0.134842E+05,0.135503E+05, QSAT2B.472
*0.136168E+05,0.136836E+05,0.137505E+05,0.138180E+05,0.138854E+05/ QSAT2B.473
DATA (ES(IES),IES=1426,1520) / QSAT2B.474
*0.139534E+05,0.140216E+05,0.140900E+05,0.141588E+05,0.142277E+05, QSAT2B.475
*0.142971E+05,0.143668E+05,0.144366E+05,0.145069E+05,0.145773E+05, QSAT2B.476
*0.146481E+05,0.147192E+05,0.147905E+05,0.148622E+05,0.149341E+05, QSAT2B.477
*0.150064E+05,0.150790E+05,0.151517E+05,0.152250E+05,0.152983E+05, QSAT2B.478
*0.153721E+05,0.154462E+05,0.155205E+05,0.155952E+05,0.156701E+05, QSAT2B.479
*0.157454E+05,0.158211E+05,0.158969E+05,0.159732E+05,0.160496E+05, QSAT2B.480
*0.161265E+05,0.162037E+05,0.162811E+05,0.163589E+05,0.164369E+05, QSAT2B.481
*0.165154E+05,0.165942E+05,0.166732E+05,0.167526E+05,0.168322E+05, QSAT2B.482
*0.169123E+05,0.169927E+05,0.170733E+05,0.171543E+05,0.172356E+05, QSAT2B.483
*0.173173E+05,0.173993E+05,0.174815E+05,0.175643E+05,0.176471E+05, QSAT2B.484
*0.177305E+05,0.178143E+05,0.178981E+05,0.179826E+05,0.180671E+05, QSAT2B.485
*0.181522E+05,0.182377E+05,0.183232E+05,0.184093E+05,0.184955E+05, QSAT2B.486
*0.185823E+05,0.186695E+05,0.187568E+05,0.188447E+05,0.189326E+05, QSAT2B.487
*0.190212E+05,0.191101E+05,0.191991E+05,0.192887E+05,0.193785E+05, QSAT2B.488
*0.194688E+05,0.195595E+05,0.196503E+05,0.197417E+05,0.198332E+05, QSAT2B.489
*0.199253E+05,0.200178E+05,0.201105E+05,0.202036E+05,0.202971E+05, QSAT2B.490
*0.203910E+05,0.204853E+05,0.205798E+05,0.206749E+05,0.207701E+05, QSAT2B.491
*0.208659E+05,0.209621E+05,0.210584E+05,0.211554E+05,0.212524E+05, QSAT2B.492
*0.213501E+05,0.214482E+05,0.215465E+05,0.216452E+05,0.217442E+05/ QSAT2B.493
DATA (ES(IES),IES=1521,1552) / QSAT2B.494
*0.218439E+05,0.219439E+05,0.220440E+05,0.221449E+05,0.222457E+05, QSAT2B.495
*0.223473E+05,0.224494E+05,0.225514E+05,0.226542E+05,0.227571E+05, QSAT2B.496
*0.228606E+05,0.229646E+05,0.230687E+05,0.231734E+05,0.232783E+05, QSAT2B.497
*0.233839E+05,0.234898E+05,0.235960E+05,0.237027E+05,0.238097E+05, QSAT2B.498
*0.239173E+05,0.240254E+05,0.241335E+05,0.242424E+05,0.243514E+05, QSAT2B.499
*0.244611E+05,0.245712E+05,0.246814E+05,0.247923E+05,0.249034E+05, QSAT2B.500
*0.250152E+05,0.250152E+05/ QSAT2B.501
C QSAT2B.502
DO 10 I=1,NPNTS QSAT2B.503
C QSAT2B.504
C COMPUTE THE FACTOR THAT CONVERTS FROM SAT VAPOUR PRESSURE QSAT2B.505
C IN A PURE WATER SYSTEM TO SAT VAPOUR PRESSURE IN AIR, FSUBW. QSAT2B.506
C THIS FORMULA IS TAKEN FROM EQUATION A4.7 OF ADRIAN GILL'S QSAT2B.507
C BOOK; ATMOSPHERE-OCEAN DYNAMICS. NOTE THAT HIS FORMULA QSAT2B.508
C WORKS IN TERMS OF PRESSURE IN MB AND TEMPERATURE IN CELSIUS, QSAT2B.509
C SO CONVERSION OF UNITS LEADS TO THE SLIGHTLY DIFFERENT QSAT2B.510
C EQUATION USED HERE. QSAT2B.511
C QSAT2B.512
FSUBW = 1.0 + 1.0E-8*P(I)*( 4.5 + QSAT2B.513
+ 6.0E-4*( T(I) - ZERODEGC )*( T(I) - ZERODEGC) ) QSAT2B.514
C QSAT2B.515
C USE THE LOOKUP TABLE TO FIND SATURATED VAPOUR PRESSURE, QSAT2B.516
C AND STORE IT IN QS. QSAT2B.517
C QSAT2B.518
TT=MAX(T_LOW,T(I)) QSAT2B.519
TT=MIN(T_HIGH,TT) QSAT2B.520
C QSAT2B.521
ATABLE = (TT - T_LOW + DELTA_T) / DELTA_T QSAT2B.522
ITABLE = ATABLE QSAT2B.523
ATABLE = ATABLE - ITABLE QSAT2B.524
C QSAT2B.525
QS(I) = (1.0 - ATABLE)*ES(ITABLE) QSAT2B.526
+ +ATABLE*ES(ITABLE+1) QSAT2B.527
C QSAT2B.528
C QSAT2B.529
C MULTIPLY BY FSUBW TO CONVERT TO SATURATED VAPOUR PRESSURE QSAT2B.530
C IN AIR (EQUATION A4.6 OF ADRIAN GILL'S BOOK). QSAT2B.531
C QSAT2B.532
QS(I) = QS(I)*FSUBW QSAT2B.533
C QSAT2B.534
C NOW FORM THE ACCURATE EXPRESSION FOR QS, WHICH IS A QSAT2B.535
C REARRANGED VERSION OF EQUATION A4.3 OF GILL'S BOOK. QSAT2B.536
C QSAT2B.537
C NOTE THAT AT VERY LOW PRESSURES WE APPLY A FIX, TO QSAT2B.538
C PREVENT A SINGULARITY. QSAT2B.539
C QSAT2B.540
QS(I) = ( EPSILON*QS(I) ) / QSAT2B.541
+ ( MAX(P(I),QS(I)) - ONE_MINUS_EPSILON*QS(I) ) QSAT2B.542
C QSAT2B.543
10 CONTINUE QSAT2B.544
C QSAT2B.545
RETURN QSAT2B.546
END QSAT2B.547
CLL SUBROUTINE QSAT--------------------------------------------------- QSAT2B.548
CLL QSAT2B.549
CLL PURPOSE : RETURNS A SATURATION MIXING RATIO GIVEN QSAT2B.550
CLL A TEMPERATURE AND PRESSURE USING SATURATION QSAT2B.551
CLL VAPOUR PRESSURES CALCULATED USING THE QSAT2B.552
CLL GOFF-GRATCH FORMULAE, ADOPTED BY THE WMO AS QSAT2B.553
CLL TAKEN FROM LANDOLT-BORNSTEIN, 1987 NUMERICAL QSAT2B.554
CLL DATA AND FUNCTIONAL RELATIONSHIPS IN SCIENCE QSAT2B.555
CLL AND TECHNOLOGY. GROUP V/VOL 4B METEOROLOGY. QSAT2B.556
CLL PHYSICAL AND CHEMICAL PROPERTIES OF AIR, P35 QSAT2B.557
CLL QSAT2B.558
CLL VALUES IN THE LOOKUP TABLE ARE OVER WATER ABOVE QSAT2B.559
CLL 0 DEG C AND OVER ICE BELOW THIS TEMPERATURE QSAT2B.560
CLL QSAT2B.561
CLL SUITABLE FOR SINGLE COLUMN MODEL USE QSAT2B.562
CLL QSAT2B.563
CLL CODE REWORKED FOR CRAY Y-MP BY D.GREGORY AUTUMN/WINTER 1989/90 QSAT2B.564
CLL QSAT2B.565
CLL MODEL MODIFICATION HISTORY: QSAT2B.566
CLL VERSION DATE QSAT2B.567
CLL 4.5 25/06/98 Correct potential failure introduced by optimising: ARR6F405.12
CLL Extend lookup array from (1:) to (0:) to cater for ARR6F405.13
CLL special case of extreme low temperatures ARR6F405.14
CLL (.LE.T_LOW) for which the array index is rounded ARR6F405.15
CLL down due to machine precision. R.Rawlins ARR6F405.16
CLL QSAT2B.568
CLL PROGRAMMING STANDARDS : QSAT2B.569
CLL QSAT2B.570
CLL LOGICAL COMPONENTS COVERED: QSAT2B.571
CLL QSAT2B.572
CLL LOGICAL COMPONENTS COVERED: P27 QSAT2B.573
CLL QSAT2B.574
CLL DOCUMENTATION : QSAT2B.575
CLL QSAT2B.576
CLLEND----------------------------------------------------------------- QSAT2B.577
C QSAT2B.578
C*L ARGUMENTS--------------------------------------------------------- QSAT2B.579
C QSAT2B.580
SUBROUTINE QSAT (QS,T,P,NPNTS) 105QSAT2B.581
C QSAT2B.582
IMPLICIT NONE QSAT2B.583
C QSAT2B.584
C---------------------------------------------------------------------- QSAT2B.585
C MODEL CONSTANTS QSAT2B.586
C---------------------------------------------------------------------- QSAT2B.587
C QSAT2B.588
*CALL C_EPSLON QSAT2B.589
*CALL C_0_DG_C QSAT2B.590
C QSAT2B.591
C---------------------------------------------------------------------- QSAT2B.592
C LOCAL CONSTANTS QSAT2B.593
C---------------------------------------------------------------------- QSAT2B.594
C QSAT2B.595
REAL T_LOW ! LOWEST TEMPERATURE FOR WHICH LOOK-UP QSAT2B.596
! TABLE OF SATURATION WATER VAPOUR QSAT2B.597
! PRESSURE IS VALID (K) QSAT2B.598
C QSAT2B.599
REAL T_HIGH ! HIGHEST TEMPERATURE FOR WHICH LOOK-UP QSAT2B.600
! TABLE OF SATURATION WATER VAPOUR QSAT2B.601
! PRESSURES IS VALID (K) QSAT2B.602
C QSAT2B.603
REAL DELTA_T ! TEMPERATURE INCREMENT OF THE LOOK-UP QSAT2B.604
! TABLE OF SATURATION VAPOUR PRESSURES QSAT2B.605
C QSAT2B.606
INTEGER N ! SIZE OF LOOK-UP TABLE OF SATURATION QSAT2B.607
! WATER VAPOUR PRESSURES QSAT2B.608
C QSAT2B.609
PARAMETER ( T_LOW = 183.15, QSAT2B.610
* T_HIGH = 338.15, QSAT2B.611
* DELTA_T = 0.1, QSAT2B.612
* N = ((T_HIGH - T_LOW + (DELTA_T*0.5))/DELTA_T) + 1.0 QSAT2B.613
* ) ! gives N=1551 ARR6F405.17
C QSAT2B.615
C QSAT2B.616
C---------------------------------------------------------------------- QSAT2B.617
C VECTOR LENGTHS AND LOOP COUNTERS QSAT2B.618
C---------------------------------------------------------------------- QSAT2B.619
C QSAT2B.620
INTEGER NPNTS ! VECTOR LENGTH QSAT2B.621
C QSAT2B.622
INTEGER I ! LOOP COUNTER QSAT2B.623
C QSAT2B.624
INTEGER IES ! LOOP COUNTER FOR DATA STATEMENT QSAT2B.625
! LOOK-UP TABLE QSAT2B.626
C QSAT2B.627
C QSAT2B.628
C---------------------------------------------------------------------- QSAT2B.629
C VARIABLES WHICH ARE INPUT QSAT2B.630
C---------------------------------------------------------------------- QSAT2B.631
C QSAT2B.632
REAL T(NPNTS) ! IN TEMPERATURE (K) QSAT2B.633
C QSAT2B.634
REAL P(NPNTS) ! IN PRESSURE (PA) QSAT2B.635
C QSAT2B.636
C QSAT2B.637
C---------------------------------------------------------------------- QSAT2B.638
C VARIABLES WHICH ARE OUTPUT QSAT2B.639
C---------------------------------------------------------------------- QSAT2B.640
C QSAT2B.641
REAL QS(NPNTS) ! OUT SATURATION MIXING RATIO AT TEMPERATURE QSAT2B.642
! T AND PRESSURE P (KG/KG) QSAT2B.643
C QSAT2B.644
C QSAT2B.645
C---------------------------------------------------------------------- QSAT2B.646
C VARIABLES WHICH ARE DEFINED LOCALLY QSAT2B.647
C---------------------------------------------------------------------- QSAT2B.648
C QSAT2B.649
REAL ES(0:N+1) ! TABLE OF SATURATION WATER VAPOUR ARR6F405.18
! PRESSURE (PA) - SET BY DATA STATEMENT QSAT2B.651
! CALCULATED FROM THE GOFF-GRATCH FORMULAE QSAT2B.652
! AS TAKEN FROM LANDOLT-BORNSTEIN, 1987 QSAT2B.653
! NUMERICAL DATA AND FUNCTIONAL RELATIONSHIPS QSAT2B.654
! IN SCIENCE AND TECHNOLOGY. GROUP V/ VOL 4B QSAT2B.655
! METEOROLOGY. PHYSICAL AND CHEMICAL QSAT2B.656
! PROPERTIES OF AIR, P35 QSAT2B.657
C QSAT2B.658
REAL ATABLE ! WORK VARIABLES QSAT2B.659
C QSAT2B.660
INTEGER ITABLE ! WORK VARIABLES QSAT2B.661
C QSAT2B.662
C VARIABLES INTRODUCED BY DLR. QSAT2B.663
C QSAT2B.664
REAL FSUBW ! FACTOR THAT CONVERTS FROM SAT VAPOUR QSAT2B.665
! PRESSURE IN A PURE WATER SYSTEM TO QSAT2B.666
! SAT VAPOUR PRESSURE IN AIR. QSAT2B.667
C QSAT2B.668
REAL ONE_MINUS_EPSILON ! ONE MINUS THE RATIO OF THE MOLECULAR QSAT2B.669
! WEIGHTS OF WATER AND DRY AIR QSAT2B.670
REAL TT QSAT2B.671
C QSAT2B.672
PARAMETER(ONE_MINUS_EPSILON = 1.0 - EPSILON) QSAT2B.673
C QSAT2B.674
C*--------------------------------------------------------------------- QSAT2B.675
CL QSAT2B.676
CL--------------------------------------------------------------------- QSAT2B.677
CL NO SIGNIFICANT STRUCTURE QSAT2B.678
CL--------------------------------------------------------------------- QSAT2B.679
CL QSAT2B.680
C QSAT2B.681
C---------------------------------------------------------------------- QSAT2B.682
C SATURATION WATER VAPOUR PRESSURE QSAT2B.683
C QSAT2B.684
C ABOVE 0 DEG C VALUES ARE OVER WATER QSAT2B.685
C QSAT2B.686
C BELOW 0 DEC C VALUES ARE OVER ICE QSAT2B.687
C---------------------------------------------------------------------- QSAT2B.688
C QSAT2B.689
C Note: 0 element is a repeat of 1st element to cater for special case ARR6F405.19
C of low temperatures (.LE.T_LOW) for which the array index is ARR6F405.20
C rounded down due to machine precision. ARR6F405.21
DATA (ES(IES),IES= 0, 95) / 0.966483E-02, ARR6F405.22
*0.966483E-02,0.984279E-02,0.100240E-01,0.102082E-01,0.103957E-01, QSAT2B.691
*0.105865E-01,0.107803E-01,0.109777E-01,0.111784E-01,0.113825E-01, QSAT2B.692
*0.115902E-01,0.118016E-01,0.120164E-01,0.122348E-01,0.124572E-01, QSAT2B.693
*0.126831E-01,0.129132E-01,0.131470E-01,0.133846E-01,0.136264E-01, QSAT2B.694
*0.138724E-01,0.141225E-01,0.143771E-01,0.146356E-01,0.148985E-01, QSAT2B.695
*0.151661E-01,0.154379E-01,0.157145E-01,0.159958E-01,0.162817E-01, QSAT2B.696
*0.165725E-01,0.168680E-01,0.171684E-01,0.174742E-01,0.177847E-01, QSAT2B.697
*0.181008E-01,0.184216E-01,0.187481E-01,0.190801E-01,0.194175E-01, QSAT2B.698
*0.197608E-01,0.201094E-01,0.204637E-01,0.208242E-01,0.211906E-01, QSAT2B.699
*0.215631E-01,0.219416E-01,0.223263E-01,0.227172E-01,0.231146E-01, QSAT2B.700
*0.235188E-01,0.239296E-01,0.243465E-01,0.247708E-01,0.252019E-01, QSAT2B.701
*0.256405E-01,0.260857E-01,0.265385E-01,0.269979E-01,0.274656E-01, QSAT2B.702
*0.279405E-01,0.284232E-01,0.289142E-01,0.294124E-01,0.299192E-01, QSAT2B.703
*0.304341E-01,0.309571E-01,0.314886E-01,0.320285E-01,0.325769E-01, QSAT2B.704
*0.331348E-01,0.337014E-01,0.342771E-01,0.348618E-01,0.354557E-01, QSAT2B.705
*0.360598E-01,0.366727E-01,0.372958E-01,0.379289E-01,0.385717E-01, QSAT2B.706
*0.392248E-01,0.398889E-01,0.405633E-01,0.412474E-01,0.419430E-01, QSAT2B.707
*0.426505E-01,0.433678E-01,0.440974E-01,0.448374E-01,0.455896E-01, QSAT2B.708
*0.463545E-01,0.471303E-01,0.479191E-01,0.487190E-01,0.495322E-01/ QSAT2B.709
DATA (ES(IES),IES= 96,190) / QSAT2B.710
*0.503591E-01,0.511977E-01,0.520490E-01,0.529145E-01,0.537931E-01, QSAT2B.711
*0.546854E-01,0.555924E-01,0.565119E-01,0.574467E-01,0.583959E-01, QSAT2B.712
*0.593592E-01,0.603387E-01,0.613316E-01,0.623409E-01,0.633655E-01, QSAT2B.713
*0.644053E-01,0.654624E-01,0.665358E-01,0.676233E-01,0.687302E-01, QSAT2B.714
*0.698524E-01,0.709929E-01,0.721490E-01,0.733238E-01,0.745180E-01, QSAT2B.715
*0.757281E-01,0.769578E-01,0.782061E-01,0.794728E-01,0.807583E-01, QSAT2B.716
*0.820647E-01,0.833905E-01,0.847358E-01,0.861028E-01,0.874882E-01, QSAT2B.717
*0.888957E-01,0.903243E-01,0.917736E-01,0.932464E-01,0.947407E-01, QSAT2B.718
*0.962571E-01,0.977955E-01,0.993584E-01,0.100942E+00,0.102551E+00, QSAT2B.719
*0.104186E+00,0.105842E+00,0.107524E+00,0.109231E+00,0.110963E+00, QSAT2B.720
*0.112722E+00,0.114506E+00,0.116317E+00,0.118153E+00,0.120019E+00, QSAT2B.721
*0.121911E+00,0.123831E+00,0.125778E+00,0.127755E+00,0.129761E+00, QSAT2B.722
*0.131796E+00,0.133863E+00,0.135956E+00,0.138082E+00,0.140241E+00, QSAT2B.723
*0.142428E+00,0.144649E+00,0.146902E+00,0.149190E+00,0.151506E+00, QSAT2B.724
*0.153859E+00,0.156245E+00,0.158669E+00,0.161126E+00,0.163618E+00, QSAT2B.725
*0.166145E+00,0.168711E+00,0.171313E+00,0.173951E+00,0.176626E+00, QSAT2B.726
*0.179342E+00,0.182096E+00,0.184893E+00,0.187724E+00,0.190600E+00, QSAT2B.727
*0.193518E+00,0.196473E+00,0.199474E+00,0.202516E+00,0.205604E+00, QSAT2B.728
*0.208730E+00,0.211905E+00,0.215127E+00,0.218389E+00,0.221701E+00/ QSAT2B.729
DATA (ES(IES),IES=191,285) / QSAT2B.730
*0.225063E+00,0.228466E+00,0.231920E+00,0.235421E+00,0.238976E+00, QSAT2B.731
*0.242580E+00,0.246232E+00,0.249933E+00,0.253691E+00,0.257499E+00, QSAT2B.732
*0.261359E+00,0.265278E+00,0.269249E+00,0.273274E+00,0.277358E+00, QSAT2B.733
*0.281498E+00,0.285694E+00,0.289952E+00,0.294268E+00,0.298641E+00, QSAT2B.734
*0.303078E+00,0.307577E+00,0.312135E+00,0.316753E+00,0.321440E+00, QSAT2B.735
*0.326196E+00,0.331009E+00,0.335893E+00,0.340842E+00,0.345863E+00, QSAT2B.736
*0.350951E+00,0.356106E+00,0.361337E+00,0.366636E+00,0.372006E+00, QSAT2B.737
*0.377447E+00,0.382966E+00,0.388567E+00,0.394233E+00,0.399981E+00, QSAT2B.738
*0.405806E+00,0.411714E+00,0.417699E+00,0.423772E+00,0.429914E+00, QSAT2B.739
*0.436145E+00,0.442468E+00,0.448862E+00,0.455359E+00,0.461930E+00, QSAT2B.740
*0.468596E+00,0.475348E+00,0.482186E+00,0.489124E+00,0.496160E+00, QSAT2B.741
*0.503278E+00,0.510497E+00,0.517808E+00,0.525224E+00,0.532737E+00, QSAT2B.742
*0.540355E+00,0.548059E+00,0.555886E+00,0.563797E+00,0.571825E+00, QSAT2B.743
*0.579952E+00,0.588198E+00,0.596545E+00,0.605000E+00,0.613572E+00, QSAT2B.744
*0.622255E+00,0.631059E+00,0.639962E+00,0.649003E+00,0.658144E+00, QSAT2B.745
*0.667414E+00,0.676815E+00,0.686317E+00,0.695956E+00,0.705728E+00, QSAT2B.746
*0.715622E+00,0.725641E+00,0.735799E+00,0.746082E+00,0.756495E+00, QSAT2B.747
*0.767052E+00,0.777741E+00,0.788576E+00,0.799549E+00,0.810656E+00, QSAT2B.748
*0.821914E+00,0.833314E+00,0.844854E+00,0.856555E+00,0.868415E+00/ QSAT2B.749
DATA (ES(IES),IES=286,380) / QSAT2B.750
*0.880404E+00,0.892575E+00,0.904877E+00,0.917350E+00,0.929974E+00, QSAT2B.751
*0.942771E+00,0.955724E+00,0.968837E+00,0.982127E+00,0.995600E+00, QSAT2B.752
*0.100921E+01,0.102304E+01,0.103700E+01,0.105116E+01,0.106549E+01, QSAT2B.753
*0.108002E+01,0.109471E+01,0.110962E+01,0.112469E+01,0.113995E+01, QSAT2B.754
*0.115542E+01,0.117107E+01,0.118693E+01,0.120298E+01,0.121923E+01, QSAT2B.755
*0.123569E+01,0.125234E+01,0.126923E+01,0.128631E+01,0.130362E+01, QSAT2B.756
*0.132114E+01,0.133887E+01,0.135683E+01,0.137500E+01,0.139342E+01, QSAT2B.757
*0.141205E+01,0.143091E+01,0.145000E+01,0.146933E+01,0.148892E+01, QSAT2B.758
*0.150874E+01,0.152881E+01,0.154912E+01,0.156970E+01,0.159049E+01, QSAT2B.759
*0.161159E+01,0.163293E+01,0.165452E+01,0.167640E+01,0.169852E+01, QSAT2B.760
*0.172091E+01,0.174359E+01,0.176653E+01,0.178977E+01,0.181332E+01, QSAT2B.761
*0.183709E+01,0.186119E+01,0.188559E+01,0.191028E+01,0.193524E+01, QSAT2B.762
*0.196054E+01,0.198616E+01,0.201208E+01,0.203829E+01,0.206485E+01, QSAT2B.763
*0.209170E+01,0.211885E+01,0.214637E+01,0.217424E+01,0.220242E+01, QSAT2B.764
*0.223092E+01,0.225979E+01,0.228899E+01,0.231855E+01,0.234845E+01, QSAT2B.765
*0.237874E+01,0.240937E+01,0.244040E+01,0.247176E+01,0.250349E+01, QSAT2B.766
*0.253560E+01,0.256814E+01,0.260099E+01,0.263431E+01,0.266800E+01, QSAT2B.767
*0.270207E+01,0.273656E+01,0.277145E+01,0.280671E+01,0.284248E+01, QSAT2B.768
*0.287859E+01,0.291516E+01,0.295219E+01,0.298962E+01,0.302746E+01/ QSAT2B.769
DATA (ES(IES),IES=381,475) / QSAT2B.770
*0.306579E+01,0.310454E+01,0.314377E+01,0.318351E+01,0.322360E+01, QSAT2B.771
*0.326427E+01,0.330538E+01,0.334694E+01,0.338894E+01,0.343155E+01, QSAT2B.772
*0.347456E+01,0.351809E+01,0.356216E+01,0.360673E+01,0.365184E+01, QSAT2B.773
*0.369744E+01,0.374352E+01,0.379018E+01,0.383743E+01,0.388518E+01, QSAT2B.774
*0.393344E+01,0.398230E+01,0.403177E+01,0.408175E+01,0.413229E+01, QSAT2B.775
*0.418343E+01,0.423514E+01,0.428746E+01,0.434034E+01,0.439389E+01, QSAT2B.776
*0.444808E+01,0.450276E+01,0.455820E+01,0.461423E+01,0.467084E+01, QSAT2B.777
*0.472816E+01,0.478607E+01,0.484468E+01,0.490393E+01,0.496389E+01, QSAT2B.778
*0.502446E+01,0.508580E+01,0.514776E+01,0.521047E+01,0.527385E+01, QSAT2B.779
*0.533798E+01,0.540279E+01,0.546838E+01,0.553466E+01,0.560173E+01, QSAT2B.780
*0.566949E+01,0.573807E+01,0.580750E+01,0.587749E+01,0.594846E+01, QSAT2B.781
*0.602017E+01,0.609260E+01,0.616591E+01,0.623995E+01,0.631490E+01, QSAT2B.782
*0.639061E+01,0.646723E+01,0.654477E+01,0.662293E+01,0.670220E+01, QSAT2B.783
*0.678227E+01,0.686313E+01,0.694495E+01,0.702777E+01,0.711142E+01, QSAT2B.784
*0.719592E+01,0.728140E+01,0.736790E+01,0.745527E+01,0.754352E+01, QSAT2B.785
*0.763298E+01,0.772316E+01,0.781442E+01,0.790676E+01,0.800001E+01, QSAT2B.786
*0.809435E+01,0.818967E+01,0.828606E+01,0.838343E+01,0.848194E+01, QSAT2B.787
*0.858144E+01,0.868207E+01,0.878392E+01,0.888673E+01,0.899060E+01, QSAT2B.788
*0.909567E+01,0.920172E+01,0.930909E+01,0.941765E+01,0.952730E+01/ QSAT2B.789
DATA (ES(IES),IES=476,570) / QSAT2B.790
*0.963821E+01,0.975022E+01,0.986352E+01,0.997793E+01,0.100937E+02, QSAT2B.791
*0.102105E+02,0.103287E+02,0.104481E+02,0.105688E+02,0.106909E+02, QSAT2B.792
*0.108143E+02,0.109387E+02,0.110647E+02,0.111921E+02,0.113207E+02, QSAT2B.793
*0.114508E+02,0.115821E+02,0.117149E+02,0.118490E+02,0.119847E+02, QSAT2B.794
*0.121216E+02,0.122601E+02,0.124002E+02,0.125416E+02,0.126846E+02, QSAT2B.795
*0.128290E+02,0.129747E+02,0.131224E+02,0.132712E+02,0.134220E+02, QSAT2B.796
*0.135742E+02,0.137278E+02,0.138831E+02,0.140403E+02,0.141989E+02, QSAT2B.797
*0.143589E+02,0.145211E+02,0.146845E+02,0.148501E+02,0.150172E+02, QSAT2B.798
*0.151858E+02,0.153564E+02,0.155288E+02,0.157029E+02,0.158786E+02, QSAT2B.799
*0.160562E+02,0.162358E+02,0.164174E+02,0.166004E+02,0.167858E+02, QSAT2B.800
*0.169728E+02,0.171620E+02,0.173528E+02,0.175455E+02,0.177406E+02, QSAT2B.801
*0.179372E+02,0.181363E+02,0.183372E+02,0.185400E+02,0.187453E+02, QSAT2B.802
*0.189523E+02,0.191613E+02,0.193728E+02,0.195866E+02,0.198024E+02, QSAT2B.803
*0.200200E+02,0.202401E+02,0.204626E+02,0.206871E+02,0.209140E+02, QSAT2B.804
*0.211430E+02,0.213744E+02,0.216085E+02,0.218446E+02,0.220828E+02, QSAT2B.805
*0.223241E+02,0.225671E+02,0.228132E+02,0.230615E+02,0.233120E+02, QSAT2B.806
*0.235651E+02,0.238211E+02,0.240794E+02,0.243404E+02,0.246042E+02, QSAT2B.807
*0.248704E+02,0.251390E+02,0.254109E+02,0.256847E+02,0.259620E+02, QSAT2B.808
*0.262418E+02,0.265240E+02,0.268092E+02,0.270975E+02,0.273883E+02/ QSAT2B.809
DATA (ES(IES),IES=571,665) / QSAT2B.810
*0.276822E+02,0.279792E+02,0.282789E+02,0.285812E+02,0.288867E+02, QSAT2B.811
*0.291954E+02,0.295075E+02,0.298222E+02,0.301398E+02,0.304606E+02, QSAT2B.812
*0.307848E+02,0.311119E+02,0.314424E+02,0.317763E+02,0.321133E+02, QSAT2B.813
*0.324536E+02,0.327971E+02,0.331440E+02,0.334940E+02,0.338475E+02, QSAT2B.814
*0.342050E+02,0.345654E+02,0.349295E+02,0.352975E+02,0.356687E+02, QSAT2B.815
*0.360430E+02,0.364221E+02,0.368042E+02,0.371896E+02,0.375790E+02, QSAT2B.816
*0.379725E+02,0.383692E+02,0.387702E+02,0.391744E+02,0.395839E+02, QSAT2B.817
*0.399958E+02,0.404118E+02,0.408325E+02,0.412574E+02,0.416858E+02, QSAT2B.818
*0.421188E+02,0.425551E+02,0.429962E+02,0.434407E+02,0.438910E+02, QSAT2B.819
*0.443439E+02,0.448024E+02,0.452648E+02,0.457308E+02,0.462018E+02, QSAT2B.820
*0.466775E+02,0.471582E+02,0.476428E+02,0.481313E+02,0.486249E+02, QSAT2B.821
*0.491235E+02,0.496272E+02,0.501349E+02,0.506479E+02,0.511652E+02, QSAT2B.822
*0.516876E+02,0.522142E+02,0.527474E+02,0.532836E+02,0.538266E+02, QSAT2B.823
*0.543737E+02,0.549254E+02,0.554839E+02,0.560456E+02,0.566142E+02, QSAT2B.824
*0.571872E+02,0.577662E+02,0.583498E+02,0.589392E+02,0.595347E+02, QSAT2B.825
*0.601346E+02,0.607410E+02,0.613519E+02,0.619689E+02,0.625922E+02, QSAT2B.826
*0.632204E+02,0.638550E+02,0.644959E+02,0.651418E+02,0.657942E+02, QSAT2B.827
*0.664516E+02,0.671158E+02,0.677864E+02,0.684624E+02,0.691451E+02, QSAT2B.828
*0.698345E+02,0.705293E+02,0.712312E+02,0.719398E+02,0.726542E+02/ QSAT2B.829
DATA (ES(IES),IES=666,760) / QSAT2B.830
*0.733754E+02,0.741022E+02,0.748363E+02,0.755777E+02,0.763247E+02, QSAT2B.831
*0.770791E+02,0.778394E+02,0.786088E+02,0.793824E+02,0.801653E+02, QSAT2B.832
*0.809542E+02,0.817509E+02,0.825536E+02,0.833643E+02,0.841828E+02, QSAT2B.833
*0.850076E+02,0.858405E+02,0.866797E+02,0.875289E+02,0.883827E+02, QSAT2B.834
*0.892467E+02,0.901172E+02,0.909962E+02,0.918818E+02,0.927760E+02, QSAT2B.835
*0.936790E+02,0.945887E+02,0.955071E+02,0.964346E+02,0.973689E+02, QSAT2B.836
*0.983123E+02,0.992648E+02,0.100224E+03,0.101193E+03,0.102169E+03, QSAT2B.837
*0.103155E+03,0.104150E+03,0.105152E+03,0.106164E+03,0.107186E+03, QSAT2B.838
*0.108217E+03,0.109256E+03,0.110303E+03,0.111362E+03,0.112429E+03, QSAT2B.839
*0.113503E+03,0.114588E+03,0.115684E+03,0.116789E+03,0.117903E+03, QSAT2B.840
*0.119028E+03,0.120160E+03,0.121306E+03,0.122460E+03,0.123623E+03, QSAT2B.841
*0.124796E+03,0.125981E+03,0.127174E+03,0.128381E+03,0.129594E+03, QSAT2B.842
*0.130822E+03,0.132058E+03,0.133306E+03,0.134563E+03,0.135828E+03, QSAT2B.843
*0.137109E+03,0.138402E+03,0.139700E+03,0.141017E+03,0.142338E+03, QSAT2B.844
*0.143676E+03,0.145025E+03,0.146382E+03,0.147753E+03,0.149133E+03, QSAT2B.845
*0.150529E+03,0.151935E+03,0.153351E+03,0.154783E+03,0.156222E+03, QSAT2B.846
*0.157678E+03,0.159148E+03,0.160624E+03,0.162117E+03,0.163621E+03, QSAT2B.847
*0.165142E+03,0.166674E+03,0.168212E+03,0.169772E+03,0.171340E+03, QSAT2B.848
*0.172921E+03,0.174522E+03,0.176129E+03,0.177755E+03,0.179388E+03/ QSAT2B.849
DATA (ES(IES),IES=761,855) / QSAT2B.850
*0.181040E+03,0.182707E+03,0.184382E+03,0.186076E+03,0.187782E+03, QSAT2B.851
*0.189503E+03,0.191240E+03,0.192989E+03,0.194758E+03,0.196535E+03, QSAT2B.852
*0.198332E+03,0.200141E+03,0.201963E+03,0.203805E+03,0.205656E+03, QSAT2B.853
*0.207532E+03,0.209416E+03,0.211317E+03,0.213236E+03,0.215167E+03, QSAT2B.854
*0.217121E+03,0.219087E+03,0.221067E+03,0.223064E+03,0.225080E+03, QSAT2B.855
*0.227113E+03,0.229160E+03,0.231221E+03,0.233305E+03,0.235403E+03, QSAT2B.856
*0.237520E+03,0.239655E+03,0.241805E+03,0.243979E+03,0.246163E+03, QSAT2B.857
*0.248365E+03,0.250593E+03,0.252830E+03,0.255093E+03,0.257364E+03, QSAT2B.858
*0.259667E+03,0.261979E+03,0.264312E+03,0.266666E+03,0.269034E+03, QSAT2B.859
*0.271430E+03,0.273841E+03,0.276268E+03,0.278722E+03,0.281185E+03, QSAT2B.860
*0.283677E+03,0.286190E+03,0.288714E+03,0.291266E+03,0.293834E+03, QSAT2B.861
*0.296431E+03,0.299045E+03,0.301676E+03,0.304329E+03,0.307006E+03, QSAT2B.862
*0.309706E+03,0.312423E+03,0.315165E+03,0.317930E+03,0.320705E+03, QSAT2B.863
*0.323519E+03,0.326350E+03,0.329199E+03,0.332073E+03,0.334973E+03, QSAT2B.864
*0.337897E+03,0.340839E+03,0.343800E+03,0.346794E+03,0.349806E+03, QSAT2B.865
*0.352845E+03,0.355918E+03,0.358994E+03,0.362112E+03,0.365242E+03, QSAT2B.866
*0.368407E+03,0.371599E+03,0.374802E+03,0.378042E+03,0.381293E+03, QSAT2B.867
*0.384588E+03,0.387904E+03,0.391239E+03,0.394604E+03,0.397988E+03, QSAT2B.868
*0.401411E+03,0.404862E+03,0.408326E+03,0.411829E+03,0.415352E+03/ QSAT2B.869
DATA (ES(IES),IES=856,950) / QSAT2B.870
*0.418906E+03,0.422490E+03,0.426095E+03,0.429740E+03,0.433398E+03, QSAT2B.871
*0.437097E+03,0.440827E+03,0.444570E+03,0.448354E+03,0.452160E+03, QSAT2B.872
*0.455999E+03,0.459870E+03,0.463765E+03,0.467702E+03,0.471652E+03, QSAT2B.873
*0.475646E+03,0.479674E+03,0.483715E+03,0.487811E+03,0.491911E+03, QSAT2B.874
*0.496065E+03,0.500244E+03,0.504448E+03,0.508698E+03,0.512961E+03, QSAT2B.875
*0.517282E+03,0.521617E+03,0.525989E+03,0.530397E+03,0.534831E+03, QSAT2B.876
*0.539313E+03,0.543821E+03,0.548355E+03,0.552938E+03,0.557549E+03, QSAT2B.877
*0.562197E+03,0.566884E+03,0.571598E+03,0.576351E+03,0.581131E+03, QSAT2B.878
*0.585963E+03,0.590835E+03,0.595722E+03,0.600663E+03,0.605631E+03, QSAT2B.879
*0.610641E+03,0.615151E+03,0.619625E+03,0.624140E+03,0.628671E+03, QSAT2B.880
*0.633243E+03,0.637845E+03,0.642465E+03,0.647126E+03,0.651806E+03, QSAT2B.881
*0.656527E+03,0.661279E+03,0.666049E+03,0.670861E+03,0.675692E+03, QSAT2B.882
*0.680566E+03,0.685471E+03,0.690396E+03,0.695363E+03,0.700350E+03, QSAT2B.883
*0.705381E+03,0.710444E+03,0.715527E+03,0.720654E+03,0.725801E+03, QSAT2B.884
*0.730994E+03,0.736219E+03,0.741465E+03,0.746756E+03,0.752068E+03, QSAT2B.885
*0.757426E+03,0.762819E+03,0.768231E+03,0.773692E+03,0.779172E+03, QSAT2B.886
*0.784701E+03,0.790265E+03,0.795849E+03,0.801483E+03,0.807137E+03, QSAT2B.887
*0.812842E+03,0.818582E+03,0.824343E+03,0.830153E+03,0.835987E+03, QSAT2B.888
*0.841871E+03,0.847791E+03,0.853733E+03,0.859727E+03,0.865743E+03/ QSAT2B.889
DATA (ES(IES),IES=951,1045) / QSAT2B.890
*0.871812E+03,0.877918E+03,0.884046E+03,0.890228E+03,0.896433E+03, QSAT2B.891
*0.902690E+03,0.908987E+03,0.915307E+03,0.921681E+03,0.928078E+03, QSAT2B.892
*0.934531E+03,0.941023E+03,0.947539E+03,0.954112E+03,0.960708E+03, QSAT2B.893
*0.967361E+03,0.974053E+03,0.980771E+03,0.987545E+03,0.994345E+03, QSAT2B.894
*0.100120E+04,0.100810E+04,0.101502E+04,0.102201E+04,0.102902E+04, QSAT2B.895
*0.103608E+04,0.104320E+04,0.105033E+04,0.105753E+04,0.106475E+04, QSAT2B.896
*0.107204E+04,0.107936E+04,0.108672E+04,0.109414E+04,0.110158E+04, QSAT2B.897
*0.110908E+04,0.111663E+04,0.112421E+04,0.113185E+04,0.113952E+04, QSAT2B.898
*0.114725E+04,0.115503E+04,0.116284E+04,0.117071E+04,0.117861E+04, QSAT2B.899
*0.118658E+04,0.119459E+04,0.120264E+04,0.121074E+04,0.121888E+04, QSAT2B.900
*0.122709E+04,0.123534E+04,0.124362E+04,0.125198E+04,0.126036E+04, QSAT2B.901
*0.126881E+04,0.127731E+04,0.128584E+04,0.129444E+04,0.130307E+04, QSAT2B.902
*0.131177E+04,0.132053E+04,0.132931E+04,0.133817E+04,0.134705E+04, QSAT2B.903
*0.135602E+04,0.136503E+04,0.137407E+04,0.138319E+04,0.139234E+04, QSAT2B.904
*0.140156E+04,0.141084E+04,0.142015E+04,0.142954E+04,0.143896E+04, QSAT2B.905
*0.144845E+04,0.145800E+04,0.146759E+04,0.147725E+04,0.148694E+04, QSAT2B.906
*0.149672E+04,0.150655E+04,0.151641E+04,0.152635E+04,0.153633E+04, QSAT2B.907
*0.154639E+04,0.155650E+04,0.156665E+04,0.157688E+04,0.158715E+04, QSAT2B.908
*0.159750E+04,0.160791E+04,0.161836E+04,0.162888E+04,0.163945E+04/ QSAT2B.909
DATA (ES(IES),IES=1046,1140) / QSAT2B.910
*0.165010E+04,0.166081E+04,0.167155E+04,0.168238E+04,0.169325E+04, QSAT2B.911
*0.170420E+04,0.171522E+04,0.172627E+04,0.173741E+04,0.174859E+04, QSAT2B.912
*0.175986E+04,0.177119E+04,0.178256E+04,0.179402E+04,0.180552E+04, QSAT2B.913
*0.181711E+04,0.182877E+04,0.184046E+04,0.185224E+04,0.186407E+04, QSAT2B.914
*0.187599E+04,0.188797E+04,0.190000E+04,0.191212E+04,0.192428E+04, QSAT2B.915
*0.193653E+04,0.194886E+04,0.196122E+04,0.197368E+04,0.198618E+04, QSAT2B.916
*0.199878E+04,0.201145E+04,0.202416E+04,0.203698E+04,0.204983E+04, QSAT2B.917
*0.206278E+04,0.207580E+04,0.208887E+04,0.210204E+04,0.211525E+04, QSAT2B.918
*0.212856E+04,0.214195E+04,0.215538E+04,0.216892E+04,0.218249E+04, QSAT2B.919
*0.219618E+04,0.220994E+04,0.222375E+04,0.223766E+04,0.225161E+04, QSAT2B.920
*0.226567E+04,0.227981E+04,0.229399E+04,0.230829E+04,0.232263E+04, QSAT2B.921
*0.233708E+04,0.235161E+04,0.236618E+04,0.238087E+04,0.239560E+04, QSAT2B.922
*0.241044E+04,0.242538E+04,0.244035E+04,0.245544E+04,0.247057E+04, QSAT2B.923
*0.248583E+04,0.250116E+04,0.251654E+04,0.253204E+04,0.254759E+04, QSAT2B.924
*0.256325E+04,0.257901E+04,0.259480E+04,0.261073E+04,0.262670E+04, QSAT2B.925
*0.264279E+04,0.265896E+04,0.267519E+04,0.269154E+04,0.270794E+04, QSAT2B.926
*0.272447E+04,0.274108E+04,0.275774E+04,0.277453E+04,0.279137E+04, QSAT2B.927
*0.280834E+04,0.282540E+04,0.284251E+04,0.285975E+04,0.287704E+04, QSAT2B.928
*0.289446E+04,0.291198E+04,0.292954E+04,0.294725E+04,0.296499E+04/ QSAT2B.929
DATA (ES(IES),IES=1141,1235) / QSAT2B.930
*0.298288E+04,0.300087E+04,0.301890E+04,0.303707E+04,0.305529E+04, QSAT2B.931
*0.307365E+04,0.309211E+04,0.311062E+04,0.312927E+04,0.314798E+04, QSAT2B.932
*0.316682E+04,0.318577E+04,0.320477E+04,0.322391E+04,0.324310E+04, QSAT2B.933
*0.326245E+04,0.328189E+04,0.330138E+04,0.332103E+04,0.334073E+04, QSAT2B.934
*0.336058E+04,0.338053E+04,0.340054E+04,0.342069E+04,0.344090E+04, QSAT2B.935
*0.346127E+04,0.348174E+04,0.350227E+04,0.352295E+04,0.354369E+04, QSAT2B.936
*0.356458E+04,0.358559E+04,0.360664E+04,0.362787E+04,0.364914E+04, QSAT2B.937
*0.367058E+04,0.369212E+04,0.371373E+04,0.373548E+04,0.375731E+04, QSAT2B.938
*0.377929E+04,0.380139E+04,0.382355E+04,0.384588E+04,0.386826E+04, QSAT2B.939
*0.389081E+04,0.391348E+04,0.393620E+04,0.395910E+04,0.398205E+04, QSAT2B.940
*0.400518E+04,0.402843E+04,0.405173E+04,0.407520E+04,0.409875E+04, QSAT2B.941
*0.412246E+04,0.414630E+04,0.417019E+04,0.419427E+04,0.421840E+04, QSAT2B.942
*0.424272E+04,0.426715E+04,0.429165E+04,0.431634E+04,0.434108E+04, QSAT2B.943
*0.436602E+04,0.439107E+04,0.441618E+04,0.444149E+04,0.446685E+04, QSAT2B.944
*0.449241E+04,0.451810E+04,0.454385E+04,0.456977E+04,0.459578E+04, QSAT2B.945
*0.462197E+04,0.464830E+04,0.467468E+04,0.470127E+04,0.472792E+04, QSAT2B.946
*0.475477E+04,0.478175E+04,0.480880E+04,0.483605E+04,0.486336E+04, QSAT2B.947
*0.489087E+04,0.491853E+04,0.494623E+04,0.497415E+04,0.500215E+04, QSAT2B.948
*0.503034E+04,0.505867E+04,0.508707E+04,0.511568E+04,0.514436E+04/ QSAT2B.949
DATA (ES(IES),IES=1236,1330) / QSAT2B.950
*0.517325E+04,0.520227E+04,0.523137E+04,0.526068E+04,0.529005E+04, QSAT2B.951
*0.531965E+04,0.534939E+04,0.537921E+04,0.540923E+04,0.543932E+04, QSAT2B.952
*0.546965E+04,0.550011E+04,0.553064E+04,0.556139E+04,0.559223E+04, QSAT2B.953
*0.562329E+04,0.565449E+04,0.568577E+04,0.571727E+04,0.574884E+04, QSAT2B.954
*0.578064E+04,0.581261E+04,0.584464E+04,0.587692E+04,0.590924E+04, QSAT2B.955
*0.594182E+04,0.597455E+04,0.600736E+04,0.604039E+04,0.607350E+04, QSAT2B.956
*0.610685E+04,0.614036E+04,0.617394E+04,0.620777E+04,0.624169E+04, QSAT2B.957
*0.627584E+04,0.631014E+04,0.634454E+04,0.637918E+04,0.641390E+04, QSAT2B.958
*0.644887E+04,0.648400E+04,0.651919E+04,0.655467E+04,0.659021E+04, QSAT2B.959
*0.662599E+04,0.666197E+04,0.669800E+04,0.673429E+04,0.677069E+04, QSAT2B.960
*0.680735E+04,0.684415E+04,0.688104E+04,0.691819E+04,0.695543E+04, QSAT2B.961
*0.699292E+04,0.703061E+04,0.706837E+04,0.710639E+04,0.714451E+04, QSAT2B.962
*0.718289E+04,0.722143E+04,0.726009E+04,0.729903E+04,0.733802E+04, QSAT2B.963
*0.737729E+04,0.741676E+04,0.745631E+04,0.749612E+04,0.753602E+04, QSAT2B.964
*0.757622E+04,0.761659E+04,0.765705E+04,0.769780E+04,0.773863E+04, QSAT2B.965
*0.777975E+04,0.782106E+04,0.786246E+04,0.790412E+04,0.794593E+04, QSAT2B.966
*0.798802E+04,0.803028E+04,0.807259E+04,0.811525E+04,0.815798E+04, QSAT2B.967
*0.820102E+04,0.824427E+04,0.828757E+04,0.833120E+04,0.837493E+04, QSAT2B.968
*0.841895E+04,0.846313E+04,0.850744E+04,0.855208E+04,0.859678E+04/ QSAT2B.969
DATA (ES(IES),IES=1331,1425) / QSAT2B.970
*0.864179E+04,0.868705E+04,0.873237E+04,0.877800E+04,0.882374E+04, QSAT2B.971
*0.886979E+04,0.891603E+04,0.896237E+04,0.900904E+04,0.905579E+04, QSAT2B.972
*0.910288E+04,0.915018E+04,0.919758E+04,0.924529E+04,0.929310E+04, QSAT2B.973
*0.934122E+04,0.938959E+04,0.943804E+04,0.948687E+04,0.953575E+04, QSAT2B.974
*0.958494E+04,0.963442E+04,0.968395E+04,0.973384E+04,0.978383E+04, QSAT2B.975
*0.983412E+04,0.988468E+04,0.993534E+04,0.998630E+04,0.100374E+05, QSAT2B.976
*0.100888E+05,0.101406E+05,0.101923E+05,0.102444E+05,0.102966E+05, QSAT2B.977
*0.103492E+05,0.104020E+05,0.104550E+05,0.105082E+05,0.105616E+05, QSAT2B.978
*0.106153E+05,0.106693E+05,0.107234E+05,0.107779E+05,0.108325E+05, QSAT2B.979
*0.108874E+05,0.109425E+05,0.109978E+05,0.110535E+05,0.111092E+05, QSAT2B.980
*0.111653E+05,0.112217E+05,0.112782E+05,0.113350E+05,0.113920E+05, QSAT2B.981
*0.114493E+05,0.115070E+05,0.115646E+05,0.116228E+05,0.116809E+05, QSAT2B.982
*0.117396E+05,0.117984E+05,0.118574E+05,0.119167E+05,0.119762E+05, QSAT2B.983
*0.120360E+05,0.120962E+05,0.121564E+05,0.122170E+05,0.122778E+05, QSAT2B.984
*0.123389E+05,0.124004E+05,0.124619E+05,0.125238E+05,0.125859E+05, QSAT2B.985
*0.126484E+05,0.127111E+05,0.127739E+05,0.128372E+05,0.129006E+05, QSAT2B.986
*0.129644E+05,0.130285E+05,0.130927E+05,0.131573E+05,0.132220E+05, QSAT2B.987
*0.132872E+05,0.133526E+05,0.134182E+05,0.134842E+05,0.135503E+05, QSAT2B.988
*0.136168E+05,0.136836E+05,0.137505E+05,0.138180E+05,0.138854E+05/ QSAT2B.989
DATA (ES(IES),IES=1426,1520) / QSAT2B.990
*0.139534E+05,0.140216E+05,0.140900E+05,0.141588E+05,0.142277E+05, QSAT2B.991
*0.142971E+05,0.143668E+05,0.144366E+05,0.145069E+05,0.145773E+05, QSAT2B.992
*0.146481E+05,0.147192E+05,0.147905E+05,0.148622E+05,0.149341E+05, QSAT2B.993
*0.150064E+05,0.150790E+05,0.151517E+05,0.152250E+05,0.152983E+05, QSAT2B.994
*0.153721E+05,0.154462E+05,0.155205E+05,0.155952E+05,0.156701E+05, QSAT2B.995
*0.157454E+05,0.158211E+05,0.158969E+05,0.159732E+05,0.160496E+05, QSAT2B.996
*0.161265E+05,0.162037E+05,0.162811E+05,0.163589E+05,0.164369E+05, QSAT2B.997
*0.165154E+05,0.165942E+05,0.166732E+05,0.167526E+05,0.168322E+05, QSAT2B.998
*0.169123E+05,0.169927E+05,0.170733E+05,0.171543E+05,0.172356E+05, QSAT2B.999
*0.173173E+05,0.173993E+05,0.174815E+05,0.175643E+05,0.176471E+05, QSAT2B.1000
*0.177305E+05,0.178143E+05,0.178981E+05,0.179826E+05,0.180671E+05, QSAT2B.1001
*0.181522E+05,0.182377E+05,0.183232E+05,0.184093E+05,0.184955E+05, QSAT2B.1002
*0.185823E+05,0.186695E+05,0.187568E+05,0.188447E+05,0.189326E+05, QSAT2B.1003
*0.190212E+05,0.191101E+05,0.191991E+05,0.192887E+05,0.193785E+05, QSAT2B.1004
*0.194688E+05,0.195595E+05,0.196503E+05,0.197417E+05,0.198332E+05, QSAT2B.1005
*0.199253E+05,0.200178E+05,0.201105E+05,0.202036E+05,0.202971E+05, QSAT2B.1006
*0.203910E+05,0.204853E+05,0.205798E+05,0.206749E+05,0.207701E+05, QSAT2B.1007
*0.208659E+05,0.209621E+05,0.210584E+05,0.211554E+05,0.212524E+05, QSAT2B.1008
*0.213501E+05,0.214482E+05,0.215465E+05,0.216452E+05,0.217442E+05/ QSAT2B.1009
DATA (ES(IES),IES=1521,1552) / QSAT2B.1010
*0.218439E+05,0.219439E+05,0.220440E+05,0.221449E+05,0.222457E+05, QSAT2B.1011
*0.223473E+05,0.224494E+05,0.225514E+05,0.226542E+05,0.227571E+05, QSAT2B.1012
*0.228606E+05,0.229646E+05,0.230687E+05,0.231734E+05,0.232783E+05, QSAT2B.1013
*0.233839E+05,0.234898E+05,0.235960E+05,0.237027E+05,0.238097E+05, QSAT2B.1014
*0.239173E+05,0.240254E+05,0.241335E+05,0.242424E+05,0.243514E+05, QSAT2B.1015
*0.244611E+05,0.245712E+05,0.246814E+05,0.247923E+05,0.249034E+05, QSAT2B.1016
*0.250152E+05,0.250152E+05/ QSAT2B.1017
C QSAT2B.1018
DO 10 I=1,NPNTS QSAT2B.1019
C COMPUTE THE FACTOR THAT CONVERTS FROM SAT VAPOUR PRESSURE QSAT2B.1020
C IN A PURE WATER SYSTEM TO SAT VAPOUR PRESSURE IN AIR, FSUBW. QSAT2B.1021
C THIS FORMULA IS TAKEN FROM EQUATION A4.7 OF ADRIAN GILL'S QSAT2B.1022
C BOOK: ATMOSPHERE-OCEAN DYNAMICS. NOTE THAT HIS FORMULA QSAT2B.1023
C WORKS IN TERMS OF PRESSURE IN MB AND TEMPERATURE IN CELSIUS, QSAT2B.1024
C SO CONVERSION OF UNITS LEADS TO THE SLIGHTLY DIFFERENT QSAT2B.1025
C EQUATION USED HERE. QSAT2B.1026
C QSAT2B.1027
FSUBW = 1.0 +1.0E-8*P(I)*( 4.5 + QSAT2B.1028
+ 6.0E-4*( T(I) - ZERODEGC )*( T(I) - ZERODEGC ) ) QSAT2B.1029
C QSAT2B.1030
C USE THE LOOKUP TABLE TO FIND SATURATED VAAPOUR PRESSURE, QSAT2B.1031
C AND STORE IT IN QS. QSAT2B.1032
C QSAT2B.1033
C QSAT2B.1034
TT=MAX(T_LOW,T(I)) QSAT2B.1035
TT=MIN(T_HIGH,TT) QSAT2B.1036
QSAT2B.1037
ATABLE = (TT - T_LOW + DELTA_T) / DELTA_T QSAT2B.1038
ITABLE = ATABLE QSAT2B.1039
ATABLE = ATABLE - ITABLE QSAT2B.1040
C QSAT2B.1041
QS(I) = (1.0 - ATABLE)*ES(ITABLE) QSAT2B.1042
+ + ATABLE*ES(ITABLE+1) QSAT2B.1043
C QSAT2B.1044
C QSAT2B.1045
C MULTIPLY BY FSUBW TO CONVERT TO SATURATED VAPOUR PRESSURE QSAT2B.1046
C IN AIR (EQUATION A4.6 OF ADRIAN GILL'S BOOK). QSAT2B.1047
C QSAT2B.1048
QS(I) = QS(I)*FSUBW QSAT2B.1049
C QSAT2B.1050
C NOW FORM THE ACCURATE EXPRESSION FOR QS, WHICH IS A QSAT2B.1051
C REARRANGED VERSION OF EQUATION A4.3 OF GILL'S BOOK. QSAT2B.1052
C QSAT2B.1053
C NOTE THAT AT VERY LOW PRESSURES WE APPLY A FIX, TO QSAT2B.1054
C PREVENT A SINGULARITY. QSAT2B.1055
C QSAT2B.1056
QS(I) = ( EPSILON*QS(I) ) / QSAT2B.1057
+ ( MAX(P(I),QS(I)) - ONE_MINUS_EPSILON*QS(I) ) QSAT2B.1058
C QSAT2B.1059
10 CONTINUE QSAT2B.1060
C QSAT2B.1061
RETURN QSAT2B.1062
END QSAT2B.1063
*ENDIF QSAT2B.1064