65 real(dp) :: g11_g(0:jmax,0:imax), g22_g(0:jmax,0:imax), &
66 g11_sgx(0:jmax,0:imax), g11_sgy(0:jmax,0:imax), &
67 g22_sgx(0:jmax,0:imax), g22_sgy(0:jmax,0:imax)
69 #if (GRID==0) /* Stereographic projection (distortion neglected) */ 81 #elif (GRID==1) /* Stereographic projection */ 84 k = (cos(0.25_dp*
pi-0.5_dp*
phi0))**2
86 k = (cos(0.25_dp*
pi+0.5_dp*
phi0))**2
88 stop
' >>> metric: PHI0 must be different from zero!' 104 g11_sgx(j,i), g22_sgx(j,i))
111 g11_sgy(j,i), g22_sgy(j,i))
115 #elif (GRID==2) /* Geographical coordinates */ 181 real(dp),
intent(in) :: x_val, y_val
182 real(dp),
intent(in) :: K
183 real(dp),
intent(out) :: g11_r, g22_r
185 g11_r = 1.0_dp / ( k**2*(1.0_dp+(x_val**2+y_val**2)/(2.0_dp*
r*k)**2)**2 )
200 real(dp),
intent(in) :: phi_val
201 real(dp),
intent(out) :: g11_r, g22_r
203 g11_r =
r**2*(cos(phi_val))**2
real(dp), dimension(0:jmax, 0:imax) insq_g22_sgy
insq_g22_sgy(j,i): Inverse square root of g22, at (i,j+1/2)
real(dp), parameter eps
eps: Small number
subroutine metric_stereo(x_val, y_val, K, g11_r, g22_r)
Components g11 and g22 of the metric tensor for the stereographical projection.
real(dp), dimension(0:jmax, 0:imax) insq_g11_sgx
insq_g11_sgx(j,i): Inverse square root of g11, at (i+1/2,j)
real(dp), dimension(0:jmax, 0:imax) insq_g22_g
insq_g22_g(j,i): Inverse square root of g22 on grid point (i,j)
real(dp), dimension(0:jmax, 0:imax) sq_g22_sgx
sq_g22_sgx(j,i): Square root of g22, at (i+1/2,j)
subroutine metric_geogr(phi_val, g11_r, g22_r)
Components g11 and g22 of the metric tensor for geographical coordinates.
real(dp), dimension(0:jmax, 0:imax) sq_g22_sgy
sq_g22_sgy(j,i): Square root of g22, at (i,j+1/2)
real(dp) r
R: Radius of the planet.
Declarations of kind types for SICOPOLIS.
real(dp), dimension(0:jmax, 0:imax) sq_g11_sgy
sq_g11_sgy(j,i): Square root of g11, at (i,j+1/2)
Definition of the components g11 and g22 of the metric tensor of the applied numerical coordinates...
real(dp) phi0
PHI0: Standard parallel of the stereographic projection.
real(dp), dimension(0:jmax, 0:imax) sq_g11_sgx
sq_g11_sgx(j,i): Square root of g11, at (i+1/2,j)
real(dp), parameter pi
pi: Constant pi
real(dp), dimension(0:jmax, 0:imax) sq_g11_g
sq_g11_g(j,i): Square root of the coefficient g11 of the metric tensor on grid point (i...
subroutine, public metric()
Main routine of module metric_m: Definition of the components g11 and g22 of the metric tensor of the...
real(dp), dimension(0:jmax) eta
eta(j): Coordinate eta (= y) of grid point j
real(dp), dimension(0:jmax, 0:imax) insq_g11_g
insq_g11_g(j,i): Inverse square root of g11 on grid point (i,j)
real(dp), dimension(0:imax) xi
xi(i): Coordinate xi (= x) of grid point i
Declarations of global variables for SICOPOLIS.
real(dp), dimension(0:jmax, 0:imax) sq_g22_g
sq_g22_g(j,i): Square root of the coefficient g22 of the metric tensor on grid point (i...