init_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 init_cpl_gcm_n (U, HPROGRAM,HINIT)
!     ########################################
!
!!****  *INIT_CPL_GCM_n* - routine to read  physical fields into  
!!                         the restart file for ARPEGE/ALADIN run
!!
!!    PURPOSE
!!    -------
!!       The purpose of this routine is to initialise some 
!!       physical fields. Indeed, when ARPEGE/ALADIN is used, 
!!       these field is not initialize at the begin of a run.
!!
!!
!!**  METHOD
!!    ------
!!      The data are read in the initial surface file :
!!        - 2D data fields
!!          
!!      It does not read the grid definition. This should have been
!!      read already.
!!
!!    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 modi_read_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
!
!
TYPE(surf_atm_t), INTENT(INOUT) :: U
!
 CHARACTER(LEN=6),  INTENT(IN)  :: HPROGRAM ! calling program
 CHARACTER(LEN=3),  INTENT(IN)  :: HINIT    ! choice of fields to initialize
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: IRESP      ! Error code after redding
! 
 CHARACTER(LEN=12) :: YRECFM     ! Name of the article to be read
!
INTEGER           :: IVERSION   ! surface version
!
LOGICAL           :: LREAD      ! work key
LOGICAL           :: GCPL_GCM   ! work key
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('INIT_CPL_GCM_N',0,zhook_handle)
!
gcpl_gcm = .false.
!
IF (hinit=='PGD') THEN
!     
   ALLOCATE(u%XRAIN (0))
   ALLOCATE(u%XSNOW (0))           
   ALLOCATE(u%XZ0   (0))           
   ALLOCATE(u%XZ0H  (0))           
   ALLOCATE(u%XQSURF(0))           
!     
ELSE
!
   ALLOCATE(u%XRAIN (u%NSIZE_FULL))
   ALLOCATE(u%XSNOW (u%NSIZE_FULL))
   ALLOCATE(u%XZ0   (u%NSIZE_FULL))
   ALLOCATE(u%XZ0H  (u%NSIZE_FULL))
   ALLOCATE(u%XQSURF(u%NSIZE_FULL))
!
   u%XRAIN (:) = 0.0
   u%XSNOW (:) = 0.0
   u%XZ0   (:) = 0.001
   u%XZ0H  (:) = 0.001
   u%XQSURF(:) = 0.0
!
ENDIF
!
yrecfm='VERSION'
 CALL read_surf(hprogram,yrecfm,iversion,iresp)
!
lread=(hinit/='PGD'.AND.hinit/='PRE'.AND.iversion>=8)
!
IF (lread) THEN
   yrecfm='LCPL_GCM'
   CALL read_surf(hprogram,yrecfm,gcpl_gcm,iresp)
ENDIF  
!
IF (lread.AND.gcpl_gcm) THEN
!
   yrecfm='RAIN_GCM'
   CALL read_surf(hprogram,yrecfm,u%XRAIN(:),iresp)
!
   yrecfm='SNOW_GCM'
   CALL read_surf(hprogram,yrecfm,u%XSNOW(:),iresp)
!
   yrecfm='Z0_GCM'
   CALL read_surf(hprogram,yrecfm,u%XZ0(:),iresp)
!
   yrecfm='Z0H_GCM'
   CALL read_surf(hprogram,yrecfm,u%XZ0H(:),iresp)
!
   yrecfm='QS_GCM'
   CALL read_surf(hprogram,yrecfm,u%XQSURF(:),iresp)
!
ENDIF        
!
IF (lhook) CALL dr_hook('INIT_CPL_GCM_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE init_cpl_gcm_n