latlonmask.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 latlonmask(HGRID,KGRID_PAR,PGRID_PAR,OLATLONMASK)
!     #####################
!
!!**** *LATLONMASK* builds the latiude and longitude mask including the grid
!!
!!    PURPOSE
!!    -------
!!
!!    AUTHOR
!!    ------
!!
!!      V. Masson       Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!      
!!      Original        03/2004
!!                      10/2007  E. Martin  IGN Grids
!!                      12/2012  P. Samuelsson SMHI Rotated lonlat
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_latlonmask_conf_proj
!
USE modi_latlonmask_ign
!
USE modi_latlonmask_lonlat_reg
!
USE modi_latlonmask_lonlatval
!
USE modi_latlonmask_lonlat_rot
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
 CHARACTER(LEN=10),           INTENT(IN)  :: HGRID       ! type of grid
INTEGER                                  :: KGRID_PAR   ! size of PGRID_PAR
REAL,    DIMENSION(:),       POINTER     :: PGRID_PAR   ! parameters defining this grid
LOGICAL, DIMENSION(720,360), INTENT(OUT) :: OLATLONMASK ! mask where domain is
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!----------------------------------------------------------------------------
IF (lhook) CALL dr_hook('LATLONMASK',0,zhook_handle)
SELECT CASE (hgrid)
  CASE('CONF PROJ ')
    CALL latlonmask_conf_proj(kgrid_par,pgrid_par,olatlonmask)

  CASE('LONLAT REG')
    CALL latlonmask_lonlat_reg(kgrid_par,pgrid_par,olatlonmask)

  CASE('IGN       ')
    CALL latlonmask_ign(kgrid_par,pgrid_par,olatlonmask)

  CASE('LONLATVAL ')
    CALL latlonmask_lonlatval(kgrid_par,pgrid_par,olatlonmask)

  CASE('LONLAT ROT')
    CALL latlonmask_lonlat_rot(kgrid_par,pgrid_par,olatlonmask)

  CASE DEFAULT
    olatlonmask(:,:) = .true.
END SELECT
IF (lhook) CALL dr_hook('LATLONMASK',1,zhook_handle)
!----------------------------------------------------------------------------
END SUBROUTINE latlonmask