init_naturen.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_nature_n (DTCO, OREAD_BUDGETC, UG, U, USS, GCP, IM, &
                                DTZ, DGO, DL, DLC, NDST, SLT, SV,         &
                                HPROGRAM,HINIT,OLAND_USE, KI,KSV,KSW,     &
                                HSV,PCO2,PRHOA,PZENITH,PAZIM,PSW_BANDS,   &
                                PDIR_ALB,PSCA_ALB,PEMIS,PTSRAD,PTSURF,    &
                                KYEAR, KMONTH,KDAY, PTIME, TPDATE_END,    &
                                HATMFILE,HATMFILETYPE,HTEST              )  
!     #############################################################
!
!!****  *INIT_NATURE_n* - routine to choose initialization of vegetation scheme
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    27/09/96 
!!       V.Masson   18/08/97 call to fmread directly with dates and strings
!!       V.Masson   15/03/99 new PGD treatment with COVER types
!        F.Solmon  06/00   adaptation for patch approach
!!      B. Decharme  04/2013 new coupling variables
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE modd_surfex_n, ONLY : isba_model_t
!
USE modd_type_date_surf, ONLY : date
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_diag_n, ONLY : diag_t, diag_options_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_sso_n, ONLY : sso_t
USE modd_grid_conf_proj_n, ONLY : grid_conf_proj_t
USE modd_data_tsz0_n, ONLY : data_tsz0_t
USE modd_dst_n, ONLY : dst_np_t
USE modd_slt_n, ONLY : slt_t
USE modd_sv_n, ONLY : sv_t
!
USE modd_csts,       ONLY : xtt
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_init_ideal_flux
!
USE modi_init_isba_n
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
TYPE(isba_model_t), INTENT(INOUT) :: IM
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
LOGICAL, INTENT(IN) :: OREAD_BUDGETC
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
TYPE(surf_atm_t), INTENT(INOUT) :: U
TYPE(sso_t), INTENT(INOUT) :: USS
TYPE(grid_conf_proj_t),INTENT(INOUT) :: GCP
TYPE(data_tsz0_t), INTENT(INOUT) :: DTZ
TYPE(diag_options_t), INTENT(INOUT) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DL
TYPE(diag_t), INTENT(INOUT) :: DLC
TYPE(dst_np_t), INTENT(INOUT) :: NDST
TYPE(slt_t), INTENT(INOUT) :: SLT
TYPE(sv_t), INTENT(INOUT) :: SV
!
 CHARACTER(LEN=6),                 INTENT(IN)  :: HPROGRAM  ! program calling surf. schemes
 CHARACTER(LEN=3),                 INTENT(IN)  :: HINIT     ! choice of fields to initialize
LOGICAL,                          INTENT(IN)  :: OLAND_USE ! 
INTEGER,                          INTENT(IN)  :: KI        ! number of points
INTEGER,                          INTENT(IN)  :: KSV       ! number of scalars
INTEGER,                          INTENT(IN)  :: KSW       ! number of short-wave spectral bands
 CHARACTER(LEN=6), DIMENSION(KSV), INTENT(IN)  :: HSV       ! name of all scalar variables
REAL,             DIMENSION(KI),  INTENT(IN)  :: PCO2      ! CO2 concentration (kg/m3)
REAL,             DIMENSION(KI),  INTENT(IN)  :: PRHOA     ! air density
REAL,             DIMENSION(KI),  INTENT(IN)  :: PZENITH   ! solar zenithal angle
REAL,             DIMENSION(KI),  INTENT(IN)  :: PAZIM     ! solar azimuthal angle (rad from N, clock)
REAL,             DIMENSION(KSW), INTENT(IN)  :: PSW_BANDS ! middle wavelength of each band
REAL,             DIMENSION(KI,KSW),INTENT(OUT) :: PDIR_ALB  ! direct albedo for each band
REAL,             DIMENSION(KI,KSW),INTENT(OUT) :: PSCA_ALB  ! diffuse albedo for each band
REAL,             DIMENSION(KI),  INTENT(OUT) :: PEMIS     ! emissivity
REAL,             DIMENSION(KI),  INTENT(OUT) :: PTSRAD    ! radiative temperature
REAL,             DIMENSION(KI),  INTENT(OUT) :: PTSURF    ! surface effective temperature         (K)
INTEGER,                          INTENT(IN)  :: KYEAR     ! current year (UTC)
INTEGER,                          INTENT(IN)  :: KMONTH    ! current month (UTC)
INTEGER,                          INTENT(IN)  :: KDAY      ! current day (UTC)
REAL,                             INTENT(IN)  :: PTIME     ! current time since
                                                          !  midnight (UTC, s)
TYPE(date), INTENT(INOUT) :: TPDATE_END
!
 CHARACTER(LEN=28),                INTENT(IN)  :: HATMFILE    ! atmospheric file name
 CHARACTER(LEN=6),                 INTENT(IN)  :: HATMFILETYPE! atmospheric file type
 CHARACTER(LEN=2),                 INTENT(IN)  :: HTEST       ! must be equal to 'OK'
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
!-------------------------------------------------------------------------------
!
!*       2.     Selection of surface scheme
!               ---------------------------
!
IF (lhook) CALL dr_hook('INIT_NATURE_N',0,zhook_handle)
IF (u%CNATURE=='NONE  ') THEN
  pdir_alb=0.
  psca_alb=0.
  pemis   =1.
  ptsrad  =xtt
  ptsurf  =xtt
ELSE IF (u%CNATURE=='FLUX  ') THEN
  CALL init_ideal_flux(dgo, dl, dlc, oread_budgetc, hprogram, hinit, &
                       ki, ksv, ksw, hsv, pdir_alb, psca_alb, pemis,   &
                       ptsrad, ptsurf, 'OK'    )  
ELSE IF (u%CNATURE=='ISBA  ' .OR. u%CNATURE=='TSZ0') THEN
  CALL init_isba_n(dtco, oread_budgetc, ug, u, uss, gcp, &
                   im, dtz, ndst, slt, sv, &
                   hprogram, hinit, oland_use, ki, ksv, ksw, hsv, &
                   pco2, prhoa, pzenith, pazim, psw_bands,        &
                   pdir_alb, psca_alb, pemis, ptsrad, ptsurf,     &
                   kyear, kmonth, kday, ptime, tpdate_end,        &
                   hatmfile, hatmfiletype, 'OK'     )  
END IF
IF (lhook) CALL dr_hook('INIT_NATURE_N',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE init_nature_n