mode_gridtype_cartesian.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_cartesian
!     ##############################
!
!############################################################################
!############################################################################
!############################################################################
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
 CONTAINS
!############################################################################
!############################################################################
!############################################################################
!     ####################################################################
      SUBROUTINE put_gridtype_cartesian(PGRID_PAR,PLAT0,PLON0,         &
                                          kimax,kjmax,px,py,pdx,pdy      )  
!     ####################################################################
!
!!****  *PUT_GRIDTYPE_CARTESIAN* - routine to store in PGRID_PAR the horizontal grid
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
IMPLICIT NONE
!
!
!*       0.1   Declarations of arguments
!              -------------------------
!
REAL,               INTENT(IN)  :: plat0    ! reference latitude
REAL,               INTENT(IN)  :: plon0    ! reference longitude
INTEGER,            INTENT(IN)  :: kimax    ! number of points in I direction
INTEGER,            INTENT(IN)  :: kjmax    ! number of points in J direction
REAL, DIMENSION(:), INTENT(IN)  :: px       ! X conformal coordinate of left boundary of grid mesh
REAL, DIMENSION(:), INTENT(IN)  :: py       ! Y conformal coordinate of bottom boundary of grid mesh
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_CARTESIAN:PUT_GRIDTYPE_CARTESIAN',0,zhook_handle)
il = SIZE(px)
ALLOCATE(pgrid_par(4+4*il))
pgrid_par(1) = plat0
pgrid_par(2) = plon0
pgrid_par(3) = float(kimax)
pgrid_par(4) = float(kjmax)
pgrid_par(4     +1:4+  il) = px(:)
pgrid_par(4+  il+1:4+2*il) = py(:)
pgrid_par(4+2*il+1:4+3*il) = pdx(:)
pgrid_par(4+3*il+1:4+4*il) = pdy(:)
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_CARTESIAN:PUT_GRIDTYPE_CARTESIAN',1,zhook_handle)
!-------------------------------------------------------------------------------
END SUBROUTINE put_gridtype_cartesian
!############################################################################
!############################################################################
!############################################################################
!     ####################################################################
      SUBROUTINE get_gridtype_cartesian(PGRID_PAR,PLAT0,PLON0,           &
                                          kimax,kjmax,px,py,pdx,pdy,kl     )  
!     ####################################################################
!
!!****  *GET_GRIDTYPE_CARTESIAN* - routine to get from PGRID_PAR the horizontal grid
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
IMPLICIT NONE
!
!
!*       0.1   Declarations of arguments
!              -------------------------
!
REAL, DIMENSION(:), INTENT(IN)            :: pgrid_par! parameters defining this grid
REAL,               INTENT(OUT), OPTIONAL :: plat0    ! reference latitude
REAL,               INTENT(OUT), OPTIONAL :: plon0    ! reference longitude
INTEGER,            INTENT(OUT), OPTIONAL :: kimax    ! number of points in I direction
INTEGER,            INTENT(OUT), OPTIONAL :: kjmax    ! number of points in J direction
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL :: px       ! X conformal coor. of grid mesh 
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL :: py       ! Y conformal coor. of grid mesh
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL :: pdx      ! X grid mesh size
REAL, DIMENSION(:), INTENT(OUT), OPTIONAL :: pdy      ! Y grid mesh size
INTEGER,            INTENT(OUT), OPTIONAL :: kl       ! number of points
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER :: il
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_CARTESIAN:GET_GRIDTYPE_CARTESIAN',0,zhook_handle)
IF (present(plat0))  plat0 = pgrid_par(1)
IF (present(plon0))  plon0 = pgrid_par(2)
IF (present(kimax))  kimax = nint(pgrid_par(3))
IF (present(kjmax))  kjmax = nint(pgrid_par(4))
!
IF (present(px)) THEN
  il = SIZE(px)
  px(:) = pgrid_par(4+1:4+il)
END IF

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

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

IF (present(pdy)) THEN
  il = SIZE(pdy)
  pdy(:)= pgrid_par(4+3*il+1:4+4*il)
END IF
!
IF (present(kl)) kl = (SIZE(pgrid_par)-4)/4
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_CARTESIAN:GET_GRIDTYPE_CARTESIAN',1,zhook_handle)
!
!-------------------------------------------------------------------------------
END SUBROUTINE get_gridtype_cartesian
!############################################################################
!############################################################################
!############################################################################
!      ###################################################
       SUBROUTINE latlon_cartesian(PLAT0,PLON0,PLAT,PLON)
!      ###################################################
!
!!****  *LATLON_CARTESIAN * - Routine to compute geographical coordinates
!!
!!     PURPOSE
!!     -------
!
!
!!**   METHOD
!!     ------
!!       
!!     AUTHOR
!!     ------
!!      V. Masson  *Meteo France*
!!
!!     MODIFICATION
!!     ------------
!!       Original  06/2004
!-------------------------------------------------------------------------------
!
!*     0.     DECLARATIONS
!             ------------
!
!
IMPLICIT NONE
!
!*     0.1    Declarations of arguments and results
!
REAL,                 INTENT(IN) :: plat0  ! Reference latitude
REAL,                 INTENT(IN) :: plon0  ! Reference longitude
REAL, DIMENSION(:),   INTENT(OUT):: plat,plon    
REAL(KIND=JPRB) :: zhook_handle
                                           ! returned geographic latitudes and 
                                           ! longitudes of the processed points 
                                           ! (degrees).
!
!*     0.2    Declarations of local variables
! 
!--------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_CARTESIAN:LATLON_CARTESIAN',0,zhook_handle)
plon(:) = plon0 
plat(:) = plat0 
IF (lhook) CALL dr_hook('MODE_GRIDTYPE_CARTESIAN:LATLON_CARTESIAN',1,zhook_handle)
!
!---------------------------------------------------------------------------------
END SUBROUTINE latlon_cartesian
!---------------------------------------------------------------------------------
!
END MODULE mode_gridtype_cartesian