47 real(dp),
intent(out) :: dxi, deta
67 #if (GRID==0 || GRID==1)
77 stop
' topography2: GRID==2 not allowed for this application!'
88 xi(i) = xi0 +
real(i,dp)*dxi
89 eta(j) = eta0 +
real(j,dp)*deta
91 call
geo_coord(phi(j,i), lambda(j,i), xi(i), eta(j))
100 dzs_dtau(j,i) = 0.0_dp
101 dzm_dtau(j,i) = 0.0_dp
102 dzb_dtau(j,i) = 0.0_dp
103 dzl_dtau(j,i) = 0.0_dp
104 dh_c_dtau(j,i) = 0.0_dp
105 dh_t_dtau(j,i) = 0.0_dp
124 area(j,i) = sq_g11_g(j,i)*sq_g22_g(j,i)*dxi*deta
subroutine topograd_2(dxi, deta, n_switch)
Calculation of topography gradients on the staggered grid and on the grid points (the latter by fourt...
Declarations of kind types for SICOPOLIS.
subroutine topograd_1(dxi, deta, n_switch)
Calculation of topography gradients on the staggered grid and on the grid points (the latter by secon...
Declarations of global variables for SICOPOLIS (for the ANT domain).
subroutine topography2(dxi, deta)
Definition of the initial surface and bedrock topography (including gradients) and of the horizontal ...
subroutine geo_coord(phi_val, lambda_val, x_val, y_val)
Computation of longitude lambda and latitude phi for position (x,y) in the numerical domain...
subroutine metric()
Definition of the components g11 and g22 of the metric tensor of the applied coordinates.
Declarations of global variables for SICOPOLIS.