read_topd_file.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!-----------------------------------------------------------------
!     #######################
      SUBROUTINE read_topd_file(HPROGRAM,HFILE,HFORM,KNPT,PTOPD_READ)
!     #######################
!
!!****  *READ_TOPD_FILE*  
!!
!!    PURPOSE
!!    -------
!     This routine aims at reading topographic files
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!    
!!    
!!
!!      
!!    REFERENCE
!!    ---------
!!
!!    
!!      
!!    AUTHOR
!!    ------
!!
!!      B. Vincendon    * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      Original   11/2006
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modi_get_luout
USE modi_open_file
USE modi_close_file
!
USE modd_topd_par, ONLY : nunit
USE modd_topodyn, ONLY : npmax
USE modd_surf_par,  ONLY : xundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
 CHARACTER(LEN=*),  INTENT(IN)  :: hprogram    !
 CHARACTER(LEN=*),  INTENT(IN)  :: hfile       ! File to be read
 CHARACTER(LEN=*),  INTENT(IN)  :: hform       ! Format of the file to be read
INTEGER,           INTENT(IN)  :: knpt        ! Number of points in the catchment
REAL, DIMENSION(:), INTENT(OUT) :: ptopd_read ! Topographic parameter read on file
!
!*      0.2    declarations of local variables
!
INTEGER                   :: jj,ji  ! loop control
INTEGER                   :: iluout ! Unit of the files
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_TOPD_FILE',0,zhook_handle)
!
!*       0.2    preparing file openning
!               ----------------------
!
 CALL get_luout(hprogram,iluout)
!
 CALL open_file(hprogram,nunit,hfile,hform,haction='READ')
!PTOPD_READ(:)=XUNDEF
!
DO jj=1,13
  READ(nunit,*) 
ENDDO
!
DO ji=1,knpt
  !
  READ(nunit,*,end=110) ptopd_read(ji)
  !
  IF (ptopd_read(ji).LT.0.0) ptopd_read(ji) = xundef
  !
ENDDO
!
110 CALL close_file(hprogram,nunit)
!
IF (lhook) CALL dr_hook('READ_TOPD_FILE',1,zhook_handle)
!
END SUBROUTINE read_topd_file