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-2013 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 netcdf

implicit none

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

! 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

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

call check( nf90_open(target_topo_dat_path//'/'//trim(target_topo_dat_name), &
            nf90_nowrite, ncid) )

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

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

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

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

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

call check( nf90_close(ncid) )

!-------- 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

contains

!-------------------------------------------------------------------------------
!> NetCDF error capturing.
!<------------------------------------------------------------------------------
   subroutine check(status)
      integer(i4b), intent (in) :: status
      if(status /= nf90_noerr) then 
         write(6,'(a)') trim(nf90_strerror(status))
         stop ' read_target_topo_nc: Stopped due to netcdf error!'
      end if
   end subroutine check  

end subroutine read_target_topo_nc
!