prep_grid_conf_proj.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_grid_conf_proj (GCP,HFILETYPE,HINTERP_TYPE,KNI)
!     ##########################################################################
!
!!****  *PREP_GRID_CONF_PROJ* - reads EXTERNALIZED Surface grid.
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!
!!      V. Masson
!!
!!    MODIFICATIONS
!!    -------------
!!      Original   06/2003
!-------------------------------------------------------------------------------
!
!*      0. DECLARATIONS
!          ------------
!
!
!
!
USE modi_read_surf
!
USE modd_grid_conf_proj_n, ONLY : grid_conf_proj_t, xx, xy
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!* 0.1. Declaration of arguments
!       ------------------------
!
TYPE(grid_conf_proj_t),INTENT(INOUT) :: GCP
!
 CHARACTER(LEN=6),  INTENT(IN)    :: HFILETYPE    ! file type
 CHARACTER(LEN=6),  INTENT(OUT)   :: HINTERP_TYPE ! Grid type
INTEGER,           INTENT(OUT)   :: KNI          ! number of points
!
!* 0.2 Declaration of local variables
!      ------------------------------
!
REAL, DIMENSION(:), ALLOCATABLE :: ZW ! work array
!
 CHARACTER(LEN=12) :: YRECFM    ! Name of the article to be read
 CHARACTER(LEN=1) :: YDIR
INTEGER           :: IRESP
INTEGER           :: JL        ! loop counter
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-----------------------------------------------------------------------
!
!*   1 Projection
!      ----------
!
IF (lhook) CALL dr_hook('PREP_GRID_CONF_PROJ',0,zhook_handle)
yrecfm = 'LAT0'
 CALL read_surf(hfiletype,yrecfm,gcp%XLAT0,iresp)
yrecfm = 'LON0'
 CALL read_surf(hfiletype,yrecfm,gcp%XLON0,iresp)
yrecfm = 'RPK'
 CALL read_surf(hfiletype,yrecfm,gcp%XRPK,iresp)
yrecfm = 'BETA'
 CALL read_surf(hfiletype,yrecfm,gcp%XBETA,iresp)
!
!-----------------------------------------------------------------------
!
!*   2 Grid
!      ----
!
yrecfm = 'LATORI'
 CALL read_surf(hfiletype,yrecfm,gcp%XLATORI,iresp)
yrecfm = 'LONORI'
 CALL read_surf(hfiletype,yrecfm,gcp%XLONORI,iresp)
!
yrecfm = 'IMAX'
 CALL read_surf(hfiletype,yrecfm,gcp%NX,iresp)
yrecfm = 'JMAX'
 CALL read_surf(hfiletype,yrecfm,gcp%NY,iresp)
!
kni = gcp%NX * gcp%NY
!
ydir = '-'
IF (hfiletype=='MESONH') ydir = 'A'
!
ALLOCATE(zw(kni))
!
IF (ALLOCATED(xx)) DEALLOCATE(xx)
ALLOCATE(xx(gcp%NX))
yrecfm = 'XX'
 CALL read_surf(hfiletype,yrecfm,zw,iresp,hdir=ydir)
xx = zw(1:gcp%NX)


IF (ALLOCATED(xy)) DEALLOCATE(xy)
ALLOCATE(xy(gcp%NY))
yrecfm = 'YY'
 CALL read_surf(hfiletype,yrecfm,zw,iresp,hdir=ydir)
DO jl=1,kni
  IF (mod(jl,gcp%NX)==0) xy(jl/gcp%NX) = zw(jl)
END DO
DEALLOCATE(zw)
!
!-----------------------------------------------------------------------
IF(kni==1)THEN
  hinterp_type = 'UNIF  '
ELSE
  hinterp_type = 'BILIN '
ENDIF
IF (lhook) CALL dr_hook('PREP_GRID_CONF_PROJ',1,zhook_handle)
!-----------------------------------------------------------------------
!
END SUBROUTINE prep_grid_conf_proj