45 character(len=100),
intent(in) :: target_topo_dat_name
46 character(len=256) :: filename_with_path
56 integer(i2b),
dimension(0:JMAX,0:IMAX) :: maske_conv
57 real(sp),
dimension(0:JMAX,0:IMAX) :: zs_conv, zb_conv, zl_conv, h_conv
59 integer(i4b) :: ncid, ncv
63 character(len=64),
parameter :: thisroutine =
'read_target_topo_nc'
68 filename_with_path = trim(target_topo_dat_path)//
'/'//trim(target_topo_dat_name)
70 call
check( nf90_open(trim(filename_with_path), nf90_nowrite, ncid), &
73 call
check( nf90_inq_varid(ncid,
'maske', ncv), thisroutine )
74 call
check( nf90_get_var(ncid, ncv, maske_conv), thisroutine )
76 call
check( nf90_inq_varid(ncid,
'zs', ncv), thisroutine )
77 call
check( nf90_get_var(ncid, ncv, zs_conv), thisroutine )
79 call
check( nf90_inq_varid(ncid,
'zb', ncv), thisroutine )
80 call
check( nf90_get_var(ncid, ncv, zb_conv), thisroutine )
82 call
check( nf90_inq_varid(ncid,
'zl', ncv), thisroutine )
83 call
check( nf90_get_var(ncid, ncv, zl_conv), thisroutine )
85 call
check( nf90_inq_varid(ncid,
'H', ncv), thisroutine )
86 call
check( nf90_get_var(ncid, ncv, h_conv), thisroutine )
88 call
check( nf90_close(ncid), thisroutine )
94 maske_target(j,i) = maske_conv(j,i)
95 zs_target(j,i) =
real(zs_conv(j,i),dp)
96 zb_target(j,i) =
real(zb_conv(j,i),dp)
97 zl_target(j,i) =
real(zl_conv(j,i),dp)
98 h_target(j,i) =
real(H_conv(j,i), dp)
Declarations of kind types for SICOPOLIS.
subroutine check(status, ch_calling_routine)
NetCDF error capturing.
Declarations of global variables for SICOPOLIS (for the ANT domain).
subroutine read_target_topo_nc(target_topo_dat_name)
Reading of the target-topography file (in NetCDF format).
Declarations of global variables for SICOPOLIS.