init_pgd_surf_atm.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 init_pgd_surf_atm (YSC, &
                              hprogram,hinit,hatmfile,hatmfiletype, &
                               kyear, kmonth, kday, ptime            )  
!     #################################################################################
!
!!****  *INIT_PGD_SURF_ATM* - Call surface initialization for PGD fields only
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     V. Masson 
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      B. Decharme  04/2013 new coupling variables
!!------------------------------------------------------------------
!
!
!
USE modd_surfex_n, ONLY : surfex_t
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_init_surf_atm_n
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
TYPE(surfex_t), INTENT(INOUT) :: ysc
!
 CHARACTER(LEN=6),   INTENT(IN)  :: hprogram  ! program calling surf. schemes
 CHARACTER(LEN=3),   INTENT(IN)  :: hinit     ! fields to initialize 'ALL', 'PRE', 'PGD'
 CHARACTER(LEN=28), INTENT(IN)   :: hatmfile    ! name of the Atmospheric file
 CHARACTER(LEN=6),  INTENT(IN)   :: hatmfiletype! type of the Atmospheric file
INTEGER,           INTENT(IN)   :: kyear       ! year
INTEGER,           INTENT(IN)   :: kmonth      ! month
INTEGER,           INTENT(IN)   :: kday        ! day
REAL,              INTENT(IN)   :: ptime       ! time
!
!
!*      0.2    declarations of local variables
!
 CHARACTER(LEN=6), DIMENSION(0)  :: ysv       ! name of all scalar variables
REAL,             DIMENSION(0)  :: zco2      ! CO2 concentration (kg/m3)
REAL,             DIMENSION(0)  :: zrhoa     ! air density (kg/m3)
REAL,             DIMENSION(0)  :: zzenith   ! solar zenithal angle
REAL,             DIMENSION(0)  :: zazim     ! solar azimuthal angle (rad from N, clock)
REAL,             DIMENSION(1)  :: zsw_bands ! middle wavelength of each band
REAL,             DIMENSION(0,1):: zdir_alb  ! direct albedo for each band
REAL,             DIMENSION(0,1):: zsca_alb  ! diffuse albedo for each band
REAL,             DIMENSION(0)  :: zemis     ! emissivity
REAL,             DIMENSION(0)  :: ztsrad    ! radiative temperature
REAL,             DIMENSION(0)  :: ztsurf    ! radiative temperature
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------------
!
!* initialization of PGD fields of output domain
!
IF (lhook) CALL dr_hook('INIT_PGD_SURF_ATM',0,zhook_handle)
 CALL init_surf_atm_n(ysc, &
                      hprogram,hinit,.false.,                     &
                      0,0,1,                                     &
                      ysv,zco2,zrhoa,                            &
                      zzenith,zazim,zsw_bands,zdir_alb,zsca_alb, &
                      zemis,ztsrad,ztsurf,                       &
                      kyear, kmonth, kday, ptime,                &
                      hatmfile,hatmfiletype, 'OK'                )  
IF (lhook) CALL dr_hook('INIT_PGD_SURF_ATM',1,zhook_handle)

!
!-------------------------------------------------------------------------------------
!
END SUBROUTINE init_pgd_surf_atm