pgd_dummy.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 pgd_dummy (DTCO, DUU, UG, U, USS, &
                            hprogram)
!     ##############################################################
!
!!**** *PGD_DUMMY* monitor for averaging and interpolations of physiographic fields
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    10/12/97
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_dummy_surf_fields_n, ONLY : dummy_surf_fields_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_surf_atm_sso_n, ONLY : surf_atm_sso_t
!
USE modd_pgd_grid,           ONLY : nl
USE modd_pgdwork,            ONLY : catype
USE modd_surf_par,           ONLY : xundef
!
USE modi_get_luout
USE modi_pgd_field
USE modi_read_nam_pgd_dummy
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: dtco
TYPE(dummy_surf_fields_t), INTENT(INOUT) :: duu
TYPE(surf_atm_grid_t), INTENT(INOUT) :: ug
TYPE(surf_atm_t), INTENT(INOUT) :: u
TYPE(surf_atm_sso_t), INTENT(INOUT) :: uss
!
 CHARACTER(LEN=6),    INTENT(IN)    :: hprogram     ! Type of program
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER                           :: iluout    ! output listing logical unit
INTEGER                           :: jnbr      ! loop counter on dummy fields
!
!*    0.3    Declaration of namelists
!            ------------------------
!
INTEGER                             :: idummy_nbr
 CHARACTER(LEN=20), DIMENSION(1000)  :: ydummy_name
 CHARACTER(LEN=3),  DIMENSION(1000)  :: ydummy_area
 CHARACTER(LEN=3),  DIMENSION(1000)  :: cdummy_atype    ! avg type for dummy pgd fields
!                                                      ! 'ARI' , 'INV'
 CHARACTER(LEN=28), DIMENSION(1000)  :: cdummy_file     ! data files
 CHARACTER(LEN=6),  DIMENSION(1000)  :: cdummy_filetype ! type of these files
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations of defaults
!             ---------------------------
!
IF (lhook) CALL dr_hook('PGD_DUMMY',0,zhook_handle)
 CALL get_luout(hprogram,iluout)
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of namelist
!             -------------------
!
 CALL read_nam_pgd_dummy(hprogram, idummy_nbr, ydummy_name, ydummy_area, &
                          cdummy_atype, cdummy_file, cdummy_filetype      )  
!
duu%NDUMMY_NBR     = idummy_nbr
ALLOCATE(duu%CDUMMY_NAME(duu%NDUMMY_NBR))
ALLOCATE(duu%CDUMMY_AREA(duu%NDUMMY_NBR))
duu%CDUMMY_NAME(:) = ydummy_name(1:duu%NDUMMY_NBR)
duu%CDUMMY_AREA(:) = ydummy_area(1:duu%NDUMMY_NBR)
!
!-------------------------------------------------------------------------------
!
!*    3.      Allocation
!             ----------
!
ALLOCATE(duu%XDUMMY_FIELDS(nl,duu%NDUMMY_NBR))
!
!-------------------------------------------------------------------------------
!
!*    4.      Computations
!             ------------
!
DO jnbr=1,duu%NDUMMY_NBR
  catype = cdummy_atype(jnbr)
  CALL pgd_field(dtco, ug, u, uss, &
                 hprogram,duu%CDUMMY_NAME(jnbr),duu%CDUMMY_AREA(jnbr),cdummy_file(jnbr), &
                   cdummy_filetype(jnbr),xundef,duu%XDUMMY_FIELDS(:,jnbr)              )  
  catype = 'ARI'
END DO
IF (lhook) CALL dr_hook('PGD_DUMMY',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE pgd_dummy