7 !! Definition of the initial surface and bedrock topography
8 !! (including gradients) and of the horizontal grid spacings dxi, deta.
9 !! For an initial topography with a thin ice layer (thickness 2*H_MIN)
10 !! everywhere on the land area.
14 !! Copyright 2009-2013 Ralf Greve
18 !! This file is part of SICOPOLIS.
20 !! SICOPOLIS is free software: you can redistribute it and/or modify
21 !! it under the terms of the GNU General Public License as published by
22 !! the Free Software Foundation, either version 3 of the License, or
23 !! (at your option) any later version.
25 !! SICOPOLIS is distributed in the hope that it will be useful,
26 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
27 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
28 !! GNU General Public License for more details.
30 !! You should have received a copy of the GNU General Public License
31 !! along with SICOPOLIS. If not, see <http://www.gnu.org/licenses/>.
33 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
36 !> Definition of the initial surface and bedrock topography
37 !! (including gradients) and of the horizontal grid spacings dxi, deta.
38 !! For an initial topography with a thin ice layer (thickness 2*H_MIN)
39 !! everywhere on the land area.
40 !<------------------------------------------------------------------------------
47 integer(i4b) :: i, j, n
48 integer(i4b) :: ios, n_dummy
58 open(24, iostat=ios, &
59 file=inpath//
'/heino/'//mask_present_file, &
60 recl=1024, status=
'old')
62 if (ios /= 0) stop
' topography1: Error when opening the mask file!'
64 do n=1, 6;
read(24,
'(a)') ch_dummy;
end do
67 read(24,2300) (maske(j,i), i=0,imax)
70 close(24, status=
'keep')
72 2300 format(imax(i1),i1)
76 #if (GRID==0 || GRID==1)
86 stop
' topography1: GRID==2 not allowed for this application!'
97 if (maske(j,i) <= 1)
then
99 zs(j,i) = zs(j,i) + 2.0_dp*h_min
102 xi(i) = xi0 +
real(i,dp)*dxi
103 eta(j) = eta0 +
real(j,dp)*deta
105 call
geo_coord(phi(j,i), lambda(j,i), xi(i), eta(j))
110 h_c(j,i) = zs(j,i)-zm(j,i)
113 dzs_dtau(j,i) = 0.0_dp
114 dzm_dtau(j,i) = 0.0_dp
115 dzb_dtau(j,i) = 0.0_dp
116 dzl_dtau(j,i) = 0.0_dp
117 dh_c_dtau(j,i) = 0.0_dp
118 dh_t_dtau(j,i) = 0.0_dp
139 area(j,i) = sq_g11_g(j,i)*sq_g22_g(j,i)*dxi*deta