*IF DEF,RECON QSATPF1A.2
C *****************************COPYRIGHT****************************** QSATPF1A.3
C (c) CROWN COPYRIGHT 1996, METEOROLOGICAL OFFICE, All Rights Reserved. QSATPF1A.4
C QSATPF1A.5
C Use, duplication or disclosure of this code is subject to the QSATPF1A.6
C restrictions as set forth in the contract. QSATPF1A.7
C QSATPF1A.8
C Meteorological Office QSATPF1A.9
C London Road QSATPF1A.10
C BRACKNELL QSATPF1A.11
C Berkshire UK QSATPF1A.12
C RG12 2SZ QSATPF1A.13
C QSATPF1A.14
C If no contract has been raised with this copy of the code, the use, QSATPF1A.15
C duplication or disclosure of it is strictly prohibited. Permission QSATPF1A.16
C to do so must first be obtained in writing from the Head of Numerical QSATPF1A.17
C Modelling at the above address. QSATPF1A.18
C ******************************COPYRIGHT****************************** QSATPF1A.19
!+ QSATPF1A.20
! QSATPF1A.21
! Subroutine Interface: QSATPF1A.22
SUBROUTINE QSAT_PF (qs,dlnesBydT,T,P,npnts) 1QSATPF1A.23
QSATPF1A.24
IMPLICIT NONE QSATPF1A.25
QSATPF1A.26
! QSATPF1A.27
! Description: RETURNS A SATURATION MIXING RATIO AND RATE OF CHANGE QSATPF1A.28
! OF NATURAL LOG OF SATURATED VAPOUR PRESSURE WRT TEMP QSATPF1A.29
! GIVEN A TEMPERATURE AND PRESSURE USING SATURATION QSATPF1A.30
! VAPOUR PRESSURES CALCULATED USING THE GOFF-GRATCH QSATPF1A.31
! FORMULAE, ADOPTED BY THE WMO AS TAKEN FROM QSATPF1A.32
! LANDOLT-BORNSTEIN, 1987 NUMERICAL DATA QSATPF1A.33
! AND FUNCTIONAL RELATIONSHIPS IN SCIENCE AND TECHNOLOGY. QSATPF1A.34
! GROUP V/VOL 4B METEOROLOGY. QSATPF1A.35
! QSATPF1A.36
! VALUES IN THE LOOKUP TABLE ARE OVER WATER ABOVE QSATPF1A.37
! 0 DEG C AND OVER ICE BELOW THIS TEMPERATURE QSATPF1A.38
! QSATPF1A.39
! QSATPF1A.40
! CODE BASED ON UM ROUTINE QSAT QSATPF1A.41
! QSATPF1A.42
! Current Code Owner: I Edmond QSATPF1A.43
! QSATPF1A.44
! History: QSATPF1A.45
! Version Date Comment QSATPF1A.46
! ------- ---- ------- QSATPF1A.47
! 4.1 15/6/96 Original code. Ian Edmond QSATPF1A.48
! vn4.4 9/4/97 Declaration change to allow code to compile UIE2F404.22
! with the Nag f90 compiler. Ian Edmond UIE2F404.23
! 4.5 29/07/98 Optimisation changes for T3E Rewrote **KAPPA UDG5F405.468
! calculations to reduce number of "**"'s and UDG5F405.469
! replaced "**"'s with vector function powr_v UDG5F405.470
! Author D.M. Goddard UDG5F405.471
! QSATPF1A.49
! Code Description: QSATPF1A.50
! Language: FORTRAN 77 + common extensions. QSATPF1A.51
! This code is written to UMDP3 v6 programming standards. QSATPF1A.52
! QSATPF1A.53
! System component covered: <appropriate code> QSATPF1A.54
! System Task: <appropriate code> QSATPF1A.55
! QSATPF1A.56
! Declarations: QSATPF1A.57
! These are of the form:- QSATPF1A.58
! INTEGER ExampleVariable !Description of variable QSATPF1A.59
! 1.0 Subroutine arguments QSATPF1A.60
! 1.1 Array arguments with intent(in): QSATPF1A.61
INTEGER UIE2F404.25
& npnts ! VECTOR LENGTH UIE2F404.26
REAL QSATPF1A.62
& T(npnts) ! IN TEMPERATURE (K) QSATPF1A.63
&,P(npnts) ! IN PRESSURE (PA) QSATPF1A.64
QSATPF1A.65
! 1.2 Array arguments with intent(out) QSATPF1A.66
REAL QSATPF1A.67
& qs(npnts) ! SATURATION MIXING RATIO AT TEMPERATURE QSATPF1A.68
! T AND PRESSURE P (KG/KG) QSATPF1A.69
&,dlnesBydT(npnts) ! RATE OF CHANGE OF NATURAL LOGARITHM OF QSATPF1A.70
! SATURATION VAPOUR PRESSURE WRT TEMPERATURE QSATPF1A.71
QSATPF1A.72
! 2.0 Local parameters: QSATPF1A.73
QSATPF1A.74
REAL epsilon ! EPSILON IS RATIO OF MOLECULAR QSATPF1A.75
PARAMETER(epsilon=0.62198)! WEIGHTS OF WATER AND DRY AIR QSATPF1A.76
QSATPF1A.77
REAL T_low ! LOWEST TEMPERATURE FOR WHICH LOOK-UP QSATPF1A.78
PARAMETER(T_low=183.15) !TABLE OF SATURATION WATER VAPOUR QSATPF1A.79
! PRESSURE IS VALID (K) QSATPF1A.80
QSATPF1A.81
REAL T_high ! HIGHEST TEMPERATURE FOR WHICH LOOK-UP QSATPF1A.82
PARAMETER(T_high = 338.15)! TABLE OF SATURATION WATER VAPOUR QSATPF1A.83
! PRESSURES IS VALID (K) QSATPF1A.84
QSATPF1A.85
REAL delta_T ! TEMPERATURE INCREMENT OF THE LOOK-UP QSATPF1A.86
PARAMETER(delta_T = 0.1) ! TABLE OF SATURATION VAPOUR PRESSURES QSATPF1A.87
! VECTOR LENGTH QSATPF1A.88
QSATPF1A.89
INTEGER n ! SIZE OF LOOK-UP TABLE OF SATURATION QSATPF1A.90
! WATER VAPOUR PRESSURES QSATPF1A.91
PARAMETER(n=((T_high - T_LOW + (delta_T*0.5))/delta_T) + 1.0) QSATPF1A.92
QSATPF1A.93
! 3.0 Local scalars: QSATPF1A.94
QSATPF1A.95
INTEGER QSATPF1A.96
& i ! LOOP COUNTER UIE2F404.24
&,ies ! LOOP COUNTER FOR DATA STATEMENT LOOK-UP TABLE QSATPF1A.99
QSATPF1A.100
INTEGER QSATPF1A.101
& itable ! WORK VARIABLES QSATPF1A.102
&,jtable ! WORK VARIABLES QSATPF1A.103
QSATPF1A.104
REAL QSATPF1A.105
& atable ! WORK VARIABLES QSATPF1A.106
&,tt UDG5F405.472
QSATPF1A.107
! 4.0 Local dynamic arrays: QSATPF1A.108
QSATPF1A.109
REAL QSATPF1A.110
& t1(2,npnts), UDG5F405.473
& es(0:n+1) ! TABLE OF SATURATION WATER VAPOUR UDG5F405.474
! PRESSURE (PA) - SET BY DATA STATEMENT QSATPF1A.112
! CALCULATED FROM THE GOFF-GRATCH FORMULAE QSATPF1A.113
! AS TAKEN FROM LANDOLT-BORNSTEIN, 1987 QSATPF1A.114
! NUMERICAL DATA AND FUNCTIONAL RELATIONSHIPS QSATPF1A.115
! IN SCIENCE AND TECHNOLOGY. GROUP V/ VOL 4B QSATPF1A.116
! METEOROLOGY. PHYSICAL AND CHEMICAL QSATPF1A.117
! PROPERTIES OF AIR, P35 QSATPF1A.118
QSATPF1A.119
!---------------------------------------------------------------------- QSATPF1A.120
! SATURATION WATER VAPOUR PRESSURE QSATPF1A.121
! QSATPF1A.122
! ABOVE 0 DEG C VALUES ARE OVER WATER QSATPF1A.123
! QSATPF1A.124
! BELOW 0 DEC C VALUES ARE OVER ICE QSATPF1A.125
!---------------------------------------------------------------------- QSATPF1A.126
QSATPF1A.127
DATA (ES(IES),IES= 0, 95) / 0.966483E-02, UDG5F405.475
*0.966483E-02,0.984279E-02,0.100240E-01,0.102082E-01,0.103957E-01, QSATPF1A.129
*0.105865E-01,0.107803E-01,0.109777E-01,0.111784E-01,0.113825E-01, QSATPF1A.130
*0.115902E-01,0.118016E-01,0.120164E-01,0.122348E-01,0.124572E-01, QSATPF1A.131
*0.126831E-01,0.129132E-01,0.131470E-01,0.133846E-01,0.136264E-01, QSATPF1A.132
*0.138724E-01,0.141225E-01,0.143771E-01,0.146356E-01,0.148985E-01, QSATPF1A.133
*0.151661E-01,0.154379E-01,0.157145E-01,0.159958E-01,0.162817E-01, QSATPF1A.134
*0.165725E-01,0.168680E-01,0.171684E-01,0.174742E-01,0.177847E-01, QSATPF1A.135
*0.181008E-01,0.184216E-01,0.187481E-01,0.190801E-01,0.194175E-01, QSATPF1A.136
*0.197608E-01,0.201094E-01,0.204637E-01,0.208242E-01,0.211906E-01, QSATPF1A.137
*0.215631E-01,0.219416E-01,0.223263E-01,0.227172E-01,0.231146E-01, QSATPF1A.138
*0.235188E-01,0.239296E-01,0.243465E-01,0.247708E-01,0.252019E-01, QSATPF1A.139
*0.256405E-01,0.260857E-01,0.265385E-01,0.269979E-01,0.274656E-01, QSATPF1A.140
*0.279405E-01,0.284232E-01,0.289142E-01,0.294124E-01,0.299192E-01, QSATPF1A.141
*0.304341E-01,0.309571E-01,0.314886E-01,0.320285E-01,0.325769E-01, QSATPF1A.142
*0.331348E-01,0.337014E-01,0.342771E-01,0.348618E-01,0.354557E-01, QSATPF1A.143
*0.360598E-01,0.366727E-01,0.372958E-01,0.379289E-01,0.385717E-01, QSATPF1A.144
*0.392248E-01,0.398889E-01,0.405633E-01,0.412474E-01,0.419430E-01, QSATPF1A.145
*0.426505E-01,0.433678E-01,0.440974E-01,0.448374E-01,0.455896E-01, QSATPF1A.146
*0.463545E-01,0.471303E-01,0.479191E-01,0.487190E-01,0.495322E-01/ QSATPF1A.147
DATA (es(ies),ies= 96,190) / QSATPF1A.148
*0.503591E-01,0.511977E-01,0.520490E-01,0.529145E-01,0.537931E-01, QSATPF1A.149
*0.546854E-01,0.555924E-01,0.565119E-01,0.574467E-01,0.583959E-01, QSATPF1A.150
*0.593592E-01,0.603387E-01,0.613316E-01,0.623409E-01,0.633655E-01, QSATPF1A.151
*0.644053E-01,0.654624E-01,0.665358E-01,0.676233E-01,0.687302E-01, QSATPF1A.152
*0.698524E-01,0.709929E-01,0.721490E-01,0.733238E-01,0.745180E-01, QSATPF1A.153
*0.757281E-01,0.769578E-01,0.782061E-01,0.794728E-01,0.807583E-01, QSATPF1A.154
*0.820647E-01,0.833905E-01,0.847358E-01,0.861028E-01,0.874882E-01, QSATPF1A.155
*0.888957E-01,0.903243E-01,0.917736E-01,0.932464E-01,0.947407E-01, QSATPF1A.156
*0.962571E-01,0.977955E-01,0.993584E-01,0.100942E+00,0.102551E+00, QSATPF1A.157
*0.104186E+00,0.105842E+00,0.107524E+00,0.109231E+00,0.110963E+00, QSATPF1A.158
*0.112722E+00,0.114506E+00,0.116317E+00,0.118153E+00,0.120019E+00, QSATPF1A.159
*0.121911E+00,0.123831E+00,0.125778E+00,0.127755E+00,0.129761E+00, QSATPF1A.160
*0.131796E+00,0.133863E+00,0.135956E+00,0.138082E+00,0.140241E+00, QSATPF1A.161
*0.142428E+00,0.144649E+00,0.146902E+00,0.149190E+00,0.151506E+00, QSATPF1A.162
*0.153859E+00,0.156245E+00,0.158669E+00,0.161126E+00,0.163618E+00, QSATPF1A.163
*0.166145E+00,0.168711E+00,0.171313E+00,0.173951E+00,0.176626E+00, QSATPF1A.164
*0.179342E+00,0.182096E+00,0.184893E+00,0.187724E+00,0.190600E+00, QSATPF1A.165
*0.193518E+00,0.196473E+00,0.199474E+00,0.202516E+00,0.205604E+00, QSATPF1A.166
*0.208730E+00,0.211905E+00,0.215127E+00,0.218389E+00,0.221701E+00/ QSATPF1A.167
DATA (es(ies),ies=191,285) / QSATPF1A.168
*0.225063E+00,0.228466E+00,0.231920E+00,0.235421E+00,0.238976E+00, QSATPF1A.169
*0.242580E+00,0.246232E+00,0.249933E+00,0.253691E+00,0.257499E+00, QSATPF1A.170
*0.261359E+00,0.265278E+00,0.269249E+00,0.273274E+00,0.277358E+00, QSATPF1A.171
*0.281498E+00,0.285694E+00,0.289952E+00,0.294268E+00,0.298641E+00, QSATPF1A.172
*0.303078E+00,0.307577E+00,0.312135E+00,0.316753E+00,0.321440E+00, QSATPF1A.173
*0.326196E+00,0.331009E+00,0.335893E+00,0.340842E+00,0.345863E+00, QSATPF1A.174
*0.350951E+00,0.356106E+00,0.361337E+00,0.366636E+00,0.372006E+00, QSATPF1A.175
*0.377447E+00,0.382966E+00,0.388567E+00,0.394233E+00,0.399981E+00, QSATPF1A.176
*0.405806E+00,0.411714E+00,0.417699E+00,0.423772E+00,0.429914E+00, QSATPF1A.177
*0.436145E+00,0.442468E+00,0.448862E+00,0.455359E+00,0.461930E+00, QSATPF1A.178
*0.468596E+00,0.475348E+00,0.482186E+00,0.489124E+00,0.496160E+00, QSATPF1A.179
*0.503278E+00,0.510497E+00,0.517808E+00,0.525224E+00,0.532737E+00, QSATPF1A.180
*0.540355E+00,0.548059E+00,0.555886E+00,0.563797E+00,0.571825E+00, QSATPF1A.181
*0.579952E+00,0.588198E+00,0.596545E+00,0.605000E+00,0.613572E+00, QSATPF1A.182
*0.622255E+00,0.631059E+00,0.639962E+00,0.649003E+00,0.658144E+00, QSATPF1A.183
*0.667414E+00,0.676815E+00,0.686317E+00,0.695956E+00,0.705728E+00, QSATPF1A.184
*0.715622E+00,0.725641E+00,0.735799E+00,0.746082E+00,0.756495E+00, QSATPF1A.185
*0.767052E+00,0.777741E+00,0.788576E+00,0.799549E+00,0.810656E+00, QSATPF1A.186
*0.821914E+00,0.833314E+00,0.844854E+00,0.856555E+00,0.868415E+00/ QSATPF1A.187
DATA (es(ies),ies=286,380) / QSATPF1A.188
*0.880404E+00,0.892575E+00,0.904877E+00,0.917350E+00,0.929974E+00, QSATPF1A.189
*0.942771E+00,0.955724E+00,0.968837E+00,0.982127E+00,0.995600E+00, QSATPF1A.190
*0.100921E+01,0.102304E+01,0.103700E+01,0.105116E+01,0.106549E+01, QSATPF1A.191
*0.108002E+01,0.109471E+01,0.110962E+01,0.112469E+01,0.113995E+01, QSATPF1A.192
*0.115542E+01,0.117107E+01,0.118693E+01,0.120298E+01,0.121923E+01, QSATPF1A.193
*0.123569E+01,0.125234E+01,0.126923E+01,0.128631E+01,0.130362E+01, QSATPF1A.194
*0.132114E+01,0.133887E+01,0.135683E+01,0.137500E+01,0.139342E+01, QSATPF1A.195
*0.141205E+01,0.143091E+01,0.145000E+01,0.146933E+01,0.148892E+01, QSATPF1A.196
*0.150874E+01,0.152881E+01,0.154912E+01,0.156970E+01,0.159049E+01, QSATPF1A.197
*0.161159E+01,0.163293E+01,0.165452E+01,0.167640E+01,0.169852E+01, QSATPF1A.198
*0.172091E+01,0.174359E+01,0.176653E+01,0.178977E+01,0.181332E+01, QSATPF1A.199
*0.183709E+01,0.186119E+01,0.188559E+01,0.191028E+01,0.193524E+01, QSATPF1A.200
*0.196054E+01,0.198616E+01,0.201208E+01,0.203829E+01,0.206485E+01, QSATPF1A.201
*0.209170E+01,0.211885E+01,0.214637E+01,0.217424E+01,0.220242E+01, QSATPF1A.202
*0.223092E+01,0.225979E+01,0.228899E+01,0.231855E+01,0.234845E+01, QSATPF1A.203
*0.237874E+01,0.240937E+01,0.244040E+01,0.247176E+01,0.250349E+01, QSATPF1A.204
*0.253560E+01,0.256814E+01,0.260099E+01,0.263431E+01,0.266800E+01, QSATPF1A.205
*0.270207E+01,0.273656E+01,0.277145E+01,0.280671E+01,0.284248E+01, QSATPF1A.206
*0.287859E+01,0.291516E+01,0.295219E+01,0.298962E+01,0.302746E+01/ QSATPF1A.207
DATA (es(ies),ies=381,475) / QSATPF1A.208
*0.306579E+01,0.310454E+01,0.314377E+01,0.318351E+01,0.322360E+01, QSATPF1A.209
*0.326427E+01,0.330538E+01,0.334694E+01,0.338894E+01,0.343155E+01, QSATPF1A.210
*0.347456E+01,0.351809E+01,0.356216E+01,0.360673E+01,0.365184E+01, QSATPF1A.211
*0.369744E+01,0.374352E+01,0.379018E+01,0.383743E+01,0.388518E+01, QSATPF1A.212
*0.393344E+01,0.398230E+01,0.403177E+01,0.408175E+01,0.413229E+01, QSATPF1A.213
*0.418343E+01,0.423514E+01,0.428746E+01,0.434034E+01,0.439389E+01, QSATPF1A.214
*0.444808E+01,0.450276E+01,0.455820E+01,0.461423E+01,0.467084E+01, QSATPF1A.215
*0.472816E+01,0.478607E+01,0.484468E+01,0.490393E+01,0.496389E+01, QSATPF1A.216
*0.502446E+01,0.508580E+01,0.514776E+01,0.521047E+01,0.527385E+01, QSATPF1A.217
*0.533798E+01,0.540279E+01,0.546838E+01,0.553466E+01,0.560173E+01, QSATPF1A.218
*0.566949E+01,0.573807E+01,0.580750E+01,0.587749E+01,0.594846E+01, QSATPF1A.219
*0.602017E+01,0.609260E+01,0.616591E+01,0.623995E+01,0.631490E+01, QSATPF1A.220
*0.639061E+01,0.646723E+01,0.654477E+01,0.662293E+01,0.670220E+01, QSATPF1A.221
*0.678227E+01,0.686313E+01,0.694495E+01,0.702777E+01,0.711142E+01, QSATPF1A.222
*0.719592E+01,0.728140E+01,0.736790E+01,0.745527E+01,0.754352E+01, QSATPF1A.223
*0.763298E+01,0.772316E+01,0.781442E+01,0.790676E+01,0.800001E+01, QSATPF1A.224
*0.809435E+01,0.818967E+01,0.828606E+01,0.838343E+01,0.848194E+01, QSATPF1A.225
*0.858144E+01,0.868207E+01,0.878392E+01,0.888673E+01,0.899060E+01, QSATPF1A.226
*0.909567E+01,0.920172E+01,0.930909E+01,0.941765E+01,0.952730E+01/ QSATPF1A.227
DATA (es(ies),ies=476,570) / QSATPF1A.228
*0.963821E+01,0.975022E+01,0.986352E+01,0.997793E+01,0.100937E+02, QSATPF1A.229
*0.102105E+02,0.103287E+02,0.104481E+02,0.105688E+02,0.106909E+02, QSATPF1A.230
*0.108143E+02,0.109387E+02,0.110647E+02,0.111921E+02,0.113207E+02, QSATPF1A.231
*0.114508E+02,0.115821E+02,0.117149E+02,0.118490E+02,0.119847E+02, QSATPF1A.232
*0.121216E+02,0.122601E+02,0.124002E+02,0.125416E+02,0.126846E+02, QSATPF1A.233
*0.128290E+02,0.129747E+02,0.131224E+02,0.132712E+02,0.134220E+02, QSATPF1A.234
*0.135742E+02,0.137278E+02,0.138831E+02,0.140403E+02,0.141989E+02, QSATPF1A.235
*0.143589E+02,0.145211E+02,0.146845E+02,0.148501E+02,0.150172E+02, QSATPF1A.236
*0.151858E+02,0.153564E+02,0.155288E+02,0.157029E+02,0.158786E+02, QSATPF1A.237
*0.160562E+02,0.162358E+02,0.164174E+02,0.166004E+02,0.167858E+02, QSATPF1A.238
*0.169728E+02,0.171620E+02,0.173528E+02,0.175455E+02,0.177406E+02, QSATPF1A.239
*0.179372E+02,0.181363E+02,0.183372E+02,0.185400E+02,0.187453E+02, QSATPF1A.240
*0.189523E+02,0.191613E+02,0.193728E+02,0.195866E+02,0.198024E+02, QSATPF1A.241
*0.200200E+02,0.202401E+02,0.204626E+02,0.206871E+02,0.209140E+02, QSATPF1A.242
*0.211430E+02,0.213744E+02,0.216085E+02,0.218446E+02,0.220828E+02, QSATPF1A.243
*0.223241E+02,0.225671E+02,0.228132E+02,0.230615E+02,0.233120E+02, QSATPF1A.244
*0.235651E+02,0.238211E+02,0.240794E+02,0.243404E+02,0.246042E+02, QSATPF1A.245
*0.248704E+02,0.251390E+02,0.254109E+02,0.256847E+02,0.259620E+02, QSATPF1A.246
*0.262418E+02,0.265240E+02,0.268092E+02,0.270975E+02,0.273883E+02/ QSATPF1A.247
DATA (es(ies),ies=571,665) / QSATPF1A.248
*0.276822E+02,0.279792E+02,0.282789E+02,0.285812E+02,0.288867E+02, QSATPF1A.249
*0.291954E+02,0.295075E+02,0.298222E+02,0.301398E+02,0.304606E+02, QSATPF1A.250
*0.307848E+02,0.311119E+02,0.314424E+02,0.317763E+02,0.321133E+02, QSATPF1A.251
*0.324536E+02,0.327971E+02,0.331440E+02,0.334940E+02,0.338475E+02, QSATPF1A.252
*0.342050E+02,0.345654E+02,0.349295E+02,0.352975E+02,0.356687E+02, QSATPF1A.253
*0.360430E+02,0.364221E+02,0.368042E+02,0.371896E+02,0.375790E+02, QSATPF1A.254
*0.379725E+02,0.383692E+02,0.387702E+02,0.391744E+02,0.395839E+02, QSATPF1A.255
*0.399958E+02,0.404118E+02,0.408325E+02,0.412574E+02,0.416858E+02, QSATPF1A.256
*0.421188E+02,0.425551E+02,0.429962E+02,0.434407E+02,0.438910E+02, QSATPF1A.257
*0.443439E+02,0.448024E+02,0.452648E+02,0.457308E+02,0.462018E+02, QSATPF1A.258
*0.466775E+02,0.471582E+02,0.476428E+02,0.481313E+02,0.486249E+02, QSATPF1A.259
*0.491235E+02,0.496272E+02,0.501349E+02,0.506479E+02,0.511652E+02, QSATPF1A.260
*0.516876E+02,0.522142E+02,0.527474E+02,0.532836E+02,0.538266E+02, QSATPF1A.261
*0.543737E+02,0.549254E+02,0.554839E+02,0.560456E+02,0.566142E+02, QSATPF1A.262
*0.571872E+02,0.577662E+02,0.583498E+02,0.589392E+02,0.595347E+02, QSATPF1A.263
*0.601346E+02,0.607410E+02,0.613519E+02,0.619689E+02,0.625922E+02, QSATPF1A.264
*0.632204E+02,0.638550E+02,0.644959E+02,0.651418E+02,0.657942E+02, QSATPF1A.265
*0.664516E+02,0.671158E+02,0.677864E+02,0.684624E+02,0.691451E+02, QSATPF1A.266
*0.698345E+02,0.705293E+02,0.712312E+02,0.719398E+02,0.726542E+02/ QSATPF1A.267
DATA (es(ies),ies=666,760) / QSATPF1A.268
*0.733754E+02,0.741022E+02,0.748363E+02,0.755777E+02,0.763247E+02, QSATPF1A.269
*0.770791E+02,0.778394E+02,0.786088E+02,0.793824E+02,0.801653E+02, QSATPF1A.270
*0.809542E+02,0.817509E+02,0.825536E+02,0.833643E+02,0.841828E+02, QSATPF1A.271
*0.850076E+02,0.858405E+02,0.866797E+02,0.875289E+02,0.883827E+02, QSATPF1A.272
*0.892467E+02,0.901172E+02,0.909962E+02,0.918818E+02,0.927760E+02, QSATPF1A.273
*0.936790E+02,0.945887E+02,0.955071E+02,0.964346E+02,0.973689E+02, QSATPF1A.274
*0.983123E+02,0.992648E+02,0.100224E+03,0.101193E+03,0.102169E+03, QSATPF1A.275
*0.103155E+03,0.104150E+03,0.105152E+03,0.106164E+03,0.107186E+03, QSATPF1A.276
*0.108217E+03,0.109256E+03,0.110303E+03,0.111362E+03,0.112429E+03, QSATPF1A.277
*0.113503E+03,0.114588E+03,0.115684E+03,0.116789E+03,0.117903E+03, QSATPF1A.278
*0.119028E+03,0.120160E+03,0.121306E+03,0.122460E+03,0.123623E+03, QSATPF1A.279
*0.124796E+03,0.125981E+03,0.127174E+03,0.128381E+03,0.129594E+03, QSATPF1A.280
*0.130822E+03,0.132058E+03,0.133306E+03,0.134563E+03,0.135828E+03, QSATPF1A.281
*0.137109E+03,0.138402E+03,0.139700E+03,0.141017E+03,0.142338E+03, QSATPF1A.282
*0.143676E+03,0.145025E+03,0.146382E+03,0.147753E+03,0.149133E+03, QSATPF1A.283
*0.150529E+03,0.151935E+03,0.153351E+03,0.154783E+03,0.156222E+03, QSATPF1A.284
*0.157678E+03,0.159148E+03,0.160624E+03,0.162117E+03,0.163621E+03, QSATPF1A.285
*0.165142E+03,0.166674E+03,0.168212E+03,0.169772E+03,0.171340E+03, QSATPF1A.286
*0.172921E+03,0.174522E+03,0.176129E+03,0.177755E+03,0.179388E+03/ QSATPF1A.287
DATA (es(ies),ies=761,855) / QSATPF1A.288
*0.181040E+03,0.182707E+03,0.184382E+03,0.186076E+03,0.187782E+03, QSATPF1A.289
*0.189503E+03,0.191240E+03,0.192989E+03,0.194758E+03,0.196535E+03, QSATPF1A.290
*0.198332E+03,0.200141E+03,0.201963E+03,0.203805E+03,0.205656E+03, QSATPF1A.291
*0.207532E+03,0.209416E+03,0.211317E+03,0.213236E+03,0.215167E+03, QSATPF1A.292
*0.217121E+03,0.219087E+03,0.221067E+03,0.223064E+03,0.225080E+03, QSATPF1A.293
*0.227113E+03,0.229160E+03,0.231221E+03,0.233305E+03,0.235403E+03, QSATPF1A.294
*0.237520E+03,0.239655E+03,0.241805E+03,0.243979E+03,0.246163E+03, QSATPF1A.295
*0.248365E+03,0.250593E+03,0.252830E+03,0.255093E+03,0.257364E+03, QSATPF1A.296
*0.259667E+03,0.261979E+03,0.264312E+03,0.266666E+03,0.269034E+03, QSATPF1A.297
*0.271430E+03,0.273841E+03,0.276268E+03,0.278722E+03,0.281185E+03, QSATPF1A.298
*0.283677E+03,0.286190E+03,0.288714E+03,0.291266E+03,0.293834E+03, QSATPF1A.299
*0.296431E+03,0.299045E+03,0.301676E+03,0.304329E+03,0.307006E+03, QSATPF1A.300
*0.309706E+03,0.312423E+03,0.315165E+03,0.317930E+03,0.320705E+03, QSATPF1A.301
*0.323519E+03,0.326350E+03,0.329199E+03,0.332073E+03,0.334973E+03, QSATPF1A.302
*0.337897E+03,0.340839E+03,0.343800E+03,0.346794E+03,0.349806E+03, QSATPF1A.303
*0.352845E+03,0.355918E+03,0.358994E+03,0.362112E+03,0.365242E+03, QSATPF1A.304
*0.368407E+03,0.371599E+03,0.374802E+03,0.378042E+03,0.381293E+03, QSATPF1A.305
*0.384588E+03,0.387904E+03,0.391239E+03,0.394604E+03,0.397988E+03, QSATPF1A.306
*0.401411E+03,0.404862E+03,0.408326E+03,0.411829E+03,0.415352E+03/ QSATPF1A.307
DATA (es(ies),ies=856,950) / QSATPF1A.308
*0.418906E+03,0.422490E+03,0.426095E+03,0.429740E+03,0.433398E+03, QSATPF1A.309
*0.437097E+03,0.440827E+03,0.444570E+03,0.448354E+03,0.452160E+03, QSATPF1A.310
*0.455999E+03,0.459870E+03,0.463765E+03,0.467702E+03,0.471652E+03, QSATPF1A.311
*0.475646E+03,0.479674E+03,0.483715E+03,0.487811E+03,0.491911E+03, QSATPF1A.312
*0.496065E+03,0.500244E+03,0.504448E+03,0.508698E+03,0.512961E+03, QSATPF1A.313
*0.517282E+03,0.521617E+03,0.525989E+03,0.530397E+03,0.534831E+03, QSATPF1A.314
*0.539313E+03,0.543821E+03,0.548355E+03,0.552938E+03,0.557549E+03, QSATPF1A.315
*0.562197E+03,0.566884E+03,0.571598E+03,0.576351E+03,0.581131E+03, QSATPF1A.316
*0.585963E+03,0.590835E+03,0.595722E+03,0.600663E+03,0.605631E+03, QSATPF1A.317
*0.610641E+03,0.615151E+03,0.619625E+03,0.624140E+03,0.628671E+03, QSATPF1A.318
*0.633243E+03,0.637845E+03,0.642465E+03,0.647126E+03,0.651806E+03, QSATPF1A.319
*0.656527E+03,0.661279E+03,0.666049E+03,0.670861E+03,0.675692E+03, QSATPF1A.320
*0.680566E+03,0.685471E+03,0.690396E+03,0.695363E+03,0.700350E+03, QSATPF1A.321
*0.705381E+03,0.710444E+03,0.715527E+03,0.720654E+03,0.725801E+03, QSATPF1A.322
*0.730994E+03,0.736219E+03,0.741465E+03,0.746756E+03,0.752068E+03, QSATPF1A.323
*0.757426E+03,0.762819E+03,0.768231E+03,0.773692E+03,0.779172E+03, QSATPF1A.324
*0.784701E+03,0.790265E+03,0.795849E+03,0.801483E+03,0.807137E+03, QSATPF1A.325
*0.812842E+03,0.818582E+03,0.824343E+03,0.830153E+03,0.835987E+03, QSATPF1A.326
*0.841871E+03,0.847791E+03,0.853733E+03,0.859727E+03,0.865743E+03/ QSATPF1A.327
DATA (es(ies),ies=951,1045) / QSATPF1A.328
*0.871812E+03,0.877918E+03,0.884046E+03,0.890228E+03,0.896433E+03, QSATPF1A.329
*0.902690E+03,0.908987E+03,0.915307E+03,0.921681E+03,0.928078E+03, QSATPF1A.330
*0.934531E+03,0.941023E+03,0.947539E+03,0.954112E+03,0.960708E+03, QSATPF1A.331
*0.967361E+03,0.974053E+03,0.980771E+03,0.987545E+03,0.994345E+03, QSATPF1A.332
*0.100120E+04,0.100810E+04,0.101502E+04,0.102201E+04,0.102902E+04, QSATPF1A.333
*0.103608E+04,0.104320E+04,0.105033E+04,0.105753E+04,0.106475E+04, QSATPF1A.334
*0.107204E+04,0.107936E+04,0.108672E+04,0.109414E+04,0.110158E+04, QSATPF1A.335
*0.110908E+04,0.111663E+04,0.112421E+04,0.113185E+04,0.113952E+04, QSATPF1A.336
*0.114725E+04,0.115503E+04,0.116284E+04,0.117071E+04,0.117861E+04, QSATPF1A.337
*0.118658E+04,0.119459E+04,0.120264E+04,0.121074E+04,0.121888E+04, QSATPF1A.338
*0.122709E+04,0.123534E+04,0.124362E+04,0.125198E+04,0.126036E+04, QSATPF1A.339
*0.126881E+04,0.127731E+04,0.128584E+04,0.129444E+04,0.130307E+04, QSATPF1A.340
*0.131177E+04,0.132053E+04,0.132931E+04,0.133817E+04,0.134705E+04, QSATPF1A.341
*0.135602E+04,0.136503E+04,0.137407E+04,0.138319E+04,0.139234E+04, QSATPF1A.342
*0.140156E+04,0.141084E+04,0.142015E+04,0.142954E+04,0.143896E+04, QSATPF1A.343
*0.144845E+04,0.145800E+04,0.146759E+04,0.147725E+04,0.148694E+04, QSATPF1A.344
*0.149672E+04,0.150655E+04,0.151641E+04,0.152635E+04,0.153633E+04, QSATPF1A.345
*0.154639E+04,0.155650E+04,0.156665E+04,0.157688E+04,0.158715E+04, QSATPF1A.346
*0.159750E+04,0.160791E+04,0.161836E+04,0.162888E+04,0.163945E+04/ QSATPF1A.347
DATA (es(ies),ies=1046,1140) / QSATPF1A.348
*0.165010E+04,0.166081E+04,0.167155E+04,0.168238E+04,0.169325E+04, QSATPF1A.349
*0.170420E+04,0.171522E+04,0.172627E+04,0.173741E+04,0.174859E+04, QSATPF1A.350
*0.175986E+04,0.177119E+04,0.178256E+04,0.179402E+04,0.180552E+04, QSATPF1A.351
*0.181711E+04,0.182877E+04,0.184046E+04,0.185224E+04,0.186407E+04, QSATPF1A.352
*0.187599E+04,0.188797E+04,0.190000E+04,0.191212E+04,0.192428E+04, QSATPF1A.353
*0.193653E+04,0.194886E+04,0.196122E+04,0.197368E+04,0.198618E+04, QSATPF1A.354
*0.199878E+04,0.201145E+04,0.202416E+04,0.203698E+04,0.204983E+04, QSATPF1A.355
*0.206278E+04,0.207580E+04,0.208887E+04,0.210204E+04,0.211525E+04, QSATPF1A.356
*0.212856E+04,0.214195E+04,0.215538E+04,0.216892E+04,0.218249E+04, QSATPF1A.357
*0.219618E+04,0.220994E+04,0.222375E+04,0.223766E+04,0.225161E+04, QSATPF1A.358
*0.226567E+04,0.227981E+04,0.229399E+04,0.230829E+04,0.232263E+04, QSATPF1A.359
*0.233708E+04,0.235161E+04,0.236618E+04,0.238087E+04,0.239560E+04, QSATPF1A.360
*0.241044E+04,0.242538E+04,0.244035E+04,0.245544E+04,0.247057E+04, QSATPF1A.361
*0.248583E+04,0.250116E+04,0.251654E+04,0.253204E+04,0.254759E+04, QSATPF1A.362
*0.256325E+04,0.257901E+04,0.259480E+04,0.261073E+04,0.262670E+04, QSATPF1A.363
*0.264279E+04,0.265896E+04,0.267519E+04,0.269154E+04,0.270794E+04, QSATPF1A.364
*0.272447E+04,0.274108E+04,0.275774E+04,0.277453E+04,0.279137E+04, QSATPF1A.365
*0.280834E+04,0.282540E+04,0.284251E+04,0.285975E+04,0.287704E+04, QSATPF1A.366
*0.289446E+04,0.291198E+04,0.292954E+04,0.294725E+04,0.296499E+04/ QSATPF1A.367
DATA (es(ies),ies=1141,1235) / QSATPF1A.368
*0.298288E+04,0.300087E+04,0.301890E+04,0.303707E+04,0.305529E+04, QSATPF1A.369
*0.307365E+04,0.309211E+04,0.311062E+04,0.312927E+04,0.314798E+04, QSATPF1A.370
*0.316682E+04,0.318577E+04,0.320477E+04,0.322391E+04,0.324310E+04, QSATPF1A.371
*0.326245E+04,0.328189E+04,0.330138E+04,0.332103E+04,0.334073E+04, QSATPF1A.372
*0.336058E+04,0.338053E+04,0.340054E+04,0.342069E+04,0.344090E+04, QSATPF1A.373
*0.346127E+04,0.348174E+04,0.350227E+04,0.352295E+04,0.354369E+04, QSATPF1A.374
*0.356458E+04,0.358559E+04,0.360664E+04,0.362787E+04,0.364914E+04, QSATPF1A.375
*0.367058E+04,0.369212E+04,0.371373E+04,0.373548E+04,0.375731E+04, QSATPF1A.376
*0.377929E+04,0.380139E+04,0.382355E+04,0.384588E+04,0.386826E+04, QSATPF1A.377
*0.389081E+04,0.391348E+04,0.393620E+04,0.395910E+04,0.398205E+04, QSATPF1A.378
*0.400518E+04,0.402843E+04,0.405173E+04,0.407520E+04,0.409875E+04, QSATPF1A.379
*0.412246E+04,0.414630E+04,0.417019E+04,0.419427E+04,0.421840E+04, QSATPF1A.380
*0.424272E+04,0.426715E+04,0.429165E+04,0.431634E+04,0.434108E+04, QSATPF1A.381
*0.436602E+04,0.439107E+04,0.441618E+04,0.444149E+04,0.446685E+04, QSATPF1A.382
*0.449241E+04,0.451810E+04,0.454385E+04,0.456977E+04,0.459578E+04, QSATPF1A.383
*0.462197E+04,0.464830E+04,0.467468E+04,0.470127E+04,0.472792E+04, QSATPF1A.384
*0.475477E+04,0.478175E+04,0.480880E+04,0.483605E+04,0.486336E+04, QSATPF1A.385
*0.489087E+04,0.491853E+04,0.494623E+04,0.497415E+04,0.500215E+04, QSATPF1A.386
*0.503034E+04,0.505867E+04,0.508707E+04,0.511568E+04,0.514436E+04/ QSATPF1A.387
DATA (es(ies),ies=1236,1330) / QSATPF1A.388
*0.517325E+04,0.520227E+04,0.523137E+04,0.526068E+04,0.529005E+04, QSATPF1A.389
*0.531965E+04,0.534939E+04,0.537921E+04,0.540923E+04,0.543932E+04, QSATPF1A.390
*0.546965E+04,0.550011E+04,0.553064E+04,0.556139E+04,0.559223E+04, QSATPF1A.391
*0.562329E+04,0.565449E+04,0.568577E+04,0.571727E+04,0.574884E+04, QSATPF1A.392
*0.578064E+04,0.581261E+04,0.584464E+04,0.587692E+04,0.590924E+04, QSATPF1A.393
*0.594182E+04,0.597455E+04,0.600736E+04,0.604039E+04,0.607350E+04, QSATPF1A.394
*0.610685E+04,0.614036E+04,0.617394E+04,0.620777E+04,0.624169E+04, QSATPF1A.395
*0.627584E+04,0.631014E+04,0.634454E+04,0.637918E+04,0.641390E+04, QSATPF1A.396
*0.644887E+04,0.648400E+04,0.651919E+04,0.655467E+04,0.659021E+04, QSATPF1A.397
*0.662599E+04,0.666197E+04,0.669800E+04,0.673429E+04,0.677069E+04, QSATPF1A.398
*0.680735E+04,0.684415E+04,0.688104E+04,0.691819E+04,0.695543E+04, QSATPF1A.399
*0.699292E+04,0.703061E+04,0.706837E+04,0.710639E+04,0.714451E+04, QSATPF1A.400
*0.718289E+04,0.722143E+04,0.726009E+04,0.729903E+04,0.733802E+04, QSATPF1A.401
*0.737729E+04,0.741676E+04,0.745631E+04,0.749612E+04,0.753602E+04, QSATPF1A.402
*0.757622E+04,0.761659E+04,0.765705E+04,0.769780E+04,0.773863E+04, QSATPF1A.403
*0.777975E+04,0.782106E+04,0.786246E+04,0.790412E+04,0.794593E+04, QSATPF1A.404
*0.798802E+04,0.803028E+04,0.807259E+04,0.811525E+04,0.815798E+04, QSATPF1A.405
*0.820102E+04,0.824427E+04,0.828757E+04,0.833120E+04,0.837493E+04, QSATPF1A.406
*0.841895E+04,0.846313E+04,0.850744E+04,0.855208E+04,0.859678E+04/ QSATPF1A.407
DATA (es(ies),ies=1331,1425) / QSATPF1A.408
*0.864179E+04,0.868705E+04,0.873237E+04,0.877800E+04,0.882374E+04, QSATPF1A.409
*0.886979E+04,0.891603E+04,0.896237E+04,0.900904E+04,0.905579E+04, QSATPF1A.410
*0.910288E+04,0.915018E+04,0.919758E+04,0.924529E+04,0.929310E+04, QSATPF1A.411
*0.934122E+04,0.938959E+04,0.943804E+04,0.948687E+04,0.953575E+04, QSATPF1A.412
*0.958494E+04,0.963442E+04,0.968395E+04,0.973384E+04,0.978383E+04, QSATPF1A.413
*0.983412E+04,0.988468E+04,0.993534E+04,0.998630E+04,0.100374E+05, QSATPF1A.414
*0.100888E+05,0.101406E+05,0.101923E+05,0.102444E+05,0.102966E+05, QSATPF1A.415
*0.103492E+05,0.104020E+05,0.104550E+05,0.105082E+05,0.105616E+05, QSATPF1A.416
*0.106153E+05,0.106693E+05,0.107234E+05,0.107779E+05,0.108325E+05, QSATPF1A.417
*0.108874E+05,0.109425E+05,0.109978E+05,0.110535E+05,0.111092E+05, QSATPF1A.418
*0.111653E+05,0.112217E+05,0.112782E+05,0.113350E+05,0.113920E+05, QSATPF1A.419
*0.114493E+05,0.115070E+05,0.115646E+05,0.116228E+05,0.116809E+05, QSATPF1A.420
*0.117396E+05,0.117984E+05,0.118574E+05,0.119167E+05,0.119762E+05, QSATPF1A.421
*0.120360E+05,0.120962E+05,0.121564E+05,0.122170E+05,0.122778E+05, QSATPF1A.422
*0.123389E+05,0.124004E+05,0.124619E+05,0.125238E+05,0.125859E+05, QSATPF1A.423
*0.126484E+05,0.127111E+05,0.127739E+05,0.128372E+05,0.129006E+05, QSATPF1A.424
*0.129644E+05,0.130285E+05,0.130927E+05,0.131573E+05,0.132220E+05, QSATPF1A.425
*0.132872E+05,0.133526E+05,0.134182E+05,0.134842E+05,0.135503E+05, QSATPF1A.426
*0.136168E+05,0.136836E+05,0.137505E+05,0.138180E+05,0.138854E+05/ QSATPF1A.427
DATA (es(ies),ies=1426,1520) / QSATPF1A.428
*0.139534E+05,0.140216E+05,0.140900E+05,0.141588E+05,0.142277E+05, QSATPF1A.429
*0.142971E+05,0.143668E+05,0.144366E+05,0.145069E+05,0.145773E+05, QSATPF1A.430
*0.146481E+05,0.147192E+05,0.147905E+05,0.148622E+05,0.149341E+05, QSATPF1A.431
*0.150064E+05,0.150790E+05,0.151517E+05,0.152250E+05,0.152983E+05, QSATPF1A.432
*0.153721E+05,0.154462E+05,0.155205E+05,0.155952E+05,0.156701E+05, QSATPF1A.433
*0.157454E+05,0.158211E+05,0.158969E+05,0.159732E+05,0.160496E+05, QSATPF1A.434
*0.161265E+05,0.162037E+05,0.162811E+05,0.163589E+05,0.164369E+05, QSATPF1A.435
*0.165154E+05,0.165942E+05,0.166732E+05,0.167526E+05,0.168322E+05, QSATPF1A.436
*0.169123E+05,0.169927E+05,0.170733E+05,0.171543E+05,0.172356E+05, QSATPF1A.437
*0.173173E+05,0.173993E+05,0.174815E+05,0.175643E+05,0.176471E+05, QSATPF1A.438
*0.177305E+05,0.178143E+05,0.178981E+05,0.179826E+05,0.180671E+05, QSATPF1A.439
*0.181522E+05,0.182377E+05,0.183232E+05,0.184093E+05,0.184955E+05, QSATPF1A.440
*0.185823E+05,0.186695E+05,0.187568E+05,0.188447E+05,0.189326E+05, QSATPF1A.441
*0.190212E+05,0.191101E+05,0.191991E+05,0.192887E+05,0.193785E+05, QSATPF1A.442
*0.194688E+05,0.195595E+05,0.196503E+05,0.197417E+05,0.198332E+05, QSATPF1A.443
*0.199253E+05,0.200178E+05,0.201105E+05,0.202036E+05,0.202971E+05, QSATPF1A.444
*0.203910E+05,0.204853E+05,0.205798E+05,0.206749E+05,0.207701E+05, QSATPF1A.445
*0.208659E+05,0.209621E+05,0.210584E+05,0.211554E+05,0.212524E+05, QSATPF1A.446
*0.213501E+05,0.214482E+05,0.215465E+05,0.216452E+05,0.217442E+05/ QSATPF1A.447
DATA (es(ies),ies=1521,1552) / UDG5F405.476
*0.218439E+05,0.219439E+05,0.220440E+05,0.221449E+05,0.222457E+05, QSATPF1A.449
*0.223473E+05,0.224494E+05,0.225514E+05,0.226542E+05,0.227571E+05, QSATPF1A.450
*0.228606E+05,0.229646E+05,0.230687E+05,0.231734E+05,0.232783E+05, QSATPF1A.451
*0.233839E+05,0.234898E+05,0.235960E+05,0.237027E+05,0.238097E+05, QSATPF1A.452
*0.239173E+05,0.240254E+05,0.241335E+05,0.242424E+05,0.243514E+05, QSATPF1A.453
*0.244611E+05,0.245712E+05,0.246814E+05,0.247923E+05,0.249034E+05, QSATPF1A.454
*0.250152E+05,0.250152E+05/ UDG5F405.477
QSATPF1A.456
!- End of header QSATPF1A.457
QSATPF1A.458
Do i=1,npnts QSATPF1A.459
QSATPF1A.460
TT=MAX(T_LOW,T(I)) UDG5F405.478
TT=MIN(T_HIGH,TT) UDG5F405.479
UDG5F405.480
atable = (TT - T_low + delta_T) / delta_T UDG5F405.481
itable = atable QSATPF1A.470
atable = atable - itable QSATPF1A.471
QSATPF1A.472
qs(I) = epsilon * ((1.-atable) * es(itable) QSATPF1A.473
* + atable * es(itable+1)) / P(I) QSATPF1A.474
QSATPF1A.475
End do UDG5F405.482
UDG5F405.483
Do i=1,npnts UDG5F405.484
QSATPF1A.477
TT=MAX(T_LOW+delta_T,T(I)) UDG5F405.485
TT=MIN(T_HIGH-delta_T,TT) UDG5F405.486
UDG5F405.487
atable = (TT - T_low ) / delta_T UDG5F405.488
itable = atable QSATPF1A.488
atable = (TT - T_low + 2*delta_T) / delta_T UDG5F405.489
jtable = atable QSATPF1A.490
atable = atable - jtable QSATPF1A.491
*IF DEF,VECTLIB PXVECTLB.124
t1(1,i)=((1.-atable) * es(jtable) + UDG5F405.491
& atable * es(jtable+1)) UDG5F405.492
t1(2,i)=((1.-atable) * es(itable) + UDG5F405.493
& atable * es(itable+1)) UDG5F405.494
*ELSE UDG5F405.495
t1(1,i)=LOG((1.-atable) * es(jtable) + UDG5F405.496
& atable * es(jtable+1)) UDG5F405.497
t1(2,i)=LOG((1.-atable) * es(itable) + UDG5F405.498
& atable * es(itable+1)) UDG5F405.499
*ENDIF UDG5F405.500
UDG5F405.501
End do UDG5F405.502
UDG5F405.503
*IF DEF,VECTLIB PXVECTLB.125
call alog_v(
2*npnts,t1,t1) UDG5F405.505
UDG5F405.506
*ENDIF UDG5F405.507
Do i=1,npnts UDG5F405.508
dlnesBydT(I) = (t1(1,i)-t1(2,i))/(2*delta_T) UDG5F405.509
End do UDG5F405.510
UDG5F405.511
UDG5F405.512
RETURN QSATPF1A.500
END QSATPF1A.501
*ENDIF QSATPF1A.502