7 !! Definition of the initial surface and bedrock topography
8 !! (including gradients) and of the horizontal grid spacings dxi, deta.
9 !! For ice-free initial topography with relaxed lithosphere.
13 !! Copyright 2009-2013 Ralf Greve
17 !! This file is part of SICOPOLIS.
19 !! SICOPOLIS is free software: you can redistribute it and/or modify
20 !! it under the terms of the GNU General Public License as published by
21 !! the Free Software Foundation, either version 3 of the License, or
22 !! (at your option) any later version.
24 !! SICOPOLIS is distributed in the hope that it will be useful,
25 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
26 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
27 !! GNU General Public License for more details.
29 !! You should have received a copy of the GNU General Public License
30 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
32 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
35 !> Definition of the initial surface and bedrock topography
36 !! (including gradients) and of the horizontal grid spacings dxi, deta.
37 !! For ice-free initial topography with relaxed lithosphere.
38 !<------------------------------------------------------------------------------
64 #if (GRID==0 || GRID==1)
74 stop
' topography2: GRID==2 not allowed for this application!'
85 xi(i) = xi0 +
real(i,dp)*dxi
86 eta(j) = eta0 +
real(j,dp)*deta
88 call
geo_coord(phi(j,i), lambda(j,i), xi(i), eta(j))
96 dzs_dtau(j,i) = 0.0_dp
97 dzm_dtau(j,i) = 0.0_dp
98 dzb_dtau(j,i) = 0.0_dp
99 dzl_dtau(j,i) = 0.0_dp
100 dh_c_dtau(j,i) = 0.0_dp
101 dh_t_dtau(j,i) = 0.0_dp
122 area(j,i) = sq_g11_g(j,i)*sq_g22_g(j,i)*dxi*deta