read_prep_teb_garden_conf.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_prep_teb_garden_conf(HPROGRAM,HVAR,HFILE,HFILETYPE,HFILEPGD,HFILEPGDTYPE,&
                                           hatmfile,hatmfiletype,hpgdfile,hpgdfiletype,kluout,ounif)
!     #######################################################
!
!!****  *READ_PREP_TEB_GARDEN_CONF* - routine to read the configuration for ISBA
!!                              fields preparation
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
USE modi_read_prep_surf_atm_conf
!
USE modn_prep_teb_garden
USE modd_prep_teb_garden, ONLY : cfile_gd, ctype, cfilepgd_gd, ctypepgd,      &
                                 cfile_hug_gd, ctype_hug,                            &
                                 cfile_hug_surf_gd, cfile_hug_root_gd, cfile_hug_deep_gd,  &
                                 xhug_surf_gd, xhug_root_gd, xhug_deep_gd,                 &
                                 xhugi_surf_gd, xhugi_root_gd, xhugi_deep_gd,              &
                                 cfile_tg_gd, ctype_tg,                              &
                                 cfile_tg_surf_gd, cfile_tg_root_gd, cfile_tg_deep_gd,     &
                                 xtg_surf_gd, xtg_root_gd, xtg_deep_gd
!
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
!              -------------------------
!
 CHARACTER(LEN=6),  INTENT(IN)  :: hprogram    ! program calling
 CHARACTER(LEN=7),  INTENT(IN)  :: hvar        ! variable treated
 CHARACTER(LEN=28), INTENT(OUT) :: hfile       ! file name
 CHARACTER(LEN=6),  INTENT(OUT) :: hfiletype   ! file type
 CHARACTER(LEN=28), INTENT(OUT) :: hfilepgd    ! file name
 CHARACTER(LEN=6),  INTENT(OUT) :: hfilepgdtype! file type
 CHARACTER(LEN=28), INTENT(IN)  :: hatmfile    ! atmospheric file name
 CHARACTER(LEN=6),  INTENT(IN)  :: hatmfiletype! atmospheric file type
 CHARACTER(LEN=28), INTENT(IN)  :: hpgdfile    ! atmospheric file name
 CHARACTER(LEN=6),  INTENT(IN)  :: hpgdfiletype! atmospheric file type
INTEGER,           INTENT(IN)  :: kluout      ! logical unit of output listing
LOGICAL,           INTENT(OUT) :: ounif       ! flag for prescribed uniform field

!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: iresp          ! IRESP  : return-code if a problem appears 
                                    ! at the open of the file in LFI  routines 
INTEGER           :: ilunam         ! Logical unit of namelist file
!
 CHARACTER(LEN=28) :: ynamelist      ! namelist file
!
LOGICAL           :: gfound         ! Return code when searching namelist
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!
IF (lhook) CALL dr_hook('READ_PREP_TEB_GARDEN_CONF',0,zhook_handle)
hfile    = '                         '
hfiletype    = '      '
!
hfilepgd = '                         '
hfilepgdtype = '      '
!
ounif    = .false.
!
!-------------------------------------------------------------------------------
!
!* choice of input file
!  --------------------
!
SELECT CASE (hvar)
  CASE ('WG     ','WGI    ')
    IF (len_trim(cfile_hug_gd)>0 .AND. len_trim(ctype_hug)>0 ) THEN
      hfile     = cfile_hug_gd
      hfiletype = ctype_hug
    END IF
  CASE ('TG     ')
    IF (len_trim(cfile_tg_gd)>0 .AND. len_trim(ctype_tg)>0 ) THEN
      hfile     = cfile_tg_gd
      hfiletype = ctype_tg
    END IF
END SELECT
!
IF (len_trim(hfile)==0 .AND. len_trim(cfile_gd)>0 .AND. len_trim(ctype)>0) THEN
  hfile     = cfile_gd
  hfiletype = ctype
END IF
!
IF (len_trim(hfilepgd)==0 .AND. len_trim(cfilepgd_gd)>0 .AND. len_trim(ctypepgd)>0) THEN
  hfilepgd     = cfilepgd_gd
  hfilepgdtype = ctypepgd
END IF
!
!! If no file name in the scheme namelist,
!! try to find a name in NAM_SURF_ATM
!
IF (len_trim(hfile)==0) THEN
!
 CALL read_prep_surf_atm_conf(hprogram,hfile,hfiletype,hfilepgd,hfilepgdtype,&
                             hatmfile,hatmfiletype,hpgdfile,hpgdfiletype,kluout)
!
END IF
!
!! If no file name in the scheme namelist,
!! nor in NAM_SURF_ATM, look if ascii input files are present
!
SELECT CASE (hvar)
  CASE ('WG     ','WGI    ')
    IF ( len_trim(ctype_hug )>0       .AND. &
           len_trim(cfile_hug_surf_gd)>0   .AND. &
           len_trim(cfile_hug_root_gd)>0   .AND. &
           len_trim(cfile_hug_deep_gd)>0         ) THEN  
       hfiletype = ctype_hug 
    END IF
    IF (hvar=='WGI    ' .AND. hfiletype=='ASCLLV') THEN
       ounif = .true.
       IF (lhook) CALL dr_hook('READ_PREP_TEB_GARDEN_CONF',1,zhook_handle)
       RETURN
    ENDIF
  CASE ('TG     ')
    IF ( len_trim(ctype_tg )>0       .AND. &
           len_trim(cfile_tg_surf_gd)>0   .AND. &
           len_trim(cfile_tg_root_gd)>0   .AND. &
           len_trim(cfile_tg_deep_gd)>0         ) THEN  
       hfiletype = ctype_tg 
    END IF
END SELECT
!
!-------------------------------------------------------------------------------
!
!* Is an uniform field prescribed?
!  ------------------------------
!
SELECT CASE (hvar)
  CASE ('WG     ')
    ounif = (xhug_surf_gd/=xundef) .OR. (xhug_root_gd/=xundef) .OR. (xhug_deep_gd/=xundef)
    IF (ounif .AND. (xhug_surf_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XHUG_SURF_GD, XHUG_ROOT_GD OR XHUG_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XHUG_SURF_GD MUST BE SET')
    END IF
    IF (ounif .AND. (xhug_root_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XHUG_SURF_GD, XHUG_ROOT_GD OR XHUG_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XHUG_ROOT_GD MUST BE SET')
    END IF
    IF (ounif .AND. (xhug_deep_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XHUG_SURF_GD, XHUG_ROOT_GD OR XHUG_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XHUG_DEEP_GD MUST BE SET')
    END IF
     !
      CASE ('WGI    ')
    ounif = (xhugi_surf_gd/=xundef) .OR. (xhugi_root_gd/=xundef) .OR. (xhugi_deep_gd/=xundef)
    IF (ounif .AND. (xhugi_surf_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XHUGI_SURF_GD, XHUGI_ROOT_GD OR XHUGI_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XHUGI_SURF_GD MUST BE SET')
    END IF
    IF (ounif .AND. (xhugi_root_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XHUGI_SURF_GD, XHUGI_ROOT_GD OR XHUGI_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XHUGI_ROOT_GD MUST BE SET')
    END IF
    IF (ounif .AND. (xhugi_deep_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XHUGI_SURF_GD, XHUGI_ROOT_GD OR XHUGI_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XHUGI_DEEP_GD MUST BE SET')
    END IF
    !
  CASE ('TG     ')
    ounif = (xtg_surf_gd/=xundef)  .OR. (xtg_root_gd/=xundef)  .OR. (xtg_deep_gd/=xundef)
    IF (ounif .AND. (xtg_surf_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XTG_SURF_GD, XTG_ROOT_GD OR XTG_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XTG_SURF_GD MUST BE SET')
    END IF
    IF (ounif .AND. (xtg_root_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XTG_SURF_GD, XTG_ROOT_GD OR XTG_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XTG_ROOT_GD MUST BE SET')
    END IF
    IF (ounif .AND. (xtg_deep_gd==xundef)) THEN
       WRITE(kluout,*)'ONE OF XTG_SURF_GD, XTG_ROOT_GD OR XTG_DEEP_GD IS GIVEN'
       CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: XTG_DEEP_GD MUST BE SET')
    END IF
    !

END SELECT
!
!-------------------------------------------------------------------------------
!
!* If no file and no uniform field is prescribed: default values used
!  ---------------------------------------------
!
IF (len_trim(hfiletype)==0 .AND. .NOT. ounif) THEN
  IF (hvar(1:2)/='TG' .AND. hvar(1:2)/='WG' .OR. hvar(1:3)=='WGI') THEN
    IF (hvar(1:2)/='ZS') WRITE(kluout,*) 'NO FILE FOR FIELD ',hvar, &
                                        ': UNIFORM DEFAULT FIELD IS PRESCRIBED'
    IF (hvar(1:3)=='WGI') THEN
      xhugi_surf_gd = 0.
      xhugi_root_gd = 0.
      xhugi_deep_gd = 0.
    ENDIF                                         
    ounif = .true.
    IF (lhook) CALL dr_hook('READ_PREP_TEB_GARDEN_CONF',1,zhook_handle)
    RETURN
  ELSE
    WRITE(kluout,*) 'AN INPUT FILE OR A UNIFORM VALUE IS REQUIRED FOR FIELD: ',hvar
    WRITE(kluout,*) 'Please complete NAM_PREP_TEB_GARDEN'
    CALL abor1_sfx('READ_PREP_TEB_GARDEN_CONF: AN INPUT FILE OR A UNIFORM VALUE IS REQUIRED FOR '//hvar)
  END IF
END IF
IF (lhook) CALL dr_hook('READ_PREP_TEB_GARDEN_CONF',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_prep_teb_garden_conf