read_connex_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_connex_file(HPROGRAM,HFILE,HFORM,KNMC,PCONN,KLINE)
!     #######################
!
!!****  *READ_CONNEX_FILE*  
!!
!!    PURPOSE
!!    -------
!     This routine aims at reading connexion file
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!    
!!    
!!
!!      
!!    REFERENCE
!!    ---------
!!
!!    
!!      
!!    AUTHOR
!!    ------
!!
!!      B. Vincendon    * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      Original   11/2006
!!                 03/2014 (B. Vincendon) format correction
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_topd_par, ONLY : nunit
USE modd_topodyn, ONLY : npmax
USE modd_surf_par,  ONLY : xundef
!
USE modi_get_luout
USE modi_open_file
USE modi_close_file
!
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)  :: knmc       ! Number of pixels in the catchment
REAL, DIMENSION(:,:),INTENT(OUT)   :: pconn    ! pixels topographic slope/length flow
INTEGER, DIMENSION(:),INTENT(OUT)  :: kline    ! second index of the pixel in the array PCONN
!
!*      0.2    declarations of local variables
!
!
INTEGER                   :: jj          ! loop control 
INTEGER                   :: iluout      ! Unit of the files
INTEGER                   :: iindex      ! index of the pixel in the topo domain
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_CONNEX_FILE',0,zhook_handle)
!
!*       0.2    preparing file openning
!               ----------------------
 CALL get_luout(hprogram,iluout)
!
 CALL open_file(hprogram,nunit,hfile,hform,haction='READ')
!
WRITE(iluout,*) 'Open ',hfile,'debut'
!
DO jj=1,7
  READ(nunit,*) 
ENDDO
!
DO jj=1,knmc
  !
  READ(nunit,*,end=120) pconn(jj,:)
  iindex = int(pconn(jj,1))
  kline(iindex) = jj
  !
ENDDO
!   
120   CALL close_file(hprogram,nunit)
!
IF (lhook) CALL dr_hook('READ_CONNEX_FILE',1,zhook_handle)
!
END SUBROUTINE read_connex_file