pack_pgd_isba.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 pack_pgd_isba (DTCO, IG, I, U, &
                                hprogram,                                    &
                                 paosip, paosim, paosjp, paosjm,              &
                                 pho2ip, pho2im, pho2jp, pho2jm,              &
                                 psso_slope                                   )  
!     ##############################################################
!
!!**** *PACK_PGD_ISBA* packs ISBA physiographic fields from all surface points to ISBA points
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    03/2004
!!    Escobar J.  08/02/2005 : bug declare ILU local variable
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_isba_grid_n, ONLY : isba_grid_t
USE modd_isba_n, ONLY : isba_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modi_pack_same_rank
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_get_surf_mask_n
!
USE modi_get_type_dim_n
!
USE modi_get_luout
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: dtco
TYPE(isba_grid_t), INTENT(INOUT) :: ig
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
 CHARACTER(LEN=6),        INTENT(IN) :: hprogram  ! Type of program
REAL,    DIMENSION(:),   INTENT(IN) :: paosip    ! A/S i+ on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: paosim    ! A/S i- on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: paosjp    ! A/S j+ on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: paosjm    ! A/S j- on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: pho2ip    ! h/2 i+ on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: pho2im    ! h/2 i- on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: pho2jp    ! h/2 j+ on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: pho2jm    ! h/2 j- on all surface points
REAL,    DIMENSION(:),   INTENT(IN) :: psso_slope! subgrid slope on all surface points
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER                        :: ilu    ! expected physical size of full surface array
INTEGER                        :: iluout ! output listing logical unit
INTEGER, DIMENSION(:), POINTER :: imask  ! mask for packing from complete field to nature field
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PACK_PGD_ISBA',0,zhook_handle)
 CALL get_luout(hprogram,iluout)
!
!*    1.      Number of points and packing
!             ----------------------------
!
 CALL get_type_dim_n(dtco, u, &
                     'NATURE',ig%NDIM)
ALLOCATE(imask(ig%NDIM))
ilu=0
 CALL get_surf_mask_n(dtco, u, &
                      'NATURE',ig%NDIM,imask,ilu,iluout)
!
!
!-------------------------------------------------------------------------------
!
!*    2.      Packing of fields
!             -----------------
!
ALLOCATE(i%XAOSIP(ig%NDIM))
ALLOCATE(i%XAOSIM(ig%NDIM))
ALLOCATE(i%XAOSJP(ig%NDIM))
ALLOCATE(i%XAOSJM(ig%NDIM))
ALLOCATE(i%XHO2IP(ig%NDIM))
ALLOCATE(i%XHO2IM(ig%NDIM))
ALLOCATE(i%XHO2JP(ig%NDIM))
ALLOCATE(i%XHO2JM(ig%NDIM))
ALLOCATE(i%XSSO_SLOPE(ig%NDIM))
 CALL pack_same_rank(imask,paosip(:),i%XAOSIP(:))
 CALL pack_same_rank(imask,paosim(:),i%XAOSIM(:))
 CALL pack_same_rank(imask,paosjp(:),i%XAOSJP(:))
 CALL pack_same_rank(imask,paosjm(:),i%XAOSJM(:))
 CALL pack_same_rank(imask,pho2ip(:),i%XHO2IP(:))
 CALL pack_same_rank(imask,pho2im(:),i%XHO2IM(:))
 CALL pack_same_rank(imask,pho2jp(:),i%XHO2JP(:))
 CALL pack_same_rank(imask,pho2jm(:),i%XHO2JM(:))
 CALL pack_same_rank(imask,psso_slope(:),i%XSSO_SLOPE(:))
IF (lhook) CALL dr_hook('PACK_PGD_ISBA',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE pack_pgd_isba