get_cpl_gcmn.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 get_cpl_gcm_n (U, &
                                hprogram,ki,prain,psnow,pz0,pz0h,pqsurf)
!     ######################################################
!
!!****  *GET_CPL_GCM_n* - routine to get physical fields   
!!                      for initialise ARPEGE/ALADIN run
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!      
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      B. Decharme   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    04/2013
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_surf_atm,   ONLY : lcpl_gcm
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_get_luout
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
 CHARACTER(LEN=6),        INTENT(IN)  :: hprogram
INTEGER,                 INTENT(IN)  :: ki       ! number of points
!
REAL, DIMENSION(KI),     INTENT(OUT) :: prain    ! total rainfall rate (kg/m2/s)
REAL, DIMENSION(KI),     INTENT(OUT) :: psnow    ! total snowfall rate (kg/m2/s)
REAL, DIMENSION(KI),     INTENT(OUT) :: pz0      ! roughness length for momentum (m)
REAL, DIMENSION(KI),     INTENT(OUT) :: pz0h     ! roughness length for heat (m)
REAL, DIMENSION(KI),     INTENT(OUT) :: pqsurf   ! specific humidity at surface (kg/kg)
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER :: iluout
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('GET_CPL_GCM_N',0,zhook_handle)
!
 CALL get_luout(hprogram,iluout)
!
IF(lcpl_gcm) THEN
!
  IF(SIZE(pz0)/=SIZE(u%XZ0H))THEN
    WRITE(iluout,*)'try to get Z0 field from atmospheric model, but size is not correct'
    WRITE(iluout,*)'size of field expected by the atmospheric model (PZ0) :', SIZE(pz0)
    WRITE(iluout,*)'size of field in SURFEX                         (XZ0) :', SIZE(u%XZ0)
    CALL abor1_sfx('GET_CPL_GCMN: PZ0 ARGUMENT SIZE /= XZ0 MODULE SIZE')
  ENDIF
!
  prain(:) = u%XRAIN
  psnow(:) = u%XSNOW
  pz0(:) = u%XZ0
  pz0h(:) = u%XZ0H
  pqsurf(:) = u%XQSURF
!
ELSE
!
  WRITE(iluout,*)'LCPL_GCM must be TRUE when you use atmospheric model'
  CALL abor1_sfx('GET_CPL_GCMN: LCPL_GCM must be TRUE')
!  
ENDIF

!
IF (lhook) CALL dr_hook('GET_CPL_GCM_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE get_cpl_gcm_n