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, G, KSIZE_FULL, KLUOUT)
!     #######################################
!!
!!    PURPOSE
!!    -------
!!    Computes variables used for interpolation
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson                   Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original     01/2004
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_sfx_grid_n, ONLY : grid_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(grid_t), INTENT(INOUT) :: UG
TYPE(grid_t), INTENT(INOUT) :: G
INTEGER, INTENT(IN) :: KSIZE_FULL
!
INTEGER,           INTENT(IN)  :: KLUOUT     ! output listing logical unit
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PREP_OUTPUT_GRID',0,zhook_handle)
!
IF (.NOT.ALLOCATED(xlat_out)) ALLOCATE(xlat_out(SIZE(g%XLAT)))
IF (.NOT.ALLOCATED(xlon_out)) ALLOCATE(xlon_out(SIZE(g%XLAT)))
IF (.NOT.ALLOCATED(xx_out)) ALLOCATE(xx_out(SIZE(g%XLAT)))
IF (.NOT.ALLOCATED(xy_out)) ALLOCATE(xy_out(SIZE(g%XLAT)))
!
IF (.NOT.ALLOCATED(linterp)) ALLOCATE(linterp(SIZE(g%XLAT)))
!
IF (SIZE(g%XLAT)==0) THEN
  IF (lhook) CALL dr_hook('PREP_OUTPUT_GRID',1,zhook_handle)
  RETURN
ENDIF
!
xlat_out = g%XLAT
xlon_out = g%XLON
linterp  = .true.
!
 CALL get_grid_coord(ug%CGRID, ug%NGRID_PAR, ug%XGRID_PAR, ksize_full, &
                     kluout,xx_out,xy_out,SIZE(g%XLAT),g%CGRID,g%XGRID_PAR)
!
IF (lhook) CALL dr_hook('PREP_OUTPUT_GRID',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE prep_output_grid