init_snow_lw.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_snow_lw(PEMISSN,TPSNOW)
!     #######################################
!
!!****  *INIT_SNOW_LW* - routine to initialize snow surf. temp and emissivity
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2003 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_type_snow
USE modd_surf_par,   ONLY : xundef
USE modd_csts,       ONLY : xtt, xlmtt
!
USE mode_snow3l
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
REAL,                 INTENT(IN)    :: PEMISSN ! snow emissivity
TYPE(surf_snow),      INTENT(INOUT) :: TPSNOW  ! snow characteristics
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('INIT_SNOW_LW',0,zhook_handle)
!
!*       1.    Emissivity
!              ----------
!
IF (tpsnow%SCHEME=='1-L' .OR. tpsnow%SCHEME=='3-L' .OR. tpsnow%SCHEME=='CRO') THEN
  WHERE(tpsnow%WSNOW(:,1)==0. .OR. tpsnow%WSNOW(:,1)==xundef )
    tpsnow%EMIS (:)= xundef
  ELSEWHERE
    tpsnow%EMIS (:)= pemissn
  END WHERE
END IF
!
!*      2.     Surface temperature
!              -------------------
!
IF (tpsnow%SCHEME=='1-L') THEN
  WHERE(tpsnow%WSNOW(:,1)==0. .OR. tpsnow%WSNOW(:,1)==xundef )
    tpsnow%TS (:)= xundef
  ELSEWHERE
    tpsnow%TS(:) = tpsnow%T(:, tpsnow%NLAYER)
  END WHERE
END IF
!
IF (tpsnow%SCHEME=='3-L' .OR. tpsnow%SCHEME=='CRO') THEN
  WHERE(tpsnow%WSNOW(:,1)==0. .OR. tpsnow%WSNOW(:,1)==xundef )
    tpsnow%TS (:)= xundef
  ELSEWHERE
    tpsnow%TS(:) = xtt + (tpsnow%HEAT(:,1) + xlmtt * tpsnow%RHO(:,1)) / snow3lscap(tpsnow%RHO(:,1))  
    tpsnow%TS(:) = min(xtt, tpsnow%TS(:))
  END WHERE
END IF
!
IF (lhook) CALL dr_hook('INIT_SNOW_LW',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE init_snow_lw