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
USE modd_surfex_omp, ONLY : xwork0, nworkb, cworkb
!
#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
!
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=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)
!
!$OMP BARRIER
!
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfx0_mnh(yrec,xwork0,nworkb,cworkb)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfx0_aro(yrec,xwork0,nworkb,cworkb)
#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
    !
!$OMP SINGLE
    !
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,xwork0,nworkb,cworkb)
#endif
    ELSEIF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,xwork0,nworkb,cworkb)
#endif
    ELSEIF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,xwork0,nworkb,cworkb)
#endif
    ELSEIF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,xwork0,nworkb,cworkb)
#endif
    ELSEIF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,xwork0,nworkb,cworkb)
#endif
    ENDIF
    !
!$OMP END SINGLE
    !
#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()
!$OMP SINGLE   
    CALL mpi_bcast(xwork0,kind(xwork0)/4,mpi_real,npio,ncomm,infompi)
!$OMP END SINGLE
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)
  ENDIF
#endif
  !
ENDIF
!    
pfield=xwork0
!
kresp = nworkb
!
IF (present(hcomment)) hcomment = cworkb
!
!$OMP BARRIER
!
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, ONLY: read_surfn_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
!
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_surfn_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, ONLY: read_surfn_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
!
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_surfn_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, ONLY: read_surfn_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_surfn_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
USE modd_surfex_omp, ONLY : nwork0, nworkb, cworkb
!
#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_surfn0_mnh
#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=16)  :: yrec
 CHARACTER(LEN=1)   :: ydir
REAL   :: xtime0
INTEGER            :: infompi
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFN0',0,zhook_handle)
!
!$OMP BARRIER
!
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfn0_mnh(yrec,nwork0,nworkb,cworkb)
#endif
ENDIF
!
IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfn0_aro(yrec,nwork0,nworkb,cworkb)
#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
    !
!$OMP SINGLE
    !    
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,nwork0,nworkb,cworkb)
#endif
    ENDIF
    !
    IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,nwork0,nworkb,cworkb)
#endif
    ENDIF
    !
    IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,nwork0,nworkb,cworkb)
#endif
    ENDIF    
    !
    IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,nwork0,nworkb,cworkb)
#endif
    ENDIF
    !
    IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,nwork0,nworkb,cworkb)
#endif
    ENDIF
    !
!$OMP END SINGLE
    !    
#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()  
!$OMP SINGLE
    CALL mpi_bcast(nwork0,kind(nwork0)/4,mpi_integer,npio,ncomm,infompi)
!$OMP END SINGLE  
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF    
#endif
  !
ENDIF
!
kfield=nwork0
!
kresp = nworkb
!
IF (present(hcomment)) hcomment = cworkb
!
!$OMP BARRIER
!
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, ONLY: read_surfn_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
!
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_surfn_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
USE modd_surfex_omp, ONLY : cwork0, nworkb, cworkb
!
#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
!
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=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)
!
!$OMP BARRIER
!
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfc0_mnh(yrec,cwork0(1:40),nworkb,cworkb)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfc0_aro(yrec,cwork0(1:40),nworkb,cworkb)
#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
    !  
!$OMP SINGLE
    !    
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,cwork0(1:40),nworkb,cworkb)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,cwork0(1:40),nworkb,cworkb)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,cwork0(1:40),nworkb,cworkb)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,cwork0(1:40),nworkb,cworkb)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,cwork0(1:40),nworkb,cworkb)
#endif
    ENDIF
    !  
!$OMP END SINGLE
    !
#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()
!$OMP SINGLE
    CALL mpi_bcast(cwork0(1:40),40,mpi_character,npio,ncomm,infompi)
!$OMP END SINGLE
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)
  ENDIF
#endif
  !
ENDIF
!
hfield = cwork0(1:len(hfield))
!
kresp = nworkb
!
IF (present(hcomment)) hcomment = cworkb
!
!$OMP BARRIER
!
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
USE modd_surfex_omp, ONLY : lwork0, nworkb, cworkb
!
#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
!
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=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)
!
!$OMP BARRIER
!
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surfl0_mnh(yrec,lwork0,nworkb,cworkb)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surfl0_aro(yrec,lwork0,nworkb,cworkb)
#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
    ! 
!$OMP SINGLE
    !
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surf0_ol(yrec,lwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surf0_lfi(yrec,lwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surf0_nc(yrec,lwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surf0_asc(&
                          yrec,lwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surf0_fa(&
                         yrec,lwork0,nworkb,cworkb)
#endif
    ENDIF
    !
!$OMP END SINGLE
    ! 
#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()
!$OMP SINGLE
    CALL mpi_bcast(lwork0,1,mpi_logical,npio,ncomm,infompi)
!$OMP END SINGLE
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF
#endif
  !
ENDIF
!
ofield = lwork0
!
kresp = nworkb
!
IF (present(hcomment)) hcomment = cworkb
!
!$OMP BARRIER
!
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, ONLY: read_surfn_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
!
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
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
  CALL read_surfl1_aro(yrec,il,ofield,kresp,ycomment,ydir)
#endif
ELSE IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
  CALL read_surfn_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
USE modd_surfex_omp, ONLY : nworkd, xwork0, nworkb, cworkb
!
#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
!
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=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)
!
!$OMP BARRIER
!
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
!$OMP SINGLE
ALLOCATE(nworkd(3))
!$OMP END SINGLE
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surft0_mnh(yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surft0_aro(yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#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
    !
!$OMP SINGLE
    !    
    IF (hprogram=='OFFLIN') THEN
#ifdef SFX_OL
      CALL read_surft_ol(yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surft_lfi(yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surft_nc(yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surft_asc(&
                          yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surft_fa(&
                         yrec,nworkd(1),nworkd(2),nworkd(3),xwork0,nworkb,cworkb)
#endif
    ENDIF
    !
!$OMP END SINGLE
    ! 
#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() 
!$OMP SINGLE    
    CALL mpi_bcast(nworkd(1),kind(nworkd)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(nworkd(2),kind(nworkd)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(nworkd(3),kind(nworkd)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(xwork0,kind(xwork0)/4,mpi_real,npio,ncomm,infompi)
!$OMP END SINGLE 
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF
#endif
  !
ENDIF
!
kresp = nworkb
!
IF (present(hcomment)) hcomment = cworkb
!
!$OMP BARRIER
!
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 = nworkd(1)
  tfield%TDATE%MONTH = nworkd(2)
  tfield%TDATE%DAY = nworkd(3)
  tfield%TIME = xwork0
END IF
!
!$OMP BARRIER
!
!$OMP SINGLE
DEALLOCATE(nworkd)
!$OMP END SINGLE
!
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
USE modd_surfex_omp, ONLY : nworkd2, xworkd, nworkb, cworkb
!
#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
!
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
!
 CHARACTER(LEN=16)  :: yrec
 CHARACTER(LEN=1)   :: ydir
INTEGER            :: iluout
INTEGER            :: infompi
!
REAL   :: xtime0
INTEGER :: il1
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT1',0,zhook_handle)
!
!$OMP BARRIER
!
il1 = SIZE(tfield,1)
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
!$OMP SINGLE
ALLOCATE(nworkd2(il1,3))
ALLOCATE(xworkd(il1))
!$OMP END SINGLE
!
IF (hprogram=='MESONH') THEN
#ifdef SFX_MNH
  CALL read_surft1_mnh(yrec,il1,nworkd2(:,1),nworkd2(:,2),nworkd2(:,3),xworkd,nworkb,cworkb)
#endif
ELSE IF (hprogram=='AROME ') THEN
#ifdef SFX_ARO
  CALL read_surft1_aro(yrec,il1,nworkd2(:,1),nworkd2(:,2),nworkd2(:,3),xworkd,nworkb,cworkb)
#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
    !
!$OMP SINGLE
    !
    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,nworkd2(:,1),nworkd2(:,2),nworkd2(:,3),xworkd,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='LFI   ') THEN
#ifdef SFX_LFI
      CALL read_surft_lfi(yrec,nworkd2(:,1),nworkd2(:,2),nworkd2(:,3),xworkd,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surft_nc(yrec,nworkd2(:,1),nworkd2(:,2),nworkd2(:,3),xworkd,nworkb,cworkb)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !
!$OMP END SINGLE
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()         
!$OMP SINGLE    
    CALL mpi_bcast(nworkd2(:,1),il1*kind(nworkd2)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(nworkd2(:,2),il1*kind(nworkd2)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(nworkd2(:,3),il1*kind(nworkd2)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(xworkd,il1*kind(xworkd)/4,mpi_real,npio,ncomm,infompi)
!$OMP END SINGLE
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)    
  ENDIF
#endif
  !
ENDIF  
! 
kresp = nworkb
IF (present(hcomment)) hcomment = cworkb
!
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  = nworkd2(:,1)
  tfield(:)%TDATE%MONTH = nworkd2(:,2)
  tfield(:)%TDATE%DAY   = nworkd2(:,3)
  tfield(:)%TIME        = xworkd(:)
END IF
!
!$OMP BARRIER
!
!$OMP SINGLE
DEALLOCATE(nworkd2)
DEALLOCATE(xworkd)
!$OMP END SINGLE
!
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_surfex_omp, ONLY : nworkd3, xworkd2, nworkb, cworkb
!
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
!
 CHARACTER(LEN=16)  :: yrec
 CHARACTER(LEN=1)   :: ydir
INTEGER            :: iluout
INTEGER            :: infompi
!
INTEGER :: il1, il2
REAL   :: xtime0
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT2',0,zhook_handle)
!
!$OMP BARRIER
!
nworkb = 0
 cworkb = ""
!
yrec = hrec
ydir = 'H'
IF (present(hdir)) ydir = hdir
!
il1 = SIZE(tfield,1)
il2 = SIZE(tfield,2)
!
!$OMP SINGLE
ALLOCATE(xworkd2(il1,il2))
ALLOCATE(nworkd3(il1,il2,3))
!$OMP END SINGLE
!
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
    !
!$OMP SINGLE
    !
    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,nworkd3(:,:,1),nworkd3(:,:,2),nworkd3(:,:,3),&
        xworkd2,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='ASCII ') THEN
#ifdef SFX_ASC
      CALL read_surft_asc(&
                          yrec,nworkd3(:,:,1),nworkd3(:,:,2),nworkd3(:,:,3),&
        xworkd2,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='FA    ') THEN
#ifdef SFX_FA
      CALL read_surft_fa(&
                         yrec,il1,il2,nworkd3(:,:,1),nworkd3(:,:,2),nworkd3(:,:,3),&
        xworkd2,nworkb,cworkb)
#endif
    ELSE IF (hprogram=='NC    ') THEN
#ifdef SFX_NC
      CALL read_surft_nc(yrec,nworkd3(:,:,1),nworkd3(:,:,2),nworkd3(:,:,3),&
        xworkd2,nworkb,cworkb)
#endif
    ENDIF
    !
#ifdef SFX_MPI
    xtime_npio_read = xtime_npio_read + (mpi_wtime() - xtime0)
#endif
    !    
!$OMP END SINGLE
    !
  ENDIF
  !
#ifdef SFX_MPI
  IF (ydir/='A' .AND. nproc>1) THEN
    xtime0 = mpi_wtime()   
!$OMP SINGLE
    CALL mpi_bcast(nworkd3(:,:,1),il1*il2*kind(nworkd3)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(nworkd3(:,:,2),il1*il2*kind(nworkd3)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(nworkd3(:,:,3),il1*il2*kind(nworkd3)/4,mpi_integer,npio,ncomm,infompi)
    CALL mpi_bcast(xworkd2,il1*il2*kind(xworkd2)/4,mpi_real,npio,ncomm,infompi)
!$OMP END SINGLE     
    xtime_comm_read = xtime_comm_read + (mpi_wtime() - xtime0)
  ENDIF
#endif
  !
ENDIF
!
kresp = nworkb
IF (present(hcomment)) hcomment = cworkb
!
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  = nworkd3(:,:,1)
  tfield(:,:)%TDATE%MONTH = nworkd3(:,:,2)
  tfield(:,:)%TDATE%DAY   = nworkd3(:,:,3)
  tfield(:,:)%TIME        = xworkd2(:,:) 
END IF
!
!$OMP BARRIER
!
!$OMP SINGLE
DEALLOCATE(nworkd3)
DEALLOCATE(xworkd2)
!$OMP END SINGLE
!
IF (lhook) CALL dr_hook('MODI_READ_SURF:READ_SURFT2',1,zhook_handle)
!
END SUBROUTINE read_surft2