psmile_range_subgrid_1d_real.F90

Go to the documentation of this file.

!-----------------------------------------------------------------------
! Copyright 2006-2010, NEC Europe Ltd., London, UK.
! All rights reserved. Use is subject to OASIS4 license terms.
!-----------------------------------------------------------------------
!BOP
!
! !ROUTINE: PSMILe_Range_Subgrid_1d_real
!
! !INTERFACE:

      subroutine psmile_range_subgrid_1d_real (          &
                    array, idlow, idhigh, nbr_corners,   &
                           ibeg,  iend,                  &
                    rinter, inter, ierror)
!
! !USES:
!
      use PRISM_constants
!
      use PSMILe, dummy_interface => PSMILe_Range_Subgrid_1d_real

      implicit none
!
! !INPUT PARAMETERS:
!
      integer, Intent (In)             :: idlow, idhigh, nbr_corners

!     Dimensions of "array"

      real, Intent (In)                :: array (idlow:idhigh, nbr_corners)

!     Fully dimensioned Array for which the range of intersection "rinter"
!     should be computed.

      integer, Intent (In)             :: ibeg, iend

!     Definintion of the subgrid

      Real, Intent (In)                :: rinter (2)

!     Specifies the intersection for which the range should be found.
!
! !OUTPUT PARAMETERS:
!
      Integer, Intent (Out)            :: inter (2)

!     Returns the range of the intersection

      integer, Intent (Out)            :: ierror

!     Returns the error code of PSMILe_Range_Subgrid_1d_real;
!             ierror = 0 : No error
!             ierror > 0 : Severe error
!
! !LOCAL VARIABLES
!
      integer            :: i
!
! !DESCRIPTION:
!
! Subroutine "PSMILe_Range_Subgrid_1d_real" computes the range of
! of grid contained in array~(ibeg:iend,~nbr_corners) which contains
! the intersection "rinter".
!
!
! !REVISION HISTORY:
!
!   Date      Programmer   Description
! ----------  ----------   -----------
! 03.06.17    H. Ritzdorf  created
!
!EOP
!----------------------------------------------------------------------
!
! $Id: psmile_range_subgrid_1d_real.F90 2687 2010-10-28 15:15:52Z coquart $
! $Author: coquart $
!
   Character(len=len_cvs_string), save :: mycvs = 
       '$Id: psmile_range_subgrid_1d_real.F90 2687 2010-10-28 15:15:52Z coquart $'
!
!----------------------------------------------------------------------
!
#ifdef VERBOSE
      print *, trim(ch_id), ': PSMILe_Range_Subgrid_1d_real'

      call psmile_flushstd
#endif /* VERBOSE */
!
!  Initialization
!
      ierror = 0
!
      inter (1) = iend + 1
      inter (2) = ibeg - 1
!
!  Get range of the intersection "rinter"
!
!cdir vector
         do i = ibeg, iend
         if (max(array (i,1), array(i,2)) >= rinter (1) .and. &
             min(array (i,1), array(i,2)) <= rinter (2)) then
            inter (1) = min (inter (1), i)
            inter (2) = max (inter (2), i)
         endif
         enddo
!
!===> All done
!
#ifdef VERBOSE
      print *, trim(ch_id), ': PSMILe_Range_Subgrid_1d_real eof', &
                            ': ierror =', ierror, &
                            '; inter  =', inter

      call psmile_flushstd
#endif /* VERBOSE */
!
      end subroutine PSMILe_Range_Subgrid_1d_real