read_topd_header_dtm.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_header_dtm(HPROGRAM,HFILE,HFORM,PX0,PY0,KNXC,KNYC,PNUL,PDXT)
!     #######################
!
!!****  *READ_TOPD_HEADER*  
!!
!!    PURPOSE
!!    -------
!     This routine aims at reading topographic files header 
!     for a given file (then for a cathment)
!
!!**  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 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
REAL,              INTENT(OUT) :: px0        ! abcissa of bottom-left pixel   
REAL,              INTENT(OUT) :: py0        ! ordinate of bottom-left pixel   
INTEGER,           INTENT(OUT) :: knxc       ! number of topographixc grid points along abcissa axis   
INTEGER,           INTENT(OUT) :: knyc       ! number of topographixc grid points along ordinate axis   
REAL,              INTENT(OUT) :: pnul       ! undifined value in topographic files
REAL,              INTENT(OUT) :: pdxt       ! catchment rid mesh size   
!
!*      0.2    declarations of local variables
!
INTEGER                   :: jj          ! loop control 
INTEGER                   :: iluout      ! Unit of the files
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_TOPD_HEADER_DTM',0,zhook_handle)
!
!*       0.2    preparing file openning
!               ----------------------
 CALL get_luout(hprogram,iluout)
!
WRITE(iluout,*) 'Open ',hfile,'header'
!
 CALL open_file(hprogram,nunit,hfile,hform,haction='READ')
!
DO jj=1,5
  READ(nunit,*) 
ENDDO
!
READ(nunit,*) px0
READ(nunit,*) py0
READ(nunit,*) knxc
READ(nunit,*) knyc
READ(nunit,*) pnul
READ(nunit,*) pdxt
READ(nunit,*)
READ(nunit,*)
!
 CALL close_file(hprogram,nunit)
!
IF (lhook) CALL dr_hook('READ_TOPD_HEADER_DTM',1,zhook_handle)
!
END SUBROUTINE read_topd_header_dtm