ol_read_atm_conf.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
SUBROUTINE ol_read_atm_conf (DTCO, U, HGRID, HSURF_FILETYPE, HFORCING_FILETYPE, ODELAYEDSTART_NC, &
                             KDATESTOP, PDURATION, PTSTEP_FORC, KNI, KYEAR, KMONTH, KDAY,  &
                             PTIME, PLAT, PLON, PZS, PZREF, PUREF, KTIMESTARTINDEX    )  
!
!==================================================================
!!****  *OL_READ_ATM_CONF* - Initialization routine
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      F. Habets   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      Modified by P. Le Moigne (04/2005): cleaning and checking
!!      Modified by P. Le Moigne (04/2006): init_io_surf for nature
!!                  with GTMSK to read dimensions.
!!      Modified by M. Lafaysse 04/2015 : option ODELAYEDSTART_NC
!==================================================================
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modi_ol_read_atm_conf_netcdf
USE modi_ol_read_atm_conf_ascii
USE modd_surf_conf,      ONLY : cprogname
!==================================================================
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(surf_atm_t), INTENT(INOUT) :: U
!
 CHARACTER(LEN=*), INTENT(IN)  :: HGRID
 CHARACTER(LEN=6), INTENT(IN)  :: HSURF_FILETYPE
 CHARACTER(LEN=6), INTENT(IN)  :: HFORCING_FILETYPE
LOGICAL, INTENT(IN)            :: ODELAYEDSTART_NC !Allow the simulation to start from a different time step than the first record of a netcdf file 
INTEGER,DIMENSION(4),INTENT(IN) :: KDATESTOP !Allow the simulation to end at a different time step than the last record of a netcdf file
INTEGER,          INTENT(OUT) :: KNI
INTEGER,          INTENT(OUT) :: KYEAR, KMONTH, KDAY
REAL,             INTENT(OUT) :: PDURATION,PTSTEP_FORC
REAL,             INTENT(OUT) :: PTIME
REAL, DIMENSION(:),  POINTER  :: PLAT, PLON
REAL, DIMENSION(:),  POINTER  :: PZS
REAL, DIMENSION(:),  POINTER  :: PZREF, PUREF
INTEGER,          INTENT(OUT) :: KTIMESTARTINDEX ! index from which we start reading FORCING.nc
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!==================================================================
!
IF (lhook) CALL dr_hook('OL_READ_ATM_CONF',0,zhook_handle)
cprogname = hsurf_filetype
!
IF (hforcing_filetype == 'NETCDF') THEN
!
 CALL ol_read_atm_conf_netcdf(dtco, u, hgrid, hsurf_filetype, odelayedstart_nc, kdatestop, &
                              pduration, ptstep_forc, kni, kyear, kmonth, kday, ptime, &
                              plat, plon, pzs, pzref, puref, ktimestartindex        )  
!
ELSE IF (hforcing_filetype == 'ASCII ' .OR. hforcing_filetype == 'BINARY') THEN
!
 CALL ol_read_atm_conf_ascii(dtco, u, hsurf_filetype,            &
                             hforcing_filetype, pduration,       &
                             ptstep_forc, kni, kyear,kmonth,     &
                             kday, ptime, plat, plon,            &
                             pzs, pzref, puref                   )  
!    
 ktimestartindex = 1
!
ENDIF
IF (lhook) CALL dr_hook('OL_READ_ATM_CONF',1,zhook_handle)
!
!==================================================================
!
END SUBROUTINE ol_read_atm_conf