45 real(dp) :: g11_g(0:jmax,0:imax), g22_g(0:jmax,0:imax), &
46 g11_sgx(0:jmax,0:imax), g11_sgy(0:jmax,0:imax), &
47 g22_sgx(0:jmax,0:imax), g22_sgy(0:jmax,0:imax)
48 real(dp) :: xhelp, yhelp
71 k = (cos(0.25_dp*pi-0.5_dp*phi0))**2
72 else if (phi0 < (-eps))
then
73 k = (cos(0.25_dp*pi+0.5_dp*phi0))**2
75 stop
' metric: PHI0 must be different from zero!'
80 g11_g(j,i) = 1.0_dp/ &
81 ( k**2*(1.0_dp+(xi(i)**2+eta(j)**2)/(2.0_dp*r*k)**2)**2 )
82 g22_g(j,i) = g11_g(j,i)
90 xhelp = 0.5_dp*(xi(i)+xi(i+1))
91 g11_sgx(j,i) = 1.0_dp/ &
92 ( k**2*(1.0_dp+(xhelp**2+eta(j)**2)/(2.0_dp*r*k)**2)**2 )
93 g22_sgx(j,i) = g11_sgx(j,i)
99 yhelp = 0.5_dp*(eta(j)+eta(j+1))
100 g22_sgy(j,i) = 1.0_dp/ &
101 ( k**2*(1.0_dp+(xi(i)**2+yhelp**2)/(2.0_dp*r*k)**2)**2 )
102 g11_sgy(j,i) = g22_sgy(j,i)
112 g11_g(j,i) = r**2*(cos(eta(j)))**2
121 g11_sgx(j,i) = r**2*(cos(eta(j)))**2
128 yhelp = 0.5_dp*(eta(j)+eta(j+1))
130 g11_sgy(j,i) = r**2*(cos(yhelp))**2
141 sq_g11_g(j,i) = sqrt(g11_g(j,i))
142 sq_g22_g(j,i) = sqrt(g22_g(j,i))
143 insq_g11_g(j,i) = 1.0_dp/sq_g11_g(j,i)
144 insq_g22_g(j,i) = 1.0_dp/sq_g22_g(j,i)
150 sq_g11_sgx(j,i) = sqrt(g11_sgx(j,i))
151 sq_g22_sgx(j,i) = sqrt(g22_sgx(j,i))
152 insq_g11_sgx(j,i) = 1.0_dp/sq_g11_sgx(j,i)
158 sq_g22_sgy(j,i) = sqrt(g22_sgy(j,i))
159 sq_g11_sgy(j,i) = sqrt(g11_sgy(j,i))
160 insq_g22_sgy(j,i) = 1.0_dp/sq_g22_sgy(j,i)