read_isba_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_isba_conf_n (CHI, DGEI, DGI, DGMI, I, &
                                   hprogram)
!     #######################################################
!
!!****  *READ_ISBA_CONF* - routine to read the configuration for ISBA
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2003 
!!      P Le Moigne 09/2005 AGS modifs of L. Jarlan
!!      P Le Moigne 09/2005 CSNOWRES option
!!      Modified by P. Le Moigne (06/2006): seeding and irrigation
!!      Modified by P. Le Moigne (05/2008): deep soil characteristics
!!      Modified by R. El Khatib 05-Apr-2012 Fix message handling
!!      Modified by C. Carmagnola (01/2013): CSNOWMETAMO = 'B92','C13','F06','T07'
!!      B. Decharme  04/2013 delete CTOPREG
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
USE modd_ch_isba_n, ONLY : ch_isba_t
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
USE modd_diag_isba_n, ONLY : diag_isba_t
USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_t
USE modd_isba_n, ONLY : isba_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_isba_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_get_surf_atm_sso_rough
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(ch_isba_t), INTENT(INOUT) :: chi
TYPE(diag_evap_isba_t), INTENT(INOUT) :: dgei
TYPE(diag_isba_t), INTENT(INOUT) :: dgi
TYPE(diag_misc_isba_t), INTENT(INOUT) :: dgmi
TYPE(isba_t), INTENT(INOUT) :: i
!
 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
 CHARACTER(LEN=4)  :: yrough         ! CROUGH from surf_atm
INTEGER           :: imi
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!* get output listing file logical unit
!
IF (lhook) CALL dr_hook('READ_ISBA_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_isban(i)
 CALL init_nam_isba_agsn(i)
 CALL init_nam_sgh_isban(i)
 CALL init_nam_diag_isban(dgei, dgi, dgmi)
 CALL init_nam_diag_surfn(dgei, dgi)
 CALL init_nam_ch_controln(chi)
 CALL init_nam_ch_isban(chi)
 CALL init_nam_spinup_carb_isban(i)
 CALL init_nam_isba_snown(i) 
ENDIF
!
IF (lnam_read) THEN
 !
 !* open namelist file
 !
 CALL open_namelist(hprogram,inam)
 !
 !* reading of namelist
 !  -------------------
 !
 CALL posnam(inam,'NAM_ISBAN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_isban)
 CALL posnam(inam,'NAM_ISBA_AGSN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_isba_agsn) 
 CALL posnam(inam,'NAM_SGH_ISBAN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_sgh_isban)
 CALL posnam(inam,'NAM_DIAG_SURFN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_diag_surfn)
 CALL posnam(inam,'NAM_DIAG_ISBAN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_diag_isban)
 CALL posnam(inam,'NAM_CH_CONTROLN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_ch_controln)
 CALL posnam(inam,'NAM_CH_ISBAN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_ch_isban)
  CALL posnam(inam,'NAM_SPINUP_CARBN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_spinup_carbn)
 CALL posnam(inam,'NAM_ISBA_SNOWN',gfound,iluout)
 IF (gfound) READ(unit=inam,nml=nam_isba_snown) 
 !
 CALL test_nam_var_surf(iluout,'CROUGH',crough,'NONE','Z01D','Z04D','BE04','UNDE')
 CALL test_nam_var_surf(iluout,'CSCOND',cscond,'NP89','PL98')
 CALL test_nam_var_surf(iluout,'CALBEDO',calbedo,'EVOL','DRY ','WET ','MEAN','USER','CM13')
 CALL test_nam_var_surf(iluout,'CC1DRY',cc1dry,'DEF ','GB93')
 CALL test_nam_var_surf(iluout,'CSOILFRZ',csoilfrz,'DEF','LWT')
 CALL test_nam_var_surf(iluout,'CDIFSFCOND',cdifsfcond,'DEF ','MLCH')
 CALL test_nam_var_surf(iluout,'CSNOWRES',csnowres,'DEF','RIL')
 CALL test_nam_var_surf(iluout,'CCPSURF',ccpsurf,'DRY','HUM')
 !
 CALL test_nam_var_surf(iluout,'CRUNOFF',crunoff,'WSAT','DT92','SGH ','TOPD')
 CALL test_nam_var_surf(iluout,'CKSAT',cksat,'DEF','SGH','EXP')
 CALL test_nam_var_surf(iluout,'CRAIN',crain,'DEF','SGH')
 CALL test_nam_var_surf(iluout,'CHORT',chort,'DEF','SGH')
 !
 CALL test_nam_var_surf(iluout,'CCH_DRY_DEP',cch_dry_dep,'      ','WES89 ','NONE  ')
 !
 CALL test_nam_var_surf(iluout,'CSNOWMETAMO',csnowmetamo,'B92','C13','F06','T07')
 CALL test_nam_var_surf(iluout,'CSNOWRAD',csnowrad,'B92','TAR','TA1','TA2')
 ! 
 !* close namelist file
 !
 CALL close_namelist(hprogram,inam)
 !
ENDIF
!
IF (imi.NE.-1) THEN
 CALL update_nam_isban(i)
 CALL update_nam_isba_agsn(i)
 CALL update_nam_sgh_isban(i)
 CALL update_nam_diag_isban(dgei, dgi, dgmi)
 CALL update_nam_diag_surfn(dgei, dgi)
 CALL update_nam_ch_controln(chi)
 CALL update_nam_ch_isban(chi)
 CALL update_nam_spinup_carb_isban(i)
 CALL update_nam_isba_snown(i) 
ENDIF
!
!-------------------------------------------------------------------------------
!
!* surface time-step forced by the atmosphere
!
!XTSTEP = XUNDEF
!
!-------------------------------------------------------------------------------
!
!* coherence check : orographic drag must be done only once
!
 CALL get_surf_atm_sso_rough(yrough)
!
IF (yrough/="NONE") THEN
  IF (lnam_read) THEN
    WRITE(iluout,*) " -------------------------------------------------------------------------- "
    WRITE(iluout,*) " Orographic roughness is done in SURF_ATM, so it should not be done in ISBA "
    WRITE(iluout,*) " CROUGH in NAM_ISBAn is set to 'NONE' "
    WRITE(iluout,*) " "
    WRITE(iluout,*) " -------------------------------------------------------------------------- "
  ENDIF
  crough = "NONE"
END IF
IF (lhook) CALL dr_hook('READ_ISBA_CONF_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_isba_conf_n