pgd_surf_atm.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 pgd_surf_atm (YSC,HPROGRAM,HFILE,HFILETYPE,OZS)
!     ###########################################################
!!
!!    PURPOSE
!!    -------
!!   This program prepares the physiographic data fields.
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson                   Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original     13/10/03
!!      A. Lemonsu      05/2009         Ajout de la clef LGARDEN pour TEB
!!      J. Escobar      11/2013         Add USE MODI_READ_NAM_PGD_CHEMISTRY
!!      B. Decharme     02/2014         Add LRM_RIVER
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE modd_surfex_n, ONLY : surfex_t
!
USE modd_surfex_mpi, ONLY : nrank, npio, nsize, nindex, nnum
!
USE modd_surf_conf,       ONLY : cprogname
USE modd_pgd_grid,        ONLY : nl, llatlonmask, ngrid_par
!
USE modi_get_size_full_n
USE modi_get_luout
USE modi_read_pgd_arrange_cover
USE modi_read_pgd_cover_garden
USE modi_ini_data_cover
USE modi_read_pgd_schemes
USE modi_read_nam_pgd_chemistry
USE modi_read_nam_write_cover_tex
USE modi_write_cover_tex_start
USE modi_write_cover_tex_cover
USE modi_latlon_grid
USE modi_put_pgd_grid
USE modi_latlonmask
USE modi_pgd_frac
USE modi_pgd_cover
USE modi_pgd_orography
USE modi_pgd_nature
USE modi_pgd_town
USE modi_pgd_inland_water
USE modi_pgd_sea
USE modi_pgd_dummy
USE modi_pgd_chemistry
USE modi_pgd_chemistry_snap
USE modi_write_cover_tex_end
USE modi_init_read_data_cover
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of dummy arguments
!            ------------------------------
!
!
TYPE(surfex_t), INTENT(INOUT) :: YSC
!
 CHARACTER(LEN=6),     INTENT(IN)  :: HPROGRAM ! program calling
 CHARACTER(LEN=28),    INTENT(IN)  :: HFILE    ! atmospheric file name
 CHARACTER(LEN=6),     INTENT(IN)  :: HFILETYPE! atmospheric file type
LOGICAL,              INTENT(IN)  :: OZS      ! .true. if orography is imposed by atm. model
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
LOGICAL :: LRM_RIVER   !delete inland river coverage. Default is false
!
INTEGER :: ISIZE_FULL, JI, IDIM_FULL
INTEGER :: ILUOUT ! logical unit of output listing file
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('PGD_SURF_ATM',0,zhook_handle)
!
lrm_river = .false.
!
CPROGNAME=HPROGRAM
!
 CALL get_luout(hprogram,iluout)
!
!*    1.      Set default constant values 
!             ---------------------------
!
!*    1.2     surface schemes
 CALL read_pgd_schemes(hprogram, ysc%U%CNATURE, ysc%U%CSEA, ysc%U%CTOWN, ysc%U%CWATER)
!
 CALL read_nam_write_cover_tex(hprogram)
!
!-------------------------------------------------------------------------------
!
!*    2.      Grid
!             ----
!
ALLOCATE(ysc%UG%G%XLAT      (ysc%U%NSIZE_FULL))
ALLOCATE(ysc%UG%G%XLON      (ysc%U%NSIZE_FULL))
ALLOCATE(ysc%UG%G%XMESH_SIZE(ysc%U%NSIZE_FULL))
ALLOCATE(ysc%UG%XJPDIR      (ysc%U%NSIZE_FULL))
 CALL latlon_grid(ysc%UG%G, ysc%U%NSIZE_FULL, ysc%UG%XJPDIR)
!
!
!*    2.3     Stores the grid in the module MODD_PGD_GRID
!
 CALL put_pgd_grid(ysc%UG%G%CGRID, ysc%U%NSIZE_FULL,ysc%UG%G%NGRID_PAR, ysc%UG%G%XGRID_PAR)
!
IF (hprogram=='MESONH') THEN
  idim_full = ysc%U%NDIM_FULL
  ysc%U%NDIM_FULL = nl
  nsize = nl
  ALLOCATE(nindex(nl))
  nindex(:) = 0
  ALLOCATE(nnum(nl))
  DO ji = 1,nl
    nnum(ji) = ji
  ENDDO
ENDIF
!
IF (.NOT.ASSOCIATED(ysc%UG%XGRID_FULL_PAR)) THEN
  ALLOCATE(ysc%UG%XGRID_FULL_PAR(SIZE(ysc%UG%G%XGRID_PAR)))
  ysc%UG%XGRID_FULL_PAR(:) = ysc%UG%G%XGRID_PAR(:)
  ysc%UG%NGRID_FULL_PAR = ngrid_par
ENDIF
!
!*    2.4     mask to limit the number of input data to read
 CALL latlonmask(ysc%UG%G%CGRID, ysc%UG%NGRID_FULL_PAR, ysc%UG%XGRID_FULL_PAR, llatlonmask)
!
!-------------------------------------------------------------------------------
!
!*    3.      surface cover
!             -------------
!
 CALL pgd_frac(ysc%DTCO, ysc%UG, ysc%U, ysc%USS, hprogram)
!
 CALL read_pgd_arrange_cover(hprogram, ysc%U%LWATER_TO_NATURE, ysc%U%LTOWN_TO_ROCK)
!
 CALL read_pgd_cover_garden(hprogram, ysc%U%LGARDEN)
!
 CALL init_read_data_cover(hprogram)
!
 CALL ini_data_cover(ysc%DTCO, ysc%U)
IF (ysc%U%LECOCLIMAP) CALL pgd_cover(ysc%DTCO, ysc%UG, ysc%U, ysc%USS, hprogram,lrm_river)
!
IF (nrank==npio) THEN
  CALL write_cover_tex_start(hprogram)
  CALL write_cover_tex_cover
ENDIF
!
!-------------------------------------------------------------------------------
!
!*    4.      Orography
!             ---------
!
 CALL pgd_orography(ysc%DTCO, ysc%UG, ysc%U, ysc%USS, hprogram, hfile, hfiletype, ozs)
!
!_______________________________________________________________________________
!
!*    5.      Additionnal fields for nature scheme
!             ------------------------------------
!
IF (ysc%U%NDIM_NATURE>0) CALL pgd_nature(ysc%DTCO, ysc%DTZ, ysc%IM, ysc%UG, ysc%U, ysc%USS, hprogram)  
!_______________________________________________________________________________
!
!*    6.      Additionnal fields for town scheme
!             ----------------------------------
!
IF (ysc%U%NDIM_TOWN>0) CALL pgd_town(ysc%DTCO, ysc%UG, ysc%U, ysc%USS, &
                                     ysc%IM%DTV, ysc%TM, ysc%GDM, ysc%GRM, hprogram)  
!_______________________________________________________________________________
!
!*    7.      Additionnal fields for inland water scheme
!             ------------------------------------------
!
IF (ysc%U%NDIM_WATER>0) CALL pgd_inland_water(ysc%DTCO, ysc%FM%G, ysc%FM%F, ysc%UG, ysc%U, &
                                              ysc%USS, ysc%WM%G, ysc%WM%W, hprogram,lrm_river)   
!_______________________________________________________________________________
!
!*    8.      Additionnal fields for sea scheme
!             ---------------------------------
!
IF (ysc%U%NDIM_SEA>0) CALL pgd_sea(ysc%DTCO, ysc%SM%DTS, ysc%SM%G, ysc%SM%S, &
                                   ysc%UG, ysc%U, ysc%USS, hprogram)  
!_______________________________________________________________________________
!
!*    9.      Dummy fields
!             ------------
!
 CALL pgd_dummy(ysc%DTCO, ysc%DUU, ysc%UG, ysc%U, ysc%USS, hprogram)
!_______________________________________________________________________________
!
!*   10.      Chemical Emission fields
!             ------------------------
!
 CALL read_nam_pgd_chemistry(hprogram,ysc%CHU%CCH_EMIS)
IF (ysc%CHU%CCH_EMIS=='SNAP') THEN
  CALL pgd_chemistry_snap(ysc%CHN, ysc%DTCO, ysc%UG, ysc%U, ysc%USS, &
                          hprogram,ysc%CHU%LCH_EMIS)
ELSE IF (ysc%CHU%CCH_EMIS=='AGGR') THEN
  CALL pgd_chemistry(ysc%CHE, ysc%DTCO, ysc%UG, ysc%U, ysc%USS, &
                     hprogram,ysc%CHU%LCH_EMIS)
ENDIF
!_______________________________________________________________________________
!
!*   11.     Writing in cover latex file
!            ---------------------------
!
IF (nrank==npio) CALL write_cover_tex_end(hprogram)
!
IF (hprogram=='MESONH') THEN
 ysc%U%NDIM_FULL = idim_full
ENDIF
!
IF (lhook) CALL dr_hook('PGD_SURF_ATM',1,zhook_handle)
!_______________________________________________________________________________
!
END SUBROUTINE pgd_surf_atm