init_surf_topd.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_surf_topd (DEC, IO, S, K, NP, NPE, UG, U, HPROGRAM,KI)
!     #############################################################
!
!!****  *INIT_SURF_TOPD* - routine to initialize variables needed for coupling with Topmodel
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!    The routine open and read the namelists NAM_COUPL_TOPD and NAM_TOPD,
!! calculates the number of catchments concerned, the different time step 
!! variables and all the variables nedded for coupling with Topmodel. 
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      B. Vincendon   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    11/2006
!!      Modif 04/2007: Arguments PTOPD_STEP,KNB_TOPD_STEP become module
!!                     variables from MODD_TOPDDYN
!!      Modif 03/2014: New organisation of routines with init_topd_ol and
!!                     displacement of init_budget
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY : isba_s_t, isba_k_t, isba_np_t, isba_npe_t
!
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_surfex_mpi, ONLY : nproc
USE modd_surfex_omp, ONLY : nblocktot
!
USE modd_topodyn, ONLY :ccat, xspeedr, xspeedh, nncat, &
                        xrtop_d2, xspeedg
USE modd_coupling_topd, ONLY :  lcoupl_topd, nnb_topd, lbudget_topd
!
!
USE modi_get_luout
USE modi_abor1_sfx
!
USE modi_init_topd_ol
USE modi_init_coupl_topd
USE modi_init_budget_coupl_rout
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEC
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(isba_s_t), INTENT(INOUT) :: S
TYPE(isba_k_t), INTENT(INOUT) :: K
TYPE(isba_np_t), INTENT(INOUT) :: NP
TYPE(isba_npe_t), INTENT(INOUT) :: NPE
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
TYPE(surf_atm_t), INTENT(INOUT) :: U
!
 CHARACTER(LEN=*),  INTENT(IN)     :: HPROGRAM      !
INTEGER,           INTENT(IN)     :: KI            ! grid dimension
!
INTEGER :: ILUOUT
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('INIT_SURF_TOPD',0,zhook_handle)
!
IF (lcoupl_topd) THEN
  IF (nproc>1) CALL abor1_sfx('INIT_SURF_TOPD: TOPD CANNOT RUN WITH MORE THAN 1 MPI TASK') 
  IF (nblocktot>1) CALL abor1_sfx("INIT_SURF_TOPD: TOPD CANNOT RUN WITH NUMEROUS OPENMP BLOCKS")
  IF (npe%AL(1)%TSNOW%SCHEME/='3-L') &
        CALL abor1_sfx("INIT_SURF_TOPD: coupling with topmodel only runs with TSNOW%SCHEME=3-L")
ENDIF
!
 CALL get_luout(hprogram,iluout)
!            
!         1.   Reads the namelists
!              --------------------
!
WRITE(iluout,*) 'Debut init_surf_topo_n'
!
IF (lcoupl_topd) THEN
  !
  !         3.   Initialises variables specific to Topmodel
  !              -------------------------------------------
  WRITE(iluout,*) 'NNCAT',nncat
  !
  CALL init_topd_ol(hprogram)
  !
  !         4.   Initialises variables nedded for coupling with Topmodel
  !              -------------------------------------------------------
  !
  CALL init_coupl_topd(dec, io, s, k, np, npe, ug, u, hprogram)
  !
  WRITE(iluout,*) 'Couplage avec TOPMODEL active'
  !
  !IF (LBUDGET_TOPD) CALL INIT_BUDGET_COUPL_ROUT(KI)
  !
ELSE
  !
  WRITE(iluout,*) 'Pas de couplage avec TOPMODEL'
  !
ENDIF
!
IF (lhook) CALL dr_hook('INIT_SURF_TOPD',1,zhook_handle)
!
END SUBROUTINE init_surf_topd