o1 = mm**2 o2 = 5.0000000000000d1*o1 o3 = -8.0000000000000d1 + o2 o4 = 1/o3 o5 = 2.00000000000000d0*q(i,j,k) o6 = cos(o5) o7 = 5.0000000000000d0*o6 o8 = 3.00000000000000d0 + o7 o9 = dnnnnngtil o10 = dnnngtil o11 = -2.60000000000000d1*o10 o12 = dngtil o13 = 2.50000000000000d1*o12 o14 = o11 + o13 + o9 o15 = -3.00000000000000d0*eta(i,j,k) o16 = exp(o15) o17 = sin(phi(i,j,k)) o18 = o17**2 o19 = mm*phi(i,j,k) o20 = sin(o19) o21 = sin(q(i,j,k)) o22 = o21**2 o23 = 2.50000000000000d1*o1 o24 = -4.0000000000000d1 + o23 o25 = 1/o24 o26 = cos(phi(i,j,k)) o27 = o26**2 o28 = -1.60000000000000d1 + o1 o29 = 3.00000000000000d0*o6 o30 = 2.00000000000000d0 + o29 o31 = 2.00000000000000d0*o28*o30*o9 o32 = -1.07000000000000d2*o1 o33 = 6.0900000000000d2*o6 o34 = -1.59000000000000d2*o1*o6 o35 = 4.0700000000000d2 + o32 + o33 + o34 o36 = o10*o35 o37 = 3.5000000000000d1*o1 o38 = 5.0000000000000d0*o1 o39 = -1.30000000000000d1 + o38 o40 = 9.0000000000000d0*o39*o6 o41 = -8.3000000000000d1 + o37 + o40 o42 = 5.0000000000000d0*o12*o41 o43 = o31 + o36 + o42 o44 = o22**2 o45 = cos(q(i,j,k)) o46 = o45**2 o47 = -3.00000000000000d0*mm*o14*o20*o22*o4*o8 o48 = -(mm*o20*o22*o25*o43) o49 = o47 + o48 o50 = -8.0000000000000d0 + o38 o51 = 1/o50 o52 = cos(o19) o53 = -9.0000000000000d0*o1 o54 = 5.0000000000000d0*o50*o6 o55 = 8.0000000000000d0 + o53 + o54 o56 = -7.5000000000000d1*o12*o55 o57 = 2.88000000000000d2*o9 o58 = -1.20000000000000d1*o1*o9 o59 = -9.4000000000000d1*o1 o60 = -4.1500000000000d2*o6 o61 = 1.75000000000000d2*o1*o6 o62 = 1.15000000000000d2 + o59 + o60 + o61 o63 = 4.0000000000000d0*o10*o62 o64 = dnnnnnnngtil o65 = 3.2000000000000d1*o22*o64 o66 = -2.00000000000000d0*o1*o22*o64 o67 = -9.2800000000000d2*o22*o9 o68 = 1.03000000000000d2*o1*o22*o9 o69 = o57 + o58 + o63 + o65 + o66 + o67 + o68 o70 = 2.00000000000000d0*o69 o71 = o56 + o70 o72 = sin(o5) o73 = dnnnnnngtil o74 = 3.2000000000000d1*o73 o75 = -2.00000000000000d0*o1*o73 o76 = -3.2000000000000d1*o6*o73 o77 = 2.00000000000000d0*o1*o6*o73 o78 = dnnnngtil o79 = -3.4000000000000d2*o78 o80 = 5.5000000000000d1*o1*o78 o81 = -8.1200000000000d2*o6*o78 o82 = 1.70000000000000d1*o1*o6*o78 o83 = 7.0000000000000d0*o1 o84 = -2.09000000000000d2*o6 o85 = 8.9000000000000d1*o1*o6 o86 = -3.10000000000000d1 + o83 + o84 + o85 o87 = dnngtil o88 = -2.00000000000000d1*o86*o87 o89 = 7.0000000000000d0*o6 o90 = 1.00000000000000d0 + o89 o91 = gtil o92 = 7.5000000000000d1*o50*o90*o91 o93 = o74 + o75 + o76 + o77 + o79 + o80 + o81 + o82 + o88 + o92 o94 = -2.40000000000000d2*byJ o95 = 1.50000000000000d2*byJ*o1 o96 = -1.15200000000000d3*o46*o52*o78 o97 = 7.2000000000000d1*o1*o46*o52*o78 o98 = 3.9000000000000d1*o1 o99 = -2.68000000000000d2*o6 o100 = 1.00000000000000d2*o1*o6 o101 = 4.0000000000000d0*q(i,j,k) o102 = cos(o101) o103 = -5.8100000000000d2*o102 o104 = 2.45000000000000d2*o1*o102 o105 = -1.11000000000000d2 + o100 + o103 + o104 + o98 + o99 o106 = -5.0000000000000d0*o105*o52*o87 o107 = -1.92000000000000d2*o22*o46*o52*o73 o108 = 1.20000000000000d1*o1*o22*o46*o52*o73 o109 = 2.88000000000000d2*o22*o52*o78 o110 = -1.80000000000000d1*o1*o22*o52*o78 o111 = 5.5680000000000d3*o22*o46*o52*o78 o112 = -6.1800000000000d2*o1*o22*o46*o52*o78 o113 = 3.2000000000000d1*o44*o52*o73 o114 = -2.00000000000000d0*o1*o44*o52*o73 o115 = -9.2800000000000d2*o44*o52*o78 o116 = 1.03000000000000d2*o1*o44*o52*o78 o117 = o106 + o107 + o108 + o109 + o110 + o111 + o112 + o113 + o & 114 + o115 + o116 + o94 + o95 + o96 + o97 o118 = 2.00000000000000d0*o117 o119 = 2.00000000000000d1*o6 o120 = 3.5000000000000d1*o102 o121 = 9.0000000000000d0 + o119 + o120 o122 = 1.87500000000000d1*o121*o50*o52*o91 o123 = o118 + o122 o124 = o21*o22 kxx(i,j,k) = mm*o16*o20*o25*o27*o43*o44 + o16*o27*o46*o49 - 1.00 & 000000000000d-1*o123*o16*o17*o22*o26*o51 - 1.00000000000000d-1*o & 16*o17*o22*o26*o46*o51*o52*o71 + 3.00000000000000d0*mm*o14*o16*o & 18*o20*o22*o4*o8 + 5.0000000000000d-2*mm*o16*o20*o21*o27*o45*o51 & *o72*o93 kxy(i,j,k) = mm*o16*o17*o20*o25*o26*o43*o44 + o16*o17*o26*o46*o4 & 9 - 5.0000000000000d-2*o123*o16*o18*o22*o51 + 5.0000000000000d-2 & *o123*o16*o22*o27*o51 - 5.0000000000000d-2*o16*o18*o22*o46*o51*o & 52*o71 + 5.0000000000000d-2*o16*o22*o27*o46*o51*o52*o71 - 3.0000 & 0000000000d0*mm*o14*o16*o17*o20*o22*o26*o4*o8 + 5.0000000000000d & -2*mm*o16*o17*o20*o21*o26*o45*o51*o72*o93 kxz(i,j,k) = mm*o124*o16*o20*o25*o26*o43*o45 - o16*o21*o26*o45*o & 49 - 5.0000000000000d-2*o123*o16*o17*o21*o45*o51 + 5.00000000000 & 00d-2*o124*o16*o17*o45*o51*o52*o71 - 2.50000000000000d-2*mm*o16* & o20*o22*o26*o51*o72*o93 + 2.50000000000000d-2*mm*o16*o20*o26*o46 & *o51*o72*o93 kyy(i,j,k) = mm*o16*o18*o20*o25*o43*o44 + o16*o18*o46*o49 + 1.00 & 000000000000d-1*o123*o16*o17*o22*o26*o51 + 1.00000000000000d-1*o & 16*o17*o22*o26*o46*o51*o52*o71 + 3.00000000000000d0*mm*o14*o16*o & 20*o22*o27*o4*o8 + 5.0000000000000d-2*mm*o16*o18*o20*o21*o45*o51 & *o72*o93 kyz(i,j,k) = mm*o124*o16*o17*o20*o25*o43*o45 - o16*o17*o21*o45*o & 49 + 5.0000000000000d-2*o123*o16*o21*o26*o45*o51 - 5.00000000000 & 00d-2*o124*o16*o26*o45*o51*o52*o71 - 2.50000000000000d-2*mm*o16* & o17*o20*o22*o51*o72*o93 + 2.50000000000000d-2*mm*o16*o17*o20*o46 & *o51*o72*o93 kzz(i,j,k) = mm*o16*o20*o22*o25*o43*o46 + o16*o22*o49 - 5.000000 & 0000000d-2*mm*o16*o20*o21*o45*o51*o72*o93