write_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.
!     #################################################################
      SUBROUTINE write_gridtype_lonlatval (DGU, U, &
                                           hprogram,klu,kgrid_par,pgrid_par,kresp)
!     #################################################################
!
!!****  *WRITE_GRIDTYPE_IGN* - routine to write the horizontal grid
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004 
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modi_write_surf
!
USE mode_gridtype_lonlatval
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
 CHARACTER(LEN=6),           INTENT(IN)  :: hprogram   ! calling program
INTEGER,                    INTENT(IN)  :: klu        ! number of points
INTEGER,                    INTENT(IN)  :: kgrid_par  ! size of PGRID_PAR
REAL, DIMENSION(KGRID_PAR), INTENT(IN)  :: pgrid_par  ! parameters defining this grid
INTEGER,                    INTENT(OUT) :: kresp      ! error return code
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(:), ALLOCATABLE   :: zx       ! X conformal coordinate of grid mesh (dim IIMAX)
REAL, DIMENSION(:), ALLOCATABLE   :: zy       ! Y conformal coordinate of grid mesh (dim IJMAX)
REAL, DIMENSION(:), ALLOCATABLE   :: zdx      ! X grid mesh size (dim IIMAX)
REAL, DIMENSION(:), ALLOCATABLE   :: zdy      ! Y grid mesh size (dim IJMAX)
!
 CHARACTER(LEN=100)                :: ycomment ! comment written in the file
REAL(KIND=JPRB) :: zhook_handle
!---------------------------------------------------------------------------
!
!*       1.    Projection and 2D grid parameters
!              ---------------------------------
!
IF (lhook) CALL dr_hook('WRITE_GRIDTYPE_LONLATVAL',0,zhook_handle)
ALLOCATE(zx(klu))
ALLOCATE(zy(klu))
ALLOCATE(zdx(klu))
ALLOCATE(zdy(klu))
!
 CALL get_gridtype_lonlatval(pgrid_par, &
                        px=zx,py=zy,pdx=zdx,pdy=zdy  )  
!
!---------------------------------------------------------------------------
!
!*       2.    Writing of the grid definition parameters
!              -----------------------------------------
!
ycomment='XX'
 CALL write_surf(dgu, u, &
                 hprogram,'XX',zx,kresp,ycomment)
!
ycomment='XY'
 CALL write_surf(dgu, u, &
                 hprogram,'XY',zy,kresp,ycomment)
!
ycomment='XDX'
 CALL write_surf(dgu, u, &
                 hprogram,'DX',zdx,kresp,ycomment)
!
ycomment='XDY'
 CALL write_surf(dgu, u, &
                 hprogram,'DY',zdy,kresp,ycomment)
!
!---------------------------------------------------------------------------
DEALLOCATE(zx)
DEALLOCATE(zy)
DEALLOCATE(zdx)
DEALLOCATE(zdy)
IF (lhook) CALL dr_hook('WRITE_GRIDTYPE_LONLATVAL',1,zhook_handle)
!---------------------------------------------------------------------------
!
END SUBROUTINE write_gridtype_lonlatval