diff options
Diffstat (limited to 'src/gij.x')
| -rw-r--r-- | src/gij.x | 49 |
1 files changed, 27 insertions, 22 deletions
@@ -1,24 +1,29 @@ - gxx(i,j,k) = 1.00000000000000d0 -c dxgxx(i,j,k) = 0 -c dygxx(i,j,k) = 0 -c dzgxx(i,j,k) = 0 - gxy(i,j,k) = 0 -c dxgxy(i,j,k) = 0 -c dygxy(i,j,k) = 0 -c dzgxy(i,j,k) = 0 + o1 = 2.d0*phi(i,j,k) + o2 = cos(o1) + o3 = -xi + o4 = cos(q(i,j,k)) + o5 = o4*xi + o6 = 1.d0 + o3 + o5 + o7 = -eta0_me + o8 = eta(i,j,k) + o7 + o9 = o8**2 + o10 = sigma_me**2 + o11 = 1/o10 + o12 = -(o11*o9) + o13 = exp(o12) + o14 = eta(i,j,k) + eta0_me + o15 = o14**2 + o16 = -(o11*o15) + o17 = exp(o16) + o18 = o13 + o17 + o19 = sin(q(i,j,k)) + o20 = o19**n_me + o21 = 2.d0*amp_me*o18*o20*o6 + o22 = exp(o21) + o23 = -1.d0 + o22 + gxx(i,j,k) = 5.d-1*(1.d0 + o22 + o2*o23) + gxy(i,j,k) = o23*cos(phi(i,j,k))*sin(phi(i,j,k)) gxz(i,j,k) = 0 -c dxgxz(i,j,k) = 0 -c dygxz(i,j,k) = 0 -c dzgxz(i,j,k) = 0 - gyy(i,j,k) = 1.00000000000000d0 -c dxgyy(i,j,k) = 0 -c dygyy(i,j,k) = 0 -c dzgyy(i,j,k) = 0 + gyy(i,j,k) = 5.d-1*(1.d0 + o22 - o2*o23) gyz(i,j,k) = 0 -c dxgyz(i,j,k) = 0 -c dygyz(i,j,k) = 0 -c dzgyz(i,j,k) = 0 - gzz(i,j,k) = 1.00000000000000d0 -c dxgzz(i,j,k) = 0 -c dygzz(i,j,k) = 0 -c dzgzz(i,j,k) = 0 + gzz(i,j,k) = o22 |