init_water_sbl.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_water_sbl(SB, PPA, PPS, PTA, PQA, PRHOA, PU, PV, PRAIN, PSNOW,     &
                              PSFTH, PSFTQ, PZREF, PUREF, PTS, PZ0 )
!     #################################################################################
!
!!****  *INIT_WATER_SBL* - inits water SBL profiles
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     S. Riette
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    03/2010
!!------------------------------------------------------------------
!
USE modd_canopy_n, ONLY : canopy_t
!
USE modd_csts,             ONLY : xcpd, xrd, xp00, xtt, xg
USE modd_canopy_turb,      ONLY : xalpsbl
!
USE modi_cls_wind
USE modi_cls_tq
USE modi_water_flux
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
TYPE(canopy_t), INTENT(INOUT) :: SB
!
REAL, DIMENSION(:), INTENT(IN)  :: PPA       ! pressure at forcing level             (Pa)
REAL, DIMENSION(:), INTENT(IN)  :: PPS       ! pressure at atmospheric model surface (Pa)
REAL, DIMENSION(:), INTENT(IN)  :: PTA       ! air temperature forcing               (K)
REAL, DIMENSION(:), INTENT(IN)  :: PQA       ! air humidity forcing                  (kg/m3)
REAL, DIMENSION(:), INTENT(IN)  :: PRHOA     ! air density                           (kg/m3)
REAL, DIMENSION(:), INTENT(IN)  :: PU        ! zonal wind                            (m/s)
REAL, DIMENSION(:), INTENT(IN)  :: PV        ! meridian wind                         (m/s)
REAL, DIMENSION(:), INTENT(IN)  :: PSNOW     ! snow precipitation                    (kg/m2/s)
REAL, DIMENSION(:), INTENT(IN)  :: PRAIN     ! liquid precipitation                  (kg/m2/s)
REAL, DIMENSION(:), INTENT(IN)  :: PZREF     ! height of T,q forcing                 (m)
REAL, DIMENSION(:), INTENT(IN)  :: PUREF     ! height of wind forcing                (m)
REAL, DIMENSION(:), INTENT(IN)  :: PTS       ! surface temperature
!
REAL, DIMENSION(:), INTENT(INOUT) :: PZ0       ! roughness length
REAL, DIMENSION(:), INTENT(OUT) :: PSFTH     ! flux of heat                          (W/m2)
REAL, DIMENSION(:), INTENT(OUT) :: PSFTQ     ! flux of water vapor                   (kg/m2/s)
!
!*      0.2    declarations of local variables
!
!* forcing variables
!
REAL, DIMENSION(SIZE(PTA))   :: ZWIND    ! lowest atmospheric level wind speed           (m/s)
REAL, DIMENSION(SIZE(PTA))   :: ZEXNA    ! Exner function at lowest SBL scheme level     (-)
REAL, DIMENSION(SIZE(PTA))   :: ZQA      ! specific humidity                             (kg/m3)
!
! SBL turbulence scheme
!
REAL, DIMENSION(SIZE(PTA))   :: ZUSTAR       ! friction velocity (m/s)
!
REAL, DIMENSION(SIZE(PTA))   :: ZEXNS
REAL, DIMENSION(SIZE(PTA))   :: ZQSAT
REAL, DIMENSION(SIZE(PTA))   :: ZCD
REAL, DIMENSION(SIZE(PTA))   :: ZCDN
REAL, DIMENSION(SIZE(PTA))   :: ZCH
REAL, DIMENSION(SIZE(PTA))   :: ZRI
REAL, DIMENSION(SIZE(PTA))   :: ZRESA_SEA
REAL, DIMENSION(SIZE(PTA))   :: ZZ0H
REAL, DIMENSION(SIZE(PTA))   :: ZCLS_WIND_ZON
REAL, DIMENSION(SIZE(PTA))   :: ZCLS_WIND_MER
REAL, DIMENSION(SIZE(PTA))   :: ZTNM
REAL, DIMENSION(SIZE(PTA))   :: ZQNM
REAL, DIMENSION(SIZE(PTA))   :: ZHUNM
REAL, DIMENSION(SIZE(PTA))   :: ZHU
INTEGER                      :: J
INTEGER                      :: JLAYER
LOGICAL                      :: GHANDLE_SIC = .false. ! no sea-ice model
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------------
!
!*     1.2     Initialisation at first time step
!              ---------------------------------
!
IF (lhook) CALL dr_hook('INIT_WATER_SBL',0,zhook_handle)
!
zexna(:) = (ppa(:)/xp00)**(xrd/xcpd)
zexns(:) = (pps(:)/xp00)**(xrd/xcpd)
zqa(:) = pqa(:) / prhoa(:)
zwind(:) = sqrt(pu**2+pv**2)
!
!We choose the case CSEA_FLUX=DIRECT to compute CD, CDN, CH, RI and ZZ0H
!Iterative computation of PZ0 / CD, CDN, CH, ZH0
DO j=1,5
  CALL water_flux(pz0, pta, zexna, prhoa, pts, zexns, zqa, prain, &
                  psnow, xtt, zwind, pzref, puref,                &
                  pps, ghandle_sic, zqsat,  psfth, psftq, zustar, &
                  zcd, zcdn, zch, zri, zresa_sea, zz0h            )
ENDDO
!
!Initialisation of T, Q, Wind and TKE on all canopy levels
zhu(:)=1.
DO jlayer=1,sb%NLVL
  !
  CALL cls_tq(pta, zqa, ppa, pps, pzref, zcd, zch, zri, pts, zhu, zz0h, &
              sb%XZ(:,jlayer), ztnm, zqnm, zhunm                        )
  !
  sb%XT(:,jlayer)=ztnm
  sb%XQ(:,jlayer)=zqnm
  !
  CALL cls_wind(pu, pv, puref, zcd, zcdn, zri, sb%XZ(:,jlayer), &
                zcls_wind_zon, zcls_wind_mer                 )
  !
  sb%XU   (:,jlayer) = sqrt( zcls_wind_zon(:)**2 + zcls_wind_mer(:)**2 )
  sb%XTKE (:,jlayer) = xalpsbl * zcd(:) * ( pu(:)**2 + pv(:)**2 )
  sb%XP   (:,jlayer) = ppa(:) + xg * prhoa(:) * ( sb%XZ(:,sb%NLVL) - sb%XZ(:,jlayer) )
  !
ENDDO
!
IF (lhook) CALL dr_hook('INIT_WATER_SBL',1,zhook_handle)
!
END SUBROUTINE init_water_sbl