47 real(dp),
intent(out) :: dxi, deta
49 integer(i4b) :: i, j, n
54 character(len= 8) :: ch_imax
55 character(len=128) :: fmt4
57 write(ch_imax, fmt=
'(i8)') imax
58 write(fmt4, fmt=
'(a)')
'('//trim(adjustl(ch_imax))//
'(i1),i1)'
62 open(23, iostat=ios, &
63 file=inpath//
'/'//trim(ch_domain_short)//
'/'//zl0_file, &
64 recl=8192, status=
'old')
66 if (ios.ne.0) stop
' topography2: Error when opening the zl0 file!'
68 open(24, iostat=ios, &
69 file=inpath//
'/'//trim(ch_domain_short)//
'/'//mask_present_file, &
70 recl=1024, status=
'old')
72 if (ios.ne.0) stop
' topography2: Error when opening the mask file!'
74 do n=1, 6;
read(23, fmt=
'(a)') ch_dummy;
end do
77 read(23, fmt=*) (zl0(j,i), i=0,imax)
80 do n=1, 6;
read(24, fmt=
'(a)') ch_dummy;
end do
83 read(24, fmt=trim(fmt4)) (maske(j,i), i=0,imax)
86 close(23, status=
'keep')
87 close(24, status=
'keep')
100 if (maske(j,i) <= 1 )
then
106 stop
' topography2: maske(j,i)>=2 not allowed for initial topography!'
109 zs(j,i) = zs(j,i) *1000.0_dp
110 zb(j,i) = zb(j,i) *1000.0_dp
111 zl(j,i) = zl(j,i) *1000.0_dp
112 zl0(j,i) = zl0(j,i)*1000.0_dp
114 xi(i) = xi0 +
real(i,dp)*dxi
115 eta(j) = eta0 +
real(j,dp)*deta
117 call
geo_coord(phi(j,i), lambda(j,i), xi(i), eta(j))
126 dzs_dtau(j,i) = 0.0_dp
127 dzm_dtau(j,i) = 0.0_dp
128 dzb_dtau(j,i) = 0.0_dp
129 dzl_dtau(j,i) = 0.0_dp
130 dh_c_dtau(j,i) = 0.0_dp
131 dh_t_dtau(j,i) = 0.0_dp
150 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.