prep_flake_unif.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 prep_flake_unif(KLUOUT,HSURF,PFIELD)
!     #################################################################################
!
!!****  *PREP_FLAKE_UNIF* - prepares FLAKE field from prescribed values
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     S. Malardel
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!    09/2010 E. Kourzeneva: Renamed the lake surface temperature 
!!                           from the namelist
!!------------------------------------------------------------------
!
!
USE modd_prep,       ONLY : cinterp_type
USE modd_prep_flake,   ONLY : xts_unif, &
                                  xunif_t_snow  , &
                                  xunif_t_ice   , &
                                  xunif_t_mnw   , &
                                  xunif_t_wml   , &
                                  xunif_t_bot   , &
                                  xunif_t_b1    , &
                                  xunif_ct      , &
                                  xunif_h_snow  , &
                                  xunif_h_ice   , &
                                  xunif_h_ml    , &
                                  xunif_h_b1         
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
INTEGER,            INTENT(IN)  :: KLUOUT    ! output listing logical unit
 CHARACTER(LEN=7),   INTENT(IN)  :: HSURF     ! type of field
REAL, POINTER, DIMENSION(:,:)   :: PFIELD    ! field to interpolate horizontally
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!*      0.2    declarations of local variables
!
!
!-------------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PREP_FLAKE_UNIF',0,zhook_handle)
SELECT CASE(hsurf)
!
!*      3.0    Orography
!
  CASE('ZS     ')
    ALLOCATE(pfield(1,1))
    pfield = 0.
!
!*      3.1    FLake variables
!
  CASE('TS     ')
    ALLOCATE(pfield(1,1))
    pfield = xts_unif
!
  CASE('T_SNOW ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_t_snow
!
  CASE('T_ICE  ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_t_ice
!
  CASE('T_MNW  ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_t_mnw
!
  CASE('T_WML  ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_t_wml
!
  CASE('T_BOT  ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_t_bot
!
  CASE('T_B1   ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_t_b1
!
  CASE('CT     ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_ct
!
  CASE('H_SNOW ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_h_snow
!
  CASE('H_ICE  ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_h_ice
!
  CASE('H_ML   ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_h_ml
!
  CASE('H_B1   ')
    ALLOCATE(pfield(1,1))
    pfield = xunif_h_b1
!
!
END SELECT
!
!*      4.     Interpolation method
!              --------------------
!
cinterp_type='UNIF  '
IF (lhook) CALL dr_hook('PREP_FLAKE_UNIF',1,zhook_handle)
!
!
!-------------------------------------------------------------------------------------
END SUBROUTINE prep_flake_unif