write_gridtype_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 write_gridtype_conf_proj (DGU, U, &
                                           hprogram,klu,kgrid_par,pgrid_par,kresp,hdir)
!     #################################################################
!
!!****  *WRITE_GRIDTYPE_CONF_PROJ* - 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_conf_proj
!
!
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
REAL                              :: zrpk     ! projection parameter 
!                                             !   K=1 : stereographic north pole
!                                             ! 0<K<1 : Lambert, north hemisphere
!                                             !   K=0 : Mercator
!                                             !-1<K<0 : Lambert, south hemisphere
!                                             !   K=-1: stereographic south pole
REAL                              :: zbeta    ! angle between grid and reference longitude
REAL                              :: zlatori  ! latitude  of point of coordinates X=0, Y=0
REAL                              :: zlonori  ! longitude of point of coordinates X=0, Y=0
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_CONF_PROJ',0,zhook_handle)
ALLOCATE(zx(klu))
ALLOCATE(zy(klu))
ALLOCATE(zdx(klu))
ALLOCATE(zdy(klu))
!
 CALL get_gridtype_conf_proj(pgrid_par,zlat0,zlon0,zrpk,zbeta,&
                              zlatori,zlonori,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,'RPK ',zrpk, kresp,ycomment)
 CALL write_surf(dgu, u, &
                 hprogram,'BETA',zbeta,kresp,ycomment)
 CALL write_surf(dgu, u, &
                 hprogram,'LATORI',zlatori,kresp,ycomment)
 CALL write_surf(dgu, u, &
                 hprogram,'LONORI',zlonori,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_CONF_PROJ',1,zhook_handle)
!---------------------------------------------------------------------------
!
END SUBROUTINE write_gridtype_conf_proj