mode_gridtype_lonlatval.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.
!     ##############################
      MODULE mode_gridtype_lonlatval
!     ##############################
!
!############################################################################
!############################################################################
!############################################################################
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
 CONTAINS
!############################################################################
!############################################################################
!############################################################################
!     ####################################################################
      SUBROUTINE put_gridtype_lonlatval(PGRID_PAR,PX,PY,PDX,PDY)
!     ####################################################################
!
!!****  *PUT_GRIDTYPE_LONLATVAL* - routine to store in PGRID_PAR the horizontal grid
!!
!!    AUTHOR
!!    ------
!!      S.Faroux   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2010 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
IMPLICIT NONE
!
!
!*       0.1   Declarations of arguments
!              -------------------------
!
REAL, DIMENSION(:), INTENT(IN)  :: px       ! X coordinate of grid mesh center
REAL, DIMENSION(:), INTENT(IN)  :: py       ! Y coordinate of grid mesh center
REAL, DIMENSION(:), INTENT(IN)  :: pdx      ! X grid mesh size
REAL, DIMENSION(:), INTENT(IN)  :: pdy      ! Y grid mesh size
REAL, DIMENSION(:), POINTER     :: pgrid_par! parameters defining this grid
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER :: il ! number of points
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_LONLATVAL:PUT_GRIDTYPE_LONLATVAL',0,zhook_handle)
il = SIZE(px)
ALLOCATE(pgrid_par(1+4*il))

pgrid_par(1) = float(il)
pgrid_par(2:il+1)        = px(:)
pgrid_par(il+2:2*il+1)   = py(:)
pgrid_par(2*il+2:3*il+1) = pdx(:)
pgrid_par(3*il+2:4*il+1) = pdy(:)
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_LONLATVAL:PUT_GRIDTYPE_LONLATVAL',1,zhook_handle)
!
!-------------------------------------------------------------------------------
END SUBROUTINE put_gridtype_lonlatval
!############################################################################
!############################################################################
!############################################################################
!     ####################################################################
      SUBROUTINE get_gridtype_lonlatval(PGRID_PAR,KL,PX,PY,PDX,PDY)
!     ####################################################################
!
!!****  *GET_GRIDTYPE_LONLATVAL* - routine to get from PGRID_PAR the horizontal grid
!!
!!    AUTHOR
!!    ------
!!      S. Faroux   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2010 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
IMPLICIT NONE
!
!
!*       0.1   Declarations of arguments
!              -------------------------
!
INTEGER,            INTENT(OUT), OPTIONAL  :: kl ! number of points
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL  :: px       ! X coordinate of grid mesh center
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL  :: py       ! Y coordinate of grid mesh center
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL  :: pdx      ! X grid mesh size
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL  :: pdy      ! Y grid mesh size
REAL, DIMENSION(:), INTENT(IN)            :: pgrid_par! parameters defining this grid
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER :: il
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_LONLATVAL:GET_GRIDTYPE_LONLATVAL',0,zhook_handle)
IF (present(kl))        kl       = nint(pgrid_par(1))
!
IF (present(px)) THEN
  il = nint(pgrid_par(1))
  px(:) = pgrid_par(2:1+il)
END IF

IF (present(py)) THEN
  il = nint(pgrid_par(1))
  py(:) = pgrid_par(1+il+1:1+2*il)
END IF

IF (present(pdx)) THEN
  il = nint(pgrid_par(1))
  pdx(:)= pgrid_par(1+2*il+1:1+3*il)
END IF

IF (present(pdy)) THEN
  il = nint(pgrid_par(1))
  pdy(:)= pgrid_par(1+3*il+1:1+4*il)
END IF
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_LONLATVAL:GET_GRIDTYPE_LONLATVAL',1,zhook_handle)
!
!-------------------------------------------------------------------------------
END SUBROUTINE get_gridtype_lonlatval
!############################################################################
!############################################################################
!############################################################################
!      ###################################################
       SUBROUTINE latlon_lonlatval(PX,PY,PLAT,PLON)
!      ###################################################
!
!!****  *LATLON_LONLATVAL * - Routine to compute geographical coordinates
!!
!!     PURPOSE
!!     -------
!        This routine computes the latitude and longitude of
!      an array given in LAMBERT coordinates
!
!
!!**   METHOD
!!     ------
!!
!!     EXTERNAL
!!     --------
!!       None
!!
!!     REFERENCE
!!     ---------
!!       
!!     AUTHOR
!!     ------
!!      S.Faroux   *Meteo-France*
!!
!!     MODIFICATION
!!     ------------
!!       Original  02/2010        
!-------------------------------------------------------------------------------
!
!*     0.     DECLARATIONS
!             ------------
!
IMPLICIT NONE
!
!*     0.1    Declarations of arguments and results
!       
REAL, DIMENSION(:),   INTENT(IN) :: px,py
                                           ! given conformal coordinates of the 
                                           ! processed points (meters);
REAL, DIMENSION(:),   INTENT(OUT):: plat,plon    
REAL(KIND=JPRB) :: zhook_handle
                                           ! returned geographic latitudes and 
                                           ! longitudes of the processed points 
                                           ! (degrees).
!
!
      IF (lhook) CALL dr_hook('MODE_GRIDTYPE_LONLATVAL:LATLON_LONLATVAL',0,zhook_handle)
      plon(:)=px(:)
!   
      plat(:)=py(:)
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_LONLATVAL:LATLON_LONLATVAL',1,zhook_handle)
!
!---------------------------------------------------------------------------------
END SUBROUTINE latlon_lonlatval
!---------------------------------------------------------------------------------
!
END MODULE mode_gridtype_lonlatval