31 REAL aa,s,tt,u,ustar,w,y
35 REAL a(32),d(31),h(31),t(31)
49 DATA a/0.0,.3917609e-1,.7841241e-1,.1177699,.1573107,.1970991,
50 + .2372021,.2776904,.3186394,.3601299,.4022501,.4450965,
51 + .4887764,.5334097,.5791322,.6260990,.6744898,.7245144,
52 + .7764218,.8305109,.8871466,.9467818,1.009990,1.077516,
53 + 1.150349,1.229859,1.318011,1.417797,1.534121,1.675940,
55 DATA d/5*0.0,.2636843,.2425085,.2255674,.2116342,.1999243,
56 + .1899108,.1812252,.1736014,.1668419,.1607967,.1553497,
57 + .1504094,.1459026,.1417700,.1379632,.1344418,.1311722,
58 + .1281260,.1252791,.1226109,.1201036,.1177417,.1155119,
59 + .1134023,.1114027,.1095039/
60 DATA t/.7673828e-3,.2306870e-2,.3860618e-2,.5438454e-2,
61 + .7050699e-2,.8708396e-2,.1042357e-1,.1220953e-1,.1408125e-1,
62 + .1605579e-1,.1815290e-1,.2039573e-1,.2281177e-1,.2543407e-1,
63 + .2830296e-1,.3146822e-1,.3499233e-1,.3895483e-1,.4345878e-1,
64 + .4864035e-1,.5468334e-1,.6184222e-1,.7047983e-1,.8113195e-1,
65 + .9462444e-1,.1123001,.1364980,.1716886,.2276241,.3304980,
67 DATA h/.3920617e-1,.3932705e-1,.3950999e-1,.3975703e-1,
68 + .4007093e-1,.4045533e-1,.4091481e-1,.4145507e-1,.4208311e-1,
69 + .4280748e-1,.4363863e-1,.4458932e-1,.4567523e-1,.4691571e-1,
70 + .4833487e-1,.4996298e-1,.5183859e-1,.5401138e-1,.5654656e-1,
71 + .5953130e-1,.6308489e-1,.6737503e-1,.7264544e-1,.7926471e-1,
72 + .8781922e-1,.9930398e-1,.1155599,.1404344,.1836142,.2790016,
88 30 ustar = u - float(i)
90 40
IF (ustar.LE.t(i)) go to 60
97 IF (s.EQ.1.0) snorm = -y
109 80
IF (ustar.GT.tt) go to 50
111 IF (ustar.GE.u) go to 70
124 IF (u.LT.1.0) go to 110
132 IF (ustar.GT.tt) go to 50
134 IF (ustar.GE.u) go to 150