real function LTXFIT (x,m) dimension x(m) dimension xmin(10),xmax(10),scale(10)!,xmean(10) dimension coeff( 3) data ncoeff/ 2/ data avdat/ -0.3687619E-02/ data xmin/ 1 -0.57780E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data xmax/ 1 0.50125E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data scale /10*0./ data coeff/ + -0.31657170E-02, 0.89058414E-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) c set up monomials functions x11 = x1 c c function c LTXFIT =avdat 1 +coeff( 1) 2 +coeff( 2)*x11 c return end real function ldelta (x,m) dimension x(m) dimension xmin(10),xmax(10),scale(10) dimension coeff( 76) data ncoeff/ 75/ data avdat/ -0.1882823E-02/ data xmin/ 1 -0.69010E+00,-0.28144E-01,-0.49814E-01,-0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data xmax/ 1 0.54634E+00, 0.24553E-01, 0.49814E-01, 0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data scale /10*0./ data coeff/ + -0.39884862E-02, 0.50622913E-01,-0.12377950E-02, 0.40097270E-02, + 0.39700102E-02,-0.11728850E-02, 0.10129562E-02, 0.12215863E-02, + 0.10103060E-03, 0.38528920E-03, 0.25860510E-03, 0.27198652E-03, + -0.12176270E-03,-0.56424102E-03, 0.86635900E-03,-0.97910120E-03, + 0.71599050E-03,-0.11752060E-02, 0.59314200E-03,-0.86932464E-04, + -0.35982830E-03,-0.83722552E-04,-0.86343750E-03,-0.23961050E-04, + 0.53753033E-05,-0.17723300E-03,-0.76363590E-03,-0.75737080E-03, + -0.91090620E-03,-0.43663560E-04,-0.71769364E-03, 0.34288471E-03, + 0.20210870E-04, 0.24764690E-03, 0.16725961E-03,-0.68604580E-03, + 0.36006620E-03,-0.79142430E-03, 0.29020872E-03,-0.37886980E-04, + 0.93915630E-04, 0.13965122E-03, 0.74772513E-04,-0.23570310E-03, + -0.60953380E-03, 0.23991301E-05,-0.47437360E-03, 0.62906343E-04, + -0.72657584E-03, 0.38419023E-03,-0.65336270E-03,-0.23461560E-03, + -0.11416872E-02,-0.30595992E-03, 0.22321100E-03,-0.33125851E-03, + 0.17155780E-03,-0.19372800E-03,-0.23453880E-03, 0.28205500E-04, + -0.50900413E-04, 0.30483310E-03, 0.58151203E-04, 0.33781474E-04, + -0.62264920E-04, 0.15605430E-03, 0.13921063E-03,-0.17766840E-03, + -0.20115260E-03, 0.14353590E-04, 0.19674450E-03, 0.28641452E-03, + -0.25983544E-03, 0.20277923E-03,-0.36867454E-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) c set up monomials functions x11 = x1 x12 = x11*x1 x13 = x12*x1 x14 = x13*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 c c function c ldelta =avdat 1 +coeff( 1) 2 +coeff( 2)*x11 3 +coeff( 3) *x21 4 +coeff( 4)*x12 5 +coeff( 5)*x11*x21 6 +coeff( 6) *x22 7 +coeff( 7) *x32 8 +coeff( 8) *x31*x41 ldelta =ldelta 9 +coeff( 9) *x42 1 +coeff( 10)*x13 2 +coeff( 11)*x12*x21 3 +coeff( 12)*x11*x22 4 +coeff( 13)*x11 *x32 5 +coeff( 14) *x21*x32 6 +coeff( 15)*x11 *x31*x41 7 +coeff( 16) *x21*x31*x41 8 +coeff( 17)*x11 *x42 ldelta =ldelta 9 +coeff( 18) *x21 *x42 1 +coeff( 19) *x24 2 +coeff( 20)*x12 *x31*x41 3 +coeff( 21)*x11*x21*x31*x41 4 +coeff( 22)*x11*x21 *x42 5 +coeff( 23) *x22 *x42 6 +coeff( 24) *x23 7 +coeff( 25)*x13*x21 8 +coeff( 26)*x11*x23 ldelta =ldelta 9 +coeff( 27) *x22*x32 1 +coeff( 28) *x22*x31*x41 2 +coeff( 29)*x11*x22*x31*x41 3 +coeff( 30)*x12*x21 *x42 4 +coeff( 31)*x11*x22 *x42 5 +coeff( 32) *x23 *x42 6 +coeff( 33)*x14*x22 7 +coeff( 34)*x11*x21*x32 8 +coeff( 35)*x11*x24 ldelta =ldelta 9 +coeff( 36) *x23*x32 1 +coeff( 37)*x12*x21*x31*x41 2 +coeff( 38) *x21*x33*x41 3 +coeff( 39) *x24 *x42 4 +coeff( 40) *x23*x34 5 +coeff( 41)*x14 6 +coeff( 42)*x12*x22 7 +coeff( 43) *x32*x42 8 +coeff( 44)*x12*x23 ldelta =ldelta 9 +coeff( 45) *x21*x34 1 +coeff( 46)*x13 *x31*x41 2 +coeff( 47) *x21*x32*x42 3 +coeff( 48)*x11 *x44 4 +coeff( 49)*x11*x23*x31*x41 5 +coeff( 50) *x24*x31*x41 6 +coeff( 51)*x11*x21*x33*x41 7 +coeff( 52)*x11*x23 *x42 8 +coeff( 53)*x11*x21*x32*x42 ldelta =ldelta 9 +coeff( 54)*x12 *x31*x43 1 +coeff( 55)*x13*x24 2 +coeff( 56)*x14*x21*x32 3 +coeff( 57)*x13*x22*x32 4 +coeff( 58)*x13 *x31*x43 5 +coeff( 59)*x14 *x33*x41 6 +coeff( 60)*x12 *x42 7 +coeff( 61)*x13 *x32 8 +coeff( 62)*x11*x22*x32 ldelta =ldelta 9 +coeff( 63) *x23*x31*x41 1 +coeff( 64)*x13 *x42 2 +coeff( 65)*x11 *x32*x42 3 +coeff( 66)*x11*x23*x32 4 +coeff( 67)*x11*x21*x34 5 +coeff( 68) *x22*x34 6 +coeff( 69)*x14 *x31*x41 7 +coeff( 70)*x13*x21*x31*x41 8 +coeff( 71)*x12*x22*x31*x41 ldelta =ldelta 9 +coeff( 72) *x22*x32*x42 1 +coeff( 73)*x12*x23 *x42 2 +coeff( 74)*x11*x24 *x42 3 +coeff( 75)*x11*x22*x32*x42 c return end real function ltheta (x,m) dimension x(m) dimension xmin(10),xmax(10),scale(10) dimension coeff( 76) data ncoeff/ 75/ data avdat/ 0.2427104E-02/ data xmin/ 1 -0.69010E+00,-0.28144E-01,-0.49814E-01,-0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data xmax/ 1 0.54634E+00, 0.24553E-01, 0.49814E-01, 0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data scale /10*0./ data coeff/ + 0.28862934E-02, 0.26194870E-03,-0.60644420E-01, 0.80566080E-02, + -0.43439270E-02, 0.42614890E-03, 0.38092610E-03, 0.74502330E-03, + -0.60916460E-03,-0.14864461E-02, 0.19036090E-02,-0.70173683E-03, + -0.67744503E-03, 0.47386530E-03,-0.14959230E-02,-0.47369170E-04, + 0.48181480E-04, 0.27729824E-03,-0.86152820E-03, 0.82670190E-03, + -0.38999461E-03, 0.75465090E-03,-0.12472630E-02,-0.17691211E-03, + -0.83196560E-03, 0.13777012E-03,-0.59001100E-03,-0.97239244E-04, + -0.66026080E-03, 0.75838404E-04, 0.93582460E-03,-0.14742370E-02, + -0.58090510E-03,-0.11790260E-03, 0.12113151E-03,-0.18503771E-03, + 0.14573951E-03, 0.11403491E-03, 0.52083670E-03, 0.35762650E-03, + 0.96164990E-04,-0.14210820E-03, 0.16858090E-04,-0.34694172E-03, + 0.19379613E-03,-0.62047560E-03, 0.68812492E-04, 0.29304810E-03, + -0.30095193E-02,-0.39137870E-04,-0.25663571E-02,-0.24104882E-02, + -0.24849441E-02, 0.90858864E-03,-0.16672342E-02,-0.32029822E-02, + 0.10667280E-02,-0.35314880E-02,-0.13071924E-03,-0.16768360E-03, + -0.13131740E-03, 0.10484060E-03, 0.20680910E-04,-0.17647470E-03, + -0.44763020E-03,-0.65231014E-03, 0.55418262E-03,-0.22222290E-02, + 0.10022510E-02,-0.56763580E-03, 0.56457793E-03,-0.54451090E-03, + 0.96533420E-03,-0.12792642E-02, 0.75331090E-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) c set up monomials functions x11 = x1 x12 = x11*x1 x13 = x12*x1 x14 = x13*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 c c function c ltheta =avdat 1 +coeff( 1) 2 +coeff( 2)*x11 3 +coeff( 3) *x21 4 +coeff( 4)*x11*x21 5 +coeff( 5) *x22 6 +coeff( 6) *x32 7 +coeff( 7) *x31*x41 8 +coeff( 8) *x42 ltheta =ltheta 9 +coeff( 9)*x13 1 +coeff( 10)*x12*x21 2 +coeff( 11)*x11*x22 3 +coeff( 12) *x23 4 +coeff( 13) *x21*x32 5 +coeff( 14)*x11 *x31*x41 6 +coeff( 15) *x21*x31*x41 7 +coeff( 16)*x11 *x42 8 +coeff( 17)*x14 ltheta =ltheta 9 +coeff( 18)*x13*x21 1 +coeff( 19)*x12*x22 2 +coeff( 20)*x11*x23 3 +coeff( 21)*x12 *x32 4 +coeff( 22)*x11*x21*x32 5 +coeff( 23) *x22*x32 6 +coeff( 24)*x12 *x31*x41 7 +coeff( 25)*x11*x21*x31*x41 8 +coeff( 26)*x12 *x42 ltheta =ltheta 9 +coeff( 27)*x11*x21 *x42 1 +coeff( 28)*x13*x22 2 +coeff( 29)*x12*x21*x32 3 +coeff( 30)*x13 *x31*x41 4 +coeff( 31)*x12*x21*x31*x41 5 +coeff( 32)*x11*x22*x31*x41 6 +coeff( 33)*x11*x22 *x42 7 +coeff( 34)*x11 *x44 8 +coeff( 35)*x14*x22 ltheta =ltheta 9 +coeff( 36)*x13*x23 1 +coeff( 37)*x12*x24 2 +coeff( 38)*x14 *x32 3 +coeff( 39)*x12*x22*x32 4 +coeff( 40)*x12 *x34 5 +coeff( 41)*x14 *x42 6 +coeff( 42)*x12*x22 *x42 7 +coeff( 43) *x23*x34 8 +coeff( 44)*x13 *x44 ltheta =ltheta 9 +coeff( 45)*x12*x21 *x44 1 +coeff( 46)*x14*x22*x32 2 +coeff( 47)*x14*x22 *x42 3 +coeff( 48) *x24 *x44 4 +coeff( 49)*x12 5 +coeff( 50) *x21 *x42 6 +coeff( 51) *x23*x32 7 +coeff( 52) *x21*x34 8 +coeff( 53) *x21*x33*x41 ltheta =ltheta 9 +coeff( 54)*x11*x21*x34 1 +coeff( 55)*x11*x23*x31*x41 2 +coeff( 56)*x11*x21*x32*x42 3 +coeff( 57) *x21*x34*x42 4 +coeff( 58)*x13*x21*x33*x41 5 +coeff( 59)*x11 *x32 6 +coeff( 60) *x24 7 +coeff( 61) *x22*x31*x41 8 +coeff( 62) *x22 *x42 ltheta =ltheta 9 +coeff( 63) *x32*x42 1 +coeff( 64)*x11*x24 2 +coeff( 65) *x23*x31*x41 3 +coeff( 66)*x13 *x42 4 +coeff( 67)*x12*x21 *x42 5 +coeff( 68) *x21*x32*x42 6 +coeff( 69)*x11*x23*x32 7 +coeff( 70) *x22*x34 8 +coeff( 71) *x24*x31*x41 ltheta =ltheta 9 +coeff( 72)*x11*x23 *x42 1 +coeff( 73) *x22*x32*x42 2 +coeff( 74)*x11*x21*x31*x43 3 +coeff( 75)*x13*x22*x32 c return end real function lphi (x,m) dimension x(m) dimension xmin(10),xmax(10),scale(10)!,xmean(10) dimension coeff( 76) data ncoeff/ 75/ data avdat/ 0.0000000E+00/ data xmin/ 1 -0.69010E+00,-0.28144E-01,-0.49814E-01,-0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data xmax/ 1 0.54634E+00, 0.24553E-01, 0.49814E-01, 0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data scale /10*0./ data coeff/ + -0.17487643E-01,-0.29655560E-01, 0.98397510E-02, 0.52609070E-02, + -0.52526821E-02, 0.72800270E-02,-0.79638370E-03, 0.86906173E-03, + -0.33054493E-02, 0.48561770E-02,-0.96981041E-03,-0.36170310E-03, + -0.10299570E-02,-0.10739100E-02, 0.11717530E-03, 0.16024122E-02, + -0.12350693E-02, 0.37981430E-03,-0.73734030E-03,-0.21286440E-02, + -0.27746080E-02, 0.73827170E-03,-0.10559620E-02,-0.50834060E-03, + -0.35608332E-03,-0.59472400E-03,-0.80524844E-03, 0.50698581E-03, + -0.87376902E-04,-0.16571274E-02, 0.29971050E-02, 0.11621890E-02, + -0.33671193E-03, 0.30749131E-03,-0.14976301E-02,-0.12771780E-02, + -0.41282970E-03,-0.74392280E-03, 0.12889981E-02,-0.15564570E-03, + 0.11562740E-02, 0.20422541E-04,-0.17378691E-03,-0.27197630E-02, + -0.18333783E-02,-0.12037600E-02,-0.78850891E-03,-0.17937870E-02, + 0.11084770E-02,-0.24426380E-02, 0.15780490E-02,-0.40722161E-03, + -0.27267434E-03,-0.54815900E-03,-0.79462080E-03,-0.20813270E-03, + 0.82824970E-04,-0.46812290E-02, 0.24091694E-02,-0.14005260E-02, + -0.15204910E-02, 0.43287270E-03,-0.82265920E-03, 0.18666980E-02, + -0.20100160E-02, 0.10665253E-02, 0.42962073E-02, 0.37922610E-03, + -0.94091950E-03, 0.68945984E-03,-0.11669840E-02, 0.36130550E-02, + 0.37542170E-03, 0.67739800E-03,-0.36846863E-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) c set up monomials functions x11 = x1 x12 = x11*x1 x13 = x12*x1 x14 = x13*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 c c function c lphi =avdat 1 +coeff( 1) *x31 2 +coeff( 2) *x41 3 +coeff( 3)*x11 *x31 4 +coeff( 4) *x21*x31 5 +coeff( 5)*x11 *x41 6 +coeff( 6) *x21 *x41 7 +coeff( 7)*x12 *x31 8 +coeff( 8)*x11*x21*x31 lphi =lphi 9 +coeff( 9) *x22*x31 1 +coeff( 10)*x12 *x41 2 +coeff( 11)*x11*x21 *x41 3 +coeff( 12) *x22 *x41 4 +coeff( 13) *x32*x41 5 +coeff( 14) *x31*x42 6 +coeff( 15) *x43 7 +coeff( 16)*x11*x22*x31 8 +coeff( 17) *x23*x31 lphi =lphi 9 +coeff( 18)*x11 *x33 1 +coeff( 19) *x21*x33 2 +coeff( 20)*x13 *x41 3 +coeff( 21)*x11*x22 *x41 4 +coeff( 22) *x23 *x41 5 +coeff( 23)*x11 *x32*x41 6 +coeff( 24)*x11 *x31*x42 7 +coeff( 25) *x21*x31*x42 8 +coeff( 26)*x11 *x43 lphi =lphi 9 +coeff( 27)*x14 *x31 1 +coeff( 28)*x13*x21*x31 2 +coeff( 29)*x12 *x33 3 +coeff( 30) *x22*x33 4 +coeff( 31)*x14 *x41 5 +coeff( 32)*x12*x22 *x41 6 +coeff( 33)*x11*x23 *x41 7 +coeff( 34) *x24 *x41 8 +coeff( 35) *x22*x32*x41 lphi =lphi 9 +coeff( 36)*x12 *x31*x42 1 +coeff( 37)*x11*x21*x31*x42 2 +coeff( 38) *x33*x42 3 +coeff( 39)*x12 *x43 4 +coeff( 40)*x11*x21 *x43 5 +coeff( 41) *x22 *x43 6 +coeff( 42) *x31*x44 7 +coeff( 43)*x12*x21*x33 8 +coeff( 44)*x14*x21 *x41 lphi =lphi 9 +coeff( 45)*x12*x23 *x41 1 +coeff( 46)*x11*x24 *x41 2 +coeff( 47)*x12*x21*x32*x41 3 +coeff( 48)*x11*x22*x32*x41 4 +coeff( 49)*x12*x21*x31*x42 5 +coeff( 50)*x11*x22*x31*x42 6 +coeff( 51) *x23*x31*x42 7 +coeff( 52)*x13 *x43 8 +coeff( 53)*x12*x21 *x43 lphi =lphi 9 +coeff( 54) *x23 *x43 1 +coeff( 55)*x12*x24*x31 2 +coeff( 56)*x14 *x33 3 +coeff( 57)*x14*x22 *x41 4 +coeff( 58)*x13*x23 *x41 5 +coeff( 59)*x12*x24 *x41 6 +coeff( 60)*x13*x21*x32*x41 7 +coeff( 61)*x12*x22*x31*x42 8 +coeff( 62)*x11*x23*x31*x42 lphi =lphi 9 +coeff( 63) *x24*x31*x42 1 +coeff( 64)*x12 *x33*x42 2 +coeff( 65)*x13*x21 *x43 3 +coeff( 66)*x14*x21*x33 4 +coeff( 67)*x14*x23 *x41 5 +coeff( 68)*x14*x21*x32*x41 6 +coeff( 69) *x33 7 +coeff( 70)*x13 *x31 8 +coeff( 71)*x12*x21*x31 lphi =lphi 9 +coeff( 72)*x12*x21 *x41 1 +coeff( 73) *x21*x32*x41 2 +coeff( 74)*x11*x23*x31 3 +coeff( 75) *x24*x31 c return end real function ly00 (x,m) dimension x(m) dimension xmin(10),xmax(10),scale(10) dimension coeff( 76) data ncoeff/ 75/ data avdat/ 0.0000000E+00/ data xmin/ 1 -0.69010E+00,-0.28144E-01,-0.49814E-01,-0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data xmax/ 1 0.54634E+00, 0.24553E-01, 0.49814E-01, 0.45070E-01, 0.00000E+00, 2 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/ data scale /10*0./ data coeff/ + -0.52208714E-01, 0.39929950E-01,-0.22094262E-01,-0.29987074E-01, + -0.31714051E-01,-0.16258660E-01, 0.13512740E-02, 0.27650652E-03, + 0.98976800E-02, 0.30010233E-02,-0.11765970E-01,-0.14765980E-01, + 0.17623210E-01, 0.39708290E-02, 0.19894510E-02, 0.12103980E-02, + -0.87466993E-03,-0.46694830E-02, 0.33423493E-02,-0.14643111E-02, + 0.11619210E-02, 0.83493990E-03, 0.21909354E-02, 0.58576660E-02, + 0.26590302E-02,-0.24663520E-03, 0.25359560E-02,-0.68681060E-03, + 0.14518880E-02, 0.38831902E-03,-0.12668310E-02, 0.14313760E-02, + 0.52571980E-02, 0.42008124E-02,-0.38940230E-02, 0.19341380E-02, + -0.10813130E-02,-0.21085320E-02, 0.33043194E-02, 0.48266010E-02, + 0.53931500E-03, 0.20872962E-02, 0.95260280E-03, 0.21073660E-02, + -0.19217660E-02, 0.14807960E-02,-0.12154962E-02,-0.14662470E-02, + -0.22175804E-03, 0.15139044E-02, 0.25094180E-02,-0.24052790E-02, + 0.62094214E-02,-0.19165310E-02, 0.19272270E-03, 0.64600270E-02, + 0.54539710E-03,-0.17769890E-03,-0.38263690E-02, 0.32160030E-02, + -0.23764302E-02, 0.10331730E-02, 0.27449650E-02,-0.26372913E-02, + 0.24685973E-02, 0.13625590E-02,-0.94646681E-03, 0.14996681E-02, + 0.87773670E-03, 0.84332610E-03, 0.11091290E-02, 0.13481281E-02, + 0.96075114E-03, 0.35147203E-02, 0.36337173E-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) c set up monomials functions x11 = x1 x12 = x11*x1 x13 = x12*x1 x14 = x13*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 c c function c ly00 =avdat 1 +coeff( 1) *x31 2 +coeff( 2) *x41 3 +coeff( 3)*x11 *x31 4 +coeff( 4) *x21*x31 5 +coeff( 5)*x11 *x41 6 +coeff( 6) *x21 *x41 7 +coeff( 7)*x12 *x31 8 +coeff( 8)*x11*x21*x31 ly00 =ly00 9 +coeff( 9) *x22*x31 1 +coeff( 10) *x33 2 +coeff( 11)*x12 *x41 3 +coeff( 12)*x11*x21 *x41 4 +coeff( 13) *x22 *x41 5 +coeff( 14) *x32*x41 6 +coeff( 15) *x31*x42 7 +coeff( 16) *x43 8 +coeff( 17)*x13 *x31 ly00 =ly00 9 +coeff( 18)*x11*x22*x31 1 +coeff( 19) *x23*x31 2 +coeff( 20)*x11 *x33 3 +coeff( 21)*x13 *x41 4 +coeff( 22)*x12*x21 *x41 5 +coeff( 23)*x11*x22 *x41 6 +coeff( 24) *x23 *x41 7 +coeff( 25)*x11 *x32*x41 8 +coeff( 26) *x21*x32*x41 ly00 =ly00 9 +coeff( 27)*x11 *x31*x42 1 +coeff( 28) *x21*x31*x42 2 +coeff( 29)*x11 *x43 3 +coeff( 30) *x21 *x43 4 +coeff( 31)*x14 *x31 5 +coeff( 32)*x13*x21*x31 6 +coeff( 33) *x22*x33 7 +coeff( 34)*x14 *x41 8 +coeff( 35)*x13*x21 *x41 ly00 =ly00 9 +coeff( 36)*x11*x23 *x41 1 +coeff( 37) *x24 *x41 2 +coeff( 38)*x12 *x32*x41 3 +coeff( 39)*x11*x21*x32*x41 4 +coeff( 40) *x22*x32*x41 5 +coeff( 41) *x34*x41 6 +coeff( 42)*x11*x21*x31*x42 7 +coeff( 43) *x22*x31*x42 8 +coeff( 44)*x12 *x43 ly00 =ly00 9 +coeff( 45)*x11*x21 *x43 1 +coeff( 46) *x22 *x43 2 +coeff( 47)*x14*x21*x31 3 +coeff( 48)*x12*x21*x33 4 +coeff( 49)*x11*x22*x33 5 +coeff( 50) *x23*x33 6 +coeff( 51)*x14*x21 *x41 7 +coeff( 52)*x13*x22 *x41 8 +coeff( 53)*x11*x22*x32*x41 ly00 =ly00 9 +coeff( 54)*x13 *x31*x42 1 +coeff( 55)*x12*x21*x31*x42 2 +coeff( 56)*x11*x22*x31*x42 3 +coeff( 57)*x13 *x43 4 +coeff( 58) *x24*x33 5 +coeff( 59)*x14*x22 *x41 6 +coeff( 60)*x13*x23 *x41 7 +coeff( 61)*x14 *x31*x42 8 +coeff( 62)*x14*x21*x33 ly00 =ly00 9 +coeff( 63)*x14*x21 *x43 1 +coeff( 64) *x23*x31*x44 2 +coeff( 65)*x12*x21*x31 3 +coeff( 66) *x21*x33 4 +coeff( 67)*x11*x23*x31 5 +coeff( 68) *x24*x31 6 +coeff( 69)*x12*x22 *x41 7 +coeff( 70)*x13*x22*x31 8 +coeff( 71)*x11*x24 *x41 ly00 =ly00 9 +coeff( 72)*x13 *x32*x41 1 +coeff( 73)*x14*x22*x31 2 +coeff( 74)*x13*x21*x31*x42 3 +coeff( 75) *x22*x33*x42 c return end