read_flake_confn.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 read_flake_conf_n (CHF, DGF, DGMF, F, &
                                    hprogram)
!     #############################################################
!
!!****  *READ_FLAKE_CONF* - reads the configuration for FLAKE
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2003 
!!      Modified    04/2013, P. Le Moigne: FLake chemistry
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
!
!
USE modd_ch_flake_n, ONLY : ch_flake_t
USE modd_diag_flake_n, ONLY : diag_flake_t
USE modd_diag_misc_flake_n, ONLY : diag_misc_flake_t
USE modd_flake_n, ONLY : flake_t
!
USE mode_modeln_surfex_handler
!
USE mode_pos_surf
!
USE modi_test_nam_var_surf
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
!
USE modn_flake_n
!
USE modd_read_namelist, ONLY : lnam_read
USE modd_surf_par,   ONLY : xundef
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(ch_flake_t), INTENT(INOUT) :: chf
TYPE(diag_flake_t), INTENT(INOUT) :: dgf
TYPE(diag_misc_flake_t), INTENT(INOUT) :: dgmf
TYPE(flake_t), INTENT(INOUT) :: f
!
 CHARACTER(LEN=6),  INTENT(IN)  :: hprogram ! program calling ISBA

!
!*       0.2   Declarations of local variables
!              -------------------------------
!
LOGICAL           :: gfound         ! Return code when searching namelist
INTEGER           :: iluout         ! logical unit of output file
INTEGER           :: inam           ! logical unit of namelist file
INTEGER           :: imi
INTEGER           :: ilu
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!* get output listing file logical unit
!
!
IF (lhook) CALL dr_hook('READ_FLAKE_CONF_N',0,zhook_handle)
 CALL get_luout(hprogram,iluout)
!
imi=get_current_model_index_surfex()
!
IF (imi.NE.-1 .AND. lnam_read) THEN
 CALL init_nam_flaken(f)
 CALL init_nam_diag_surfn(dgf)
 CALL init_nam_diag_flaken(dgmf)
 CALL init_nam_ch_flaken(chf)
ENDIF
!
IF (lnam_read) THEN
        
 !* open namelist file
 !
 CALL open_namelist(hprogram,inam)
 !
 !* reading of namelist
 !  -------------------
 !
 CALL posnam(inam,'NAM_FLAKEN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_flaken)
 !
 CALL test_nam_var_surf(iluout,'CSNOW_FLK',csnow_flk,'DEF')
 CALL test_nam_var_surf(iluout,'CFLK_FLUX',cflk_flux,'DEF  ','FLAKE')
 CALL test_nam_var_surf(iluout,'CFLK_ALB',cflk_alb,'UNIF','TA96','MK10')
 !
 CALL posnam(inam,'NAM_DIAG_SURFN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_diag_surfn)
 CALL posnam(inam,'NAM_DIAG_FLAKEN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_diag_flaken)
 IF (lwater_profile .AND. count(xzwat_profile /= xundef) == 0) &
        CALL abor1_sfx("XZWAT_PROFILE MUST BE DEFINED IN NAMELIST NAM_DIAG_FLAKEN IF LWATER_PROFILE=T")
 CALL posnam(inam,'NAM_CH_FLAKEN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_ch_flaken)
!
 CALL test_nam_var_surf(iluout,'CCH_DRY_DEP',cch_dry_dep,'      ','WES89 ','NONE  ')
 !
 !* close namelist file
 !
 CALL close_namelist(hprogram,inam)
 !
ENDIF
!
IF (imi.NE.-1) THEN
 CALL update_nam_flaken(f)
 CALL update_nam_diag_surfn(dgf)
 CALL update_nam_diag_flaken(dgmf)
 CALL update_nam_ch_flaken(chf)
ENDIF
!
!-------------------------------------------------------------------------------
!
!* surface time-step forced by the atmosphere
!
xtstep = xundef
IF (lhook) CALL dr_hook('READ_FLAKE_CONF_N',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_flake_conf_n