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 (DGU, &
                                        dti, dtz, ig, i, u, &
                                       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
!!
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
!
!
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
!
USE modd_data_isba_n, ONLY : data_isba_t
USE modd_data_tsz0_n, ONLY : data_tsz0_t
USE modd_isba_grid_n, ONLY : isba_grid_t
USE modd_isba_n, ONLY : isba_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
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 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(data_isba_t), INTENT(INOUT) :: dti
TYPE(data_tsz0_t), INTENT(INOUT) :: dtz
TYPE(isba_grid_t), INTENT(INOUT) :: ig
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
 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(dgu, u, &
                 hprogram,yrecfm,i%CISBA,iresp,hcomment=ycomment)
!
!* Pedo-transfert function
!
yrecfm='PEDOTF'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%CPEDOTF,iresp,hcomment=ycomment)
!
!* type of photosynthesis
!
yrecfm='PHOTO'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%CPHOTO,iresp,hcomment=ycomment)
!
!* new radiative transfert
!
yrecfm='TR_ML'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LTR_ML,iresp,hcomment=ycomment)
!
!* threshold to remove little fractions of patches
!
yrecfm='RM_PATCH'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XRM_PATCH,iresp,hcomment=ycomment)

!* number of soil layers
!
yrecfm='GROUND_LAYER'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%NGROUND_LAYER,iresp,hcomment=ycomment)
!
!* Reference grid for DIF
!
IF(i%CISBA=='DIF') THEN
  DO jlayer=1,i%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(dgu, u, &
                 hprogram,yrecfm,i%XSOILGRID(jlayer),iresp,hcomment=ycomment)
  END DO 
ENDIF
!
!* number of biomass pools
!
yrecfm='NBIOMASS'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%NNBIOMASS,iresp,hcomment=ycomment)
!
!* number of tiles
!
yrecfm='PATCH_NUMBER'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%NPATCH,iresp,hcomment=ycomment)
!
!* flag indicating if fields are computed from ecoclimap or not
!
yrecfm='ECOCLIMAP'
ycomment=yrecfm
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LECOCLIMAP,iresp,hcomment=ycomment)
!
!* logical vector indicating for which patches MEB should be applied
!
yrecfm='MEB_PATCH'
ycomment='(LOGICAL LIST)'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LMEB_PATCH(:),iresp,hcomment=ycomment,hdir='-')
!
isize_lmeb_patch = count(i%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(dgu, u, &
                 hprogram,yrecfm,i%LFORC_MEASURE,iresp,hcomment=ycomment)
!
!* flag indicating if litter layer is used or not
!
   yrecfm='MEB_LITTER'
   ycomment=yrecfm
   CALL write_surf(dgu, u, &
                hprogram,yrecfm,i%LMEB_LITTER,iresp,hcomment=ycomment)
!
!* flag indicating if ground resistance is used or not
!
   yrecfm='MEB_GNDRES'
   ycomment=yrecfm
   CALL write_surf(dgu, u, &
                hprogram,yrecfm,i%LMEB_GNDRES,iresp,hcomment=ycomment)
!
ENDIF
!
!*       2.     Physiographic data fields:
!               -------------------------
!
!* cover classes
!
yrecfm='COVER_LIST'
ycomment='(LOGICAL LIST)'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LCOVER(:),iresp,hcomment=ycomment,hdir='-')
!
ycomment='COVER FIELDS'
 CALL write_surf_cov(dgu, u, &
                     hprogram,'COVER',i%XCOVER(:,:),i%LCOVER,iresp,hcomment=ycomment)
!
!* orography
!
yrecfm='ZS'
ycomment='ZS'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XZS(:),iresp,hcomment=ycomment)
!
!* latitude, longitude
!
 CALL write_grid(dgu, u, &
                 hprogram,ig%CGRID,ig%XGRID_PAR,ig%XLAT,ig%XLON,ig%XMESH_SIZE,iresp,i%XZ0EFFJPDIR)
!
!
!* clay fraction
!
!
yrecfm='CLAY'
ycomment='X_Y_CLAY'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XCLAY(:,1),iresp,hcomment=ycomment)
!
!* sand fraction
!
yrecfm='SAND'
ycomment='X_Y_SAND'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XSAND(:,1),iresp,hcomment=ycomment)
!
!* soil organic carbon
!
yrecfm='SOCP'
ycomment=''
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LSOCP,iresp,hcomment=ycomment)
!
IF(i%LSOCP)THEN
  !        
  ycomment='X_Y_SOC'
  yrecfm='SOC_TOP'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XSOC(:,1),iresp,hcomment=ycomment)
  yrecfm='SOC_SUB'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XSOC(:,2),iresp,hcomment=ycomment)
  !
ENDIF
!
!* permafrost distribution
!
yrecfm='PERMAFROST'
ycomment=''
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LPERM,iresp,hcomment=ycomment)
!
IF(i%LPERM)THEN
  ycomment='X_Y_PERM'
  yrecfm='PERM'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XPERM(:),iresp,hcomment=ycomment)
ENDIF
!
!* groundwater distribution
!
yrecfm='GWKEY'
ycomment=''
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LGW,iresp,hcomment=ycomment)
!
IF(i%LGW)THEN
  ycomment='X_Y_GWFRAC'
  yrecfm='GWFRAC'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XGW(:),iresp,hcomment=ycomment)
ENDIF
!
!SOILNOX
!
yrecfm='NO'
ycomment=''
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LNOF,iresp,hcomment=ycomment)
!
IF (i%LNOF) THEN
  !
  yrecfm='PH'
  ycomment='X_Y_PH'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XPH(:),iresp,hcomment=ycomment)
  !
  yrecfm='FERT'
  ycomment='X_Y_FERT'
  CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XFERT(:),iresp,hcomment=ycomment)
  !
ENDIF
!
!* subgrid-scale orography parameters to compute dynamical roughness length
!
yrecfm='AOSIP'
ycomment='X_Y_AOSIP'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XAOSIP,iresp,hcomment=ycomment)
!
yrecfm='AOSIM'
ycomment='X_Y_AOSIM'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XAOSIM,iresp,hcomment=ycomment)
!
yrecfm='AOSJP'
ycomment='X_Y_AOSJP'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XAOSJP,iresp,hcomment=ycomment)
!
yrecfm='AOSJM'
ycomment='X_Y_AOSJM'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XAOSJM,iresp,hcomment=ycomment)
!
yrecfm='HO2IP'
ycomment='X_Y_HO2IP'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XHO2IP,iresp,hcomment=ycomment)
!
yrecfm='HO2IM'
ycomment='X_Y_HO2IM'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XHO2IM,iresp,hcomment=ycomment)
!
yrecfm='HO2JP'
ycomment='X_Y_HO2JP'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XHO2JP,iresp,hcomment=ycomment)
!
yrecfm='HO2JM'
ycomment='X_Y_HO2JM'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XHO2JM,iresp,hcomment=ycomment)
!
yrecfm='SSO_SLOPE'
ycomment='X_Y_SSO_SLOPE (-)'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XSSO_SLOPE,iresp,hcomment=ycomment)
!
!* orographic runoff coefficient
!
yrecfm='RUNOFFB'
ycomment='X_Y_RUNOFFB'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XRUNOFFB,iresp,hcomment=ycomment)
!
!* subgrid drainage coefficient
!
yrecfm='WDRAIN'
ycomment='X_Y_WDRAIN'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XWDRAIN,iresp,hcomment=ycomment)
!
!* topographic index statistics
!
yrecfm='CTI'
ycomment=''
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%LCTI,iresp,hcomment=ycomment)
!
IF(i%LCTI)THEN
!
yrecfm='TI_MIN'
ycomment='X_Y_TI_MIN'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XTI_MIN,iresp,hcomment=ycomment)
!
yrecfm='TI_MAX'
ycomment='X_Y_TI_MAX'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XTI_MAX,iresp,hcomment=ycomment)
!
yrecfm='TI_MEAN'
ycomment='X_Y_TI_MEAN'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XTI_MEAN,iresp,hcomment=ycomment)
!
yrecfm='TI_STD'
ycomment='X_Y_TI_STD'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XTI_STD,iresp,hcomment=ycomment)
!
yrecfm='TI_SKEW'
ycomment='X_Y_TI_SKEW'
 CALL write_surf(dgu, u, &
                 hprogram,yrecfm,i%XTI_SKEW,iresp,hcomment=ycomment)
!
ENDIF
!
!-------------------------------------------------------------------------------
 CALL writesurf_pgd_isba_par_n(dgu, u, &
                               dti, &
                               hprogram)
IF (u%CNATURE=='TSZ0') CALL writesurf_pgd_tsz0_par_n(dgu, u, &
                                                     dtz, &
                                                     hprogram)
!
IF (lhook) CALL dr_hook('WRITESURF_PGD_ISBA_N',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE writesurf_pgd_isba_n