pgd_seaflux.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_seaflux (DTCO, DTS, SG, S, UG, U, USS, &
                              HPROGRAM)
!     ##############################################################
!
!!**** *PGD_SEAFLUX* monitor for averaging and interpolations of SEAFLUX physiographic fields
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    03/2004
!!    Lebeaupin-B C. 01/2008 : include bathymetry
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_data_seaflux_n, ONLY : data_seaflux_t
USE modd_sfx_grid_n, ONLY : grid_t
USE modd_seaflux_n, ONLY : seaflux_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_sso_n, ONLY : sso_t
!
USE modd_pgd_grid,       ONLY : nl
USE modd_data_cover_par,  ONLY : jpcover
!
USE modi_read_nam_pgd_seabathy
USE modi_pgd_bathyfield
!
USE modi_get_surf_size_n
USE modi_pack_pgd
USE modi_pack_pgd_seaflux
USE modi_pgd_seaflux_par
!
!
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(data_seaflux_t), INTENT(INOUT) :: DTS
TYPE(grid_t), INTENT(INOUT) :: SG
TYPE(seaflux_t), INTENT(INOUT) :: S
TYPE(surf_atm_grid_t), INTENT(INOUT) :: UG
TYPE(surf_atm_t), INTENT(INOUT) :: U
TYPE(sso_t), INTENT(INOUT) :: USS
!
 CHARACTER(LEN=6),    INTENT(IN)    :: HPROGRAM     ! Type of program
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL, DIMENSION(NL)               :: ZSEABATHY ! bathymetry on all surface points
!
!*    0.3    Declaration of namelists
!            ------------------------
!
 CHARACTER(LEN=28)        :: YSEABATHY         ! file name for bathymetrie
 CHARACTER(LEN=6)         :: YSEABATHYFILETYPE ! bathymetry data file type
 CHARACTER(LEN=28)        :: YNCVARNAME        ! variable to read in netcdf
                                              ! file
REAL                     :: XUNIF_SEABATHY    ! uniform value of bathymetry
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations of defaults
!             ---------------------------
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of namelist
!             -------------------
!
IF (lhook) CALL dr_hook('PGD_SEAFLUX',0,zhook_handle)
 CALL read_nam_pgd_seabathy(hprogram,yseabathy,yseabathyfiletype,yncvarname,&
       xunif_seabathy)  
!
!-------------------------------------------------------------------------------
!
!*    3.      Coherence of options
!             --------------------
!
!-------------------------------------------------------------------------------
!
!*    4.      Bathymetry
!             ----------
!
 CALL pgd_bathyfield(ug, u, uss, hprogram,'bathymetry','SEA',yseabathy,yseabathyfiletype,&
       yncvarname,xunif_seabathy,zseabathy(:))  
!-------------------------------------------------------------------------------
!
!*    5.      Number of points and packing
!             ----------------------------
!
 CALL get_surf_size_n(dtco, u, 'SEA   ',sg%NDIM)
!
ALLOCATE(s%LCOVER     (jpcover))
ALLOCATE(s%XZS        (sg%NDIM))
ALLOCATE(sg%XLAT       (sg%NDIM))
ALLOCATE(sg%XLON       (sg%NDIM))
ALLOCATE(sg%XMESH_SIZE (sg%NDIM))
!
 CALL pack_pgd(dtco, u, hprogram, 'SEA   ', sg, s%LCOVER, s%XCOVER, s%XZS )  
!
 CALL pack_pgd_seaflux(dtco, sg%NDIM, s, u, hprogram, zseabathy)
!
 CALL pgd_seaflux_par(dtco, dts, sg%NDIM, ug, u, uss, hprogram)
 !
IF (lhook) CALL dr_hook('PGD_SEAFLUX',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE pgd_seaflux