49 real(dp),
intent(out) :: dxi, deta
51 integer(i4b) :: i, j, n
52 integer(i4b) :: ios, n_dummy
57 character(len= 8) :: ch_imax
58 character(len=128) :: fmt4
60 write(ch_imax, fmt=
'(i8)') imax
61 write(fmt4, fmt=
'(a)')
'('//trim(adjustl(ch_imax))//
'(i1),i1)'
67 open(24, iostat=ios, &
68 file=inpath//
'/'//trim(ch_domain_short)//
'/'//mask_present_file, &
69 recl=1024, status=
'old')
71 if (ios /= 0) stop
' topography1: Error when opening the mask file!'
73 do n=1, 6;
read(24, fmt=
'(a)') ch_dummy;
end do
76 read(24, fmt=trim(fmt4)) (maske(j,i), i=0,imax)
79 close(24, status=
'keep')
83 #if (GRID==0 || GRID==1)
93 stop
' topography1: GRID==2 not allowed for this application!'
104 if (maske(j,i) <= 1)
then
106 zs(j,i) = zs(j,i) + 2.0_dp*h_min
109 xi(i) = xi0 +
real(i,dp)*dxi
110 eta(j) = eta0 +
real(j,dp)*deta
112 call
geo_coord(phi(j,i), lambda(j,i), xi(i), eta(j))
118 h_c(j,i) = zs(j,i)-zm(j,i)
121 dzs_dtau(j,i) = 0.0_dp
122 dzm_dtau(j,i) = 0.0_dp
123 dzb_dtau(j,i) = 0.0_dp
124 dzl_dtau(j,i) = 0.0_dp
125 dh_c_dtau(j,i) = 0.0_dp
126 dh_t_dtau(j,i) = 0.0_dp
145 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...
subroutine topography1(dxi, deta)
Definition of the initial surface and bedrock topography (including gradients) and of the horizontal ...
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 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.