read_surf.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.
!##################
MODULE modi_read_surf
!##################
!
!!****  *READX0* - routine to read a real scalar
!!
!!    PURPOSE
!!    -------
!
!       The purpose of READX0 is
!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!     
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!
!!      S.Malardel      *METEO-FRANCE*
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      original                                                     01/08/03
!!      J.Escobar      10/06/2013: replace DOUBLE PRECISION by REAL to handle problem for promotion of real on IBM SP
!----------------------------------------------------------------------------
!
  INTERFACE read_surf
!
     SUBROUTINE read_surfx0 (&
                             HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
REAL, INTENT(OUT) :: PFIELD            ! real scalar to be read
INTEGER,INTENT(OUT) :: KRESP             ! KRESP  : return-code if a problem appears 
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT  ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surfx0
!
     SUBROUTINE read_surfx1 (&
                             HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
REAL, DIMENSION(:), INTENT(OUT) ::PFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP            ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR       ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
END SUBROUTINE read_surfx1
!
     SUBROUTINE read_surfx2 (&
                             HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM    ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC        ! name of the article to be read
REAL, DIMENSION(:,:), INTENT(OUT) :: PFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP               ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
END SUBROUTINE read_surfx2
!
!RJ: interface to READ_SURFX2COV moved out
!
     SUBROUTINE read_surfx3 (&
                             HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM      ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC          ! name of the article to be read
REAL, DIMENSION(:,:,:), INTENT(OUT) :: PFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP                 ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
END SUBROUTINE read_surfx3
!
     SUBROUTINE read_surfn0 (&
                             HPROGRAM,HREC,KFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
INTEGER, INTENT(OUT) :: KFIELD           ! integer to be read  
INTEGER, INTENT(OUT) :: KRESP            ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surfn0
!
     SUBROUTINE read_surfn1 (&
                             HPROGRAM,HREC,KFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM     ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC         ! name of the article to be read
INTEGER, DIMENSION(:), INTENT(OUT) :: KFIELD ! integer to be read  
INTEGER, INTENT(OUT) :: KRESP                ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
END SUBROUTINE read_surfn1
!
     SUBROUTINE read_surfc0 (&
                             HPROGRAM,HREC,HFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM   ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC       ! name of the article to be read
 CHARACTER(LEN=*), INTENT(OUT) :: HFIELD    ! caracter to be read  
INTEGER, INTENT(OUT) :: KRESP              ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surfc0
!
      SUBROUTINE read_surfl0 (&
                              HPROGRAM,HREC,OFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
LOGICAL, INTENT(OUT)         :: OFIELD   ! array containing the data field
INTEGER, INTENT(OUT) :: KRESP            ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surfl0
!
      SUBROUTINE read_surfl1 (&
                              HPROGRAM,HREC,OFIELD,KRESP,HCOMMENT,HDIR)
!
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM     ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC         ! name of the article to be read
LOGICAL, DIMENSION(:), INTENT(OUT) :: OFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP                ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
END SUBROUTINE read_surfl1
!
      SUBROUTINE read_surft0 (&
                              HPROGRAM,HREC,TFIELD,KRESP,HCOMMENT,HDIR)
!
!
USE modd_type_date_surf
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM  ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC      ! name of the article to be read
TYPE(date_time), INTENT(INOUT) ::TFIELD   ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP             ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surft0
!
      SUBROUTINE read_surft1 (&
                              HPROGRAM,HREC,TFIELD,KRESP,HCOMMENT,HDIR)
!
!
USE modd_type_date_surf
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
type(date_time), DIMENSION(:), INTENT(INOUT) :: tfield ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surft1
!
      SUBROUTINE read_surft2 (&
                              HPROGRAM,HREC,TFIELD,KRESP,HCOMMENT,HDIR)
!
!
USE modd_type_date_surf
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) ::HREC      ! name of the article to be read
type(date_time), DIMENSION(:,:), INTENT(INOUT)::tfield ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
END SUBROUTINE read_surft2
!
END INTERFACE
!
END MODULE modi_read_surf
!
!     #############################################################
      SUBROUTINE read_surfx0 (&
                             HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
#ifdef SFX_OL
USE mode_read_surf_ol, ONLY: read_surf0_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surf0_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surf0_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surf0_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surf0_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfx0_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfx0_aro
#endif
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN)  :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN)  :: HREC     ! name of the article to be read
REAL, INTENT(OUT) :: PFIELD               ! the real scalar to be read  
INTEGER, INTENT(OUT) :: KRESP             ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
REAL   :: XTIME0
INTEGER            :: INFOMPI
REAL               :: ZWORK
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX0',0,zhook_handle)
!
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfx0_mnh(yrec,pfield,kresp,ycomment)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfx0_aro(yrec,pfield,kresp,ycomment)
#endif
ENDIF
!
IF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' .OR. &
    hprogram=='NC    ') THEN 
  !
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    !
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,pfield,kresp,ycomment)
#endif
    ELSEIF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,pfield,kresp,ycomment)
#endif
    ELSEIF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,pfield,kresp,ycomment)
#endif
    ELSEIF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,pfield,kresp,ycomment)
#endif
    ELSEIF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,pfield,kresp,ycomment)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()
    CALL mpi_bcast(pfield,kind(pfield)/4,mpi_real,npio,ncomm,infompi)
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)
  ENDIF
#endif
  !
ENDIF
!    
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX0',1,zhook_handle)
!
END SUBROUTINE read_surfx0
!
!     #############################################################
      SUBROUTINE read_surfx1 (HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surf_par, ONLY : xundef
USE modd_assim, ONLY : lassim, lread_all, cassim_isba
!
#ifdef SFX_OL
USE mode_read_surf_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surfn_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surfn_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surfn_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surfx_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfx1_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfx1_aro
#endif
!
USE modi_get_iok_assim
!
IMPLICIT NONE
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM  ! calling program
!
!
 CHARACTER(LEN=*), INTENT(IN) :: HREC      ! name of the article to be read
REAL, DIMENSION(:), INTENT(OUT) :: PFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP             ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: IL, IOK
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX1',0,zhook_handle)
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
il = SIZE(pfield)
!
IF (lassim .AND. cassim_isba/="OI   " .AND. .NOT.lread_all) THEN
  !
  CALL get_iok_assim(yrec,iok)
  !
  IF (iok==0) THEN
    pfield(:) = xundef
    IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX1',1,zhook_handle)
    RETURN
  ENDIF
  !
ENDIF
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfx1_mnh(yrec,il,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfx1_aro(yrec,il,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
  CALL read_surfx1n1_ol(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
  CALL read_surfn_lfi(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
  CALL read_surfn_nc(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
  CALL read_surfn_asc(&
                      yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
  CALL read_surfx_fa(yrec,il,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX1',1,zhook_handle)
!
END SUBROUTINE read_surfx1
!
!     #############################################################
      SUBROUTINE read_surfx2 (HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surf_par, ONLY : xundef
USE modd_assim, ONLY : lassim, lread_all, cassim_isba
!
#ifdef SFX_OL
USE mode_read_surf_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surfn_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surfn_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surfn_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surfx_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfx2_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfx2_aro
#endif
!
USE modi_get_iok_assim
!
IMPLICIT NONE
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM    ! calling program
!
!
 CHARACTER(LEN=*), INTENT(IN) :: HREC        ! name of the article to be read
REAL, DIMENSION(:,:), INTENT(OUT) :: PFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP               ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: IL1, IL2, IOK
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX2',0,zhook_handle)
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
il1 = SIZE(pfield,1)
il2 = SIZE(pfield,2)
!
IF (lassim .AND. cassim_isba/="OI   "  .AND. .NOT.lread_all) THEN
  !
  CALL get_iok_assim(yrec,iok)
  !
  IF (iok==0) THEN
    pfield(:,:) = xundef
    IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX2',1,zhook_handle)
    RETURN
  ENDIF
  !
ENDIF
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfx2_mnh(yrec,il1,il2,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfx2_aro(yrec,il1,il2,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
  CALL read_surfl1x2_ol(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
  CALL read_surfn_lfi(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
  CALL read_surfn_nc(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
  CALL read_surfn_asc(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
  CALL read_surfx_fa(yrec,il1,il2,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX2',1,zhook_handle)
!
END SUBROUTINE read_surfx2
!
!     #############################################################
      SUBROUTINE read_surfx3 (&
                             HPROGRAM,HREC,PFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
#ifdef SFX_OL
USE mode_read_surf_ol
#endif
!
IMPLICIT NONE
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM      ! calling program
!
!
 CHARACTER(LEN=*), INTENT(IN) :: HREC          ! name of the article to be read
REAL, DIMENSION(:,:,:), INTENT(OUT) :: PFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP                 ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: IL1, IL2, IL3
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX3',0,zhook_handle)
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
il1 = SIZE(pfield,1)
il2 = SIZE(pfield,2)
il3 = SIZE(pfield,3)
!
!plmIF (HPROGRAM=='MESONH') THEN
!plm  CALL READ_SURFX3_MNH(YREC,IL1,IL2,PFIELD,KRESP,YCOMMENT,YDIR)
!plmENDIF
!
IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
  CALL read_surfl1x2_ol(yrec,pfield,kresp,ycomment,ydir)
#endif
ENDIF
!
!plmIF (HPROGRAM=='ASCII ') THEN
!plm  CALL READ_SURFX3_ASC(YREC,IL1,IL2,PFIELD,KRESP,YCOMMENT,YDIR)
!plmENDIF
!
!plmIF (HPROGRAM=='AROME ') THEN
!plm  CALL READ_SURFX3_ARO(YREC,IL1,IL2,PFIELD,KRESP,YCOMMENT,YDIR)
!plmENDIF
!
!plmIF (HPROGRAM=='FA    ') THEN
!plm  CALL READ_SURFX3_FA(YREC,IL1,IL2,PFIELD,KRESP,YCOMMENT,YDIR)
!plmENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFX3',1,zhook_handle)
!
END SUBROUTINE read_surfx3
!
!     #############################################################
      SUBROUTINE read_surfn0 (&
                             HPROGRAM,HREC,KFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
#ifdef SFX_OL
USE mode_read_surf_ol, ONLY: read_surf0_ol
#endif
#ifdef SFX_LFI
USE modd_io_surf_lfi, ONLY : lmnh_compatible, niu,nib,nie,nju,njb,nje
USE mode_read_surf_lfi, ONLY: read_surf0_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surf0_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surf0_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surf0_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfn0_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfn0_aro
#endif
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
INTEGER, INTENT(OUT) :: KFIELD           ! the integer to be read  
INTEGER, INTENT(OUT) :: KRESP            ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
REAL   :: XTIME0
INTEGER, DIMENSION(6) :: IDIMS
INTEGER            :: INFOMPI
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFN0',0,zhook_handle)
!
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfn0_mnh(yrec,kfield,kresp,ycomment)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfn0_aro(yrec,kfield,kresp,ycomment)
#endif
ENDIF
!
IF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' .OR. &
    hprogram=='NC    ' ) THEN 
  !
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    !
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,kfield,kresp,ycomment)
#endif
    ENDIF
    !
    IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,kfield,kresp,ycomment)
#endif
    ENDIF
    !
    IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,kfield,kresp,ycomment)
#endif
    ENDIF    
    !
    IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,kfield,kresp,ycomment)
#endif
    ENDIF
    !
    IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,kfield,kresp,ycomment)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN          
    xtime0 = mpi_wtime()  
    CALL mpi_bcast(kfield,kind(kfield)/4,mpi_integer,npio,ncomm,infompi)
#ifdef SFX_LFI
    IF (hprogram=='LFI   '.AND.trim(hrec)=="DIM_FULL") THEN
      CALL mpi_bcast(lmnh_compatible,1,mpi_logical,npio,ncomm,infompi)
      IF (nrank==npio) THEN
        idims(1) = niu
        idims(2) = nib
        idims(3) = nie
        idims(4) = nju
        idims(5) = njb
        idims(6) = nje
      ENDIF
      CALL mpi_bcast(idims,SIZE(idims)*kind(idims)/4,mpi_integer,npio,ncomm,infompi)
      IF (nrank/=npio) THEN
        niu = idims(1)
        nib = idims(2)
        nie = idims(3)
        nju = idims(4)
        njb = idims(5)
        nje = idims(6)
      ENDIF
    ENDIF
#endif    
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF    
#endif
  !
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
!
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFN0',1,zhook_handle)
!
END SUBROUTINE read_surfn0
!
!     #############################################################
      SUBROUTINE read_surfn1 (&
                             HPROGRAM,HREC,KFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
#ifdef SFX_OL
USE mode_read_surf_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surfn_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surfn_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surfn_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surfn_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfn1_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfn1_aro
#endif
!
IMPLICIT NONE
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM     ! calling program
!
!
 CHARACTER(LEN=*), INTENT(IN) :: HREC         ! name of the article to be read
INTEGER, DIMENSION(:), INTENT(OUT) :: KFIELD ! the integer to be read  
INTEGER, INTENT(OUT) :: KRESP                ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: IL
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFN1',0,zhook_handle)
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
il = SIZE(kfield,1)
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfn1_mnh(yrec,il,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfn1_aro(yrec,il,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
  CALL read_surfx1n1_ol(yrec,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
  CALL read_surfn_lfi(yrec,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
  CALL read_surfn_nc(yrec,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
  CALL read_surfn_asc(&
                      yrec,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
  CALL read_surfn_fa(&
                     yrec,il,kfield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFN1',1,zhook_handle)
!
END SUBROUTINE read_surfn1

!     #############################################################
      SUBROUTINE read_surfc0 (&
                             HPROGRAM,HREC,HFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
#ifdef SFX_OL
USE mode_read_surf_ol, ONLY: read_surf0_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surf0_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surf0_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surf0_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surf0_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfc0_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfc0_aro
#endif
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN)  :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN)  :: HREC     ! name of the article to be read
 CHARACTER(LEN=*), INTENT(OUT) :: HFIELD   ! the integer to be read  
INTEGER, INTENT(OUT) :: KRESP             ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL,INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL,INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=40) :: YFIELD
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
REAL   :: XTIME0
INTEGER            :: INFOMPI
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFC0',0,zhook_handle)
!
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfc0_mnh(yrec,yfield,kresp,ycomment)
  hfield = yfield
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfc0_aro(yrec,yfield,kresp,ycomment)
  hfield = yfield
#endif
ELSEIF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' .OR. &
    hprogram=='NC    ' ) THEN 
  !
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    !  
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,yfield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,yfield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,yfield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,yfield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,yfield,kresp,ycomment)
#endif
    ENDIF
    !  
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)  
#endif
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()
    CALL mpi_bcast(yfield,40,mpi_character,npio,ncomm,infompi)
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)
  ENDIF
#endif
  !
  hfield = yfield(1:len(hfield))
  !
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFC0',1,zhook_handle)
!
END SUBROUTINE read_surfc0
!
!     #############################################################
      SUBROUTINE read_surfl0 (&
                              HPROGRAM,HREC,OFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
#ifdef SFX_OL
USE mode_read_surf_ol, ONLY: read_surf0_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surf0_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surf0_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surf0_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surf0_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfl0_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfl0_aro
#endif
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
LOGICAL, INTENT(OUT) :: OFIELD           ! array containing the data field
INTEGER, INTENT(OUT) :: KRESP           ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT  ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
REAL   :: XTIME0
INTEGER            :: INFOMPI
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFL0',0,zhook_handle)
!
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfl0_mnh(yrec,ofield,kresp,ycomment)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfl0_aro(yrec,ofield,kresp,ycomment)
#endif
ELSEIF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' .OR. &
    hprogram=='NC    ' ) THEN 
  !  
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    ! 
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,ofield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,ofield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,ofield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,ofield,kresp,ycomment)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,ofield,kresp,ycomment)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()
    CALL mpi_bcast(ofield,1,mpi_logical,npio,ncomm,infompi)
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF
#endif
  !
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFL0',1,zhook_handle)
!
END SUBROUTINE read_surfl0
!
!     #############################################################
      SUBROUTINE read_surfl1 (HPROGRAM,HREC,OFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
#ifdef SFX_OL
USE mode_read_surf_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surfn_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surfn_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surfn_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surfn_fa
#endif
#ifdef SFX_MNH
USE modi_read_surfl1_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surfl1_aro
#endif

IMPLICIT NONE
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM     ! calling program
!
!
 CHARACTER(LEN=*), INTENT(IN) :: HREC         ! name of the article to be read
LOGICAL, DIMENSION(:), INTENT(OUT) :: OFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP                ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR     ! type of field :
!                                                   ! 'H' : field with
!                                                   !       horizontal spatial dim.
!                                                   ! '-' : no horizontal dim.
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: IL
LOGICAL :: GDATA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFL1',0,zhook_handle)
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
il = SIZE(ofield)
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfl1_mnh(yrec,il,ofield,kresp,ycomment,ydir)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  ! ascendant reproducibility
  IF (yrec(1:2)=="L_".AND.il==1) THEN
    CALL read_surfl0_aro(yrec,gdata,kresp,ycomment)
    ofield = gdata
  ELSE
    CALL read_surfl1_aro(yrec,il,ofield,kresp,ycomment,ydir)
  ENDIF
#endif
ELSE IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
  CALL read_surfl1x2_ol(yrec,ofield,kresp,ycomment,ydir)
#endif
ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
  CALL read_surfn_lfi(yrec,ofield,kresp,ycomment,ydir)
#endif
ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
  CALL read_surfn_nc(yrec,ofield,kresp,ycomment,ydir)
#endif
ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
  CALL read_surfn_asc(&
                      yrec,ofield,kresp,ycomment,ydir)
#endif
ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
  CALL read_surfn_fa(&
                     yrec,il,ofield,kresp,ycomment,ydir)
#endif
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFL1',1,zhook_handle)
!
END SUBROUTINE read_surfl1
!
!     #############################################################
      SUBROUTINE read_surft0 (&
                              HPROGRAM,HREC,TFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_type_date_surf
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
#ifdef SFX_OL
USE mode_read_surf_ol, ONLY: read_surft_ol
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surft_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surft_nc
#endif
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surft_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surft_fa
#endif
#ifdef SFX_MNH
USE modi_read_surft0_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surft0_aro
#endif
!
USE modi_get_luout
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
!RJ: to match actual interface used above, can silently not update values in TFIELD
!RJ TYPE(DATE_TIME), INTENT(OUT) :: TFIELD   ! array containing the data field
TYPE(date_time), INTENT(INOUT) :: TFIELD   ! array containing the data field
INTEGER, INTENT(OUT) :: KRESP            ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
 CHARACTER(LEN=100) :: YCOMMENT
REAL :: ZWORK
 INTEGER, DIMENSION(3) :: IWORK
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
!
REAL   :: XTIME0
INTEGER :: ILUOUT
INTEGER :: INFOMPI
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT0',0,zhook_handle)
!
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surft0_mnh(yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surft0_aro(yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
ELSEIF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' .OR. &
    hprogram=='NC    ' ) THEN 
  !  
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    !
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surft_ol(yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surft_lfi(yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surft_nc(yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surft_asc(&
                          yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surft_fa(&
                         yrec,iwork(1),iwork(2),iwork(3),zwork,kresp,ycomment)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime() 
    CALL mpi_bcast(iwork(1),kind(iwork)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(iwork(2),kind(iwork)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(iwork(3),kind(iwork)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(zwork,kind(zwork)/4,mpi_real,npio,ncomm,infompi)
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF
#endif
  !
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (kresp==-2) THEN
  CALL get_luout(hprogram,iluout)
  WRITE(iluout,*) '-------'
  WRITE(iluout,*) 'WARNING'
  WRITE(iluout,*) '-------'
  WRITE(iluout,*) ' '
  WRITE(iluout,*) 'Date is not present file'
  WRITE(iluout,*) 'Forcing value is kept'
  WRITE(iluout,*) ' '
ELSE
  tfield%TDATE%YEAR = iwork(1)
  tfield%TDATE%MONTH = iwork(2)
  tfield%TDATE%DAY = iwork(3)
  tfield%TIME = zwork
END IF
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT0',1,zhook_handle)
!
END SUBROUTINE read_surft0
!
!     #############################################################
      SUBROUTINE read_surft1 (&
                              HPROGRAM,HREC,TFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_type_date_surf
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surft_asc
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surft_lfi
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surft_nc
#endif
#ifdef SFX_MNH
USE modi_read_surft1_mnh
#endif
#ifdef SFX_ARO
USE modi_read_surft1_aro
#endif
!
USE modi_abor1_sfx
USE modi_get_luout
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM   ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC       ! name of the article to be read
TYPE(date_time), DIMENSION(:), INTENT(INOUT)::TFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP              ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
INTEGER, DIMENSION(:,:), ALLOCATABLE :: IWORK2
REAL, DIMENSION(:), ALLOCATABLE :: ZWORK2
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: ILUOUT
INTEGER            :: INFOMPI
!
 CHARACTER(LEN=100) :: YCOMMENT
REAL   :: XTIME0
INTEGER :: IL1
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT1',0,zhook_handle)
!
il1 = SIZE(tfield,1)
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
ALLOCATE(iwork2(il1,3))
ALLOCATE(zwork2(il1))
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surft1_mnh(yrec,il1,iwork2(:,1),iwork2(:,2),iwork2(:,3),zwork2,kresp,ycomment)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surft1_aro(yrec,il1,iwork2(:,1),iwork2(:,2),iwork2(:,3),zwork2,kresp,ycomment)
#endif
ELSEIF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' .OR. &
    hprogram=='NC    ') THEN 
  !
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    !
    IF (hprogram=='OFFLIN') THEN
      CALL abor1_sfx('READ_SURFT1: NOT AVAILABLE FOR OFFLIN')
    ELSE IF (hprogram=='FA    ') THEN
      CALL abor1_sfx('READ_SURFT1: NOT AVAILABLE FOR FA')      
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surft_asc(&
                          yrec,iwork2(:,1),iwork2(:,2),iwork2(:,3),zwork2,kresp,ycomment)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surft_lfi(yrec,iwork2(:,1),iwork2(:,2),iwork2(:,3),zwork2,kresp,ycomment)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surft_nc(yrec,iwork2(:,1),iwork2(:,2),iwork2(:,3),zwork2,kresp,ycomment)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()         
    CALL mpi_bcast(iwork2(:,1),il1*kind(iwork2)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(iwork2(:,2),il1*kind(iwork2)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(iwork2(:,3),il1*kind(iwork2)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(zwork2,il1*kind(zwork2)/4,mpi_real,npio,ncomm,infompi)
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF
#endif
  !
ENDIF  
! 
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (kresp==-2) THEN
  CALL get_luout(hprogram,iluout)
  WRITE(iluout,*) '-------'
  WRITE(iluout,*) 'WARNING'
  WRITE(iluout,*) '-------'
  WRITE(iluout,*) ' '
  WRITE(iluout,*) 'Date is not present file'
  WRITE(iluout,*) 'Forcing value is kept'
  WRITE(iluout,*) ' '
ELSE
  tfield(:)%TDATE%YEAR  = iwork2(:,1)
  tfield(:)%TDATE%MONTH = iwork2(:,2)
  tfield(:)%TDATE%DAY   = iwork2(:,3)
  tfield(:)%TIME        = zwork2(:)
END IF
!
DEALLOCATE(iwork2)
DEALLOCATE(zwork2)
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT1',1,zhook_handle)
!
END SUBROUTINE read_surft1
!
!     #############################################################
      SUBROUTINE read_surft2 (&
                              HPROGRAM,HREC,TFIELD,KRESP,HCOMMENT,HDIR)
!     #############################################################
!
!
!
!
USE yomhook ,ONLY : lhook, dr_hook
USE parkind1 ,ONLY : jprb
!
USE modd_surfex_mpi, ONLY : nrank, npio, ncomm, nproc, xtime_npio_read, xtime_comm_read
!
USE modd_type_date_surf
!
USE modi_abor1_sfx
USE modi_get_luout
!
#ifdef SFX_ASC
USE mode_read_surf_asc, ONLY: read_surft_asc
#endif
#ifdef SFX_FA
USE mode_read_surf_fa, ONLY: read_surft_fa
#endif
#ifdef SFX_NC
USE mode_read_surf_nc, ONLY: read_surft_nc
#endif
#ifdef SFX_LFI
USE mode_read_surf_lfi, ONLY: read_surft_lfi
#endif
!
IMPLICIT NONE
!
#ifdef SFX_MPI
include "mpif.h"
#endif
!
!*      0.1   Declarations of arguments
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! calling program
 CHARACTER(LEN=*), INTENT(IN) :: HREC     ! name of the article to be read
TYPE(date_time), DIMENSION(:,:), INTENT(INOUT) :: TFIELD ! array containing the data field  
INTEGER, INTENT(OUT) :: KRESP            ! KRESP  : return-code if a problem appears
 CHARACTER(LEN=*), OPTIONAL, INTENT(OUT) :: HCOMMENT   ! name of the article to be read
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: HDIR
!
!*      0.2   Declarations of local variables
!
INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: IWORK3
REAL, DIMENSION(:,:), ALLOCATABLE :: ZWORK3
 CHARACTER(LEN=16)  :: YREC
 CHARACTER(LEN=1)   :: YDIR
INTEGER            :: ILUOUT
INTEGER            :: INFOMPI
!
 CHARACTER(LEN=100) :: YCOMMENT
INTEGER :: IL1, IL2
REAL   :: XTIME0
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT2',0,zhook_handle)
!
kresp = 0
ycomment = ""
!
yrec = hrec
ydir = 'H'
IF (PRESENT(hdir)) ydir = hdir
!
il1 = SIZE(tfield,1)
il2 = SIZE(tfield,2)
!
ALLOCATE(zwork3(il1,il2))
ALLOCATE(iwork3(il1,il2,3))
!
IF (hprogram=='MESONH') THEN
  CALL abor1_sfx('READ_SURFT2: NOT AVAILABLE FOR MESONH')
ELSE IF (hprogram=='AROME ') THEN
  CALL abor1_sfx('READ_SURFT2: NOT AVAILABLE FOR AROME')  
ELSEIF (hprogram=='OFFLIN' .OR. hprogram=='ASCII ' .OR. &
    hprogram=='FA    ' .OR. hprogram=='LFI   ' ) THEN 
  !  
  IF (nrank==npio) THEN
    !
#ifdef SFX_MPI
    xtime0 = mpi_wtime()
#endif
    !
    IF (hprogram=='OFFLIN') THEN
      CALL abor1_sfx('READ_SURFT2: NOT AVAILABLE FOR OFFLIN')
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surft_lfi(yrec,iwork3(:,:,1),iwork3(:,:,2),iwork3(:,:,3),&
        zwork3,kresp,ycomment)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surft_asc(&
                          yrec,iwork3(:,:,1),iwork3(:,:,2),iwork3(:,:,3),&
        zwork3,kresp,ycomment)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surft_fa(&
                         yrec,il1,il2,iwork3(:,:,1),iwork3(:,:,2),iwork3(:,:,3),&
        zwork3,kresp,ycomment)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surft_nc(yrec,iwork3(:,:,1),iwork3(:,:,2),iwork3(:,:,3),&
        zwork3,kresp,ycomment)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !    
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()   
    CALL mpi_bcast(iwork3(:,:,1),il1*il2*kind(iwork3)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(iwork3(:,:,2),il1*il2*kind(iwork3)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(iwork3(:,:,3),il1*il2*kind(iwork3)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(zwork3,il1*il2*kind(zwork3)/4,mpi_real,npio,ncomm,infompi)
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)
  ENDIF
#endif
  !
ENDIF
!
IF (PRESENT(hcomment)) hcomment = ycomment
!
IF (kresp==-2) THEN
  CALL get_luout(hprogram,iluout)
  WRITE(iluout,*) '-------'
  WRITE(iluout,*) 'WARNING'
  WRITE(iluout,*) '-------'
  WRITE(iluout,*) ' '
  WRITE(iluout,*) 'Date is not present file'
  WRITE(iluout,*) 'Forcing value is kept'
  WRITE(iluout,*) ' '
ELSE
  tfield(:,:)%TDATE%YEAR  = iwork3(:,:,1)
  tfield(:,:)%TDATE%MONTH = iwork3(:,:,2)
  tfield(:,:)%TDATE%DAY   = iwork3(:,:,3)
  tfield(:,:)%TIME        = zwork3(:,:) 
END IF
!
DEALLOCATE(iwork3)
DEALLOCATE(zwork3)
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT2',1,zhook_handle)
!
END SUBROUTINE read_surft2