85     real t[3],zk;  in wl90c_energy()  local
94     t[1] = rhob+rhoa;  in wl90c_energy()
95     t[2] = 1/pow(t[1],.3333333333333333);  in wl90c_energy()
96 …t[2]*sqrt(pow(gradb,2.0)+2.0*gradab+pow(grada,2.0))-0.7486*t[1])*sqrt(1.0-1.0*pow(rhoa-1.0*rhob,2.…  in wl90c_energy()
107     real t[20];  in wl90c_first()  local
117     t[1] = rhoa-1.0*rhob;  in wl90c_first()
118     t[2] = pow(t[1],2.0);  in wl90c_first()
119     t[3] = rhob+rhoa;  in wl90c_first()
120     t[4] = 1/pow(t[3],2.0);  in wl90c_first()
121     t[5] = sqrt(1.0-1.0*t[2]*t[4]);  in wl90c_first()
122     t[6] = 1/pow(t[3],1.333333333333333);  in wl90c_first()
123     t[7] = -.2067834969664666*t[6];  in wl90c_first()
124     t[8] = 1/pow(rhoa,1.333333333333333);  in wl90c_first()
125     t[9] = 1/pow(rhob,1.333333333333333);  in wl90c_first()
126     t[10] = 1/pow(t[3],.3333333333333333);  in wl90c_first()
127     t[11] = 0.6203504908994*t[10]+0.9*(gradb*t[9]+grada*t[8])+3.60073;  in wl90c_first()
128     t[12] = 1/pow(t[11],2.0);  in wl90c_first()
129     t[13] = sqrt(pow(gradb,2.0)+2.0*gradab+pow(grada,2.0));  in wl90c_first()
130     t[14] = 0.06001*t[13]*t[10]-0.7486*t[3];  in wl90c_first()
131     t[15] = 1/t[5];  in wl90c_first()
132     t[16] = 2.0*t[2]/pow(t[3],3.0);  in wl90c_first()
133     t[17] = 1/t[11];  in wl90c_first()
134     t[18] = t[5]*(-.02000333333333333*t[13]*t[6]-0.7486)*t[17];  in wl90c_first()
135     t[19] = 1/t[13];  in wl90c_first()
136 …dfdra = -1.0*t[12]*t[14]*t[5]*(t[7]-1.2*grada/pow(rhoa,2.333333333333333))+0.5*t[14]*t[15]*t[17]*(…  in wl90c_first()
137 …dfdrb = -1.0*t[12]*t[14]*t[5]*(t[7]-1.2*gradb/pow(rhob,2.333333333333333))+0.5*t[14]*t[15]*t[17]*(…  in wl90c_first()
138     dfdga = 0.06001*grada*t[19]*t[10]*t[5]*t[17]-0.9*t[8]*t[5]*t[12]*t[14];  in wl90c_first()
139     dfdgb = 0.06001*gradb*t[19]*t[10]*t[5]*t[17]-0.9*t[9]*t[5]*t[12]*t[14];  in wl90c_first()
140     dfdab = 0.06001*t[19]*t[10]*t[5]*t[17];  in wl90c_first()
152     real t[44];  in wl90c_second()  local
166     t[1] = rhoa-1.0*rhob;  in wl90c_second()
167     t[2] = pow(t[1],2.0);  in wl90c_second()
168     t[3] = rhob+rhoa;  in wl90c_second()
169     t[4] = 1/pow(t[3],2.0);  in wl90c_second()
170     t[5] = sqrt(1.0-1.0*t[2]*t[4]);  in wl90c_second()
171     t[6] = 1/pow(rhoa,2.333333333333333);  in wl90c_second()
172     t[7] = 1/pow(t[3],1.333333333333333);  in wl90c_second()
173     t[8] = -.2067834969664666*t[7];  in wl90c_second()
174     t[9] = t[8]-1.2*grada*t[6];  in wl90c_second()
175     t[10] = 1/pow(rhoa,1.333333333333333);  in wl90c_second()
176     t[11] = 1/pow(rhob,1.333333333333333);  in wl90c_second()
177     t[12] = 1/pow(t[3],.3333333333333333);  in wl90c_second()
178     t[13] = 0.6203504908994*t[12]+0.9*(gradb*t[11]+grada*t[10])+3.60073;  in wl90c_second()
179     t[14] = 1/pow(t[13],2.0);  in wl90c_second()
180     t[15] = pow(grada,2.0);  in wl90c_second()
181     t[16] = pow(gradb,2.0);  in wl90c_second()
182     t[17] = sqrt(2.0*gradab+t[16]+t[15]);  in wl90c_second()
183     t[18] = 0.06001*t[17]*t[12]-0.7486*t[3];  in wl90c_second()
184     t[19] = 1/t[5];  in wl90c_second()
185     t[20] = 1/pow(t[3],3.0);  in wl90c_second()
186     t[21] = 2.0*t[2]*t[20];  in wl90c_second()
187     t[22] = t[21]-2.0*t[1]*t[4];  in wl90c_second()
188     t[23] = 1/t[13];  in wl90c_second()
189     t[24] = -.02000333333333333*t[17]*t[7]-0.7486;  in wl90c_second()
190     t[25] = t[5]*t[24]*t[23];  in wl90c_second()
191     t[26] = 1/pow(rhob,2.333333333333333);  in wl90c_second()
192     t[27] = t[8]-1.2*gradb*t[26];  in wl90c_second()
193     t[28] = 2.0*t[1]*t[4]+t[21];  in wl90c_second()
194     t[29] = 1/t[17];  in wl90c_second()
195     t[30] = 0.06001*t[29]*t[12]*t[5]*t[23];  in wl90c_second()
196     t[31] = 1/pow(t[13],3.0);  in wl90c_second()
197     t[32] = 1/pow(t[3],2.333333333333333);  in wl90c_second()
198     t[33] = .2757113292886222*t[32];  in wl90c_second()
199     t[34] = -6.0*t[2]/pow(t[3],4.0);  in wl90c_second()
200     t[35] = -2.0*t[4];  in wl90c_second()
201     t[36] = 1/pow(t[5],3.0);  in wl90c_second()
202     t[37] = .02667111111111111*t[17]*t[32]*t[5]*t[23];  in wl90c_second()
203     t[38] = -0.9*t[10]*t[5]*t[24]*t[14];  in wl90c_second()
204     t[39] = -.02000333333333333*grada*t[29]*t[7]*t[5]*t[23];  in wl90c_second()
205     t[40] = -0.9*t[11]*t[5]*t[24]*t[14];  in wl90c_second()
206     t[41] = -.02000333333333333*gradb*t[29]*t[7]*t[5]*t[23];  in wl90c_second()
207     t[42] = -.02000333333333333*t[29]*t[7]*t[5]*t[23];  in wl90c_second()
208     t[43] = 1/pow(t[17],3.0);  in wl90c_second()
209     dfdra =-1.0*t[14]*t[18]*t[5]*t[9]+t[25]+0.5*t[18]*t[19]*t[22]*t[23];  in wl90c_second()
210     dfdrb = -1.0*t[14]*t[18]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[28]+t[25];  in wl90c_second()
211     dfdga = 0.06001*grada*t[29]*t[12]*t[5]*t[23]-0.9*t[10]*t[5]*t[14]*t[18];  in wl90c_second()
212     dfdgb = 0.06001*gradb*t[29]*t[12]*t[5]*t[23]-0.9*t[11]*t[5]*t[14]*t[18];  in wl90c_second()
213     dfdab = t[30];  in wl90c_second()
214 …t[14]*t[18]*t[5]*(2.8*grada/pow(rhoa,3.333333333333333)+t[33])+2.0*t[18]*t[31]*t[5]*pow(t[9],2.0)-…  in wl90c_second()
215 …t[18]*t[27]*t[31]*t[5]*t[9]-1.0*t[14]*t[24]*t[5]*t[9]-0.5*t[14]*t[18]*t[19]*t[28]*t[9]-1.0*t[14]*t…  in wl90c_second()
216 …t[39]+0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[23]+t[38]-0.06001*grada*t[29]*t[12]*t[5]*t[9]*t[14]…  in wl90c_second()
217 …t[41]+0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[23]+t[40]-0.06001*gradb*t[29]*t[12]*t[5]*t[9]*t[14]…  in wl90c_second()
218     d2fdraab = t[42]+0.030005*t[29]*t[12]*t[19]*t[22]*t[23]-0.06001*t[29]*t[12]*t[5]*t[9]*t[14];  in wl90c_second()
219 …t[14]*t[18]*t[5]*(2.8*gradb/pow(rhob,3.333333333333333)+t[33])+2.0*t[18]*pow(t[27],2.0)*t[31]*t[5]…  in wl90c_second()
220 …t[39]+0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[23]+t[38]-0.06001*grada*t[29]*t[12]*t[5]*t[27]*t[14…  in wl90c_second()
221 …t[41]+0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[23]+t[40]-0.06001*gradb*t[29]*t[12]*t[5]*t[27]*t[14…  in wl90c_second()
222     d2fdrbab = t[42]+0.030005*t[29]*t[12]*t[19]*t[28]*t[23]-0.06001*t[29]*t[12]*t[5]*t[27]*t[14];  in wl90c_second()
223 …dgaga = 1.62*t[18]*t[31]*t[5]/pow(rhoa,2.666666666666667)+t[30]-0.06001*t[15]*t[43]*t[12]*t[5]*t[2…  in wl90c_second()
224 …gradb*t[43]*t[12]*t[5]*t[23]-0.054009*grada*t[29]*t[11]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[10…  in wl90c_second()
225     d2fdgaab = -0.06001*grada*t[43]*t[12]*t[5]*t[23]-0.054009*t[29]*t[10]*t[12]*t[5]*t[14];  in wl90c_second()
226 …dgbgb = 1.62*t[18]*t[31]*t[5]/pow(rhob,2.666666666666667)+t[30]-0.06001*t[16]*t[43]*t[12]*t[5]*t[2…  in wl90c_second()
227     d2fdgbab = -0.06001*gradb*t[43]*t[12]*t[5]*t[23]-0.054009*t[29]*t[11]*t[12]*t[5]*t[14];  in wl90c_second()
228     d2fdabab = -0.06001*t[43]*t[12]*t[5]*t[23];  in wl90c_second()
255     real t[104];  in wl90c_third()  local
278     t[1] = rhoa-1.0*rhob;  in wl90c_third()
279     t[2] = pow(t[1],2.0);  in wl90c_third()
280     t[3] = rhob+rhoa;  in wl90c_third()
281     t[4] = 1/pow(t[3],2.0);  in wl90c_third()
282     t[5] = sqrt(1.0-1.0*t[2]*t[4]);  in wl90c_third()
283     t[6] = 1/pow(rhoa,2.333333333333333);  in wl90c_third()
284     t[7] = 1/pow(t[3],1.333333333333333);  in wl90c_third()
285     t[8] = -.2067834969664666*t[7];  in wl90c_third()
286     t[9] = t[8]-1.2*grada*t[6];  in wl90c_third()
287     t[10] = 1/pow(rhoa,1.333333333333333);  in wl90c_third()
288     t[11] = 1/pow(rhob,1.333333333333333);  in wl90c_third()
289     t[12] = 1/pow(t[3],.3333333333333333);  in wl90c_third()
290     t[13] = 0.6203504908994*t[12]+0.9*(gradb*t[11]+grada*t[10])+3.60073;  in wl90c_third()
291     t[14] = 1/pow(t[13],2.0);  in wl90c_third()
292     t[15] = pow(grada,2.0);  in wl90c_third()
293     t[16] = pow(gradb,2.0);  in wl90c_third()
294     t[17] = sqrt(2.0*gradab+t[16]+t[15]);  in wl90c_third()
295     t[18] = 0.06001*t[17]*t[12]-0.7486*t[3];  in wl90c_third()
296     t[19] = 1/t[5];  in wl90c_third()
297     t[20] = 1/pow(t[3],3.0);  in wl90c_third()
298     t[21] = 2.0*t[2]*t[20];  in wl90c_third()
299     t[22] = t[21]-2.0*t[1]*t[4];  in wl90c_third()
300     t[23] = 1/t[13];  in wl90c_third()
301     t[24] = -.02000333333333333*t[17]*t[7]-0.7486;  in wl90c_third()
302     t[25] = t[5]*t[24]*t[23];  in wl90c_third()
303     t[26] = 1/pow(rhob,2.333333333333333);  in wl90c_third()
304     t[27] = t[8]-1.2*gradb*t[26];  in wl90c_third()
305     t[28] = 2.0*t[1]*t[4]+t[21];  in wl90c_third()
306     t[29] = 1/t[17];  in wl90c_third()
307     t[30] = 0.06001*t[29]*t[12]*t[5]*t[23];  in wl90c_third()
308     t[31] = pow(t[9],2.0);  in wl90c_third()
309     t[32] = 1/pow(t[13],3.0);  in wl90c_third()
310     t[33] = 1/pow(rhoa,3.333333333333333);  in wl90c_third()
311     t[34] = 1/pow(t[3],2.333333333333333);  in wl90c_third()
312     t[35] = .2757113292886222*t[34];  in wl90c_third()
313     t[36] = t[35]+2.8*grada*t[33];  in wl90c_third()
314     t[37] = 1/pow(t[3],4.0);  in wl90c_third()
315     t[38] = -6.0*t[2]*t[37];  in wl90c_third()
316     t[39] = -2.0*t[4];  in wl90c_third()
317     t[40] = t[39]+t[38]+8.0*t[1]*t[20];  in wl90c_third()
318     t[41] = 1/pow(t[5],3.0);  in wl90c_third()
319     t[42] = pow(t[22],2.0);  in wl90c_third()
320     t[43] = .02667111111111111*t[17]*t[34]*t[5]*t[23];  in wl90c_third()
321     t[44] = 2.0*t[4]+t[38];  in wl90c_third()
322     t[45] = -0.9*t[10]*t[5]*t[24]*t[14];  in wl90c_third()
323     t[46] = -.02000333333333333*grada*t[29]*t[7]*t[5]*t[23];  in wl90c_third()
324     t[47] = -0.9*t[11]*t[5]*t[24]*t[14];  in wl90c_third()
325     t[48] = -.02000333333333333*gradb*t[29]*t[7]*t[5]*t[23];  in wl90c_third()
326     t[49] = -0.06001*t[29]*t[12]*t[5]*t[9]*t[14];  in wl90c_third()
327     t[50] = 0.030005*t[29]*t[12]*t[19]*t[22]*t[23];  in wl90c_third()
328     t[51] = -.02000333333333333*t[29]*t[7]*t[5]*t[23];  in wl90c_third()
329     t[52] = pow(t[27],2.0);  in wl90c_third()
330     t[53] = 1/pow(rhob,3.333333333333333);  in wl90c_third()
331     t[54] = t[35]+2.8*gradb*t[53];  in wl90c_third()
332     t[55] = t[39]+t[38]-8.0*t[1]*t[20];  in wl90c_third()
333     t[56] = pow(t[28],2.0);  in wl90c_third()
334     t[57] = -0.06001*t[29]*t[12]*t[5]*t[27]*t[14];  in wl90c_third()
335     t[58] = 0.030005*t[29]*t[12]*t[19]*t[28]*t[23];  in wl90c_third()
336     t[59] = 1/pow(rhoa,2.666666666666667);  in wl90c_third()
337     t[60] = 1/pow(t[17],3.0);  in wl90c_third()
338     t[61] = -0.054009*t[29]*t[10]*t[12]*t[5]*t[14];  in wl90c_third()
339     t[62] = -0.06001*grada*t[60]*t[12]*t[5]*t[23];  in wl90c_third()
340     t[63] = 1/pow(rhob,2.666666666666667);  in wl90c_third()
341     t[64] = -0.054009*t[29]*t[11]*t[12]*t[5]*t[14];  in wl90c_third()
342     t[65] = -0.06001*gradb*t[60]*t[12]*t[5]*t[23];  in wl90c_third()
343     t[66] = -0.06001*t[60]*t[12]*t[5]*t[23];  in wl90c_third()
344     t[67] = 1/pow(t[13],4.0);  in wl90c_third()
345     t[68] = 1/pow(t[3],3.333333333333333);  in wl90c_third()
346     t[69] = -.6433264350067851*t[68];  in wl90c_third()
347     t[70] = 24.0*t[2]/pow(t[3],5.0);  in wl90c_third()
348     t[71] = 12.0*t[20];  in wl90c_third()
349     t[72] = 1/pow(t[5],5.0);  in wl90c_third()
350     t[73] = -.06223259259259259*t[17]*t[68]*t[5]*t[23];  in wl90c_third()
351     t[74] = .6433264350067851*t[68]*t[5]*t[14]*t[18];  in wl90c_third()
352     t[75] = -4.0*t[20];  in wl90c_third()
353     t[76] = 4.0*t[24]*t[27]*t[32]*t[5]*t[9];  in wl90c_third()
354     t[77] = -1.0*t[14]*t[19]*t[24]*t[28]*t[9];  in wl90c_third()
355     t[78] = -1.0*t[14]*t[19]*t[22]*t[24]*t[27];  in wl90c_third()
356     t[79] = -.5514226585772444*t[34]*t[5]*t[24]*t[14];  in wl90c_third()
357     t[80] = -0.5*t[22]*t[23]*t[24]*t[28]*t[41];  in wl90c_third()
358     t[81] = t[44]*t[19]*t[24]*t[23];  in wl90c_third()
359     t[82] = -0.024004*t[17]*t[10]*t[34]*t[5]*t[14];  in wl90c_third()
360     t[83] = .02667111111111111*grada*t[29]*t[34]*t[5]*t[23];  in wl90c_third()
361     t[84] = -0.024004*t[17]*t[11]*t[34]*t[5]*t[14];  in wl90c_third()
362     t[85] = .02667111111111111*gradb*t[29]*t[34]*t[5]*t[23];  in wl90c_third()
363     t[86] = .02667111111111111*t[29]*t[34]*t[5]*t[23];  in wl90c_third()
364     t[87] = 1/pow(t[3],2.666666666666667);  in wl90c_third()
365     t[88] = 1.62*t[59]*t[5]*t[24]*t[32];  in wl90c_third()
366     t[89] = 0.036006*grada*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_third()
367     t[90] = .02000333333333333*t[15]*t[60]*t[7]*t[5]*t[23];  in wl90c_third()
368     t[91] = 1.62*t[10]*t[11]*t[5]*t[24]*t[32];  in wl90c_third()
369     t[92] = 0.018003*gradb*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_third()
370     t[93] = 0.018003*grada*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_third()
371     t[94] = .02000333333333333*grada*gradb*t[60]*t[7]*t[5]*t[23];  in wl90c_third()
372     t[95] = 0.018003*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_third()
373     t[96] = .02000333333333333*grada*t[60]*t[7]*t[5]*t[23];  in wl90c_third()
374     t[97] = 1.62*t[63]*t[5]*t[24]*t[32];  in wl90c_third()
375     t[98] = 0.036006*gradb*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_third()
376     t[99] = .02000333333333333*t[16]*t[60]*t[7]*t[5]*t[23];  in wl90c_third()
377     t[100] = 0.018003*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_third()
378     t[101] = .02000333333333333*gradb*t[60]*t[7]*t[5]*t[23];  in wl90c_third()
379     t[102] = .02000333333333333*t[60]*t[7]*t[5]*t[23];  in wl90c_third()
380     t[103] = 1/pow(t[17],5.0);  in wl90c_third()
381     dfdra = -1.0*t[14]*t[18]*t[5]*t[9]+t[25]+0.5*t[18]*t[19]*t[22]*t[23];  in wl90c_third()
382     dfdrb = -1.0*t[14]*t[18]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[28]+t[25];  in wl90c_third()
383     dfdga = 0.06001*grada*t[29]*t[12]*t[5]*t[23]-0.9*t[10]*t[5]*t[14]*t[18];  in wl90c_third()
384     dfdgb = 0.06001*gradb*t[29]*t[12]*t[5]*t[23]-0.9*t[11]*t[5]*t[14]*t[18];  in wl90c_third()
385     dfdab = t[30];  in wl90c_third()
386 …t[14]*t[24]*t[5]*t[9]-1.0*t[14]*t[18]*t[19]*t[22]*t[9]-1.0*t[14]*t[18]*t[36]*t[5]+2.0*t[18]*t[31]*…  in wl90c_third()
387 …t[18]*t[27]*t[32]*t[5]*t[9]-1.0*t[14]*t[24]*t[5]*t[9]-0.5*t[14]*t[18]*t[19]*t[28]*t[9]-1.0*t[14]*t…  in wl90c_third()
388 …t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[9]*t[14]…  in wl90c_third()
389 …t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[9]*t[14]…  in wl90c_third()
390     d2fdraab = t[51]+t[50]+t[49];  in wl90c_third()
391 …t[18]*t[23]*t[41]*t[56]+0.5*t[18]*t[19]*t[23]*t[55]-1.0*t[14]*t[18]*t[5]*t[54]+2.0*t[18]*t[32]*t[5…  in wl90c_third()
392 …t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[27]*t[14…  in wl90c_third()
393 …t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[27]*t[14…  in wl90c_third()
394     d2fdrbab = t[51]+t[58]+t[57];  in wl90c_third()
395 …d2fdgaga = t[30]-0.06001*t[15]*t[60]*t[12]*t[5]*t[23]-0.108018*grada*t[29]*t[10]*t[12]*t[5]*t[14]+…  in wl90c_third()
396 …gradb*t[60]*t[12]*t[5]*t[23]-0.054009*grada*t[29]*t[11]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[10…  in wl90c_third()
397     d2fdgaab = t[62]+t[61];  in wl90c_third()
398 …d2fdgbgb = t[30]-0.06001*t[16]*t[60]*t[12]*t[5]*t[23]-0.108018*gradb*t[29]*t[11]*t[12]*t[5]*t[14]+…  in wl90c_third()
399     d2fdgbab = t[65]+t[64];  in wl90c_third()
400     d2fdabab = t[66];  in wl90c_third()
401 …t[14]*t[18]*t[5]*(t[69]-9.333333333333334*grada/pow(rhoa,4.333333333333333))-6.0*t[18]*t[5]*t[67]*…  in wl90c_third()
402 …t[14]*t[18]*t[19]*t[44]*t[9]+0.5*t[14]*t[18]*t[22]*t[28]*t[41]*t[9]+2.0*t[18]*t[19]*t[22]*t[27]*t[…  in wl90c_third()
403 …t[29]*t[12]*t[41]*t[42]*t[23]+t[83]+0.030005*grada*t[29]*t[12]*t[40]*t[19]*t[23]-.0200033333333333…  in wl90c_third()
404 …t[29]*t[12]*t[41]*t[42]*t[23]+t[85]+0.030005*gradb*t[29]*t[12]*t[40]*t[19]*t[23]-.0200033333333333…  in wl90c_third()
405 …t[29]*t[12]*t[41]*t[42]*t[23]+t[86]+0.030005*t[29]*t[12]*t[40]*t[19]*t[23]-.02000333333333333*t[29…  in wl90c_third()
406 …t[18]*t[5]*t[52]*t[67]*t[9]+0.25*t[14]*t[18]*t[41]*t[56]*t[9]-0.5*t[14]*t[18]*t[19]*t[55]*t[9]+2.0…  in wl90c_third()
407 …t[83]+0.030005*grada*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*grada*t[29]*t[7]*t[19]*t[28]…  in wl90c_third()
408 …t[85]+0.030005*gradb*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*gradb*t[29]*t[7]*t[19]*t[28]…  in wl90c_third()
409 …t[86]+0.030005*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*t[29]*t[7]*t[19]*t[28]*t[23]-0.015…  in wl90c_third()
410 …t[18]*t[32]*t[5]/pow(rhoa,3.666666666666667)+t[90]+t[89]+t[88]+t[51]+t[50]+t[49]+0.216036*grada*t[…  in wl90c_third()
411 …t[94]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+.07201199999999999*gradb*t[29]*t[6]*t[12]…  in wl90c_third()
412 …t[96]-0.030005*grada*t[60]*t[12]*t[19]*t[22]*t[23]+.07201199999999999*t[29]*t[6]*t[12]*t[5]*t[14]+…  in wl90c_third()
413 …t[51]+t[99]+t[50]-0.030005*t[16]*t[60]*t[12]*t[19]*t[22]*t[23]+t[98]+t[49]+0.06001*t[16]*t[60]*t[1…  in wl90c_third()
414 …t[101]-0.030005*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+t[100]+0.06001*gradb*t[60]*t[12]*t[5]*t[9]*t[1…  in wl90c_third()
415     d3fdraabab = t[102]-0.030005*t[60]*t[12]*t[19]*t[22]*t[23]+0.06001*t[60]*t[12]*t[5]*t[9]*t[14];  in wl90c_third()
416 …t[14]*t[18]*t[5]*(t[69]-9.333333333333334*gradb/pow(rhob,4.333333333333333))+t[73]+0.375*t[18]*t[2…  in wl90c_third()
417 …t[29]*t[12]*t[41]*t[56]*t[23]+t[83]+0.030005*grada*t[29]*t[12]*t[55]*t[19]*t[23]-.0200033333333333…  in wl90c_third()
418 …t[29]*t[12]*t[41]*t[56]*t[23]+t[85]+0.030005*gradb*t[29]*t[12]*t[55]*t[19]*t[23]-.0200033333333333…  in wl90c_third()
419 …t[29]*t[12]*t[41]*t[56]*t[23]+t[86]+0.030005*t[29]*t[12]*t[55]*t[19]*t[23]-.02000333333333333*t[29…  in wl90c_third()
420 …t[51]+t[90]+t[58]-0.030005*t[15]*t[60]*t[12]*t[19]*t[28]*t[23]+t[89]+t[57]+0.06001*t[15]*t[60]*t[1…  in wl90c_third()
421 …t[94]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+.07201199999999999*grada*t[29]*t[26]*t[12…  in wl90c_third()
422 …t[96]-0.030005*grada*t[60]*t[12]*t[19]*t[28]*t[23]+t[95]+0.06001*grada*t[60]*t[12]*t[5]*t[27]*t[14…  in wl90c_third()
423 …t[18]*t[32]*t[5]/pow(rhob,3.666666666666667)+t[99]+t[98]+t[97]+t[58]+t[57]+t[51]+0.216036*gradb*t[…  in wl90c_third()
424 …t[101]-0.030005*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+.07201199999999999*t[29]*t[26]*t[12]*t[5]*t[14…  in wl90c_third()
425     d3fdrbabab = t[102]-0.030005*t[60]*t[12]*t[19]*t[28]*t[23]+0.06001*t[60]*t[12]*t[5]*t[27]*t[14];  in wl90c_third()
426 …t[18]*t[5]*t[67]/pow(rhoa,4.0)+0.18003*t[103]*t[12]*t[23]*t[5]*pow(grada,3.0)+0.2916486*grada*t[29…  in wl90c_third()
427 …t[65]+0.18003*t[15]*gradb*t[103]*t[12]*t[5]*t[23]+t[64]+0.054009*t[15]*t[60]*t[11]*t[12]*t[5]*t[14…  in wl90c_third()
428 …d3fdgagaab = t[66]+0.18003*t[15]*t[103]*t[12]*t[5]*t[23]+0.108018*grada*t[60]*t[10]*t[12]*t[5]*t[1…  in wl90c_third()
429 …t[62]+0.18003*grada*t[16]*t[103]*t[12]*t[5]*t[23]+0.108018*grada*gradb*t[60]*t[11]*t[12]*t[5]*t[14…  in wl90c_third()
430 …t[103]*t[12]*t[5]*t[23]+0.054009*grada*t[60]*t[11]*t[12]*t[5]*t[14]+0.054009*gradb*t[60]*t[10]*t[1…  in wl90c_third()
431     d3fdgaabab = 0.18003*grada*t[103]*t[12]*t[5]*t[23]+0.054009*t[60]*t[10]*t[12]*t[5]*t[14];  in wl90c_third()
432 …t[18]*t[5]*t[67]/pow(rhob,4.0)+0.18003*t[103]*t[12]*t[23]*t[5]*pow(gradb,3.0)+0.2916486*gradb*t[29…  in wl90c_third()
433 …d3fdgbgbab= t[66]+0.18003*t[16]*t[103]*t[12]*t[5]*t[23]+0.108018*gradb*t[60]*t[11]*t[12]*t[5]*t[14…  in wl90c_third()
434     d3fdgbabab = 0.18003*gradb*t[103]*t[12]*t[5]*t[23]+0.054009*t[60]*t[11]*t[12]*t[5]*t[14];  in wl90c_third()
435     d3fdababab = 0.18003*t[103]*t[12]*t[5]*t[23];  in wl90c_third()
497     real t[300];  in wl90c_fourth()  local
538     t[1] = rhoa-1.0*rhob;  in wl90c_fourth()
539     t[2] = pow(t[1],2.0);  in wl90c_fourth()
540     t[3] = rhob+rhoa;  in wl90c_fourth()
541     t[4] = 1/pow(t[3],2.0);  in wl90c_fourth()
542     t[5] = sqrt(1.0-1.0*t[2]*t[4]);  in wl90c_fourth()
543     t[6] = 1/pow(rhoa,2.333333333333333);  in wl90c_fourth()
544     t[7] = 1/pow(t[3],1.333333333333333);  in wl90c_fourth()
545     t[8] = -.2067834969664666*t[7];  in wl90c_fourth()
546     t[9] = t[8]-1.2*grada*t[6];  in wl90c_fourth()
547     t[10] = 1/pow(rhoa,1.333333333333333);  in wl90c_fourth()
548     t[11] = 1/pow(rhob,1.333333333333333);  in wl90c_fourth()
549     t[12] = 1/pow(t[3],.3333333333333333);  in wl90c_fourth()
550     t[13] = 0.6203504908994*t[12]+0.9*(gradb*t[11]+grada*t[10])+3.60073;  in wl90c_fourth()
551     t[14] = 1/pow(t[13],2.0);  in wl90c_fourth()
552     t[15] = pow(grada,2.0);  in wl90c_fourth()
553     t[16] = pow(gradb,2.0);  in wl90c_fourth()
554     t[17] = sqrt(2.0*gradab+t[16]+t[15]);  in wl90c_fourth()
555     t[18] = 0.06001*t[17]*t[12]-0.7486*t[3];  in wl90c_fourth()
556     t[19] = 1/t[5];  in wl90c_fourth()
557     t[20] = 1/pow(t[3],3.0);  in wl90c_fourth()
558     t[21] = 2.0*t[2]*t[20];  in wl90c_fourth()
559     t[22] = t[21]-2.0*t[1]*t[4];  in wl90c_fourth()
560     t[23] = 1/t[13];  in wl90c_fourth()
561     t[24] = -.02000333333333333*t[17]*t[7]-0.7486;  in wl90c_fourth()
562     t[25] = t[5]*t[24]*t[23];  in wl90c_fourth()
563     t[26] = 1/pow(rhob,2.333333333333333);  in wl90c_fourth()
564     t[27] = t[8]-1.2*gradb*t[26];  in wl90c_fourth()
565     t[28] = 2.0*t[1]*t[4]+t[21];  in wl90c_fourth()
566     t[29] = 1/t[17];  in wl90c_fourth()
567     t[30] = 0.06001*t[29]*t[12]*t[5]*t[23];  in wl90c_fourth()
568     t[31] = pow(t[9],2.0);  in wl90c_fourth()
569     t[32] = 1/pow(t[13],3.0);  in wl90c_fourth()
570     t[33] = 1/pow(rhoa,3.333333333333333);  in wl90c_fourth()
571     t[34] = 1/pow(t[3],2.333333333333333);  in wl90c_fourth()
572     t[35] = .2757113292886222*t[34];  in wl90c_fourth()
573     t[36] = t[35]+2.8*grada*t[33];  in wl90c_fourth()
574     t[37] = 1/pow(t[3],4.0);  in wl90c_fourth()
575     t[38] = -6.0*t[2]*t[37];  in wl90c_fourth()
576     t[39] = -2.0*t[4];  in wl90c_fourth()
577     t[40] = t[39]+t[38]+8.0*t[1]*t[20];  in wl90c_fourth()
578     t[41] = 1/pow(t[5],3.0);  in wl90c_fourth()
579     t[42] = pow(t[22],2.0);  in wl90c_fourth()
580     t[43] = .02667111111111111*t[17]*t[34]*t[5]*t[23];  in wl90c_fourth()
581     t[44] = 2.0*t[4]+t[38];  in wl90c_fourth()
582     t[45] = -0.9*t[10]*t[5]*t[24]*t[14];  in wl90c_fourth()
583     t[46] = -.02000333333333333*grada*t[29]*t[7]*t[5]*t[23];  in wl90c_fourth()
584     t[47] = -0.9*t[11]*t[5]*t[24]*t[14];  in wl90c_fourth()
585     t[48] = -.02000333333333333*gradb*t[29]*t[7]*t[5]*t[23];  in wl90c_fourth()
586     t[49] = -0.06001*t[29]*t[12]*t[5]*t[9]*t[14];  in wl90c_fourth()
587     t[50] = 0.030005*t[29]*t[12]*t[19]*t[22]*t[23];  in wl90c_fourth()
588     t[51] = -.02000333333333333*t[29]*t[7]*t[5]*t[23];  in wl90c_fourth()
589     t[52] = pow(t[27],2.0);  in wl90c_fourth()
590     t[53] = 1/pow(rhob,3.333333333333333);  in wl90c_fourth()
591     t[54] = t[35]+2.8*gradb*t[53];  in wl90c_fourth()
592     t[55] = t[39]+t[38]-8.0*t[1]*t[20];  in wl90c_fourth()
593     t[56] = pow(t[28],2.0);  in wl90c_fourth()
594     t[57] = -0.06001*t[29]*t[12]*t[5]*t[27]*t[14];  in wl90c_fourth()
595     t[58] = 0.030005*t[29]*t[12]*t[19]*t[28]*t[23];  in wl90c_fourth()
596     t[59] = 1/pow(rhoa,2.666666666666667);  in wl90c_fourth()
597     t[60] = 1/pow(t[17],3.0);  in wl90c_fourth()
598     t[61] = -0.054009*t[29]*t[10]*t[12]*t[5]*t[14];  in wl90c_fourth()
599     t[62] = -0.06001*grada*t[60]*t[12]*t[5]*t[23];  in wl90c_fourth()
600     t[63] = 1/pow(rhob,2.666666666666667);  in wl90c_fourth()
601     t[64] = -0.054009*t[29]*t[11]*t[12]*t[5]*t[14];  in wl90c_fourth()
602     t[65] = -0.06001*gradb*t[60]*t[12]*t[5]*t[23];  in wl90c_fourth()
603     t[66] = -0.06001*t[60]*t[12]*t[5]*t[23];  in wl90c_fourth()
604     t[67] = pow(t[9],3.0);  in wl90c_fourth()
605     t[68] = 1/pow(t[13],4.0);  in wl90c_fourth()
606     t[69] = 1/pow(rhoa,4.333333333333333);  in wl90c_fourth()
607     t[70] = 1/pow(t[3],3.333333333333333);  in wl90c_fourth()
608     t[71] = -.6433264350067851*t[70];  in wl90c_fourth()
609     t[72] = t[71]-9.333333333333334*grada*t[69];  in wl90c_fourth()
610     t[73] = 1/pow(t[3],5.0);  in wl90c_fourth()
611     t[74] = 24.0*t[2]*t[73];  in wl90c_fourth()
612     t[75] = 12.0*t[20];  in wl90c_fourth()
613     t[76] = t[75]+t[74]-36.0*t[1]*t[37];  in wl90c_fourth()
614     t[77] = 1/pow(t[5],5.0);  in wl90c_fourth()
615     t[78] = pow(t[22],3.0);  in wl90c_fourth()
616     t[79] = -.06223259259259259*t[17]*t[70]*t[5]*t[23];  in wl90c_fourth()
617     t[80] = .6433264350067851*t[70]*t[5]*t[14]*t[18];  in wl90c_fourth()
618     t[81] = -4.0*t[20];  in wl90c_fourth()
619     t[82] = t[81]+t[74]-12.0*t[1]*t[37];  in wl90c_fourth()
620     t[83] = 4.0*t[24]*t[27]*t[32]*t[5]*t[9];  in wl90c_fourth()
621     t[84] = -1.0*t[14]*t[19]*t[24]*t[28]*t[9];  in wl90c_fourth()
622     t[85] = -1.0*t[14]*t[19]*t[22]*t[24]*t[27];  in wl90c_fourth()
623     t[86] = -.5514226585772444*t[34]*t[5]*t[24]*t[14];  in wl90c_fourth()
624     t[87] = -0.5*t[22]*t[23]*t[24]*t[28]*t[41];  in wl90c_fourth()
625     t[88] = t[44]*t[19]*t[24]*t[23];  in wl90c_fourth()
626     t[89] = -0.024004*t[17]*t[10]*t[34]*t[5]*t[14];  in wl90c_fourth()
627     t[90] = .02667111111111111*grada*t[29]*t[34]*t[5]*t[23];  in wl90c_fourth()
628     t[91] = -0.024004*t[17]*t[11]*t[34]*t[5]*t[14];  in wl90c_fourth()
629     t[92] = .02667111111111111*gradb*t[29]*t[34]*t[5]*t[23];  in wl90c_fourth()
630     t[93] = 0.12002*t[29]*t[12]*t[5]*t[31]*t[32];  in wl90c_fourth()
631     t[94] = -0.06001*t[29]*t[12]*t[19]*t[22]*t[9]*t[14];  in wl90c_fourth()
632     t[95] = .04000666666666666*t[29]*t[7]*t[5]*t[9]*t[14];  in wl90c_fourth()
633     t[96] = -0.06001*t[29]*t[12]*t[36]*t[5]*t[14];  in wl90c_fourth()
634     t[97] = -.02000333333333333*t[29]*t[7]*t[19]*t[22]*t[23];  in wl90c_fourth()
635     t[98] = 0.030005*t[29]*t[12]*t[40]*t[19]*t[23];  in wl90c_fourth()
636     t[99] = .02667111111111111*t[29]*t[34]*t[5]*t[23];  in wl90c_fourth()
637     t[100] = -0.0150025*t[29]*t[12]*t[41]*t[42]*t[23];  in wl90c_fourth()
638     t[101] = t[81]+t[74]+12.0*t[1]*t[37];  in wl90c_fourth()
639     t[102] = 1/pow(t[3],2.666666666666667);  in wl90c_fourth()
640     t[103] = 0.12002*t[29]*t[12]*t[5]*t[9]*t[27]*t[32];  in wl90c_fourth()
641     t[104] = -0.030005*t[29]*t[12]*t[19]*t[28]*t[9]*t[14];  in wl90c_fourth()
642     t[105] = .02000333333333333*t[29]*t[7]*t[5]*t[9]*t[14];  in wl90c_fourth()
643     t[106] = -0.030005*t[29]*t[12]*t[19]*t[22]*t[27]*t[14];  in wl90c_fourth()
644     t[107] = .02000333333333333*t[29]*t[7]*t[5]*t[27]*t[14];  in wl90c_fourth()
645     t[108] = -.01654543687061022*t[29]*t[102]*t[5]*t[14];  in wl90c_fourth()
646     t[109] = -.01000166666666667*t[29]*t[7]*t[19]*t[22]*t[23];  in wl90c_fourth()
647     t[110] = -0.0150025*t[29]*t[12]*t[41]*t[22]*t[28]*t[23];  in wl90c_fourth()
648     t[111] = -.01000166666666667*t[29]*t[7]*t[19]*t[28]*t[23];  in wl90c_fourth()
649     t[112] = 0.030005*t[29]*t[12]*t[44]*t[19]*t[23];  in wl90c_fourth()
650     t[113] = 1/pow(rhoa,3.666666666666667);  in wl90c_fourth()
651     t[114] = 1.62*t[59]*t[5]*t[24]*t[32];  in wl90c_fourth()
652     t[115] = 0.036006*grada*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
653     t[116] = .02000333333333333*t[15]*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
654     t[117] = 1.62*t[10]*t[11]*t[5]*t[24]*t[32];  in wl90c_fourth()
655     t[118] = 0.018003*gradb*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
656     t[119] = 0.018003*grada*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
657     t[120] = .02000333333333333*grada*gradb*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
658     t[121] = 0.108018*t[29]*t[10]*t[12]*t[5]*t[9]*t[32];  in wl90c_fourth()
659     t[122] = -0.0270045*t[29]*t[10]*t[12]*t[19]*t[22]*t[14];  in wl90c_fourth()
660     t[123] = 0.06001*grada*t[60]*t[12]*t[5]*t[9]*t[14];  in wl90c_fourth()
661     t[124] = 0.018003*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
662     t[125] = .07201199999999999*t[29]*t[6]*t[12]*t[5]*t[14];  in wl90c_fourth()
663     t[126] = -0.030005*grada*t[60]*t[12]*t[19]*t[22]*t[23];  in wl90c_fourth()
664     t[127] = .02000333333333333*grada*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
665     t[128] = 1.62*t[63]*t[5]*t[24]*t[32];  in wl90c_fourth()
666     t[129] = 0.036006*gradb*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
667     t[130] = .02000333333333333*t[16]*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
668     t[131] = 0.108018*t[29]*t[11]*t[12]*t[5]*t[9]*t[32];  in wl90c_fourth()
669     t[132] = -0.0270045*t[29]*t[11]*t[12]*t[19]*t[22]*t[14];  in wl90c_fourth()
670     t[133] = 0.06001*gradb*t[60]*t[12]*t[5]*t[9]*t[14];  in wl90c_fourth()
671     t[134] = 0.018003*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
672     t[135] = -0.030005*gradb*t[60]*t[12]*t[19]*t[22]*t[23];  in wl90c_fourth()
673     t[136] = .02000333333333333*gradb*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
674     t[137] = 0.06001*t[60]*t[12]*t[5]*t[9]*t[14];  in wl90c_fourth()
675     t[138] = -0.030005*t[60]*t[12]*t[19]*t[22]*t[23];  in wl90c_fourth()
676     t[139] = .02000333333333333*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
677     t[140] = pow(t[27],3.0);  in wl90c_fourth()
678     t[141] = 1/pow(rhob,4.333333333333333);  in wl90c_fourth()
679     t[142] = t[71]-9.333333333333334*gradb*t[141];  in wl90c_fourth()
680     t[143] = t[75]+t[74]+36.0*t[1]*t[37];  in wl90c_fourth()
681     t[144] = pow(t[28],3.0);  in wl90c_fourth()
682     t[145] = 0.12002*t[29]*t[12]*t[5]*t[52]*t[32];  in wl90c_fourth()
683     t[146] = -0.06001*t[29]*t[12]*t[19]*t[28]*t[27]*t[14];  in wl90c_fourth()
684     t[147] = .04000666666666666*t[29]*t[7]*t[5]*t[27]*t[14];  in wl90c_fourth()
685     t[148] = -0.06001*t[29]*t[12]*t[54]*t[5]*t[14];  in wl90c_fourth()
686     t[149] = -.02000333333333333*t[29]*t[7]*t[19]*t[28]*t[23];  in wl90c_fourth()
687     t[150] = 0.030005*t[29]*t[12]*t[55]*t[19]*t[23];  in wl90c_fourth()
688     t[151] = -0.0150025*t[29]*t[12]*t[41]*t[56]*t[23];  in wl90c_fourth()
689     t[152] = 0.108018*t[29]*t[10]*t[12]*t[5]*t[27]*t[32];  in wl90c_fourth()
690     t[153] = -0.0270045*t[29]*t[10]*t[12]*t[19]*t[28]*t[14];  in wl90c_fourth()
691     t[154] = 0.06001*grada*t[60]*t[12]*t[5]*t[27]*t[14];  in wl90c_fourth()
692     t[155] = -0.030005*grada*t[60]*t[12]*t[19]*t[28]*t[23];  in wl90c_fourth()
693     t[156] = 1/pow(rhob,3.666666666666667);  in wl90c_fourth()
694     t[157] = 0.108018*t[29]*t[11]*t[12]*t[5]*t[27]*t[32];  in wl90c_fourth()
695     t[158] = -0.0270045*t[29]*t[11]*t[12]*t[19]*t[28]*t[14];  in wl90c_fourth()
696     t[159] = 0.06001*gradb*t[60]*t[12]*t[5]*t[27]*t[14];  in wl90c_fourth()
697     t[160] = .07201199999999999*t[29]*t[26]*t[12]*t[5]*t[14];  in wl90c_fourth()
698     t[161] = -0.030005*gradb*t[60]*t[12]*t[19]*t[28]*t[23];  in wl90c_fourth()
699     t[162] = 0.06001*t[60]*t[12]*t[5]*t[27]*t[14];  in wl90c_fourth()
700     t[163] = -0.030005*t[60]*t[12]*t[19]*t[28]*t[23];  in wl90c_fourth()
701     t[164] = 1/pow(rhoa,4.0);  in wl90c_fourth()
702     t[165] = pow(grada,3.0);  in wl90c_fourth()
703     t[166] = 1/pow(t[17],5.0);  in wl90c_fourth()
704     t[167] = 0.0972162*t[29]*t[59]*t[12]*t[5]*t[32];  in wl90c_fourth()
705     t[168] = 0.108018*grada*t[60]*t[10]*t[12]*t[5]*t[14];  in wl90c_fourth()
706     t[169] = 0.18003*t[15]*t[166]*t[12]*t[5]*t[23];  in wl90c_fourth()
707     t[170] = 0.054009*t[60]*t[10]*t[12]*t[5]*t[14];  in wl90c_fourth()
708     t[171] = 0.18003*grada*t[166]*t[12]*t[5]*t[23];  in wl90c_fourth()
709     t[172] = 1/pow(rhob,4.0);  in wl90c_fourth()
710     t[173] = pow(gradb,3.0);  in wl90c_fourth()
711     t[174] = 0.0972162*t[29]*t[63]*t[12]*t[5]*t[32];  in wl90c_fourth()
712     t[175] = 0.108018*gradb*t[60]*t[11]*t[12]*t[5]*t[14];  in wl90c_fourth()
713     t[176] = 0.18003*t[16]*t[166]*t[12]*t[5]*t[23];  in wl90c_fourth()
714     t[177] = 0.054009*t[60]*t[11]*t[12]*t[5]*t[14];  in wl90c_fourth()
715     t[178] = 0.18003*gradb*t[166]*t[12]*t[5]*t[23];  in wl90c_fourth()
716     t[179] = 0.18003*t[166]*t[12]*t[5]*t[23];  in wl90c_fourth()
717     t[180] = 1/pow(t[13],5.0);  in wl90c_fourth()
718     t[181] = 1/pow(rhoa,5.333333333333333);  in wl90c_fourth()
719     t[182] = 1/pow(t[3],4.333333333333333);  in wl90c_fourth()
720     t[183] = 2.144421450022617*t[182];  in wl90c_fourth()
721     t[184] = -120.0*t[2]/pow(t[3],6.0);  in wl90c_fourth()
722     t[185] = -72.0*t[37];  in wl90c_fourth()
723     t[186] = 1/pow(t[5],7.0);  in wl90c_fourth()
724     t[187] = 0.207441975308642*t[17]*t[182]*t[5]*t[23];  in wl90c_fourth()
725     t[188] = -2.144421450022617*t[182]*t[5]*t[14]*t[18];  in wl90c_fourth()
726     t[189] = .1600266666666667*t[17]*t[34]*t[5]*t[9]*t[27]*t[32];  in wl90c_fourth()
727     t[190] = -.04000666666666666*t[17]*t[34]*t[19]*t[28]*t[9]*t[14];  in wl90c_fourth()
728     t[191] = -.04000666666666666*t[17]*t[34]*t[19]*t[22]*t[27]*t[14];  in wl90c_fourth()
729     t[192] = 1.929979305020355*t[70]*t[5]*t[24]*t[14];  in wl90c_fourth()
730     t[193] = 1/pow(t[3],4.666666666666667);  in wl90c_fourth()
731     t[194] = -.02206058249414695*t[17]*t[193]*t[5]*t[14];  in wl90c_fourth()
732     t[195] = -.02000333333333333*t[17]*t[34]*t[41]*t[22]*t[28]*t[23];  in wl90c_fourth()
733     t[196] = .04000666666666666*t[17]*t[34]*t[44]*t[19]*t[23];  in wl90c_fourth()
734     t[197] = .05600933333333333*t[17]*t[10]*t[70]*t[5]*t[14];  in wl90c_fourth()
735     t[198] = -.06223259259259259*grada*t[29]*t[70]*t[5]*t[23];  in wl90c_fourth()
736     t[199] = .05600933333333333*t[17]*t[11]*t[70]*t[5]*t[14];  in wl90c_fourth()
737     t[200] = -.06223259259259259*gradb*t[29]*t[70]*t[5]*t[23];  in wl90c_fourth()
738     t[201] = -.06223259259259259*t[29]*t[70]*t[5]*t[23];  in wl90c_fourth()
739     t[202] = -1.157987583012213*t[10]*t[70]*t[5]*t[32]*t[18];  in wl90c_fourth()
740     t[203] = -10.8*t[10]*t[5]*t[9]*t[27]*t[24]*t[68];  in wl90c_fourth()
741     t[204] = -.08001333333333333*grada*t[29]*t[7]*t[5]*t[9]*t[27]*t[32];  in wl90c_fourth()
742     t[205] = 1.8*t[10]*t[19]*t[28]*t[9]*t[24]*t[32];  in wl90c_fourth()
743     t[206] = 1.8*t[10]*t[19]*t[22]*t[27]*t[24]*t[32];  in wl90c_fourth()
744     t[207] = -4.8*t[6]*t[5]*t[27]*t[24]*t[32];  in wl90c_fourth()
745     t[208] = .9925607854390399*t[10]*t[34]*t[5]*t[24]*t[32];  in wl90c_fourth()
746     t[209] = .02000333333333333*grada*t[29]*t[7]*t[19]*t[28]*t[9]*t[14];  in wl90c_fourth()
747     t[210] = .02000333333333333*grada*t[29]*t[7]*t[19]*t[22]*t[27]*t[14];  in wl90c_fourth()
748     t[211] = 0.45*t[10]*t[41]*t[22]*t[28]*t[24]*t[14];  in wl90c_fourth()
749     t[212] = 1.2*t[6]*t[19]*t[28]*t[24]*t[14];  in wl90c_fourth()
750     t[213] = -0.9*t[10]*t[44]*t[19]*t[24]*t[14];  in wl90c_fourth()
751     t[214] = 1/pow(t[3],3.666666666666667);  in wl90c_fourth()
752     t[215] = .04963631061183065*grada*t[29]*t[214]*t[5]*t[14];  in wl90c_fourth()
753     t[216] = .01000166666666667*grada*t[29]*t[7]*t[41]*t[22]*t[28]*t[23];  in wl90c_fourth()
754     t[217] = -.02000333333333333*grada*t[29]*t[7]*t[44]*t[19]*t[23];  in wl90c_fourth()
755     t[218] = -1.157987583012213*t[11]*t[70]*t[5]*t[32]*t[18];  in wl90c_fourth()
756     t[219] = -10.8*t[11]*t[5]*t[9]*t[27]*t[24]*t[68];  in wl90c_fourth()
757     t[220] = -.08001333333333333*gradb*t[29]*t[7]*t[5]*t[9]*t[27]*t[32];  in wl90c_fourth()
758     t[221] = 1.8*t[11]*t[19]*t[28]*t[9]*t[24]*t[32];  in wl90c_fourth()
759     t[222] = -4.8*t[26]*t[5]*t[9]*t[24]*t[32];  in wl90c_fourth()
760     t[223] = 1.8*t[11]*t[19]*t[22]*t[27]*t[24]*t[32];  in wl90c_fourth()
761     t[224] = .9925607854390399*t[11]*t[34]*t[5]*t[24]*t[32];  in wl90c_fourth()
762     t[225] = .02000333333333333*gradb*t[29]*t[7]*t[19]*t[28]*t[9]*t[14];  in wl90c_fourth()
763     t[226] = .02000333333333333*gradb*t[29]*t[7]*t[19]*t[22]*t[27]*t[14];  in wl90c_fourth()
764     t[227] = 1.2*t[26]*t[19]*t[22]*t[24]*t[14];  in wl90c_fourth()
765     t[228] = 0.45*t[11]*t[41]*t[22]*t[28]*t[24]*t[14];  in wl90c_fourth()
766     t[229] = -0.9*t[11]*t[44]*t[19]*t[24]*t[14];  in wl90c_fourth()
767     t[230] = .04963631061183065*gradb*t[29]*t[214]*t[5]*t[14];  in wl90c_fourth()
768     t[231] = .01000166666666667*gradb*t[29]*t[7]*t[41]*t[22]*t[28]*t[23];  in wl90c_fourth()
769     t[232] = -.02000333333333333*gradb*t[29]*t[7]*t[44]*t[19]*t[23];  in wl90c_fourth()
770     t[233] = -.08001333333333333*t[29]*t[7]*t[5]*t[9]*t[27]*t[32];  in wl90c_fourth()
771     t[234] = .02000333333333333*t[29]*t[7]*t[19]*t[28]*t[9]*t[14];  in wl90c_fourth()
772     t[235] = .02000333333333333*t[29]*t[7]*t[19]*t[22]*t[27]*t[14];  in wl90c_fourth()
773     t[236] = .04963631061183065*t[29]*t[214]*t[5]*t[14];  in wl90c_fourth()
774     t[237] = .01000166666666667*t[29]*t[7]*t[41]*t[22]*t[28]*t[23];  in wl90c_fourth()
775     t[238] = -.02000333333333333*t[29]*t[7]*t[44]*t[19]*t[23];  in wl90c_fourth()
776     t[239] = .04320719999999999*t[17]*t[59]*t[34]*t[5]*t[32];  in wl90c_fourth()
777     t[240] = -0.048008*grada*t[29]*t[10]*t[34]*t[5]*t[14];  in wl90c_fourth()
778     t[241] = -.02667111111111111*t[15]*t[60]*t[34]*t[5]*t[23];  in wl90c_fourth()
779     t[242] = .04320719999999999*t[17]*t[10]*t[11]*t[34]*t[5]*t[32];  in wl90c_fourth()
780     t[243] = -0.024004*gradb*t[29]*t[10]*t[34]*t[5]*t[14];  in wl90c_fourth()
781     t[244] = -0.024004*grada*t[29]*t[11]*t[34]*t[5]*t[14];  in wl90c_fourth()
782     t[245] = -.02667111111111111*grada*gradb*t[60]*t[34]*t[5]*t[23];  in wl90c_fourth()
783     t[246] = -0.024004*t[29]*t[10]*t[34]*t[5]*t[14];  in wl90c_fourth()
784     t[247] = -.02667111111111111*grada*t[60]*t[34]*t[5]*t[23];  in wl90c_fourth()
785     t[248] = .04320719999999999*t[17]*t[63]*t[34]*t[5]*t[32];  in wl90c_fourth()
786     t[249] = -0.048008*gradb*t[29]*t[11]*t[34]*t[5]*t[14];  in wl90c_fourth()
787     t[250] = -.02667111111111111*t[16]*t[60]*t[34]*t[5]*t[23];  in wl90c_fourth()
788     t[251] = -0.024004*t[29]*t[11]*t[34]*t[5]*t[14];  in wl90c_fourth()
789     t[252] = -.02667111111111111*gradb*t[60]*t[34]*t[5]*t[23];  in wl90c_fourth()
790     t[253] = -.02667111111111111*t[60]*t[34]*t[5]*t[23];  in wl90c_fourth()
791     t[254] = -4.374*t[164]*t[5]*t[24]*t[68];  in wl90c_fourth()
792     t[255] = -0.0972162*grada*t[29]*t[59]*t[7]*t[5]*t[32];  in wl90c_fourth()
793     t[256] = -.05400899999999999*t[15]*t[60]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
794     t[257] = .05400899999999999*t[29]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
795     t[258] = -.06000999999999999*t[165]*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
796     t[259] = .06000999999999999*grada*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
797     t[260] = -4.374*t[59]*t[11]*t[5]*t[24]*t[68];  in wl90c_fourth()
798     t[261] = -0.0324054*gradb*t[29]*t[59]*t[7]*t[5]*t[32];  in wl90c_fourth()
799     t[262] = -0.0648108*grada*t[29]*t[10]*t[11]*t[7]*t[5]*t[32];  in wl90c_fourth()
800     t[263] = -0.036006*grada*gradb*t[60]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
801     t[264] = -0.018003*t[15]*t[60]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
802     t[265] = -.06000999999999999*t[15]*gradb*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
803     t[266] = -0.0324054*t[29]*t[59]*t[7]*t[5]*t[32];  in wl90c_fourth()
804     t[267] = -0.036006*grada*t[60]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
805     t[268] = -.06000999999999999*t[15]*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
806     t[269] = -4.374*t[10]*t[63]*t[5]*t[24]*t[68];  in wl90c_fourth()
807     t[270] = -0.0324054*grada*t[29]*t[63]*t[7]*t[5]*t[32];  in wl90c_fourth()
808     t[271] = -0.0648108*gradb*t[29]*t[10]*t[11]*t[7]*t[5]*t[32];  in wl90c_fourth()
809     t[272] = -0.018003*t[16]*t[60]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
810     t[273] = -0.036006*grada*gradb*t[60]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
811     t[274] = -.06000999999999999*grada*t[16]*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
812     t[275] = -0.0324054*t[29]*t[10]*t[11]*t[7]*t[5]*t[32];  in wl90c_fourth()
813     t[276] = -0.018003*gradb*t[60]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
814     t[277] = -0.018003*grada*t[60]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
815     t[278] = -.06000999999999999*grada*gradb*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
816     t[279] = -0.018003*t[60]*t[10]*t[7]*t[5]*t[14];  in wl90c_fourth()
817     t[280] = -.06000999999999999*grada*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
818     t[281] = -4.374*t[172]*t[5]*t[24]*t[68];  in wl90c_fourth()
819     t[282] = -0.0972162*gradb*t[29]*t[63]*t[7]*t[5]*t[32];  in wl90c_fourth()
820     t[283] = -.05400899999999999*t[16]*t[60]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
821     t[284] = .05400899999999999*t[29]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
822     t[285] = -.06000999999999999*t[173]*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
823     t[286] = .06000999999999999*gradb*t[60]*t[7]*t[5]*t[23];  in wl90c_fourth()
824     t[287] = -0.0324054*t[29]*t[63]*t[7]*t[5]*t[32];  in wl90c_fourth()
825     t[288] = -0.036006*gradb*t[60]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
826     t[289] = -.06000999999999999*t[16]*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
827     t[290] = -0.018003*t[60]*t[11]*t[7]*t[5]*t[14];  in wl90c_fourth()
828     t[291] = -.06000999999999999*gradb*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
829     t[292] = -.06000999999999999*t[166]*t[7]*t[5]*t[23];  in wl90c_fourth()
830     t[293] = 1/pow(rhob,5.333333333333333);  in wl90c_fourth()
831     t[294] = 1/pow(t[17],7.0);  in wl90c_fourth()
832     t[295] = -0.18003*t[60]*t[12]*t[5]*t[23];  in wl90c_fourth()
833     t[296] = 0.2916486*t[29]*t[10]*t[11]*t[12]*t[5]*t[32];  in wl90c_fourth()
834     t[297] = 0.162027*gradb*t[60]*t[10]*t[12]*t[5]*t[14];  in wl90c_fourth()
835     t[298] = 0.162027*grada*t[60]*t[11]*t[12]*t[5]*t[14];  in wl90c_fourth()
836     t[299] = 0.54009*grada*gradb*t[166]*t[12]*t[5]*t[23];  in wl90c_fourth()
837     dfdra = -1.0*t[14]*t[18]*t[5]*t[9]+t[25]+0.5*t[18]*t[19]*t[22]*t[23];  in wl90c_fourth()
838     dfdrb = -1.0*t[14]*t[18]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[28]+t[25];  in wl90c_fourth()
839     dfdga = 0.06001*grada*t[29]*t[12]*t[5]*t[23]-0.9*t[10]*t[5]*t[14]*t[18];  in wl90c_fourth()
840     dfdgb = 0.06001*gradb*t[29]*t[12]*t[5]*t[23]-0.9*t[11]*t[5]*t[14]*t[18];  in wl90c_fourth()
841     dfdab = t[30];  in wl90c_fourth()
842 …t[14]*t[24]*t[5]*t[9]-1.0*t[14]*t[18]*t[19]*t[22]*t[9]-1.0*t[14]*t[18]*t[36]*t[5]+2.0*t[18]*t[31]*…  in wl90c_fourth()
843 …t[18]*t[27]*t[32]*t[5]*t[9]-1.0*t[14]*t[24]*t[5]*t[9]-0.5*t[14]*t[18]*t[19]*t[28]*t[9]-1.0*t[14]*t…  in wl90c_fourth()
844 …t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[9]*t[14]…  in wl90c_fourth()
845 …t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[9]*t[14]…  in wl90c_fourth()
846     d2fdraab = t[51]+t[50]+t[49];  in wl90c_fourth()
847 …t[18]*t[23]*t[41]*t[56]+0.5*t[18]*t[19]*t[23]*t[55]-1.0*t[14]*t[18]*t[5]*t[54]+2.0*t[18]*t[32]*t[5…  in wl90c_fourth()
848 …t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[27]*t[14…  in wl90c_fourth()
849 …t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[27]*t[14…  in wl90c_fourth()
850     d2fdrbab = t[51]+t[58]+t[57];  in wl90c_fourth()
851 …d2fdgaga = t[30]-0.06001*t[15]*t[60]*t[12]*t[5]*t[23]-0.108018*grada*t[29]*t[10]*t[12]*t[5]*t[14]+…  in wl90c_fourth()
852 …gradb*t[60]*t[12]*t[5]*t[23]-0.054009*grada*t[29]*t[11]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[10…  in wl90c_fourth()
853     d2fdgaab = t[62]+t[61];  in wl90c_fourth()
854 …d2fdgbgb = t[30]-0.06001*t[16]*t[60]*t[12]*t[5]*t[23]-0.108018*gradb*t[29]*t[11]*t[12]*t[5]*t[14]+…  in wl90c_fourth()
855     d2fdgbab = t[65]+t[64];  in wl90c_fourth()
856     d2fdabab = t[66];  in wl90c_fourth()
857 …t[18]*t[32]*t[36]*t[5]*t[9]+0.75*t[14]*t[18]*t[41]*t[42]*t[9]-1.5*t[14]*t[18]*t[19]*t[40]*t[9]-3.0…  in wl90c_fourth()
858 …t[14]*t[18]*t[19]*t[44]*t[9]+0.5*t[14]*t[18]*t[22]*t[28]*t[41]*t[9]+2.0*t[18]*t[19]*t[22]*t[27]*t[…  in wl90c_fourth()
859 …t[29]*t[12]*t[41]*t[42]*t[23]+t[90]+0.030005*grada*t[29]*t[12]*t[40]*t[19]*t[23]-.0200033333333333…  in wl90c_fourth()
860 …t[29]*t[12]*t[41]*t[42]*t[23]+t[92]+0.030005*gradb*t[29]*t[12]*t[40]*t[19]*t[23]-.0200033333333333…  in wl90c_fourth()
861     d3fdraraab = t[100]+t[99]+t[98]+t[97]+t[96]+t[95]+t[94]+t[93];  in wl90c_fourth()
862 …t[18]*t[5]*t[52]*t[68]*t[9]+0.25*t[14]*t[18]*t[41]*t[56]*t[9]-0.5*t[14]*t[18]*t[19]*t[55]*t[9]+2.0…  in wl90c_fourth()
863 …t[90]+0.030005*grada*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*grada*t[29]*t[7]*t[19]*t[28]…  in wl90c_fourth()
864 …t[92]+0.030005*gradb*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*gradb*t[29]*t[7]*t[19]*t[28]…  in wl90c_fourth()
865     d3fdrarbab = t[99]+t[112]+t[111]+t[110]+t[109]+t[108]+t[107]+t[106]+t[105]+t[104]+t[103];  in wl90c_fourth()
866 …t[51]+t[116]+t[50]-0.030005*t[15]*t[60]*t[12]*t[19]*t[22]*t[23]+0.144024*grada*t[29]*t[6]*t[12]*t[…  in wl90c_fourth()
867 …t[120]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+.07201199999999999*gradb*t[29]*t[6]*t[12…  in wl90c_fourth()
868     d3fdragaab = t[127]+t[126]+t[125]+t[124]+t[123]+t[122]+t[121];  in wl90c_fourth()
869 …t[51]+t[130]+t[50]-0.030005*t[16]*t[60]*t[12]*t[19]*t[22]*t[23]+t[129]+t[49]+0.06001*t[16]*t[60]*t…  in wl90c_fourth()
870     d3fdragbab = t[136]+t[135]+t[134]+t[133]+t[132]+t[131];  in wl90c_fourth()
871     d3fdraabab = t[139]+t[138]+t[137];  in wl90c_fourth()
872 …t[79]+0.375*t[144]*t[18]*t[23]*t[77]-6.0*t[140]*t[18]*t[5]*t[68]+0.75*t[14]*t[18]*t[27]*t[41]*t[56…  in wl90c_fourth()
873 …t[29]*t[12]*t[41]*t[56]*t[23]+t[90]+0.030005*grada*t[29]*t[12]*t[55]*t[19]*t[23]-.0200033333333333…  in wl90c_fourth()
874 …t[29]*t[12]*t[41]*t[56]*t[23]+t[92]+0.030005*gradb*t[29]*t[12]*t[55]*t[19]*t[23]-.0200033333333333…  in wl90c_fourth()
875     d3fdrbrbab = t[151]+t[99]+t[150]+t[149]+t[148]+t[147]+t[146]+t[145];  in wl90c_fourth()
876 …t[51]+t[116]+t[58]-0.030005*t[15]*t[60]*t[12]*t[19]*t[28]*t[23]+t[115]+t[57]+0.06001*t[15]*t[60]*t…  in wl90c_fourth()
877 …t[120]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+.07201199999999999*grada*t[29]*t[26]*t[1…  in wl90c_fourth()
878     d3fdrbgaab = t[127]+t[155]+t[124]+t[154]+t[153]+t[152];  in wl90c_fourth()
879 …t[51]+t[130]+t[58]-0.030005*t[16]*t[60]*t[12]*t[19]*t[28]*t[23]+0.144024*gradb*t[29]*t[26]*t[12]*t…  in wl90c_fourth()
880     d3fdrbgbab = t[136]+t[161]+t[160]+t[134]+t[159]+t[158]+t[157];  in wl90c_fourth()
881     d3fdrbabab = t[139]+t[163]+t[162];  in wl90c_fourth()
882 …t[60]*t[12]*t[5]*t[23]+0.18003*t[165]*t[166]*t[12]*t[5]*t[23]-0.162027*t[29]*t[10]*t[12]*t[5]*t[14…  in wl90c_fourth()
883 …t[65]+0.18003*t[15]*gradb*t[166]*t[12]*t[5]*t[23]+t[64]+0.054009*t[15]*t[60]*t[11]*t[12]*t[5]*t[14…  in wl90c_fourth()
884     d3fdgagaab = t[66]+t[169]+t[168]+t[167];  in wl90c_fourth()
885 …t[62]+0.18003*grada*t[16]*t[166]*t[12]*t[5]*t[23]+0.108018*grada*gradb*t[60]*t[11]*t[12]*t[5]*t[14…  in wl90c_fourth()
886 …t[166]*t[12]*t[5]*t[23]+0.054009*grada*t[60]*t[11]*t[12]*t[5]*t[14]+0.054009*gradb*t[60]*t[10]*t[1…  in wl90c_fourth()
887     d3fdgaabab = t[171]+t[170];  in wl90c_fourth()
888 …t[60]*t[12]*t[5]*t[23]+0.18003*t[173]*t[166]*t[12]*t[5]*t[23]-0.162027*t[29]*t[11]*t[12]*t[5]*t[14…  in wl90c_fourth()
889     d3fdgbgbab = t[66]+t[176]+t[175]+t[174];  in wl90c_fourth()
890     d3fdgbabab = t[178]+t[177];  in wl90c_fourth()
891     d3fdababab = t[179];  in wl90c_fourth()
892 …t[18]*t[180]*t[5]*pow(t[9],4.0)-1.5*t[14]*t[18]*t[77]*t[78]*t[9]-2.0*t[14]*t[18]*t[19]*t[76]*t[9]+…  in wl90c_fourth()
893 …t[14]*t[18]*t[19]*t[82]*t[9]-1.125*t[14]*t[18]*t[28]*t[42]*t[77]*t[9]-18.0*t[18]*t[27]*t[36]*t[5]*…  in wl90c_fourth()
894 …t[29]*t[12]*t[77]*t[78]*t[23]+0.0150025*grada*t[29]*t[7]*t[41]*t[42]*t[23]+t[198]-0.030005*grada*t…  in wl90c_fourth()
895 …t[29]*t[12]*t[77]*t[78]*t[23]+0.0150025*gradb*t[29]*t[7]*t[41]*t[42]*t[23]+t[200]-0.030005*gradb*t…  in wl90c_fourth()
896 …t[29]*t[12]*t[77]*t[78]*t[23]+0.0150025*t[29]*t[7]*t[41]*t[42]*t[23]+t[201]-0.030005*t[29]*t[7]*t[…  in wl90c_fourth()
897 …t[14]*t[18]*t[22]*t[56]*t[77]*t[9]-12.0*t[24]*t[5]*t[52]*t[68]*t[9]-6.0*t[18]*t[19]*t[22]*t[52]*t[…  in wl90c_fourth()
898 …t[29]*t[7]*t[41]*t[42]*t[23]+t[198]+t[217]-.01000166666666667*grada*t[29]*t[7]*t[40]*t[19]*t[23]+0…  in wl90c_fourth()
899 …t[29]*t[7]*t[41]*t[42]*t[23]+t[200]+t[232]-.01000166666666667*gradb*t[29]*t[7]*t[40]*t[19]*t[23]+0…  in wl90c_fourth()
900 …t[29]*t[7]*t[41]*t[42]*t[23]+t[201]+t[238]-.01000166666666667*t[29]*t[7]*t[40]*t[19]*t[23]+0.03000…  in wl90c_fourth()
901 …t[18]*t[32]*t[5]/pow(rhoa,4.666666666666667)+t[99]+t[98]+t[97]+t[96]+t[95]+t[94]+t[93]-0.648108*gr…  in wl90c_fourth()
902 …t[60]*t[12]*t[41]*t[42]*t[23]+t[245]-0.030005*grada*gradb*t[60]*t[12]*t[40]*t[19]*t[23]+.020003333…  in wl90c_fourth()
903 …t[60]*t[12]*t[41]*t[42]*t[23]+t[247]-0.030005*grada*t[60]*t[12]*t[40]*t[19]*t[23]+.020003333333333…  in wl90c_fourth()
904 …t[100]+0.0150025*t[16]*t[60]*t[12]*t[41]*t[42]*t[23]+t[99]+t[250]+t[98]-0.030005*t[16]*t[60]*t[12]…  in wl90c_fourth()
905 …t[60]*t[12]*t[41]*t[42]*t[23]+t[252]-0.030005*gradb*t[60]*t[12]*t[40]*t[19]*t[23]+.020003333333333…  in wl90c_fourth()
906 …t[60]*t[12]*t[41]*t[42]*t[23]+t[253]-0.030005*t[60]*t[12]*t[40]*t[19]*t[23]+.02000333333333333*t[6…  in wl90c_fourth()
907 …t[14]*t[144]*t[18]*t[77]*t[9]-18.0*t[18]*t[27]*t[5]*t[54]*t[68]*t[9]-18.0*t[24]*t[5]*t[52]*t[68]*t…  in wl90c_fourth()
908 …t[29]*t[7]*t[41]*t[56]*t[23]+0.02250375*grada*t[29]*t[12]*t[77]*t[22]*t[56]*t[23]+t[198]+t[217]-.0…  in wl90c_fourth()
909 …t[29]*t[7]*t[41]*t[56]*t[23]+0.02250375*gradb*t[29]*t[12]*t[77]*t[22]*t[56]*t[23]+t[200]+t[232]-.0…  in wl90c_fourth()
910 …t[29]*t[7]*t[41]*t[56]*t[23]+0.02250375*t[29]*t[12]*t[77]*t[22]*t[56]*t[23]+t[201]+t[238]-.0100016…  in wl90c_fourth()
911 …t[99]+t[241]+t[112]-0.030005*t[15]*t[60]*t[12]*t[44]*t[19]*t[23]+t[111]+.01000166666666667*t[15]*t…  in wl90c_fourth()
912 …t[245]-0.030005*grada*gradb*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*grada*gradb*t[60]*t[7…  in wl90c_fourth()
913 …t[247]-0.030005*grada*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*grada*t[60]*t[7]*t[19]*t[28…  in wl90c_fourth()
914 …t[99]+t[250]+t[112]-0.030005*t[16]*t[60]*t[12]*t[44]*t[19]*t[23]+t[111]+.01000166666666667*t[16]*t…  in wl90c_fourth()
915 …t[252]-0.030005*gradb*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*gradb*t[60]*t[7]*t[19]*t[28…  in wl90c_fourth()
916 …t[253]-0.030005*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*t[60]*t[7]*t[19]*t[28]*t[23]+0.01…  in wl90c_fourth()
917 …t[18]*t[5]*t[68]/pow(rhoa,5.0)-.8749458000000001*grada*t[29]*t[59]*t[12]*t[5]*t[9]*t[68]-0.324054*…  in wl90c_fourth()
918 …t[136]+t[265]+t[135]+0.090015*t[15]*gradb*t[166]*t[12]*t[19]*t[22]*t[23]-0.144024*grada*gradb*t[60…  in wl90c_fourth()
919 …t[139]+t[268]+t[138]+0.090015*t[15]*t[166]*t[12]*t[19]*t[22]*t[23]-0.144024*grada*t[60]*t[6]*t[12]…  in wl90c_fourth()
920 …t[127]+t[274]+t[126]+0.090015*grada*t[16]*t[166]*t[12]*t[19]*t[22]*t[23]+t[125]-.07201199999999999…  in wl90c_fourth()
921 …t[278]+0.090015*grada*gradb*t[166]*t[12]*t[19]*t[22]*t[23]-.07201199999999999*gradb*t[60]*t[6]*t[1…  in wl90c_fourth()
922 …t[280]+0.090015*grada*t[166]*t[12]*t[19]*t[22]*t[23]-.07201199999999999*t[60]*t[6]*t[12]*t[5]*t[14…  in wl90c_fourth()
923 …t[286]+t[285]-0.090015*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+0.090015*t[173]*t[166]*t[12]*t[19]*t[22…  in wl90c_fourth()
924 …t[139]+t[289]+t[138]+0.090015*t[16]*t[166]*t[12]*t[19]*t[22]*t[23]+t[288]+t[137]-0.18003*t[16]*t[1…  in wl90c_fourth()
925 …t[291]+0.090015*gradb*t[166]*t[12]*t[19]*t[22]*t[23]+t[290]-0.18003*gradb*t[166]*t[12]*t[5]*t[9]*t…  in wl90c_fourth()
926 …d4fdraababab = t[292]+0.090015*t[166]*t[12]*t[19]*t[22]*t[23]-0.18003*t[166]*t[12]*t[5]*t[9]*t[14];  in wl90c_fourth()
927 …t[18]*t[23]*t[55]*t[56]*t[77]-1.5*t[14]*t[144]*t[18]*t[27]*t[77]+1.5*t[144]*t[23]*t[24]*t[77]+0.5*…  in wl90c_fourth()
928 …t[29]*t[12]*t[77]*t[144]*t[23]+0.0150025*grada*t[29]*t[7]*t[41]*t[56]*t[23]+t[198]-0.030005*grada*…  in wl90c_fourth()
929 …t[29]*t[12]*t[77]*t[144]*t[23]+0.0150025*gradb*t[29]*t[7]*t[41]*t[56]*t[23]+t[200]-0.030005*gradb*…  in wl90c_fourth()
930 …t[29]*t[12]*t[77]*t[144]*t[23]+0.0150025*t[29]*t[7]*t[41]*t[56]*t[23]+t[201]-0.030005*t[29]*t[7]*t…  in wl90c_fourth()
931 …t[151]+0.0150025*t[15]*t[60]*t[12]*t[41]*t[56]*t[23]+t[99]+t[241]+t[150]-0.030005*t[15]*t[60]*t[12…  in wl90c_fourth()
932 …t[60]*t[12]*t[41]*t[56]*t[23]+t[245]-0.030005*grada*gradb*t[60]*t[12]*t[55]*t[19]*t[23]+.020003333…  in wl90c_fourth()
933 …t[60]*t[12]*t[41]*t[56]*t[23]+t[247]-0.030005*grada*t[60]*t[12]*t[55]*t[19]*t[23]+.020003333333333…  in wl90c_fourth()
934 …t[18]*t[32]*t[5]/pow(rhob,4.666666666666667)+t[99]-0.648108*gradb*t[29]*t[11]*t[12]*t[5]*t[52]*t[6…  in wl90c_fourth()
935 …t[60]*t[12]*t[41]*t[56]*t[23]+t[252]-0.030005*gradb*t[60]*t[12]*t[55]*t[19]*t[23]+.020003333333333…  in wl90c_fourth()
936 …t[60]*t[12]*t[41]*t[56]*t[23]+t[253]-0.030005*t[60]*t[12]*t[55]*t[19]*t[23]+.02000333333333333*t[6…  in wl90c_fourth()
937 …t[259]+t[258]-0.090015*grada*t[60]*t[12]*t[19]*t[28]*t[23]+0.090015*t[165]*t[166]*t[12]*t[19]*t[28…  in wl90c_fourth()
938 …t[136]+t[265]+t[161]+0.090015*t[15]*gradb*t[166]*t[12]*t[19]*t[28]*t[23]+t[160]-.07201199999999999…  in wl90c_fourth()
939 …t[139]+t[268]+t[163]+0.090015*t[15]*t[166]*t[12]*t[19]*t[28]*t[23]+t[267]+t[162]-0.18003*t[15]*t[1…  in wl90c_fourth()
940 …t[127]+t[274]+t[155]+0.090015*grada*t[16]*t[166]*t[12]*t[19]*t[28]*t[23]-0.144024*grada*gradb*t[60…  in wl90c_fourth()
941 …t[278]+0.090015*grada*gradb*t[166]*t[12]*t[19]*t[28]*t[23]-.07201199999999999*grada*t[60]*t[26]*t[…  in wl90c_fourth()
942 …t[280]+0.090015*grada*t[166]*t[12]*t[19]*t[28]*t[23]+t[279]-0.18003*grada*t[166]*t[12]*t[5]*t[27]*…  in wl90c_fourth()
943 …t[18]*t[5]*t[68]/pow(rhob,5.0)-.8749458000000001*gradb*t[29]*t[63]*t[12]*t[5]*t[27]*t[68]-0.777729…  in wl90c_fourth()
944 …t[139]+t[289]+t[163]+0.090015*t[16]*t[166]*t[12]*t[19]*t[28]*t[23]-0.144024*gradb*t[60]*t[26]*t[12…  in wl90c_fourth()
945 …t[291]+0.090015*gradb*t[166]*t[12]*t[19]*t[28]*t[23]-.07201199999999999*t[60]*t[26]*t[12]*t[5]*t[1…  in wl90c_fourth()
946 …d4fdrbababab = t[292]+0.090015*t[166]*t[12]*t[19]*t[28]*t[23]-0.18003*t[166]*t[12]*t[5]*t[27]*t[14…  in wl90c_fourth()
947 …t[12]*t[23]*t[294]*t[5]*pow(grada,4.0)-1.04993496*grada*t[29]*t[164]*t[12]*t[5]*t[68]-.58329720000…  in wl90c_fourth()
948 …t[299]-0.90015*t[165]*gradb*t[294]*t[12]*t[5]*t[23]+t[298]-0.162027*t[165]*t[166]*t[11]*t[12]*t[5]…  in wl90c_fourth()
949 …t[166]*t[12]*t[5]*t[23]-0.90015*t[165]*t[294]*t[12]*t[5]*t[23]+0.162027*t[60]*t[10]*t[12]*t[5]*t[1…  in wl90c_fourth()
950 …t[66]+t[176]+t[169]-0.90015*t[15]*t[16]*t[294]*t[12]*t[5]*t[23]+t[175]-0.324054*t[15]*gradb*t[166]…  in wl90c_fourth()
951 …t[178]-0.90015*t[15]*gradb*t[294]*t[12]*t[5]*t[23]+t[177]-0.162027*t[15]*t[166]*t[11]*t[12]*t[5]*t…  in wl90c_fourth()
952 …4fdgagaabab = t[179]-0.90015*t[15]*t[294]*t[12]*t[5]*t[23]-0.324054*grada*t[166]*t[10]*t[12]*t[5]*…  in wl90c_fourth()
953 …t[299]-0.90015*grada*t[173]*t[294]*t[12]*t[5]*t[23]+t[298]-0.486081*grada*t[16]*t[166]*t[11]*t[12]…  in wl90c_fourth()
954 …t[171]-0.90015*grada*t[16]*t[294]*t[12]*t[5]*t[23]-0.324054*grada*gradb*t[166]*t[11]*t[12]*t[5]*t[…  in wl90c_fourth()
955 …t[294]*t[12]*t[5]*t[23]-0.162027*grada*t[166]*t[11]*t[12]*t[5]*t[14]-0.162027*gradb*t[166]*t[10]*t…  in wl90c_fourth()
956     d4fdgaababab = -0.90015*grada*t[294]*t[12]*t[5]*t[23]-0.162027*t[166]*t[10]*t[12]*t[5]*t[14];  in wl90c_fourth()
957 …t[12]*t[23]*t[294]*t[5]*pow(gradb,4.0)-1.04993496*gradb*t[29]*t[172]*t[12]*t[5]*t[68]-.58329720000…  in wl90c_fourth()
958 …t[166]*t[12]*t[5]*t[23]-0.90015*t[173]*t[294]*t[12]*t[5]*t[23]+0.162027*t[60]*t[11]*t[12]*t[5]*t[1…  in wl90c_fourth()
959 …4fdgbgbabab = t[179]-0.90015*t[16]*t[294]*t[12]*t[5]*t[23]-0.324054*gradb*t[166]*t[11]*t[12]*t[5]*…  in wl90c_fourth()
960     d4fdgbababab = -0.90015*gradb*t[294]*t[12]*t[5]*t[23]-0.162027*t[166]*t[11]*t[12]*t[5]*t[14];  in wl90c_fourth()
961     d4fdabababab = -0.90015*t[294]*t[12]*t[5]*t[23];  in wl90c_fourth()