writesurf_pgd_isban.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 writesurf_pgd_isba_n (HSELECT, HNATURE, DTV, DTZ, G, ISS, IO, S, K, HPROGRAM)
!     ################################################
!
!!****  *WRITESURF_PGD_ISBA_n* - writes ISBA physiographic fields
!!                        
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2003 
!!      P. Le Moigne 12/2004 : add type of photosynthesis 
!!      B. Decharme  06/2009 : add topographic index statistics
!!      A.L. Gibelin 04/2009 : dimension NBIOMASS for ISBA-A-gs
!!      B. Decharme  07/2011 : delete argument HWRITE
!!      B. Decharme  07/2012 : files of data for permafrost area and for SOC top and sub soil
!!                   11/2013 : same for groundwater distribution
!!                   11/2014 : Write XSOILGRID as a series of real 
!!      P. Samuelsson 10/2014 : MEB
!!      M. Moge      02/2015 parallelization using WRITE_LCOVER
!!    10/2016 B. Decharme : bug surface/groundwater coupling   
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
#ifdef SFX_OL
USE modn_io_offline, ONLY : lwr_vegtype
#endif
!
USE modd_data_isba_n, ONLY : data_isba_t
USE modd_data_tsz0_n, ONLY : data_tsz0_t
USE modd_sfx_grid_n, ONLY : grid_t
USE modd_sso_n, ONLY : sso_t
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY : isba_s_t, isba_k_t
!
USE modd_surf_par, ONLY : xundef
USE modd_data_cover_par, ONLY : jpcover
!
USE mode_write_surf_cov, ONLY : write_surf_cov
!
USE modi_write_surf
USE modi_write_grid
USE modi_writesurf_pgd_isba_par_n
USE modi_writesurf_pgd_tsz0_par_n
USE modi_write_lcover
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
 CHARACTER(LEN=*), DIMENSION(:), INTENT(IN) :: HSELECT
 CHARACTER(LEN=*), INTENT(IN) :: HNATURE
!
TYPE(data_isba_t), INTENT(INOUT) :: DTV
TYPE(data_tsz0_t), INTENT(INOUT) :: DTZ
TYPE(grid_t), INTENT(INOUT) :: G
TYPE(sso_t), INTENT(INOUT) :: ISS
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(isba_s_t), INTENT(INOUT) :: S
TYPE(isba_k_t), INTENT(INOUT) :: K
!
 CHARACTER(LEN=6),  INTENT(IN)  :: HPROGRAM ! program calling
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: IRESP          ! IRESP  : return-code if a problem appears
 CHARACTER(LEN=12) :: YRECFM         ! Name of the article to be read
 CHARACTER(LEN=100):: YCOMMENT       ! Comment string
 CHARACTER(LEN=4 ) :: YLVL
!
INTEGER :: JJ, JLAYER
INTEGER :: ISIZE_LMEB_PATCH  ! Number of patches with MEB=true
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!
!-------------------------------------------------------------------------------
!
!
!* soil scheme option
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_ISBA_N',0,zhook_handle)
yrecfm='ISBA'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%CISBA,iresp,hcomment=ycomment)
!
!* Pedo-transfert function
!
yrecfm='PEDOTF'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%CPEDOTF,iresp,hcomment=ycomment)
!
!* type of photosynthesis
!
yrecfm='PHOTO'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%CPHOTO,iresp,hcomment=ycomment)
!
!* new radiative transfert
!
yrecfm='TR_ML'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%LTR_ML,iresp,hcomment=ycomment)
!
!* calbedo
!
yrecfm='ALBEDO'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%CALBEDO,iresp,hcomment=ycomment)
!
!* threshold to remove little fractions of patches
!
yrecfm='RM_PATCH'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%XRM_PATCH,iresp,hcomment=ycomment)

!* number of soil layers
!
yrecfm='GROUND_LAYER'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%NGROUND_LAYER,iresp,hcomment=ycomment)
!
!* Reference grid for DIF
!
IF(io%CISBA=='DIF') THEN
  DO jlayer=1,io%NGROUND_LAYER
    WRITE(ylvl,'(I4)') jlayer     
    yrecfm='SOILGRID'//adjustl(ylvl(:len_trim(ylvl)))
    ycomment='Depth of ISBA soilgrid layer '//adjustl(ylvl(:len_trim(ylvl)))
    CALL write_surf(hselect, hprogram,yrecfm,io%XSOILGRID(jlayer),iresp,hcomment=ycomment)
  END DO
ENDIF
!
!* number of biomass pools
!
yrecfm='NBIOMASS'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%NNBIOMASS,iresp,hcomment=ycomment)
!
!* number of tiles
!
yrecfm='PATCH_NUMBER'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%NPATCH,iresp,hcomment=ycomment)
!
!* flag indicating if fields are computed from ecoclimap or not
!
yrecfm='ECOCLIMAP'
ycomment=yrecfm
 CALL write_surf(hselect, hprogram,yrecfm,io%LECOCLIMAP,iresp,hcomment=ycomment)
!
!* logical vector indicating for which patches MEB should be applied
!
yrecfm='MEB_PATCH'
ycomment='(LOGICAL LIST)'
 CALL write_surf(hselect, hprogram,yrecfm,io%LMEB_PATCH(:),iresp,hcomment=ycomment,hdir='-')
!
isize_lmeb_patch = count(io%LMEB_PATCH(:))
!
IF (isize_lmeb_patch>0)THEN
!
!* flag indicating if forcing is from observed measurements or not
!
   yrecfm='FORC_MEASURE'
   ycomment=yrecfm
   CALL write_surf(hselect, hprogram,yrecfm,io%LFORC_MEASURE,iresp,hcomment=ycomment)
!
!* flag indicating if litter layer is used or not
!
   yrecfm='MEB_LITTER'
   ycomment=yrecfm
   CALL write_surf(hselect, hprogram,yrecfm,io%LMEB_LITTER,iresp,hcomment=ycomment)
!
!* flag indicating if ground resistance is used or not
!
   yrecfm='MEB_GNDRES'
   ycomment=yrecfm
   CALL write_surf(hselect, hprogram,yrecfm,io%LMEB_GNDRES,iresp,hcomment=ycomment)
!
ENDIF
!
!*       2.     Physiographic data fields:
!               -------------------------
!
!* cover classes
!
 CALL write_lcover(hselect,hprogram,s%LCOVER)
!
#ifdef SFX_OL
IF (lwr_vegtype) THEN
  yrecfm='VEGTYPE'
  ycomment='(X_Y_VEGTYPE)'
  CALL write_surf(hselect, hprogram,yrecfm,k%XVEGTYPE,iresp,hcomment=ycomment)
ENDIF
#endif
!
!* orography
!
yrecfm='ZS'
ycomment='ZS'
 CALL write_surf(hselect, hprogram,yrecfm,s%XZS(:),iresp,hcomment=ycomment)
!
!* latitude, longitude
!
 CALL write_grid(hselect, hprogram,g%CGRID,g%XGRID_PAR,g%XLAT,g%XLON,g%XMESH_SIZE,iresp)
!
!
!* clay fraction
!
!
yrecfm='CLAY'
ycomment='X_Y_CLAY'
 CALL write_surf(hselect, hprogram,yrecfm,k%XCLAY(:,1),iresp,hcomment=ycomment)
!
!* sand fraction
!
yrecfm='SAND'
ycomment='X_Y_SAND'
 CALL write_surf(hselect, hprogram,yrecfm,k%XSAND(:,1),iresp,hcomment=ycomment)
!
!* soil organic carbon
!
yrecfm='SOCP'
ycomment=''
 CALL write_surf(hselect, hprogram,yrecfm,io%LSOCP,iresp,hcomment=ycomment)
!
IF(io%LSOCP)THEN
  !        
  ycomment='X_Y_SOC'
  yrecfm='SOC_TOP'
  CALL write_surf(hselect, hprogram,yrecfm,s%XSOC(:,1),iresp,hcomment=ycomment)
  yrecfm='SOC_SUB'
  CALL write_surf(hselect, hprogram,yrecfm,s%XSOC(:,2),iresp,hcomment=ycomment)
  !
ENDIF
!
!* permafrost distribution
!
yrecfm='PERMAFROST'
ycomment=''
 CALL write_surf(hselect, hprogram,yrecfm,io%LPERM,iresp,hcomment=ycomment)
!
IF(io%LPERM)THEN
  ycomment='X_Y_PERM'
  yrecfm='PERM'
  CALL write_surf(hselect, hprogram,yrecfm,k%XPERM(:),iresp,hcomment=ycomment)
ENDIF
!
!SOILNOX
!
yrecfm='NO'
ycomment=''
 CALL write_surf(hselect, hprogram,yrecfm,io%LNOF,iresp,hcomment=ycomment)
!
IF (io%LNOF) THEN
  !
  yrecfm='PH'
  ycomment='X_Y_PH'
  CALL write_surf(hselect, hprogram,yrecfm,s%XPH(:),iresp,hcomment=ycomment)
  !
  yrecfm='FERT'
  ycomment='X_Y_FERT'
  CALL write_surf(hselect, hprogram,yrecfm,s%XFERT(:),iresp,hcomment=ycomment)
  !
ENDIF
!
!* subgrid-scale orography parameters to compute dynamical roughness length
!
yrecfm='AOSIP'
ycomment='X_Y_AOSIP'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XAOSIP,iresp,hcomment=ycomment)
!
yrecfm='AOSIM'
ycomment='X_Y_AOSIM'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XAOSIM,iresp,hcomment=ycomment)
!
yrecfm='AOSJP'
ycomment='X_Y_AOSJP'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XAOSJP,iresp,hcomment=ycomment)
!
yrecfm='AOSJM'
ycomment='X_Y_AOSJM'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XAOSJM,iresp,hcomment=ycomment)
!
yrecfm='HO2IP'
ycomment='X_Y_HO2IP'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XHO2IP,iresp,hcomment=ycomment)
!
yrecfm='HO2IM'
ycomment='X_Y_HO2IM'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XHO2IM,iresp,hcomment=ycomment)
!
yrecfm='HO2JP'
ycomment='X_Y_HO2JP'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XHO2JP,iresp,hcomment=ycomment)
!
yrecfm='HO2JM'
ycomment='X_Y_HO2JM'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XHO2JM,iresp,hcomment=ycomment)
!
yrecfm='SSO_SLOPE'
ycomment='X_Y_SSO_SLOPE (-)'
 CALL write_surf(hselect, hprogram,yrecfm,iss%XSSO_SLOPE,iresp,hcomment=ycomment)
!
!* orographic runoff coefficient
!
yrecfm='RUNOFFB'
ycomment='X_Y_RUNOFFB'
 CALL write_surf(hselect, hprogram,yrecfm,k%XRUNOFFB,iresp,hcomment=ycomment)
!
!* subgrid drainage coefficient
!
yrecfm='WDRAIN'
ycomment='X_Y_WDRAIN'
 CALL write_surf(hselect, hprogram,yrecfm,k%XWDRAIN,iresp,hcomment=ycomment)
!
!* topographic index statistics
!
yrecfm='CTI'
ycomment=''
 CALL write_surf(hselect, hprogram,yrecfm,io%LCTI,iresp,hcomment=ycomment)
!
IF(io%LCTI)THEN
!
yrecfm='TI_MIN'
ycomment='X_Y_TI_MIN'
 CALL write_surf(hselect, hprogram,yrecfm,s%XTI_MIN,iresp,hcomment=ycomment)
!
yrecfm='TI_MAX'
ycomment='X_Y_TI_MAX'
 CALL write_surf(hselect, hprogram,yrecfm,s%XTI_MAX,iresp,hcomment=ycomment)
!
yrecfm='TI_MEAN'
ycomment='X_Y_TI_MEAN'
 CALL write_surf(hselect, hprogram,yrecfm,s%XTI_MEAN,iresp,hcomment=ycomment)
!
yrecfm='TI_STD'
ycomment='X_Y_TI_STD'
 CALL write_surf(hselect, hprogram,yrecfm,s%XTI_STD,iresp,hcomment=ycomment)
!
yrecfm='TI_SKEW'
ycomment='X_Y_TI_SKEW'
 CALL write_surf(hselect, hprogram,yrecfm,s%XTI_SKEW,iresp,hcomment=ycomment)
!
ENDIF
!
!-------------------------------------------------------------------------------
 CALL writesurf_pgd_isba_par_n(hselect, dtv, hprogram)
IF (hnature=='TSZ0') CALL writesurf_pgd_tsz0_par_n(hselect, dtz, hprogram)
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_ISBA_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE writesurf_pgd_isba_n