default_ideal_flux.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 default_ideal_flux(KFORCF, KFORCT, PTIMEF, PTIMET,   &
                                    psfth, psftq, psfco2,             &
                                    hustartype, pustar, pz0m, palb,   &
                                    pemis, ptsrad)
!     ########################################################################
!
!!****  *DEFAULT_IDEAL_FLUX* - routine to set default values for the configuration for ISBA FLUX assimilation scheme
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      L. Jarlan  *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2005
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_surf_par,   ONLY : xundef
USE modd_csts,       ONLY : xtt
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
INTEGER, INTENT(OUT) :: kforcf
INTEGER, INTENT(OUT) :: kforct
REAL, DIMENSION(:), INTENT(OUT)   :: ptimef
REAL, DIMENSION(:), INTENT(OUT)   :: ptimet
REAL, DIMENSION(:), INTENT(OUT)   :: psfth      ! hourly data of heat surface flux        (W/m2)
REAL, DIMENSION(:), INTENT(OUT)   :: psftq      ! hourly data of water vapor surface flux (kg/m2/s) or (W/m2)
REAL, DIMENSION(:), INTENT(OUT)   :: psfco2     ! hourly data of CO2 surface flux         (kg/m2/s)
 CHARACTER(LEN=5), INTENT(OUT)     :: hustartype ! type of computation for friction
                                                ! 'USTAR'
                                                ! 'Z0   '
REAL, DIMENSION(:), INTENT(OUT)   :: pustar     ! hourly data of friction                 (m2/s2)
REAL, INTENT(OUT)                 :: pz0m       ! roughness length (m)
REAL, INTENT(OUT)                 :: palb       ! albedo (-)
REAL, INTENT(OUT)                 :: pemis      ! emissivity (-)
REAL, DIMENSION(:),INTENT(OUT)    :: ptsrad     ! radiative temperature (K)
!
!*       0.2   declarations of local variables
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('DEFAULT_IDEAL_FLUX',0,zhook_handle)
!
kforcf = 2.
kforct = 2.
!
ptimef(1) = 0.
ptimet(1) = 0.
ptimef(2) = xundef
ptimet(2) = xundef

!
!----------------------------------------------------------------------------------
!
!*       1.    HOURLY surface theta flux (NFORC+1 values from 00UTC to 24UTC)
!               -------------------------
!
!* unit: W/m2
!
psfth(:) = 0.
!
!----------------------------------------------------------------------------------
!
!*       2.    HOURLY surface vapor mixing ratio flux (NFORC+1 values from 00UTC to 24UTC)
!              --------------------------------------
!
!* unit: kg/m2/s
!
psftq(:) = 0.
!
!----------------------------------------------------------------------------------
!
!*       4.    HOURLY surface CO2 flux (NFORC+1 values from 00UTC to 24UTC)
!              -----------------------
!
!* unit: kg/m2/s
!
psfco2(:) = 0.
!
!----------------------------------------------------------------------------------
!
!*       5.    Type of definition for friction fluxes
!              --------------------------------------
!
!* HUSTARTYPE = 'Z0   '  ! friction is defined using a roughness length formulation
!*            = 'USTAR'  ! friction is prescribed via the friction velocity u*
!
hustartype = 'Z0   '
!
!----------------------------------------------------------------------------------
!
!*       6.    Roughness length (used if XUSTARTYPE = 'Z0   ')
!              ----------------
!
pz0m = 0.01              ! unit in meters
!      
!----------------------------------------------------------------------------------
!
!*       6.    Friction (used if XUSTARTYPE = 'USTAR')
!              --------
!
!* unit: m2/s2
!
pustar = 0.
!
!
!----------------------------------------------------------------------------------
!
!*       7.    HOURLY surface radiative temperature (NFORC+1 values from 00UTC to 24UTC)
!               -------------------------
!

ptsrad(:) = xtt            ! radiative surface temperature, (unit is K)
!
!----------------------------------------------------------------------------------
!
!*       8.    Radiative fields constant and uniform values
!              --------------------------------------------
!
palb      = 0.             ! albedo, no unit
pemis     = 1.             ! emissivity, no unit
!      
IF (lhook) CALL dr_hook('DEFAULT_IDEAL_FLUX',1,zhook_handle)
!
END SUBROUTINE default_ideal_flux