      function x_rm_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.2959515E-01/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.31312305E-01,-0.14796564E-01, 0.59280229E+00, 0.49214821E-01,
     + -0.46108749E-01, 0.14570918E-01,-0.22722213E-01,-0.49254913E-02,
     + -0.92012724E-02, 0.51151239E-02,-0.21827884E-02,-0.13570552E-02,
     +  0.54984605E-02,-0.87119229E-02, 0.12744780E-01, 0.14404791E-02,
     + -0.99429546E-03, 0.10645495E-01,-0.23127585E-02, 0.48288560E-03,
     + -0.57611644E-03, 0.10934970E-02,-0.15130508E-02,-0.64102118E-02,
     + -0.11280115E-02, 0.19355290E-03,-0.24248223E-03,-0.39885826E-02,
     + -0.30448174E-02,-0.79223473E-03, 0.10026763E-04,-0.26283744E-02,
     + -0.16864746E-03,-0.37479401E-02, 0.36058445E-02, 0.51963019E-04,
     + -0.25135151E-03,-0.20509037E-03,-0.34382399E-02,-0.79507596E-03,
     + -0.15037507E-02,-0.43960316E-02, 0.39490382E-02, 0.17908786E-03,
     +  0.40420677E-03, 0.18966637E-03, 0.54642282E-04,-0.12367414E-02,
     + -0.43591217E-03, 0.23623266E-03, 0.11139095E-03,-0.87888085E-03,
     + -0.25211426E-02,-0.23621719E-02,-0.40090407E-03, 0.27069639E-03,
     + -0.89822296E-03,-0.12260140E-02,-0.98276534E-04, 0.89903752E-03,
     + -0.16274583E-02,-0.49447739E-02, 0.48669247E-03,-0.45318613E-02,
     + -0.80988638E-03,-0.10011581E-02, 0.69831776E-05, 0.14393464E-03,
     +  0.73078752E-03,-0.17296587E-03,-0.10199330E-03, 0.14223545E-03,
     +  0.76541946E-04,-0.32006649E-03,-0.69520646E-03,-0.23075807E-03,
     + -0.19138094E-03, 0.26871811E-03,-0.34230825E-03, 0.84449811E-03,
     +  0.25456012E-03, 0.44468430E-04, 0.25204054E-03,-0.87252981E-03,
     +  0.42386522E-03,-0.10856831E-02,-0.12587398E-03, 0.89869282E-04,
     +  0.60616661E-04,-0.22865795E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_rm_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)                *x52
     6  +coeff(  6)    *x22            
     7  +coeff(  7)    *x21    *x42    
     8  +coeff(  8)        *x31*x41    
      x_rm_fp     =x_rm_fp     
     9  +coeff(  9)            *x42    
     1  +coeff( 10)                *x53
     2  +coeff( 11)*x11                
     3  +coeff( 12)    *x21*x31*x41    
     4  +coeff( 13)    *x22        *x51
     5  +coeff( 14)    *x24            
     6  +coeff( 15)    *x22    *x42    
     7  +coeff( 16)    *x23            
     8  +coeff( 17)    *x21        *x52
      x_rm_fp     =x_rm_fp     
     9  +coeff( 18)    *x22    *x42*x51
     1  +coeff( 19)    *x23    *x42*x51
     2  +coeff( 20)    *x21    *x41    
     3  +coeff( 21)    *x21*x32        
     4  +coeff( 22)        *x31*x41*x51
     5  +coeff( 23)    *x21*x31*x41*x51
     6  +coeff( 24)    *x24        *x51
     7  +coeff( 25)    *x24        *x52
     8  +coeff( 26)            *x41    
      x_rm_fp     =x_rm_fp     
     9  +coeff( 27)        *x32        
     1  +coeff( 28)    *x23        *x51
     2  +coeff( 29)    *x21    *x42*x51
     3  +coeff( 30)            *x42*x52
     4  +coeff( 31)*x11*x21            
     5  +coeff( 32)    *x23*x31*x41    
     6  +coeff( 33)*x11            *x51
     7  +coeff( 34)    *x23        *x52
     8  +coeff( 35)    *x21    *x42*x52
      x_rm_fp     =x_rm_fp     
     9  +coeff( 36)        *x31        
     1  +coeff( 37)    *x22    *x41    
     2  +coeff( 38)        *x32    *x51
     3  +coeff( 39)    *x22        *x52
     4  +coeff( 40)                *x54
     5  +coeff( 41)    *x22        *x53
     6  +coeff( 42)    *x24    *x42    
     7  +coeff( 43)    *x22    *x42*x52
     8  +coeff( 44)    *x23        *x53
      x_rm_fp     =x_rm_fp     
     9  +coeff( 45)*x11*x23            
     1  +coeff( 46)*x11        *x44    
     2  +coeff( 47)    *x21*x31        
     3  +coeff( 48)    *x22*x31*x41    
     4  +coeff( 49)            *x44    
     5  +coeff( 50)    *x21*x32    *x51
     6  +coeff( 51)        *x32    *x52
     7  +coeff( 52)    *x21        *x53
     8  +coeff( 53)    *x21*x31*x43    
      x_rm_fp     =x_rm_fp     
     9  +coeff( 54)    *x21    *x44    
     1  +coeff( 55)    *x22*x31*x41*x51
     2  +coeff( 56)        *x31*x43*x51
     3  +coeff( 57)            *x44*x51
     4  +coeff( 58)            *x42*x53
     5  +coeff( 59)*x11*x22            
     6  +coeff( 60)    *x21    *x42*x53
     7  +coeff( 61)    *x23*x32*x42    
     8  +coeff( 62)    *x24    *x42*x51
      x_rm_fp     =x_rm_fp     
     9  +coeff( 63)    *x22    *x44*x51
     1  +coeff( 64)    *x23    *x42*x52
     2  +coeff( 65)    *x22*x33*x43    
     3  +coeff( 66)    *x22    *x44*x52
     4  +coeff( 67)            *x41*x51
     5  +coeff( 68)    *x21    *x41*x51
     6  +coeff( 69)            *x42*x51
     7  +coeff( 70)    *x22    *x41*x51
     8  +coeff( 71)    *x21    *x41*x52
      x_rm_fp     =x_rm_fp     
     9  +coeff( 72)        *x31*x41*x52
     1  +coeff( 73)    *x23*x32        
     2  +coeff( 74)    *x21*x33*x41    
     3  +coeff( 75)    *x21*x32*x42    
     4  +coeff( 76)    *x21*x32    *x52
     5  +coeff( 77)    *x21    *x41*x53
     6  +coeff( 78)        *x31*x41*x53
     7  +coeff( 79)    *x21        *x54
     8  +coeff( 80)    *x24*x31*x41    
      x_rm_fp     =x_rm_fp     
     9  +coeff( 81)    *x22*x33*x41    
     1  +coeff( 82)*x11            *x52
     2  +coeff( 83)    *x22*x32    *x52
     3  +coeff( 84)            *x44*x52
     4  +coeff( 85)    *x22        *x54
     5  +coeff( 86)            *x42*x54
     6  +coeff( 87)*x11*x22    *x41    
     7  +coeff( 88)*x11*x21*x31*x41    
     8  +coeff( 89)*x11    *x32*x41    
      x_rm_fp     =x_rm_fp     
     9  +coeff( 90)    *x24*x31*x42    
c
      return
      end
      function t_rm_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 14)
      data ncoeff/ 13/
      data avdat/ -0.6260162E-02/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.61264923E-02,-0.29467598E-01, 0.97875349E-01, 0.41515026E-02,
     + -0.92660673E-02,-0.30143671E-02,-0.84616500E-03,-0.12821936E-02,
     +  0.93618815E-03, 0.17289305E-03,-0.61316736E-03,-0.11568249E-02,
     +  0.19218810E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
c
c                  function
c
      t_rm_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)                *x52
     6  +coeff(  6)    *x21    *x42    
     7  +coeff(  7)        *x31*x41    
     8  +coeff(  8)            *x42    
      t_rm_fp     =t_rm_fp     
     9  +coeff(  9)                *x53
     1  +coeff( 10)    *x22            
     2  +coeff( 11)    *x21        *x52
     3  +coeff( 12)    *x24            
     4  +coeff( 13)    *x22    *x42    
c
      return
      end
      function y_rm_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 42)
      data ncoeff/ 41/
      data avdat/ -0.3494358E-03/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.57060906E-03,-0.12975497E-01,-0.22609044E-01, 0.29703972E-06,
     + -0.34047859E-02,-0.12480460E-01,-0.13689178E-02, 0.26399836E-01,
     + -0.49356645E-03, 0.28158800E-04,-0.47017331E-03,-0.15985806E-01,
     + -0.21583430E-03,-0.11071013E-01, 0.20102144E-02, 0.37116366E-02,
     +  0.32082906E-02,-0.65779299E-02,-0.28430850E-02,-0.21790290E-03,
     + -0.19601511E-02,-0.19839173E-02,-0.13503424E-02, 0.10175610E-03,
     +  0.29074709E-03,-0.92324388E-03, 0.31127317E-02, 0.19189958E-02,
     + -0.21744177E-02,-0.31717853E-04, 0.12544284E-03, 0.10639494E-03,
     + -0.23551812E-03, 0.53717021E-03, 0.16194966E-02, 0.99165447E-03,
     +  0.13900723E-02,-0.23219418E-02, 0.77330571E-03, 0.84419025E-03,
     +  0.11442991E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      y_rm_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)        *x32        
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      y_rm_fp     =y_rm_fp     
     9  +coeff(  9)        *x31*x42    
     1  +coeff( 10)            *x42*x51
     2  +coeff( 11)    *x22*x31        
     3  +coeff( 12)    *x22    *x41    
     4  +coeff( 13)    *x21*x31    *x51
     5  +coeff( 14)    *x21    *x41*x51
     6  +coeff( 15)        *x31    *x52
     7  +coeff( 16)            *x41*x52
     8  +coeff( 17)    *x23    *x41    
      y_rm_fp     =y_rm_fp     
     9  +coeff( 18)    *x22    *x41*x51
     1  +coeff( 19)            *x41*x53
     2  +coeff( 20)                *x51
     3  +coeff( 21)            *x43    
     4  +coeff( 22)            *x43*x51
     5  +coeff( 23)    *x22*x31    *x51
     6  +coeff( 24)    *x21            
     7  +coeff( 25)*x11*x21    *x41    
     8  +coeff( 26)        *x31    *x53
      y_rm_fp     =y_rm_fp     
     9  +coeff( 27)    *x23    *x41*x51
     1  +coeff( 28)    *x22    *x41*x52
     2  +coeff( 29)            *x41*x54
     3  +coeff( 30)    *x23    *x41*x53
     4  +coeff( 31)    *x22            
     5  +coeff( 32)    *x21        *x51
     6  +coeff( 33)        *x32*x41    
     7  +coeff( 34)    *x21    *x43    
     8  +coeff( 35)    *x21    *x43*x51
      y_rm_fp     =y_rm_fp     
     9  +coeff( 36)        *x31*x42*x52
     1  +coeff( 37)    *x24    *x41    
     2  +coeff( 38)    *x21    *x41*x53
     3  +coeff( 39)        *x31    *x54
     4  +coeff( 40)    *x21    *x43*x52
     5  +coeff( 41)    *x21*x31*x44*x51
c
      return
      end
      function p_rm_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 26)
      data ncoeff/ 25/
      data avdat/  0.1300000E-03/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.24185822E-03,-0.17643868E-03, 0.56287809E-02,-0.40407673E-01,
     +  0.19915290E-02, 0.19825526E-01,-0.38024609E-02,-0.13798978E-01,
     +  0.91076409E-02, 0.17996276E-02, 0.18962107E-02, 0.17234576E-02,
     +  0.14439640E-02, 0.88304380E-03, 0.13848371E-03,-0.17799588E-03,
     +  0.55863068E-03,-0.15386121E-02,-0.12337682E-02, 0.22008267E-03,
     +  0.26394514E-03,-0.10419282E-03,-0.88530622E-03, 0.13750857E-02,
     + -0.88997610E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
c
c                  function
c
      p_rm_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)            *x41    
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      p_rm_fp     =p_rm_fp     
     9  +coeff(  9)    *x22    *x41    
     1  +coeff( 10)    *x22*x31        
     2  +coeff( 11)            *x43    
     3  +coeff( 12)            *x41*x52
     4  +coeff( 13)        *x31*x42    
     5  +coeff( 14)        *x31    *x52
     6  +coeff( 15)                *x51
     7  +coeff( 16)*x11        *x41    
     8  +coeff( 17)    *x21    *x41*x51
      p_rm_fp     =p_rm_fp     
     9  +coeff( 18)    *x23    *x41    
     1  +coeff( 19)            *x43*x53
     2  +coeff( 20)        *x32*x41    
     3  +coeff( 21)    *x21*x31    *x51
     4  +coeff( 22)    *x24            
     5  +coeff( 23)    *x21    *x43    
     6  +coeff( 24)    *x22    *x43    
     7  +coeff( 25)    *x23*x31*x42    
c
      return
      end
      function l_rm_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 10)
      data ncoeff/  9/
      data avdat/ -0.2479135E-01/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.18681193E-01,-0.32364637E+00,-0.30632297E-01,-0.36204338E-01,
     + -0.20303734E-01, 0.31383445E-02,-0.85230926E-02,-0.62010852E-02,
     +  0.13474983E-01,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      l_rm_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)                *x52
     6  +coeff(  6)*x11                
     7  +coeff(  7)            *x42    
     8  +coeff(  8)    *x21        *x51
      l_rm_fp     =l_rm_fp     
     9  +coeff(  9)    *x21    *x42    
c
      return
      end
      function x_rm_dent   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 77)
      data ncoeff/ 76/
      data avdat/ -0.1971463E-02/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.12158405E-02,-0.12687169E+00,-0.77855722E-02, 0.38910429E-02,
     +  0.10818441E-02, 0.75110787E-03, 0.90342871E-03, 0.32449004E-03,
     +  0.49296312E-03, 0.62724540E-03,-0.73074727E-04,-0.18368258E-03,
     +  0.34786157E-04, 0.46320641E-02,-0.38764574E-04, 0.49187503E-04,
     + -0.82101265E-04, 0.19797643E-02, 0.28823747E-02,-0.10506563E-03,
     + -0.15100003E-03,-0.36481106E-04,-0.13183696E-03,-0.24674000E-04,
     +  0.98501472E-03, 0.26302726E-02,-0.27694530E-05,-0.44630884E-04,
     + -0.49712762E-05,-0.57002810E-04, 0.22973648E-03, 0.10960809E-03,
     + -0.14700892E-03,-0.28773281E-03,-0.24444944E-04,-0.41279160E-04,
     + -0.58720976E-04, 0.26910613E-04, 0.29217332E-04,-0.54220367E-04,
     +  0.20748215E-04,-0.69547517E-04, 0.31978349E-03, 0.38880335E-05,
     + -0.14583477E-04, 0.32153759E-04,-0.68234243E-04, 0.10560735E-04,
     +  0.14566386E-04,-0.11627312E-03, 0.65673936E-04,-0.57999758E-04,
     + -0.10112215E-03, 0.47725234E-05,-0.32161275E-04,-0.22397295E-03,
     + -0.68008678E-03,-0.31476517E-03,-0.64382073E-03, 0.30184776E-04,
     + -0.97719667E-05,-0.18235618E-04,-0.33449160E-04, 0.14197029E-03,
     +  0.32643558E-04, 0.17701446E-04,-0.20392497E-04,-0.15862394E-03,
     +  0.68087468E-03,-0.30165378E-03,-0.48138758E-04,-0.47609385E-04,
     +  0.47448504E-03, 0.16158589E-03,-0.23437505E-03,-0.89149387E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x34 = x33*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_rm_dent   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)    *x21        *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)*x11                
     6  +coeff(  6)    *x21    *x42    
     7  +coeff(  7)    *x23            
     8  +coeff(  8)    *x21*x31*x41    
      x_rm_dent   =x_rm_dent   
     9  +coeff(  9)    *x22        *x51
     1  +coeff( 10)    *x21        *x52
     2  +coeff( 11)                *x51
     3  +coeff( 12)            *x42    
     4  +coeff( 13)    *x21*x32        
     5  +coeff( 14)    *x23    *x42    
     6  +coeff( 15)    *x21    *x41    
     7  +coeff( 16)    *x21    *x42*x51
     8  +coeff( 17)*x11*x21            
      x_rm_dent   =x_rm_dent   
     9  +coeff( 18)    *x23*x31*x41    
     1  +coeff( 19)    *x21    *x44    
     2  +coeff( 20)*x11            *x51
     3  +coeff( 21)    *x23*x31*x43    
     4  +coeff( 22)        *x31*x41    
     5  +coeff( 23)    *x23        *x51
     6  +coeff( 24)    *x21*x31*x41*x51
     7  +coeff( 25)    *x21*x32*x42    
     8  +coeff( 26)    *x21*x31*x43    
      x_rm_dent   =x_rm_dent   
     9  +coeff( 27)    *x21*x31        
     1  +coeff( 28)    *x22    *x42    
     2  +coeff( 29)    *x21*x32    *x51
     3  +coeff( 30)    *x21        *x53
     4  +coeff( 31)    *x23*x32        
     5  +coeff( 32)    *x21*x33*x41    
     6  +coeff( 33)    *x24*x31*x41    
     7  +coeff( 34)    *x24    *x42    
     8  +coeff( 35)            *x42*x51
      x_rm_dent   =x_rm_dent   
     9  +coeff( 36)    *x23    *x41    
     1  +coeff( 37)    *x21    *x43    
     2  +coeff( 38)    *x22        *x52
     3  +coeff( 39)*x11*x22            
     4  +coeff( 40)    *x23*x33        
     5  +coeff( 41)*x11        *x42    
     6  +coeff( 42)    *x24        *x52
     7  +coeff( 43)    *x23*x33*x41*x52
     8  +coeff( 44)            *x41    
      x_rm_dent   =x_rm_dent   
     9  +coeff( 45)        *x31*x41*x51
     1  +coeff( 46)    *x24            
     2  +coeff( 47)    *x21*x31*x42    
     3  +coeff( 48)    *x22*x33        
     4  +coeff( 49)    *x22*x32*x41    
     5  +coeff( 50)    *x21*x31*x41*x52
     6  +coeff( 51)    *x21    *x42*x52
     7  +coeff( 52)    *x23*x32*x41    
     8  +coeff( 53)    *x23    *x43    
      x_rm_dent   =x_rm_dent   
     9  +coeff( 54)*x11    *x31    *x51
     1  +coeff( 55)    *x21*x34    *x51
     2  +coeff( 56)    *x23*x31*x41*x51
     3  +coeff( 57)    *x23    *x42*x51
     4  +coeff( 58)    *x21*x31*x43*x51
     5  +coeff( 59)    *x21    *x44*x51
     6  +coeff( 60)*x11*x22    *x41    
     7  +coeff( 61)*x11    *x32*x41    
     8  +coeff( 62)*x11*x22        *x51
      x_rm_dent   =x_rm_dent   
     9  +coeff( 63)        *x32*x44*x51
     1  +coeff( 64)    *x21*x31*x43*x52
     2  +coeff( 65)    *x24        *x53
     3  +coeff( 66)*x11*x23*x31        
     4  +coeff( 67)*x11*x22*x31    *x51
     5  +coeff( 68)    *x23*x32*x42*x51
     6  +coeff( 69)    *x23    *x44*x51
     7  +coeff( 70)    *x23*x32*x44    
     8  +coeff( 71)*x11*x24*x31*x42    
      x_rm_dent   =x_rm_dent   
     9  +coeff( 72)*x11*x24    *x41*x52
     1  +coeff( 73)*x11*x24*x31*x42*x51
     2  +coeff( 74)*x11*x24    *x43*x51
     3  +coeff( 75)    *x22*x33*x44*x54
     4  +coeff( 76)*x11*x24*x33*x44*x53
c
      return
      end
      function t_rm_dent   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 17)
      data ncoeff/ 16/
      data avdat/  0.3064638E-02/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.27066949E-02, 0.58385540E-01, 0.47084666E-02, 0.35225335E-02,
     +  0.60269312E-03,-0.65487711E-03, 0.25465244E-02,-0.29436031E-02,
     +  0.44147638E-03,-0.89815434E-03,-0.11182391E-02, 0.20586849E-03,
     +  0.38328912E-03,-0.35607090E-03,-0.80105057E-03,-0.84323058E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      t_rm_dent   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)    *x22            
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)                *x51
     6  +coeff(  6)*x11                
     7  +coeff(  7)            *x42    
     8  +coeff(  8)    *x21    *x42    
      t_rm_dent   =t_rm_dent   
     9  +coeff(  9)        *x31*x41    
     1  +coeff( 10)    *x23            
     2  +coeff( 11)    *x21*x31*x41    
     3  +coeff( 12)    *x22        *x51
     4  +coeff( 13)            *x42*x51
     5  +coeff( 14)    *x21        *x52
     6  +coeff( 15)    *x24            
     7  +coeff( 16)    *x22    *x42    
c
      return
      end
      function y_rm_dent   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 16)
      data ncoeff/ 15/
      data avdat/  0.5548850E-03/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.19114884E-02, 0.14763227E-01, 0.14345217E+00, 0.42251404E-02,
     +  0.16082872E-01,-0.26748103E-04, 0.15344961E-02, 0.26446288E-02,
     + -0.52206521E-02,-0.24205893E-02,-0.11257696E-02,-0.11236201E-02,
     + -0.15872379E-03,-0.12699275E-02,-0.72668318E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x21 = x2
      x22 = x21*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      y_rm_dent   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)    *x21    *x41    
     5  +coeff(  5)            *x41*x51
     6  +coeff(  6)        *x33    *x51
     7  +coeff(  7)    *x21*x31        
     8  +coeff(  8)        *x31    *x51
      y_rm_dent   =y_rm_dent   
     9  +coeff(  9)    *x22    *x41    
     1  +coeff( 10)            *x43    
     2  +coeff( 11)    *x22*x31        
     3  +coeff( 12)            *x41*x52
     4  +coeff( 13)                *x51
     5  +coeff( 14)        *x31*x42    
     6  +coeff( 15)    *x21    *x41*x51
c
      return
      end
      function p_rm_dent   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 27)
      data ncoeff/ 26/
      data avdat/ -0.6493008E-03/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.12091404E-02,-0.22762461E-01,-0.56177851E-01,-0.29315110E-01,
     +  0.15938440E-01,-0.20578433E-01,-0.16071276E-03,-0.28317407E-03,
     + -0.10228791E-04,-0.64682309E-02, 0.31550433E-02,-0.33254414E-02,
     +  0.27665874E-03,-0.28241018E-02,-0.11157492E-02, 0.16156933E-04,
     + -0.15003741E-03, 0.27615027E-03,-0.15073685E-02,-0.14240200E-02,
     +  0.30296287E-02, 0.16564394E-03,-0.25714090E-03, 0.46311985E-03,
     + -0.24178743E-03, 0.11901254E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      p_rm_dent   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)    *x21    *x41    
     5  +coeff(  5)            *x41*x51
     6  +coeff(  6)    *x22    *x41    
     7  +coeff(  7)    *x23*x31        
     8  +coeff(  8)    *x22*x31    *x51
      p_rm_dent   =p_rm_dent   
     9  +coeff(  9)        *x33    *x51
     1  +coeff( 10)    *x21*x31        
     2  +coeff( 11)        *x31    *x51
     3  +coeff( 12)    *x22*x31        
     4  +coeff( 13)    *x21            
     5  +coeff( 14)            *x43    
     6  +coeff( 15)            *x41*x52
     7  +coeff( 16)        *x31*x44    
     8  +coeff( 17)                *x51
      p_rm_dent   =p_rm_dent   
     9  +coeff( 18)*x11        *x41    
     1  +coeff( 19)        *x31*x42    
     2  +coeff( 20)    *x22    *x41*x51
     3  +coeff( 21)    *x24    *x41    
     4  +coeff( 22)    *x22            
     5  +coeff( 23)        *x32*x41    
     6  +coeff( 24)    *x21    *x41*x51
     7  +coeff( 25)        *x31    *x52
     8  +coeff( 26)    *x23    *x41    
c
      return
      end
      function l_rm_dent   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff(  9)
      data ncoeff/  8/
      data avdat/ -0.1111059E-02/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.23778188E-02, 0.24361008E+00,-0.18663691E-01, 0.15685597E-01,
     + -0.20876296E-02,-0.70757973E-02,-0.21220695E-02,-0.59965081E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x31 = x3
      x41 = x4
      x42 = x41*x4
      x51 = x5
c
c                  function
c
      l_rm_dent   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)    *x22            
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)*x11                
     6  +coeff(  6)            *x42    
     7  +coeff(  7)        *x31*x41    
     8  +coeff(  8)    *x21    *x42    
c
      return
      end
      function x_rm_dext   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 84)
      data ncoeff/ 83/
      data avdat/  0.7234355E-02/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.49796188E-02, 0.35662815E+00, 0.13335755E+00, 0.48381791E-01,
     +  0.13044784E-01,-0.39394516E-02,-0.53608087E-02,-0.34038462E-02,
     + -0.61285030E-02,-0.26163252E-02,-0.90141490E-03,-0.39973724E-03,
     + -0.24161434E-03,-0.13557655E-02,-0.95651287E-03,-0.22837657E-02,
     + -0.15485235E-01, 0.10072922E-03, 0.25593414E-03,-0.64334478E-02,
     + -0.88519724E-02,-0.10010676E-01, 0.32519342E-03, 0.49466282E-04,
     + -0.10866014E-03, 0.61519362E-03,-0.24511379E-04,-0.21084878E-03,
     + -0.35551335E-02, 0.99582430E-04, 0.12490938E-03,-0.35617242E-03,
     + -0.78598497E-03,-0.37739371E-03,-0.35397182E-03, 0.11924547E-02,
     + -0.14290291E-03, 0.50355960E-03, 0.33599898E-03, 0.15481484E-04,
     +  0.13284417E-03,-0.12293318E-03,-0.79433594E-04, 0.10320821E-02,
     +  0.23767064E-03,-0.17211083E-03, 0.22859449E-03, 0.24696515E-03,
     + -0.10177018E-02,-0.14179158E-03,-0.13350659E-02,-0.13686162E-02,
     +  0.90771455E-05,-0.74218449E-04,-0.39391947E-04, 0.12937740E-03,
     +  0.21699254E-04,-0.11726027E-03,-0.74579366E-04, 0.22944018E-03,
     +  0.36689109E-03,-0.39269085E-04, 0.41631731E-04, 0.76620386E-03,
     +  0.98505919E-03, 0.13095065E-03, 0.10900787E-03,-0.11066099E-03,
     +  0.37017158E-04, 0.51889865E-03, 0.46006589E-04,-0.96272130E-03,
     +  0.41526335E-03, 0.21156073E-03,-0.27913938E-03,-0.12528799E-03,
     + -0.70844253E-04, 0.81277103E-04, 0.16173604E-02, 0.20278932E-03,
     + -0.87210891E-03,-0.20231835E-02,-0.10250336E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_rm_dext   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)    *x22            
     6  +coeff(  6)*x11                
     7  +coeff(  7)            *x42    
     8  +coeff(  8)                *x52
      x_rm_dext   =x_rm_dext   
     9  +coeff(  9)    *x21    *x42    
     1  +coeff( 10)    *x21*x31*x41    
     2  +coeff( 11)        *x31*x41    
     3  +coeff( 12)    *x23            
     4  +coeff( 13)    *x21*x32        
     5  +coeff( 14)            *x42*x51
     6  +coeff( 15)    *x21        *x52
     7  +coeff( 16)    *x24            
     8  +coeff( 17)    *x23    *x42    
      x_rm_dext   =x_rm_dext   
     9  +coeff( 18)            *x41    
     1  +coeff( 19)    *x21    *x41    
     2  +coeff( 20)    *x23*x31*x41    
     3  +coeff( 21)    *x21*x31*x43    
     4  +coeff( 22)    *x21    *x44    
     5  +coeff( 23)    *x23*x32*x42    
     6  +coeff( 24)    *x21*x31        
     7  +coeff( 25)        *x32        
     8  +coeff( 26)    *x23        *x51
      x_rm_dext   =x_rm_dext   
     9  +coeff( 27)    *x21    *x42*x51
     1  +coeff( 28)*x11*x21            
     2  +coeff( 29)    *x21*x32*x42    
     3  +coeff( 30)*x11            *x51
     4  +coeff( 31)                *x53
     5  +coeff( 32)    *x22*x31*x41    
     6  +coeff( 33)    *x23*x32        
     7  +coeff( 34)    *x21*x33*x41    
     8  +coeff( 35)    *x24        *x51
      x_rm_dext   =x_rm_dext   
     9  +coeff( 36)    *x22    *x42*x51
     1  +coeff( 37)*x11*x22            
     2  +coeff( 38)    *x23    *x43    
     3  +coeff( 39)    *x23*x31*x41*x51
     4  +coeff( 40)            *x41*x51
     5  +coeff( 41)    *x22        *x51
     6  +coeff( 42)        *x31*x41*x51
     7  +coeff( 43)    *x22*x32        
     8  +coeff( 44)    *x22    *x42    
      x_rm_dext   =x_rm_dext   
     9  +coeff( 45)    *x21*x31*x42    
     1  +coeff( 46)    *x22        *x52
     2  +coeff( 47)    *x23*x33        
     3  +coeff( 48)    *x23*x32*x41    
     4  +coeff( 49)    *x24    *x42    
     5  +coeff( 50)            *x42*x54
     6  +coeff( 51)    *x23*x33*x43*x52
     7  +coeff( 52)*x11*x24*x31*x42*x51
     8  +coeff( 53)        *x31        
      x_rm_dext   =x_rm_dext   
     9  +coeff( 54)    *x21*x33        
     1  +coeff( 55)        *x32*x42    
     2  +coeff( 56)    *x21    *x43    
     3  +coeff( 57)    *x21*x32    *x51
     4  +coeff( 58)    *x21*x31*x41*x51
     5  +coeff( 59)        *x31*x41*x52
     6  +coeff( 60)    *x22*x31*x41*x51
     7  +coeff( 61)    *x21*x31*x41*x52
     8  +coeff( 62)*x11        *x42    
      x_rm_dext   =x_rm_dext   
     9  +coeff( 63)*x11*x21        *x51
     1  +coeff( 64)    *x21*x31*x43*x51
     2  +coeff( 65)    *x21    *x44*x51
     3  +coeff( 66)    *x23        *x53
     4  +coeff( 67)*x11*x23            
     5  +coeff( 68)*x11*x22    *x41    
     6  +coeff( 69)*x11    *x32*x41    
     7  +coeff( 70)    *x21*x32*x44    
     8  +coeff( 71)    *x22*x33*x41*x51
      x_rm_dext   =x_rm_dext   
     9  +coeff( 72)    *x24    *x42*x51
     1  +coeff( 73)    *x22*x32*x42*x51
     2  +coeff( 74)    *x24    *x41*x52
     3  +coeff( 75)    *x22*x32*x41*x52
     4  +coeff( 76)    *x24        *x53
     5  +coeff( 77)*x11*x23*x31        
     6  +coeff( 78)*x11*x22    *x41*x51
     7  +coeff( 79)    *x23*x32*x42*x51
     8  +coeff( 80)    *x22*x32*x41*x53
      x_rm_dext   =x_rm_dext   
     9  +coeff( 81)    *x23*x33*x41*x52
     1  +coeff( 82)    *x23*x32*x44*x51
     2  +coeff( 83)*x11*x24    *x43*x51
c
      return
      end
      function t_rm_dext   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 15)
      data ncoeff/ 14/
      data avdat/ -0.5381034E-02/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.48889928E-02,-0.79465523E-01, 0.22764638E-01,-0.39165663E-02,
     + -0.44434885E-02, 0.68279909E-03, 0.56192471E-03,-0.20161346E-02,
     + -0.12257159E-02, 0.40552607E-02, 0.10902986E-02,-0.14241661E-02,
     +  0.44854204E-03, 0.11052612E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x31 = x3
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      t_rm_dext   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)*x11                
     7  +coeff(  7)        *x31*x41    
     8  +coeff(  8)            *x42    
      t_rm_dext   =t_rm_dext   
     9  +coeff(  9)                *x52
     1  +coeff( 10)    *x21    *x42    
     2  +coeff( 11)    *x21*x31*x41    
     3  +coeff( 12)            *x42*x51
     4  +coeff( 13)    *x21        *x52
     5  +coeff( 14)    *x22    *x42    
c
      return
      end
      function y_rm_dext   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 31)
      data ncoeff/ 30/
      data avdat/ -0.5112524E-03/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.69001969E-03,-0.19158196E-01, 0.12410419E+00, 0.28031823E-05,
     +  0.17512483E-03,-0.79481760E-02,-0.60732678E-01, 0.92906039E-02,
     +  0.52859880E-01,-0.33738330E-01,-0.59362054E-02,-0.84468198E-03,
     + -0.29962439E-04, 0.55420160E-03,-0.50696323E-03,-0.71826747E-04,
     + -0.43305755E-02,-0.59221964E-02,-0.36693225E-02,-0.59890072E-02,
     +  0.52600593E-03,-0.10966142E-02, 0.39071781E-02,-0.37391405E-02,
     +  0.33053465E-03,-0.78330288E-03,-0.77967223E-03, 0.29755482E-02,
     + -0.41175345E-02, 0.27023375E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x45 = x44*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      y_rm_dext   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)        *x32        
     5  +coeff(  5)            *x42    
     6  +coeff(  6)    *x21*x31        
     7  +coeff(  7)    *x21    *x41    
     8  +coeff(  8)        *x31    *x51
      y_rm_dext   =y_rm_dext   
     9  +coeff(  9)            *x41*x51
     1  +coeff( 10)    *x22    *x41    
     2  +coeff( 11)    *x21    *x41*x51
     3  +coeff( 12)            *x45    
     4  +coeff( 13)    *x22*x33        
     5  +coeff( 14)    *x21            
     6  +coeff( 15)                *x51
     7  +coeff( 16)        *x33        
     8  +coeff( 17)        *x31*x42    
      y_rm_dext   =y_rm_dext   
     9  +coeff( 18)    *x22*x31        
     1  +coeff( 19)            *x41*x52
     2  +coeff( 20)            *x43    
     3  +coeff( 21)*x11        *x41    
     4  +coeff( 22)    *x21*x31    *x51
     5  +coeff( 23)    *x23    *x41    
     6  +coeff( 24)    *x22    *x41*x51
     7  +coeff( 25)    *x22            
     8  +coeff( 26)        *x32*x41    
      y_rm_dext   =y_rm_dext   
     9  +coeff( 27)        *x31    *x52
     1  +coeff( 28)    *x21    *x43    
     2  +coeff( 29)    *x22    *x43    
     3  +coeff( 30)    *x23*x31*x42    
c
      return
      end
      function p_rm_dext   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 39)
      data ncoeff/ 38/
      data avdat/ -0.2026863E-03/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.12157404E-03, 0.12844315E-03,-0.74268105E-02, 0.83093941E-02,
     + -0.20211480E-05,-0.17022170E-02,-0.28337327E-05,-0.14357582E-01,
     +  0.17929034E-02, 0.12018704E-01,-0.83294958E-02,-0.11571046E-02,
     + -0.35148123E-02,-0.14070494E-03,-0.11283813E-03,-0.84036647E-03,
     + -0.49530010E-03, 0.21464094E-02,-0.14251973E-02,-0.76550740E-03,
     +  0.10250550E-03, 0.13379084E-02, 0.82220235E-04,-0.16297630E-03,
     + -0.13679752E-03, 0.12499507E-03, 0.52097032E-03,-0.17502163E-03,
     + -0.12548169E-02,-0.10197319E-03, 0.54043479E-03, 0.31057541E-04,
     +  0.31873144E-04, 0.26934850E-03, 0.41968212E-03, 0.22194002E-03,
     + -0.20225342E-03,-0.55421243E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      p_rm_dext   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)            *x41    
     5  +coeff(  5)*x11                
     6  +coeff(  6)    *x21*x31        
     7  +coeff(  7)        *x32        
     8  +coeff(  8)    *x21    *x41    
      p_rm_dext   =p_rm_dext   
     9  +coeff(  9)        *x31    *x51
     1  +coeff( 10)            *x41*x51
     2  +coeff( 11)    *x22    *x41    
     3  +coeff( 12)            *x43    
     4  +coeff( 13)    *x21    *x41*x51
     5  +coeff( 14)    *x24*x31        
     6  +coeff( 15)                *x51
     7  +coeff( 16)        *x31*x42    
     8  +coeff( 17)            *x41*x52
      p_rm_dext   =p_rm_dext   
     9  +coeff( 18)    *x23    *x41    
     1  +coeff( 19)    *x22    *x41*x51
     2  +coeff( 20)    *x22*x31        
     3  +coeff( 21)*x11        *x41    
     4  +coeff( 22)    *x23    *x41*x51
     5  +coeff( 23)    *x22            
     6  +coeff( 24)        *x32*x41    
     7  +coeff( 25)    *x21*x31    *x51
     8  +coeff( 26)*x11*x21    *x41    
      p_rm_dext   =p_rm_dext   
     9  +coeff( 27)    *x21    *x43    
     1  +coeff( 28)    *x22*x31    *x51
     2  +coeff( 29)    *x22    *x43    
     3  +coeff( 30)    *x23*x31*x42    
     4  +coeff( 31)    *x24    *x41*x52
     5  +coeff( 32)            *x42    
     6  +coeff( 33)    *x21        *x51
     7  +coeff( 34)    *x23*x31        
     8  +coeff( 35)    *x21*x31*x42    
      p_rm_dext   =p_rm_dext   
     9  +coeff( 36)            *x43*x51
     1  +coeff( 37)    *x21    *x41*x52
     2  +coeff( 38)    *x22*x31*x42    
c
      return
      end
      function l_rm_dext   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff(  8)
      data ncoeff/  7/
      data avdat/ -0.2204845E-01/
      data xmin/
     1 -0.99954E-03,-0.63963E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.18472731E-01,-0.50298321E+00,-0.98946452E-01,-0.33539932E-01,
     + -0.50107773E-01, 0.52032378E-02, 0.17350024E-01,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x41 = x4
      x42 = x41*x4
      x51 = x5
c
c                  function
c
      l_rm_dext   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)*x11                
     7  +coeff(  7)    *x21    *x42    
c
      return
      end
      function x_rm_q3en   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.4311523E-03/
      data xmin/
     1 -0.99954E-03,-0.62955E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.37071069E-02, 0.22661927E+00, 0.13900360E+00, 0.20298548E-01,
     +  0.37430096E-01,-0.58411960E-02,-0.27363829E-02,-0.47753593E-02,
     + -0.55103097E-02,-0.83406561E-03,-0.16667477E-02, 0.31641915E-02,
     + -0.10121265E-01,-0.24058043E-03,-0.12430379E-02,-0.23209283E-02,
     + -0.33555520E-03,-0.46171038E-02, 0.91404945E-04,-0.11681682E-03,
     +  0.12887787E-03,-0.27201517E-03, 0.35305441E-03,-0.85058511E-03,
     + -0.60264184E-02,-0.67297756E-02,-0.36045909E-03, 0.20443420E-04,
     +  0.89093152E-03, 0.53482858E-03,-0.20293887E-02,-0.62706869E-03,
     +  0.13194871E-02,-0.46287110E-03,-0.36096777E-03, 0.28528983E-04,
     + -0.12704886E-03,-0.14893970E-03,-0.30889161E-03, 0.19005844E-02,
     +  0.10371022E-04, 0.20249076E-04, 0.16992624E-03, 0.23089153E-03,
     + -0.47893950E-03,-0.68329959E-04,-0.28120706E-03, 0.11775459E-03,
     + -0.86062275E-04, 0.12894980E-03, 0.76653436E-04,-0.18133247E-04,
     + -0.14367489E-03,-0.15815790E-02,-0.78393688E-03, 0.15969483E-04,
     +  0.11893174E-03,-0.10669955E-03,-0.86154090E-04, 0.68544832E-04,
     + -0.58488724E-04,-0.71734025E-05, 0.20957587E-03,-0.86190021E-05,
     +  0.21394128E-04, 0.13207966E-03, 0.55377539E-04, 0.78409800E-03,
     +  0.20773386E-04, 0.41642104E-03, 0.42607120E-03, 0.46197057E-03,
     +  0.13218988E-02,-0.45281832E-03,-0.51693502E-03, 0.13245686E-03,
     + -0.60005405E-04, 0.33984624E-03,-0.17921887E-03,-0.15289700E-02,
     + -0.74020890E-03,-0.15678506E-02,-0.20317784E-03, 0.39086698E-03,
     + -0.10673234E-02, 0.37554681E-03, 0.93851355E-04,-0.53080841E-03,
     +  0.37305121E-03, 0.35723412E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x34 = x33*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_rm_q3en   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)                *x52
     7  +coeff(  7)*x11                
     8  +coeff(  8)            *x42    
      x_rm_q3en   =x_rm_q3en   
     9  +coeff(  9)    *x21    *x42    
     1  +coeff( 10)        *x31*x41    
     2  +coeff( 11)    *x21*x31*x41    
     3  +coeff( 12)    *x22        *x51
     4  +coeff( 13)    *x23    *x42    
     5  +coeff( 14)    *x21*x32        
     6  +coeff( 15)            *x42*x51
     7  +coeff( 16)    *x24            
     8  +coeff( 17)*x11*x21            
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 18)    *x23*x31*x41    
     1  +coeff( 19)            *x41    
     2  +coeff( 20)        *x32        
     3  +coeff( 21)    *x21    *x41    
     4  +coeff( 22)    *x21        *x52
     5  +coeff( 23)                *x53
     6  +coeff( 24)    *x22*x31*x41    
     7  +coeff( 25)    *x21*x31*x43    
     8  +coeff( 26)    *x21    *x44    
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 27)    *x23*x32*x42    
     1  +coeff( 28)    *x21*x31        
     2  +coeff( 29)    *x23            
     3  +coeff( 30)    *x23        *x51
     4  +coeff( 31)    *x21*x32*x42    
     5  +coeff( 32)    *x24        *x51
     6  +coeff( 33)    *x22    *x42*x51
     7  +coeff( 34)    *x23*x32        
     8  +coeff( 35)    *x21*x33*x41    
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 36)*x11            *x51
     1  +coeff( 37)*x11*x22            
     2  +coeff( 38)    *x24*x32        
     3  +coeff( 39)            *x42*x54
     4  +coeff( 40)    *x23*x32*x44*x52
     5  +coeff( 41)        *x31        
     6  +coeff( 42)            *x41*x51
     7  +coeff( 43)    *x23    *x41    
     8  +coeff( 44)    *x21    *x43    
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 45)    *x22        *x52
     1  +coeff( 46)        *x31*x41*x52
     2  +coeff( 47)    *x23        *x52
     3  +coeff( 48)    *x23*x33        
     4  +coeff( 49)    *x23*x32*x41    
     5  +coeff( 50)    *x23*x31*x42    
     6  +coeff( 51)*x11*x23            
     7  +coeff( 52)    *x23*x32    *x52
     8  +coeff( 53)*x11*x24    *x41    
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 54)*x11*x24*x31*x42*x51
     1  +coeff( 55)*x11*x24    *x43*x51
     2  +coeff( 56)        *x31*x41*x51
     3  +coeff( 57)    *x21*x31*x42    
     4  +coeff( 58)        *x32*x42    
     5  +coeff( 59)        *x31*x43    
     6  +coeff( 60)    *x21*x32    *x51
     7  +coeff( 61)    *x21*x31*x41*x51
     8  +coeff( 62)*x11    *x31        
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 63)    *x22*x31*x41*x51
     1  +coeff( 64)            *x42*x53
     2  +coeff( 65)    *x24    *x42    
     3  +coeff( 66)    *x22*x32*x42    
     4  +coeff( 67)        *x34*x42    
     5  +coeff( 68)    *x22    *x44    
     6  +coeff( 69)        *x32*x44    
     7  +coeff( 70)    *x21*x31*x43*x51
     8  +coeff( 71)    *x21    *x44*x51
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 72)    *x24        *x52
     1  +coeff( 73)    *x22    *x42*x52
     2  +coeff( 74)            *x44*x52
     3  +coeff( 75)    *x22*x31*x43*x51
     4  +coeff( 76)    *x23        *x54
     5  +coeff( 77)*x11*x23*x31        
     6  +coeff( 78)    *x23*x34*x41    
     7  +coeff( 79)    *x24*x32*x42    
     8  +coeff( 80)    *x24    *x44    
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 81)    *x23    *x44*x51
     1  +coeff( 82)    *x24    *x42*x52
     2  +coeff( 83)    *x23*x31*x42*x52
     3  +coeff( 84)            *x44*x54
     4  +coeff( 85)    *x24    *x44*x51
     5  +coeff( 86)    *x22    *x44*x53
     6  +coeff( 87)    *x24    *x41*x54
     7  +coeff( 88)    *x21*x32*x42*x54
     8  +coeff( 89)    *x21    *x44*x54
      x_rm_q3en   =x_rm_q3en   
     9  +coeff( 90)*x11*x22*x31*x42*x51
c
      return
      end
      function t_rm_q3en   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 17)
      data ncoeff/ 16/
      data avdat/ -0.4582219E-02/
      data xmin/
     1 -0.99954E-03,-0.62955E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.34629928E-02,-0.79486005E-01, 0.22677328E-01,-0.51691062E-02,
     + -0.58313254E-02, 0.68446755E-03,-0.13894757E-02,-0.25861390E-03,
     +  0.78996830E-03, 0.10637202E-02, 0.60188817E-03, 0.55687519E-03,
     +  0.29148525E-04, 0.10867831E-02, 0.63292385E-03, 0.22663670E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x31 = x3
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      t_rm_q3en   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)*x11                
     7  +coeff(  7)                *x52
     8  +coeff(  8)        *x31*x41    
      t_rm_q3en   =t_rm_q3en   
     9  +coeff(  9)            *x42    
     1  +coeff( 10)    *x21    *x42    
     2  +coeff( 11)    *x21        *x52
     3  +coeff( 12)    *x23*x31*x41    
     4  +coeff( 13)    *x23            
     5  +coeff( 14)    *x21*x31*x41    
     6  +coeff( 15)            *x42*x51
     7  +coeff( 16)    *x23    *x42    
c
      return
      end
      function y_rm_q3en   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 33)
      data ncoeff/ 32/
      data avdat/ -0.6854549E-03/
      data xmin/
     1 -0.99954E-03,-0.62955E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.62492944E-03,-0.26635434E-01, 0.13582335E+00,-0.13228284E-05,
     +  0.21789431E-03,-0.87099792E-02,-0.77386186E-01, 0.11712763E-01,
     +  0.65113142E-01,-0.40081237E-01,-0.96082306E-02,-0.50757225E-02,
     + -0.79012575E-03,-0.15358298E-03, 0.68391621E-03,-0.62887254E-03,
     + -0.75485434E-04,-0.52858596E-02,-0.64172600E-02, 0.70126499E-02,
     + -0.94705634E-03,-0.74896608E-02, 0.62810624E-03,-0.40127854E-02,
     + -0.82572835E-04,-0.99290744E-03, 0.38523498E-03,-0.11126071E-02,
     +  0.39883684E-02, 0.51329809E-03,-0.49974974E-02, 0.35373492E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x45 = x44*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      y_rm_q3en   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)        *x32        
     5  +coeff(  5)            *x42    
     6  +coeff(  6)    *x21*x31        
     7  +coeff(  7)    *x21    *x41    
     8  +coeff(  8)        *x31    *x51
      y_rm_q3en   =y_rm_q3en   
     9  +coeff(  9)            *x41*x51
     1  +coeff( 10)    *x22    *x41    
     2  +coeff( 11)    *x21    *x41*x51
     3  +coeff( 12)            *x41*x52
     4  +coeff( 13)            *x45    
     5  +coeff( 14)    *x22*x33        
     6  +coeff( 15)    *x21            
     7  +coeff( 16)                *x51
     8  +coeff( 17)        *x33        
      y_rm_q3en   =y_rm_q3en   
     9  +coeff( 18)        *x31*x42    
     1  +coeff( 19)    *x22*x31        
     2  +coeff( 20)    *x23    *x41    
     3  +coeff( 21)        *x32*x41    
     4  +coeff( 22)            *x43    
     5  +coeff( 23)*x11        *x41    
     6  +coeff( 24)    *x22    *x41*x51
     7  +coeff( 25)    *x21*x33    *x51
     8  +coeff( 26)        *x31    *x54
      y_rm_q3en   =y_rm_q3en   
     9  +coeff( 27)    *x22            
     1  +coeff( 28)    *x21*x31    *x51
     2  +coeff( 29)    *x21    *x43    
     3  +coeff( 30)*x11*x21    *x41    
     4  +coeff( 31)    *x22    *x43    
     5  +coeff( 32)    *x23*x31*x42    
c
      return
      end
      function p_rm_q3en   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 38)
      data ncoeff/ 37/
      data avdat/ -0.2050469E-03/
      data xmin/
     1 -0.99954E-03,-0.62955E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.85091749E-04, 0.13220677E-03,-0.74262372E-02, 0.13292775E-01,
     + -0.26657275E-06,-0.24366082E-05,-0.14518440E-01, 0.20092814E-02,
     +  0.13601062E-01,-0.96007567E-02,-0.15204048E-02,-0.24994400E-02,
     + -0.76610770E-04,-0.29488378E-04, 0.39744253E-04,-0.90694253E-03,
     + -0.12990909E-03,-0.22051390E-02,-0.11733577E-02,-0.88393723E-03,
     +  0.16142093E-02,-0.55409886E-03,-0.10196324E-02, 0.89103036E-04,
     +  0.10245007E-03,-0.19454979E-03,-0.17018689E-03,-0.38304125E-03,
     + -0.11648545E-02, 0.10154011E-02, 0.84695488E-03, 0.70741895E-03,
     +  0.26120892E-04, 0.13726339E-03, 0.45826996E-03, 0.56729332E-03,
     +  0.55235164E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
c
c                  function
c
      p_rm_q3en   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)            *x41    
     5  +coeff(  5)*x11                
     6  +coeff(  6)        *x32        
     7  +coeff(  7)    *x21    *x41    
     8  +coeff(  8)        *x31    *x51
      p_rm_q3en   =p_rm_q3en   
     9  +coeff(  9)            *x41*x51
     1  +coeff( 10)    *x22    *x41    
     2  +coeff( 11)            *x43    
     3  +coeff( 12)    *x21    *x41*x51
     4  +coeff( 13)    *x23*x31        
     5  +coeff( 14)    *x21*x33        
     6  +coeff( 15)    *x23*x33        
     7  +coeff( 16)    *x24*x31*x42    
     8  +coeff( 17)                *x51
      p_rm_q3en   =p_rm_q3en   
     9  +coeff( 18)    *x21*x31        
     1  +coeff( 19)    *x22*x31        
     2  +coeff( 20)        *x31*x42    
     3  +coeff( 21)    *x23    *x41    
     4  +coeff( 22)            *x41*x52
     5  +coeff( 23)    *x22    *x41*x51
     6  +coeff( 24)    *x22            
     7  +coeff( 25)*x11        *x41    
     8  +coeff( 26)        *x32*x41    
      p_rm_q3en   =p_rm_q3en   
     9  +coeff( 27)    *x21*x31    *x51
     1  +coeff( 28)    *x22*x31    *x51
     2  +coeff( 29)    *x22    *x43    
     3  +coeff( 30)    *x23    *x41*x51
     4  +coeff( 31)    *x22    *x41*x52
     5  +coeff( 32)    *x23    *x41*x52
     6  +coeff( 33)    *x21        *x51
     7  +coeff( 34)*x11*x21    *x41    
     8  +coeff( 35)    *x21    *x43    
      p_rm_q3en   =p_rm_q3en   
     9  +coeff( 36)    *x23*x31*x42    
     1  +coeff( 37)    *x23    *x41*x53
c
      return
      end
      function l_rm_q3en   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff(  9)
      data ncoeff/  8/
      data avdat/ -0.1510160E-01/
      data xmin/
     1 -0.99954E-03,-0.62955E-01,-0.19979E-01,-0.32934E-01,-0.49893E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49990E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.10176132E-01,-0.32301199E+00,-0.32675438E-01,-0.30403778E-01,
     + -0.23216445E-01, 0.31609200E-02,-0.55845478E-02, 0.10939431E-01,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x41 = x4
      x42 = x41*x4
      x51 = x5
c
c                  function
c
      l_rm_q3en   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)*x11                
     7  +coeff(  7)            *x42    
     8  +coeff(  8)    *x21    *x42    
c
      return
      end
      function x_rm_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.1157085E-01/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.13848802E-01, 0.70274472E-01, 0.31127551E+00, 0.35810959E-01,
     + -0.20400297E-01, 0.13229701E-01,-0.59033968E-02,-0.11484828E-01,
     + -0.19109958E-02,-0.24746181E-02,-0.12280558E-02,-0.71138871E-03,
     + -0.45938487E-02, 0.67913057E-02, 0.23004657E-02, 0.50693131E-02,
     +  0.11888699E-03, 0.23938922E-03,-0.33376072E-03, 0.37227606E-02,
     +  0.44959955E-03,-0.15995982E-02,-0.15192435E-02,-0.30846912E-02,
     + -0.17503517E-02,-0.47919934E-03,-0.15636432E-03,-0.79915207E-03,
     + -0.72794827E-03,-0.27149410E-02,-0.44972422E-02,-0.38748381E-02,
     + -0.98236895E-04,-0.29960773E-02, 0.24594739E-04, 0.46549321E-04,
     +  0.91227214E-03, 0.35140172E-04,-0.12614741E-02,-0.25478099E-03,
     + -0.32231377E-02,-0.62350584E-04, 0.15118547E-02,-0.72543509E-04,
     +  0.17051271E-03,-0.94384147E-03,-0.36673184E-03, 0.44097360E-04,
     + -0.15973669E-03,-0.32723113E-03,-0.10145549E-03,-0.19352895E-03,
     + -0.32342540E-03,-0.26643727E-03,-0.13070402E-02, 0.49587688E-03,
     +  0.20815988E-03, 0.49638795E-04,-0.98451134E-03, 0.11420117E-02,
     + -0.50929823E-03,-0.27177934E-03, 0.10877625E-02,-0.22286796E-02,
     + -0.15492557E-02, 0.40434365E-03,-0.76160533E-03, 0.15390204E-03,
     +  0.18008514E-04, 0.45771262E-03,-0.69817383E-03,-0.19468551E-04,
     +  0.36787460E-04,-0.53047253E-04,-0.53242678E-04, 0.12605012E-03,
     + -0.62544415E-04,-0.11490897E-03, 0.62537511E-04,-0.17671879E-03,
     +  0.17178849E-03,-0.45834412E-03, 0.17245789E-03, 0.16570285E-03,
     +  0.31739098E-03, 0.26564418E-04,-0.92158305E-04,-0.52301595E-04,
     +  0.53682965E-04, 0.44706443E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x34 = x33*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_rm_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)                *x52
     6  +coeff(  6)    *x22            
     7  +coeff(  7)            *x42    
     8  +coeff(  8)    *x21    *x42    
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff(  9)*x11                
     1  +coeff( 10)        *x31*x41    
     2  +coeff( 11)    *x21*x31*x41    
     3  +coeff( 12)    *x22        *x53
     4  +coeff( 13)    *x24            
     5  +coeff( 14)    *x22    *x42    
     6  +coeff( 15)                *x53
     7  +coeff( 16)    *x22    *x42*x51
     8  +coeff( 17)            *x41    
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 18)    *x21    *x41    
     1  +coeff( 19)    *x21*x32        
     2  +coeff( 20)    *x22        *x51
     3  +coeff( 21)        *x31*x41*x51
     4  +coeff( 22)    *x23        *x51
     5  +coeff( 23)    *x21    *x42*x51
     6  +coeff( 24)    *x24        *x51
     7  +coeff( 25)    *x23        *x52
     8  +coeff( 26)    *x24        *x52
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 27)        *x32        
     1  +coeff( 28)    *x22*x31*x41    
     2  +coeff( 29)    *x21*x31*x41*x51
     3  +coeff( 30)    *x23*x31*x41    
     4  +coeff( 31)    *x23    *x42    
     5  +coeff( 32)    *x21    *x44    
     6  +coeff( 33)*x11            *x51
     7  +coeff( 34)    *x24    *x42    
     8  +coeff( 35)        *x31        
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 36)    *x21*x31        
     1  +coeff( 37)    *x23            
     2  +coeff( 38)    *x21        *x52
     3  +coeff( 39)    *x22        *x52
     4  +coeff( 40)                *x54
     5  +coeff( 41)    *x21*x31*x43    
     6  +coeff( 42)        *x34    *x51
     7  +coeff( 43)    *x21    *x42*x52
     8  +coeff( 44)*x11*x22            
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 45)    *x21*x34    *x51
     1  +coeff( 46)    *x23    *x42*x51
     2  +coeff( 47)            *x44*x52
     3  +coeff( 48)    *x23        *x53
     4  +coeff( 49)    *x23*x32*x42    
     5  +coeff( 50)            *x44*x53
     6  +coeff( 51)    *x22    *x41    
     7  +coeff( 52)            *x42*x51
     8  +coeff( 53)    *x21        *x53
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 54)    *x21*x33*x41    
     1  +coeff( 55)    *x21*x32*x42    
     2  +coeff( 56)    *x24*x31*x41    
     3  +coeff( 57)    *x23*x32*x41    
     4  +coeff( 58)*x11        *x42    
     5  +coeff( 59)    *x22    *x44    
     6  +coeff( 60)    *x22    *x42*x52
     7  +coeff( 61)            *x42*x54
     8  +coeff( 62)*x11*x21    *x42    
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 63)    *x23    *x44    
     1  +coeff( 64)    *x24    *x42*x51
     2  +coeff( 65)    *x23    *x42*x52
     3  +coeff( 66)    *x23*x33*x42    
     4  +coeff( 67)    *x24*x31*x43    
     5  +coeff( 68)*x11*x23*x32        
     6  +coeff( 69)*x11*x21*x34        
     7  +coeff( 70)*x11*x23    *x42    
     8  +coeff( 71)    *x24*x33*x42    
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 72)        *x32*x41    
     1  +coeff( 73)    *x21    *x41*x51
     2  +coeff( 74)    *x23*x31        
     3  +coeff( 75)    *x21*x33        
     4  +coeff( 76)    *x21    *x43    
     5  +coeff( 77)        *x31*x43    
     6  +coeff( 78)    *x22    *x41*x51
     7  +coeff( 79)        *x32    *x52
     8  +coeff( 80)    *x23*x32        
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 81)    *x22*x31*x42    
     1  +coeff( 82)    *x22*x31*x41*x51
     2  +coeff( 83)    *x23*x33        
     3  +coeff( 84)    *x22*x32*x41*x51
     4  +coeff( 85)    *x21*x31*x43*x51
     5  +coeff( 86)*x11            *x52
     6  +coeff( 87)*x11*x21*x32        
     7  +coeff( 88)*x11*x22    *x41    
     8  +coeff( 89)*x11*x22        *x51
      x_rm_q3ex   =x_rm_q3ex   
     9  +coeff( 90)    *x24*x31*x41*x51
c
      return
      end
      function t_rm_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 14)
      data ncoeff/ 13/
      data avdat/ -0.6260102E-02/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.61263000E-02,-0.29467762E-01, 0.97874917E-01, 0.41515664E-02,
     + -0.92655942E-02,-0.30138511E-02,-0.84603863E-03,-0.12819642E-02,
     +  0.93882711E-03, 0.17303065E-03,-0.61184086E-03,-0.11568579E-02,
     +  0.19215203E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
c
c                  function
c
      t_rm_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x21        *x51
     5  +coeff(  5)                *x52
     6  +coeff(  6)    *x21    *x42    
     7  +coeff(  7)        *x31*x41    
     8  +coeff(  8)            *x42    
      t_rm_q3ex   =t_rm_q3ex   
     9  +coeff(  9)                *x53
     1  +coeff( 10)    *x22            
     2  +coeff( 11)    *x21        *x52
     3  +coeff( 12)    *x24            
     4  +coeff( 13)    *x22    *x42    
c
      return
      end
      function y_rm_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 34)
      data ncoeff/ 33/
      data avdat/ -0.7236771E-03/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.20925411E-03,-0.29320149E-01, 0.93478397E-01, 0.42835768E-05,
     +  0.21208262E-03,-0.90395678E-02,-0.69108874E-01, 0.91163926E-02,
     +  0.64450100E-01, 0.65700922E-04,-0.40117275E-01,-0.13032256E-01,
     +  0.71424101E-02,-0.46550017E-03,-0.12775412E-03,-0.29086710E-02,
     +  0.61017857E-03,-0.62218815E-03,-0.40116832E-02,-0.67059672E-02,
     + -0.50435974E-02,-0.22955372E-02,-0.62170331E-02,-0.96669170E-03,
     +  0.51887892E-03,-0.10145019E-02, 0.40288866E-03, 0.26684057E-02,
     +  0.56494196E-03,-0.12917406E-02,-0.60282582E-02, 0.26683060E-02,
     +  0.14761499E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x45 = x44*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      y_rm_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)        *x32        
     5  +coeff(  5)            *x42    
     6  +coeff(  6)    *x21*x31        
     7  +coeff(  7)    *x21    *x41    
     8  +coeff(  8)        *x31    *x51
      y_rm_q3ex   =y_rm_q3ex   
     9  +coeff(  9)            *x41*x51
     1  +coeff( 10)            *x42*x51
     2  +coeff( 11)    *x22    *x41    
     3  +coeff( 12)    *x21    *x41*x51
     4  +coeff( 13)    *x23    *x41    
     5  +coeff( 14)            *x45    
     6  +coeff( 15)    *x22*x33        
     7  +coeff( 16)    *x22*x31*x42    
     8  +coeff( 17)    *x21            
      y_rm_q3ex   =y_rm_q3ex   
     9  +coeff( 18)                *x51
     1  +coeff( 19)        *x31*x42    
     2  +coeff( 20)            *x43    
     3  +coeff( 21)    *x22*x31        
     4  +coeff( 22)            *x41*x52
     5  +coeff( 23)    *x22    *x41*x51
     6  +coeff( 24)        *x32*x41    
     7  +coeff( 25)*x11        *x41    
     8  +coeff( 26)    *x21*x31    *x51
      y_rm_q3ex   =y_rm_q3ex   
     9  +coeff( 27)    *x22            
     1  +coeff( 28)    *x21    *x43    
     2  +coeff( 29)*x11*x21    *x41    
     3  +coeff( 30)    *x22*x31    *x51
     4  +coeff( 31)    *x22    *x43    
     5  +coeff( 32)    *x23    *x41*x51
     6  +coeff( 33)    *x21*x33*x42    
c
      return
      end
      function p_rm_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 26)
      data ncoeff/ 25/
      data avdat/  0.1299790E-03/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.24184867E-03,-0.17650127E-03, 0.56292210E-02,-0.40408503E-01,
     +  0.19912715E-02, 0.19825120E-01,-0.38016809E-02,-0.13797985E-01,
     +  0.91066975E-02, 0.17996674E-02, 0.18971288E-02, 0.17311068E-02,
     +  0.14438445E-02, 0.87936921E-03, 0.13845197E-03,-0.17789770E-03,
     +  0.56071556E-03,-0.15386158E-02,-0.12273728E-02, 0.22020182E-03,
     +  0.26408283E-03,-0.10422838E-03,-0.88568253E-03, 0.13739286E-02,
     + -0.89020299E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
c
c                  function
c
      p_rm_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)            *x41    
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      p_rm_q3ex   =p_rm_q3ex   
     9  +coeff(  9)    *x22    *x41    
     1  +coeff( 10)    *x22*x31        
     2  +coeff( 11)            *x43    
     3  +coeff( 12)            *x41*x52
     4  +coeff( 13)        *x31*x42    
     5  +coeff( 14)        *x31    *x52
     6  +coeff( 15)                *x51
     7  +coeff( 16)*x11        *x41    
     8  +coeff( 17)    *x21    *x41*x51
      p_rm_q3ex   =p_rm_q3ex   
     9  +coeff( 18)    *x23    *x41    
     1  +coeff( 19)            *x43*x53
     2  +coeff( 20)        *x32*x41    
     3  +coeff( 21)    *x21*x31    *x51
     4  +coeff( 22)    *x24            
     5  +coeff( 23)    *x21    *x43    
     6  +coeff( 24)    *x22    *x43    
     7  +coeff( 25)    *x23*x31*x42    
c
      return
      end
      function l_rm_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff(  9)
      data ncoeff/  8/
      data avdat/ -0.1960552E-01/
      data xmin/
     1 -0.99954E-03,-0.62518E-01,-0.19967E-01,-0.32934E-01,-0.49651E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.99924E-03, 0.64732E-01, 0.20000E-01, 0.32298E-01, 0.49986E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.11299579E-01,-0.32281932E+00,-0.32077547E-01,-0.34536079E-01,
     + -0.14496966E-01, 0.32102165E-02,-0.18219585E-01, 0.11492895E-01,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x21 = x2
      x22 = x21*x2
      x41 = x4
      x42 = x41*x4
      x51 = x5
      x52 = x51*x5
c
c                  function
c
      l_rm_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)    *x22            
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)*x11                
     7  +coeff(  7)            *x42*x52
     8  +coeff(  8)    *x21    *x42    
c
      return
      end
