ol_read_atm_binary.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_binary(KFORC_STEP,                                 &
                              pta,pqa,pwind,pdir_sw,psca_sw,plw,psnow,    &
                              prain,pps,pco2,pdir                         )  
!**************************************************************************
!
!!    PURPOSE
!!    -------
!         Read in the ascii file the atmospheric forcing for the actual time
!         step KFORC_STEP, and for the next one.
!         The two time step are needed for the time interpolation of the
!         forcing.
!         If the end of the file  is reached, set the two step to the last
!         values.
!         Return undef value if the variable is not present
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      A. Lemonsu  *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original     03/2008        
!
USE modd_io_surf_ol, ONLY : xcount
USE modi_read_surf_atm
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
! global variables
REAL, DIMENSION(:,:),INTENT(INOUT) :: pta
REAL, DIMENSION(:,:),INTENT(INOUT) :: pqa
REAL, DIMENSION(:,:),INTENT(INOUT) :: pwind
REAL, DIMENSION(:,:),INTENT(INOUT) :: pdir_sw
REAL, DIMENSION(:,:),INTENT(INOUT) :: psca_sw
REAL, DIMENSION(:,:),INTENT(INOUT) :: plw
REAL, DIMENSION(:,:),INTENT(INOUT) :: psnow
REAL, DIMENSION(:,:),INTENT(INOUT) :: prain
REAL, DIMENSION(:,:),INTENT(INOUT) :: pps
REAL, DIMENSION(:,:),INTENT(INOUT) :: pco2
REAL, DIMENSION(:,:),INTENT(INOUT) :: pdir
INTEGER,INTENT(IN)               :: kforc_step
! local variables
REAL(KIND=JPRB) :: zhook_handle
! 
! read data
IF (lhook) CALL dr_hook('OL_READ_ATM_BINARY',0,zhook_handle)
 CALL read_surf_atm('BINARY',pta(:,1:xcount),kforc_step,xcount,22)
 CALL read_surf_atm('BINARY',pqa(:,1:xcount),kforc_step,xcount,23)
 CALL read_surf_atm('BINARY',pwind(:,1:xcount),kforc_step,xcount,24)
 CALL read_surf_atm('BINARY',plw(:,1:xcount),kforc_step,xcount,25)
 CALL read_surf_atm('BINARY',pdir_sw(:,1:xcount),kforc_step,xcount,26)
 CALL read_surf_atm('BINARY',psca_sw(:,1:xcount),kforc_step,xcount,27)
 CALL read_surf_atm('BINARY',prain(:,1:xcount),kforc_step,xcount,28)
 CALL read_surf_atm('BINARY',psnow(:,1:xcount),kforc_step,xcount,29)
 CALL read_surf_atm('BINARY',pps(:,1:xcount),kforc_step,xcount,30)
 CALL read_surf_atm('BINARY',pdir(:,1:xcount),kforc_step,xcount,31)
 CALL read_surf_atm('BINARY',pco2(:,1:xcount),kforc_step,xcount,32)
!
IF (lhook) CALL dr_hook('OL_READ_ATM_BINARY',1,zhook_handle)
!
END SUBROUTINE ol_read_atm_binary