write_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.
!     #################################################################
      SUBROUTINE write_gridtype_cartesian (DGU, U, &
                                           hprogram,klu,kgrid_par,pgrid_par,kresp,hdir)
!     #################################################################
!
!!****  *WRITE_GRIDTYPE_CARTESIAN* - 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_cartesian
!
!
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
 CHARACTER(LEN=1), OPTIONAL, INTENT(IN)  :: hdir       ! type of field :
                                            ! 'H' : field with
                                            !       horizontal spatial dim.
                                            ! 'A' : (complete) field with
                                            !       horizontal spatial dim.
                                            ! '-' : no horizontal dim.
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL                              :: zlat0    ! reference latitude
REAL                              :: zlon0    ! reference longitude
INTEGER                           :: iimax    ! number of points in I direction
INTEGER                           :: ijmax    ! number of points in J direction
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_CARTESIAN',0,zhook_handle)
ALLOCATE(zx(klu))
ALLOCATE(zy(klu))
ALLOCATE(zdx(klu))
ALLOCATE(zdy(klu))
!
 CALL get_gridtype_cartesian(pgrid_par,zlat0,zlon0,           &
                              iimax,ijmax,zx,zy,zdx,zdy        )  
!
!---------------------------------------------------------------------------
!
!*       2.    Writing of the grid definition parameters
!              -----------------------------------------
!
ycomment=' '
 CALL write_surf(dgu, u, &
                 hprogram,'LAT0',zlat0,kresp,ycomment)
 CALL write_surf(dgu, u, &
                 hprogram,'LON0',zlon0,kresp,ycomment)
 CALL write_surf(dgu, u, &
                 hprogram,'IMAX ',iimax, kresp,ycomment)
 CALL write_surf(dgu, u, &
                 hprogram,'JMAX ',ijmax, kresp,ycomment)
IF (present(hdir)) THEN
  CALL write_surf(dgu, u, &
                 hprogram,'XX',zx,kresp,ycomment,hdir)
  CALL write_surf(dgu, u, &
                 hprogram,'YY',zy,kresp,ycomment,hdir)
  CALL write_surf(dgu, u, &
                 hprogram,'DX',zdx,kresp,ycomment,hdir)
  CALL write_surf(dgu, u, &
                 hprogram,'DY',zdy,kresp,ycomment,hdir)
ELSE
  CALL write_surf(dgu, u, &
                 hprogram,'XX',zx,kresp,ycomment)
  CALL write_surf(dgu, u, &
                 hprogram,'YY',zy,kresp,ycomment)
  CALL write_surf(dgu, u, &
                 hprogram,'DX',zdx,kresp,ycomment)
  CALL write_surf(dgu, u, &
                 hprogram,'DY',zdy,kresp,ycomment)
END IF
!
!---------------------------------------------------------------------------
DEALLOCATE(zx)
DEALLOCATE(zy)
DEALLOCATE(zdx)
DEALLOCATE(zdy)
IF (lhook) CALL dr_hook('WRITE_GRIDTYPE_CARTESIAN',1,zhook_handle)
!---------------------------------------------------------------------------
!
END SUBROUTINE write_gridtype_cartesian