read_from_surfex_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_from_surfex_file (DTCO, U, &
                                  HFTYP,HFNAM,HMASK,HSCHEME,PFIELD,HNAM)
!
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modi_init_io_surf_n
USE modi_end_io_surf_n
USE modi_read_surf
!
#ifdef SFX_ASC
USE modi_set_surfex_file_name_asc
#endif
#ifdef SFX_FA
USE modi_set_surfex_file_name_fa
#endif
#ifdef SFX_LFI
USE modi_set_surfex_file_name_lfi
#endif
#ifdef SFX_NC
USE modi_set_surfex_file_name_nc
#endif
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(surf_atm_t), INTENT(INOUT) :: U

!
 CHARACTER(LEN=6), INTENT(IN) :: HFTYP
 CHARACTER(LEN=28), INTENT(IN) :: HFNAM
 CHARACTER(LEN=6), INTENT(IN) :: HMASK
 CHARACTER(LEN=6), INTENT(IN) :: HSCHEME
REAL, DIMENSION(:), INTENT(OUT) :: PFIELD
 CHARACTER(LEN=*), OPTIONAL, INTENT(IN) :: HNAM
!
 CHARACTER(LEN=28)     :: HNAME
 CHARACTER(LEN=12)     :: YRECFM
INTEGER               :: IRESP, I0
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE

IF (lhook) CALL dr_hook('READ_FROM_SURFEX_FILE',0,zhook_handle)

hname=adjustl(adjustr(hfnam))

IF (PRESENT(hnam)) THEN
  yrecfm=hnam
ELSE
  i0=index(hname,'.',back=.true.)
  yrecfm=hname(1:i0-1)
ENDIF

IF (hftyp=='ASCII ') THEN
#ifdef SFX_ASC
  CALL set_surfex_file_name_asc(hname_in=hname)
#endif
ELSEIF (hftyp=='FA    ') THEN
#ifdef SFX_FA
  CALL set_surfex_file_name_fa(hname_in=hname)
#endif
ELSEIF (hftyp=='LFI   ') THEN
#ifdef SFX_LFI
  CALL set_surfex_file_name_lfi(hname_in=hname)
#endif
ELSEIF (hftyp=='NC    ') THEN
#ifdef SFX_NC
  CALL set_surfex_file_name_nc(hname_in=hname)
#endif
ENDIF
 CALL init_io_surf_n(dtco, u, &
                     hftyp,hmask,hscheme,'READ ')
 CALL read_surf(&
                hftyp,yrecfm,pfield,iresp)
 CALL end_io_surf_n(hftyp)

IF (lhook) CALL dr_hook('READ_FROM_SURFEX_FILE',1,zhook_handle)

END SUBROUTINE read_from_surfex_file