read_namelists_isban.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_namelists_isba_n (IM, HPROGRAM, HINIT)
!     #######################################################
!
!---------------------------    
!
!
USE modd_surfex_n, ONLY : isba_model_t
!
USE modd_surf_par, ONLY : xundef
!
USE modn_isba_n                      
!    
USE modn_soiltemp_arp
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
USE mode_pos_surf
!
USE modi_default_isba
USE modi_default_ch_dep
USE modi_default_ch_bio_flux
USE modi_default_diag_isba
USE modi_default_crocus
!
USE modi_read_default_isba_n
USE modi_read_isba_conf_n
!
USE modi_read_nam_prep_isba_n
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!
TYPE(isba_model_t), INTENT(INOUT) :: IM
!
 CHARACTER(LEN=6),   INTENT(IN)  :: HPROGRAM  ! program calling surf. schemes
 CHARACTER(LEN=3),   INTENT(IN)  :: HINIT     ! choice of fields to initialize
!
INTEGER :: ILUNAM         ! logical unit of namelist file
INTEGER :: ILUOUT
LOGICAL :: GFOUND         ! Return code when searching namelist
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!------------------------------------
!
IF (lhook) CALL dr_hook('READ_NAMELISTS_ISBA_N',0,zhook_handle)
 CALL default_isba(xtstep, xout_tstep,                           &
                     crunoff, cscond,                            &
                     cc1dry, csoilfrz, cdifsfcond, csnowres,     &
                     ccpsurf, xcgmax, xcdrag, cksat, lsoc,       &
                     crain, chort, lglacier, lcanopy_drag,       &
                     lvegupd, lspinupcarbs, lspinupcarbw,        &
                     xspinmaxs, xspinmaxw, xco2_start, xco2_end, &
                     nnbyearspins, nnbyearspinw, lnitro_dilu     )
!
 CALL default_ch_dep(cch_dry_dep)
 CALL default_ch_bio_flux(lch_bio_flux) 
!           
 CALL default_diag_isba(n2m,lsurf_budget,l2m_min_zs,lrad_budget,  &
                       lcoef,lsurf_vars,lsurf_evap_budget,        &
                       lsurf_misc_budget,lsurf_diag_albedo,       &
                       lsurf_budgetc,lsurf_misc_dif,lpatch_budget,&
                       lpgd,lreset_budgetc,lwater_budget,lprosnow,&
                       lvolumetric_snowliq,xdiag_tstep           ) 
!      
 CALL default_crocus(lsnowdrift,lsnowdrift_sublim,lsnow_abs_zenith,&
                     csnowmetamo,csnowrad)
!
 CALL read_default_isba_n(im%CHI, im%ID%DE, im%ID%O, im%ID%DM, im%O, hprogram)
!
 CALL read_isba_conf_n(im%CHI, im%ID%DE, im%ID%O, im%ID%DM, im%O, hprogram)
!
!
sodelx(:)      = xundef
ltemp_arp      = .false.
ntemplayer_arp = 4
!
 CALL get_luout(hprogram,iluout)
!
 CALL open_namelist(hprogram,ilunam)
 CALL posnam(ilunam,'NAM_SOILTEMP_ARP',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_soiltemp_arp)
 CALL close_namelist(hprogram,ilunam)
!
IF (hinit=='PRE') CALL read_nam_prep_isba_n(hprogram)      
IF (lhook) CALL dr_hook('READ_NAMELISTS_ISBA_N',1,zhook_handle)
!
!
END SUBROUTINE read_namelists_isba_n