read_target_topo_nc.F90

Go to the documentation of this file.

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
!
!  Subroutine :  r e a d _ t a r g e t _ t o p o _ n c
!
!> @file
!!
!! Reading of the target-topography file (in NetCDF format).
!!
!! @section Copyright
!!
!! Copyright 2009-2016 Ralf Greve
!!
!! @section License
!!
!! This file is part of SICOPOLIS.
!!
!! SICOPOLIS is free software: you can redistribute it and/or modify
!! it under the terms of the GNU General Public License as published by
!! the Free Software Foundation, either version 3 of the License, or
!! (at your option) any later version.
!!
!! SICOPOLIS is distributed in the hope that it will be useful,
!! but WITHOUT ANY WARRANTY; without even the implied warranty of
!! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!! GNU General Public License for more details.
!!
!! You should have received a copy of the GNU General Public License
!! along with SICOPOLIS.  If not, see <http://www.gnu.org/licenses/>.
!<
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

!-------------------------------------------------------------------------------
!> Reading of the target-topography file (in NetCDF format).
!<------------------------------------------------------------------------------
subroutine read_target_topo_nc(target_topo_dat_name)

use sico_types
use sico_variables
use sico_vars
use netcdf
use nc_check

implicit none

character(len=100), intent(in) :: target_topo_dat_name
character(len=256) :: filename_with_path

! Return variables
! (defined as global variables in module sico_variables):
!
!    maske_target(j,i),
!    zs_target(j,i), zb_target(j,i), zl_target(j,i), H_target(j,i)

integer(i4b) :: i, j

integer(i2b), dimension(0:JMAX,0:IMAX) :: maske_conv
real(sp), dimension(0:JMAX,0:IMAX) :: zs_conv, zb_conv, zl_conv, h_conv

integer(i4b) :: ncid, ncv
!     ncid:      ID of the output file
!     ncv:       Variable ID

character(len=64), parameter :: thisroutine = 'read_target_topo_nc'

!-------- Read target-topography data
!         (from a time-slice file of a previous simulation) --------

filename_with_path = trim(target_topo_dat_path)//'/'//trim(target_topo_dat_name)

call check( nf90_open(trim(filename_with_path), nf90_nowrite, ncid), &
            thisroutine )

call check( nf90_inq_varid(ncid, 'maske', ncv),  thisroutine )
call check( nf90_get_var(ncid, ncv, maske_conv), thisroutine )

call check( nf90_inq_varid(ncid, 'zs', ncv),  thisroutine )
call check( nf90_get_var(ncid, ncv, zs_conv), thisroutine )

call check( nf90_inq_varid(ncid, 'zb', ncv),  thisroutine )
call check( nf90_get_var(ncid, ncv, zb_conv), thisroutine )

call check( nf90_inq_varid(ncid, 'zl', ncv),  thisroutine )
call check( nf90_get_var(ncid, ncv, zl_conv), thisroutine )

call check( nf90_inq_varid(ncid, 'H', ncv),  thisroutine )
call check( nf90_get_var(ncid, ncv, h_conv), thisroutine )

call check( nf90_close(ncid), thisroutine )

!-------- Convert data to double precision --------

do i=0, imax
do j=0, jmax
   maske_target(j,i) = maske_conv(j,i)
   zs_target(j,i)    = real(zs_conv(j,i),dp)
   zb_target(j,i)    = real(zb_conv(j,i),dp)
   zl_target(j,i)    = real(zl_conv(j,i),dp)
   h_target(j,i)     = real(H_conv(j,i), dp)
end do
end do

end subroutine read_target_topo_nc
!