prep_output_grid.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 prep_output_grid (UG, U, &
                                   kluout,hgrid,pgrid_par,plat,plon)
!     #######################################
!!
!!    PURPOSE
!!    -------
!!    Computes variables used for interpolation
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson                   Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original     01/2004
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modi_get_grid_coord
!
USE modd_prep, ONLY : xlat_out, xlon_out, xx_out, xy_out, linterp
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of dummy arguments
!            ------------------------------
!
!
TYPE(surf_atm_grid_t), INTENT(INOUT) :: ug
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
INTEGER,           INTENT(IN)  :: kluout     ! output listing logical unit
 CHARACTER(LEN=10), INTENT(IN)  :: hgrid      ! grid type
REAL, DIMENSION(:), POINTER    :: pgrid_par  ! parameters defining this grid
REAL, DIMENSION(:),INTENT(IN)  :: plat       ! latitudes
REAL, DIMENSION(:),INTENT(IN)  :: plon       ! longitudes
REAL(KIND=JPRB) :: zhook_handle
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
!
!------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PREP_OUTPUT_GRID',0,zhook_handle)
IF (.NOT.ALLOCATED(xlat_out)) ALLOCATE(xlat_out(SIZE(plat)))
IF (.NOT.ALLOCATED(xlon_out)) ALLOCATE(xlon_out(SIZE(plat)))
IF (.NOT.ALLOCATED(xx_out)) ALLOCATE(xx_out(SIZE(plat)))
IF (.NOT.ALLOCATED(xy_out)) ALLOCATE(xy_out(SIZE(plat)))
!
IF (.NOT.ALLOCATED(linterp)) ALLOCATE(linterp(SIZE(plat)))

xlat_out = plat
xlon_out = plon
linterp  = .true.
!
 CALL get_grid_coord(ug, u, &
                     kluout,xx_out,xy_out,SIZE(plat),hgrid,pgrid_par)
IF (lhook) CALL dr_hook('PREP_OUTPUT_GRID',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE prep_output_grid