diff options
Diffstat (limited to 'src/kij_odd.x')
-rw-r--r-- | src/kij_odd.x | 153 |
1 files changed, 153 insertions, 0 deletions
diff --git a/src/kij_odd.x b/src/kij_odd.x new file mode 100644 index 0000000..f9036f2 --- /dev/null +++ b/src/kij_odd.x @@ -0,0 +1,153 @@ + 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 |