write_cover_tex_isba_par.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_cover_tex_isba_par (DTCO, HALBEDO, OTR_ML, &
                                           KPATCH,KLAYER,HISBA,HPHOTO,PSOILGRID)
!     ##########################
!
!!**** *WRITE_COVER_TEX* writes the ISBA data arrays into a tex file
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    08/01/98
!!
!!      P Le Moigne 09/2005 AGS modifs of L. Jarlan
!!      B. Decharme    2008 Bug if ZDMAX = XUNDEF
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
!
USE mode_write_cover_tex

USE modi_convert_cover_isba

USE modd_type_date_surf
USE modd_write_cover_tex,ONLY : ntex, cname, clang, nlines
USE modd_surf_par,       ONLY : xundef
USE modd_csts,           ONLY : xday

USE modd_data_cover_par, ONLY : jpcover, nvegtype
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
!
 CHARACTER(LEN=*), INTENT(IN) :: HALBEDO
LOGICAL, INTENT(IN) :: OTR_ML
!
INTEGER,          INTENT(IN) :: KPATCH! number of patch
INTEGER,          INTENT(IN) :: KLAYER! number of soil layers
 CHARACTER(LEN=*), INTENT(IN) :: HISBA ! type of soil (Force-Restore OR Diffusion)
 CHARACTER(LEN=*), INTENT(IN) :: HPHOTO! type of photosynthesis
REAL, DIMENSION(:),INTENT(IN) :: PSOILGRID ! reference grid for DIF
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
!
INTEGER :: JI,JJ,JIP
!
 CHARACTER(LEN=6), DIMENSION(12      ) :: YDATA_VEGPARAM! vegetation parameters
 CHARACTER(LEN=6)                      :: YSTRING6
 CHARACTER(LEN=24)                     :: YDATA_CV      ! Cv
 CHARACTER(LEN=24)                     :: YDATA_RE25    ! Re25
 CHARACTER(LEN=4)                      :: YDATA_STRESS  ! Stress def.
 CHARACTER(LEN=6), DIMENSION(12      ) :: YDATA_MONTH   ! monthly veg. parameters
 CHARACTER(LEN=6), DIMENSION(10      ) :: YDATA_LAYER   ! soil layer parameters
!
 CHARACTER(LEN=20) :: YFMT  ! fortran format
!
INTEGER           :: JPAGE ! current page when the number of classes
                           ! is too long to stand on one page only
INTEGER           :: JPATCH   ! loop counter
!
!
LOGICAL           :: GLINE ! flag to write an additional horizontal line
!
REAL, DIMENSION(JPCOVER,JPCOVER     ) :: ZCOVER
LOGICAL, DIMENSION(JPCOVER          ) :: GCOVER
REAL, DIMENSION(JPCOVER,KPATCH, 12  ) :: ZVEG, ZLAI, ZZ0VEG, ZEMIS_ECO, ZF2I
REAL, DIMENSION(JPCOVER,KPATCH      ) :: ZRSMIN,ZGAMMA,ZRGL,ZCV,             &
                                           ZALBNIR_VEG,ZALBVIS_VEG,ZALBUV_VEG, &
                                           ZGMES,ZBSLAI,ZLAIMIN,ZSEFOLD,       &
                                           ZH_TREE, ZGC, ZZ0_O_Z0H,            &
                                           ZWRMAX_CF, ZDMAX, ZRE25  
REAL, DIMENSION(JPCOVER,KLAYER,KPATCH):: ZDG, ZROOTFRAC
REAL, DIMENSION(JPCOVER,KPATCH)       :: ZDROOT
REAL, DIMENSION(JPCOVER,KPATCH)       :: ZDG2
INTEGER, DIMENSION(JPCOVER,KPATCH)    :: IWG_LAYER
LOGICAL, DIMENSION(JPCOVER,KPATCH, 12   ) :: GSTRESS
REAL(KIND=JPRB) :: ZHOOK_HANDLE

!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('WRITE_COVER_TEX_ISBA_PAR',0,zhook_handle)
IF (ntex==0 .AND. lhook) CALL dr_hook('WRITE_COVER_TEX_ISBA_PAR',1,zhook_handle)
IF (ntex==0) RETURN
!
gline=.false.
!
zcover(:,:) = 0.
DO ji=1,jpcover
  zcover(ji,ji) = 1.
END DO
!
gcover(:) = .true.
!
!ocl scalar
!
DO jj=1,12
  CALL convert_cover_isba(dtco, halbedo, &
                          hisba,otr_ml,3*jj-1,zcover,gcover,hphoto, 'NAT',         &
                            pveg=zveg(:,:,jj), plai=zlai(:,:,jj),            &
                            pz0=zz0veg(:,:,jj), pemis_eco=zemis_eco(:,:,jj), &
                            pf2i=zf2i(:,:,jj),ostress=gstress(:,:,jj)        )  
END DO

 CALL convert_cover_isba(dtco, halbedo, &
                          hisba,otr_ml,2,zcover,gcover,hphoto, 'NAT',            &
                        prsmin=zrsmin,pgamma=zgamma,pwrmax_cf=zwrmax_cf, &
                        prgl=zrgl,pcv=zcv,psoilgrid=psoilgrid,           &
                        pdg=zdg,kwg_layer=iwg_layer,pdroot=zdroot,       &
                        pdg2=zdg2,pz0_o_z0h=zz0_o_z0h,                   &
                        palbnir_veg=zalbnir_veg,palbvis_veg=zalbvis_veg, &
                        palbuv_veg=zalbuv_veg,                           &
                        prootfrac=zrootfrac,                             &
                        pgmes=zgmes,pbslai=zbslai,plaimin=zlaimin,       &
                        psefold=zsefold,pgc=zgc,pdmax=zdmax,             &
                        ph_tree=zh_tree,pre25=zre25                       )  
!
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
DO jpatch=1,kpatch
!
 IF (       (kpatch>=2  .AND. jpatch==1) &
         .OR. (kpatch>=7  .AND. jpatch==2) &
         .OR. (kpatch>=10 .AND. jpatch==3) ) cycle  
!
!
!
 ji=0
 DO 

  IF (ji==jpcover) EXIT

  IF (clang=='EN') THEN
    WRITE(ntex,*) '{\bf averaged leaf area index} (patch ',jpatch,'/',kpatch,') \\'
  ELSE
    WRITE(ntex,*) "{\bf indice foliaire moyenn\'e} (partition ",jpatch,'/',kpatch,') \\'
  END IF
  WRITE(ntex,*) '\medskip\'
  WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c|c|c|c||c|c||c||}'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '&&01&02&03&04&05&06&07&08&09&10&11&12&$d_2$&$d_3$&$h$'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  JIP=0
  DO
    IF (JI==JPCOVER) EXIT
    JI=JI+1
    IF (ZLAI(JI,JPATCH,1)/=XUNDEF) THEN
      JIP=JIP+1
      DO JJ=1,12
.OR.        IF (ZLAI(JI,JPATCH,JJ)==0.  ZLAI(JI,JPATCH,JJ)==XUNDEF) THEN
          YDATA_MONTH(JJ) = ' -  '
        ELSE
          WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NBT(ZLAI(JI,JPATCH,JJ),1),'.',DEC(ZLAI(JI,JPATCH,JJ),1),')'
          WRITE(YSTRING6, FMT=YFMT) ZLAI(JI,JPATCH,JJ)
          YDATA_MONTH(JJ) = YSTRING6
        END IF
      END DO

      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZDG(JI,2,JPATCH),1),'.',DEC(ZDG(JI,2,JPATCH),1),')'
      WRITE(YSTRING6, FMT=YFMT) ZDG(JI,2,JPATCH)
      YDATA_VEGPARAM(1) = YSTRING6
      IF (KLAYER>=3) THEN
        WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZDG(JI,3,JPATCH),1),'.',DEC(ZDG(JI,3,JPATCH),1),')'
        WRITE(YSTRING6, FMT=YFMT) ZDG(JI,3,JPATCH)
        YDATA_VEGPARAM(2) = YSTRING6
      ELSE
        YDATA_VEGPARAM(2) = ' -    '
      END IF

.OR.      IF (ZH_TREE(JI,JPATCH)==XUNDEF  ZH_TREE(JI,JPATCH)==0.) THEN
        YDATA_VEGPARAM(3) = ' -    '
      ELSE
        WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZH_TREE(JI,JPATCH)),'.',DEC(ZH_TREE(JI,JPATCH)),')'
        WRITE(YSTRING6, FMT=YFMT) ZH_TREE(JI,JPATCH)
        YDATA_VEGPARAM(3) = YSTRING6
      END IF

      WRITE(NTEX, FMT=*) &
          JI,' & ',CNAME(JI),' & ',YDATA_MONTH(1), ' & ',YDATA_MONTH(2), ' & ', &
            YDATA_MONTH(3),' & ',YDATA_MONTH(4), ' & ',YDATA_MONTH(5), ' & ', &
            YDATA_MONTH(6),' & ',YDATA_MONTH(7), ' & ',YDATA_MONTH(8), ' & ', &
            YDATA_MONTH(9),' & ',YDATA_MONTH(10),' & ',YDATA_MONTH(11),' & ', &
            YDATA_MONTH(12),' & ',YDATA_VEGPARAM(1),' & ',YDATA_VEGPARAM(2),' & ',YDATA_VEGPARAM(3),' '  

      WRITE(NTEX,*) '\hline'
      GLINE=.TRUE.
    END IF
    CALL HLINE(NTEX,GLINE,JI)
    IF (JIP==NLINES) EXIT  
  END DO
  WRITE(NTEX,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(NTEX,*) '\clearpage'
!
ENDDO
!-------------------------------------------------------------------------------
!
!
JI=0
DO 
  IF (JI==JPCOVER) EXIT

  IF (CLANG=='en') THEN
    WRITE(NTEX,*) '{\bf ground layer depth} (from surface, patch ',JPATCH,'/',KPATCH,') '
  ELSE
    WRITE(NTEX,*) "{\bf profondeur des couche de sol} (depuis la surface, partition ",JPATCH,'/',KPATCH,') '
  END IF
  WRITE(NTEX,*) '\medskip
  WRITE(ntex,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c|c||}'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '&&01&02&03&04&05&06&07&08&09&10\\'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  jip=0
  DO
    IF (ji==jpcover) EXIT
    ji=ji+1
    ydata_layer(:) = ' -  '
    IF (zdg(ji,1,jpatch)/=xundef) THEN
      jip=jip+1
      DO jj=1,min(klayer,10)
        WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zdg(ji,jj,jpatch),1),'.',dec(zdg(ji,jj,jpatch),1),')'
        WRITE(ystring6, fmt=yfmt) zdg(ji,jj,jpatch)
        ydata_layer(jj) = ystring6
      END DO

      WRITE(ntex, fmt=*) &
          ji,' & ',cname(ji),' & ',ydata_layer(1), ' & ',ydata_layer(2), ' & ', &
            ydata_layer(3),' & ',ydata_layer(4), ' & ',ydata_layer(5), ' & ', &
            ydata_layer(6),' & ',ydata_layer(7), ' & ',ydata_layer(8), ' & ', &
            ydata_layer(9),' & ',ydata_layer(10),' \\'  

      WRITE(ntex,*) '\hline'
      gline=.true.
    END IF
    CALL hline(ntex,gline,ji)
    IF (jip==nlines) EXIT
  END DO
  WRITE(ntex,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(ntex,*) '\clearpage'
!
ENDDO
!-------------------------------------------------------------------------------
!
IF (hisba=='DIF') THEN
!
ji=0
DO 
  IF (ji==jpcover) EXIT  

  IF (clang=='EN') THEN
    WRITE(ntex,*) '{\bf cumulative root fraction} (patch ',jpatch,'/',kpatch,') \\'
  ELSE
    WRITE(ntex,*) "{\bf fraction de racines cumul\'ee} (partition ",jpatch,'/',kpatch,') \\'
  END IF
  WRITE(ntex,*) '\medskip\'
  WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c|c||}'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '&&01&02&03&04&05&06&07&08&09&10'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  JIP=0
  DO
    IF (JI==JPCOVER) EXIT
    JI=JI+1
    YDATA_LAYER(:) = ' -  '
    IF (ZROOTFRAC(JI,1,JPATCH)/=XUNDEF) THEN
      JIP=JIP+1
      DO JJ=1,MIN(KLAYER,10)
        WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NBT(ZROOTFRAC(JI,JJ,JPATCH),1),'.',DEC(ZROOTFRAC(JI,JJ,JPATCH),1),')'
        WRITE(YSTRING6, FMT=YFMT) ZROOTFRAC(JI,JJ,JPATCH)
        YDATA_LAYER(JJ) = YSTRING6
      END DO

      WRITE(NTEX, FMT=*) &
          JI,' & ',CNAME(JI),' & ',YDATA_LAYER(1), ' & ',YDATA_LAYER(2), ' & ', &
            YDATA_LAYER(3),' & ',YDATA_LAYER(4), ' & ',YDATA_LAYER(5), ' & ', &
            YDATA_LAYER(6),' & ',YDATA_LAYER(7), ' & ',YDATA_LAYER(8), ' & ', &
            YDATA_LAYER(9),' & ',YDATA_LAYER(10),' '  

      WRITE(NTEX,*) '\hline'
      GLINE=.TRUE.
    END IF
    CALL HLINE(NTEX,GLINE,JI)
    IF (JIP==NLINES) EXIT
  END DO
  WRITE(NTEX,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(NTEX,*) '\clearpage'

ENDDO
END IF
!
!-------------------------------------------------------------------------------
!
JI=0
DO 
  IF (JI==JPCOVER) EXIT

  IF (CLANG=='en') THEN
    WRITE(NTEX,*) '{\bf vegetation fraction(over natural or agricultural areas)} (patch ',JPATCH,'/',KPATCH,') '
  ELSE
    WRITE(NTEX,*) "{\bf fraction de v\'eg\'etation (sur les surfaces naturelles ou cultiv\'ees)} (partition ", &
  JPATCH,'/',KPATCH,') '  
  END IF
  WRITE(NTEX,*) '\medskip
  WRITE(ntex,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c|c|c|c||}'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '&&01&02&03&04&05&06&07&08&09&10&11&12\\'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  jip=0
  DO
    IF (ji==jpcover) EXIT
    ji=ji+1
    IF (zlai(ji,jpatch,1)/=xundef) THEN
      jip=jip+1
      DO jj=1,12
        IF (zveg(ji,jpatch,jj)==0.) THEN
          ydata_month(jj) = ' -  '
        ELSE
          WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zveg(ji,jpatch,jj),2),'.',dec(zveg(ji,jpatch,jj),2),')'
          WRITE(ystring6, fmt=yfmt) zveg(ji,jpatch,jj)

          ydata_month(jj) = ystring6
        END IF
      END DO

      WRITE(ntex, fmt=*) &
          ji,' & ',cname(ji),' & ',ydata_month(1), ' & ',ydata_month(2), ' & ', &
            ydata_month(3),' & ',ydata_month(4), ' & ',ydata_month(5), ' & ', &
            ydata_month(6),' & ',ydata_month(7), ' & ',ydata_month(8), ' & ', &
            ydata_month(9),' & ',ydata_month(10),' & ',ydata_month(11),' & ', &
            ydata_month(12),' \\'  

      WRITE(ntex,*) '\hline'
      gline=.true.
    END IF
    CALL hline(ntex,gline,ji)
    IF (jip==nlines) EXIT
  END DO
  WRITE(ntex,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(ntex,*) '\clearpage'
!
ENDDO
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
!
ji=0
DO 
  IF (ji==jpcover) EXIT

  IF (clang=='EN') THEN
    WRITE(ntex,*) '{\bf surface roughness length for momentum: $z_{0}$} (patch ',jpatch,'/',kpatch,') \\'
  ELSE
    WRITE(ntex,*) "{\bf longueur de rugosit\'e de la surface (qdm): $z_{0}$} (partition ",jpatch,'/',kpatch,') \\'
  END IF
  WRITE(ntex,*) '\medskip\'
  WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c|c|c|c||}'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '&&01&02&03&04&05&06&07&08&09&10&11&12'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  JIP=0
  DO
    IF (JI==JPCOVER) EXIT
    JI=JI+1

    IF (ZLAI(JI,JPATCH,1)/=XUNDEF) THEN
      JIP=JIP+1
      DO JJ=1,12
        WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NBT(ZZ0VEG(JI,JPATCH,JJ)),'.',DEC(ZZ0VEG(JI,JPATCH,JJ)),')'
        WRITE(YSTRING6, FMT=YFMT) ZZ0VEG(JI,JPATCH,JJ)
        YDATA_MONTH(JJ) = YSTRING6
      END DO

      WRITE(NTEX, FMT=*) &
          JI,' & ',CNAME(JI),' & ',YDATA_MONTH(1), ' & ',YDATA_MONTH(2), ' & ', &
            YDATA_MONTH(3),' & ',YDATA_MONTH(4), ' & ',YDATA_MONTH(5), ' & ', &
            YDATA_MONTH(6),' & ',YDATA_MONTH(7), ' & ',YDATA_MONTH(8), ' & ', &
            YDATA_MONTH(9),' & ',YDATA_MONTH(10),' & ',YDATA_MONTH(11),' & ', &
            YDATA_MONTH(12),' '  

      WRITE(NTEX,*) '\hline'
      GLINE=.TRUE.
    END IF
    CALL HLINE(NTEX,GLINE,JI)
    IF (JIP==NLINES) EXIT
   END DO

  WRITE(NTEX,*) '\end{tabular}'

!
!-------------------------------------------------------------------------------
!
  WRITE(NTEX,*) '\clearpage'
!
ENDDO
!-------------------------------------------------------------------------------
!
!
JI=0
DO 
  IF (JI==JPCOVER) EXIT

  IF (CLANG=='en') THEN
    WRITE(NTEX,*) '{\bf emissivity of natural continental surfaces} (patch ',JPATCH,'/',KPATCH,') '
  ELSE
    WRITE(NTEX,*) "{\bf \'emissivit\'e des surfaces continentales naturelles} (partition ",JPATCH,'/',KPATCH,') '
  END IF
  WRITE(NTEX,*) '\medskip
  WRITE(ntex,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c|c|c|c||}'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '&&01&02&03&04&05&06&07&08&09&10&11&12\\'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  jip=0
  DO
    IF (ji==jpcover) EXIT
    ji=ji+1
    IF (zlai(ji,jpatch,1)/=xundef) THEN
      jip=jip+1
      DO jj=1,12
        IF (zemis_eco(ji,jpatch,jj)==1.) THEN
          ydata_month(jj) = ' 1.  '
        ELSE
          WRITE(ystring6, fmt='(F3.2)') zemis_eco(ji,jpatch,jj)
          ydata_month(jj) = ystring6
        END IF
      END DO

      WRITE(ntex, fmt=*) &
          ji,' & ',cname(ji),' & ',ydata_month(1), ' & ',ydata_month(2), ' & ', &
            ydata_month(3),' & ',ydata_month(4), ' & ',ydata_month(5), ' & ', &
            ydata_month(6),' & ',ydata_month(7), ' & ',ydata_month(8), ' & ', &
            ydata_month(9),' & ',ydata_month(10),' & ',ydata_month(11),' & ', &
            ydata_month(12),' \\'  

      WRITE(ntex,*) '\hline'
      gline=.true.
    END IF
    CALL hline(ntex,gline,ji)
    IF (jip==nlines) EXIT
  END DO
  WRITE(ntex,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(ntex,*) '\clearpage'
!
ENDDO
!-------------------------------------------------------------------------------
!
!
ji=0
DO 
  IF (ji==jpcover) EXIT

  IF (clang=='EN') THEN
    WRITE(ntex,*) '{\bf other vegetation parameters} (1) (patch ',jpatch,'/',kpatch,') \\'
  ELSE
    WRITE(ntex,*) "{\bf autres param\`etres de v\'eg\'etation} (1) (partition ",jpatch,'/',kpatch,') \\'
  END IF
  WRITE(ntex,*) '\medskip\'
  WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c||}'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) &
    '&&$\alpha_{nir}$&$\alpha_{vis}$&$\alpha_{uv}$&$r_{s_{min}}$&$\gamma$&$rgl$&$c_{w_{r_{max}}}$&$z_0$/$z_{0_h}$&$c_v$'  
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  JIP=0
  DO
    IF (JI==JPCOVER) EXIT
    JI=JI+1
.AND..AND.    IF (ZLAI(JI,JPATCH,1)>0.  ZVEG(JI,JPATCH,1)>0.  ZVEG(JI,JPATCH,1)/=XUNDEF) THEN
      JIP=JIP+1
      WRITE(YSTRING6, FMT='(f3.2)') ZALBNIR_VEG(JI,JPATCH)
      YDATA_VEGPARAM(4) = YSTRING6
      WRITE(YSTRING6, FMT='(f3.2)') ZALBVIS_VEG(JI,JPATCH)
      YDATA_VEGPARAM(5) = YSTRING6
      WRITE(YSTRING6, FMT='(f3.2)') ZALBUV_VEG (JI,JPATCH)
      YDATA_VEGPARAM(6) = YSTRING6
      ZRSMIN(JI,JPATCH) = NINT(ZRSMIN(JI,JPATCH))
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZRSMIN(JI,JPATCH)),'.',DEC(ZRSMIN(JI,JPATCH)),')'
      WRITE(YSTRING6, FMT=YFMT) ZRSMIN(JI,JPATCH)
      YDATA_VEGPARAM(7) = YSTRING6
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NBT(ZGAMMA(JI,JPATCH)),'.',DEC(ZGAMMA(JI,JPATCH)),')'
      WRITE(YSTRING6, FMT=YFMT) ZGAMMA(JI,JPATCH)
      YDATA_VEGPARAM(8) = YSTRING6
      ZRGL(JI,JPATCH) = NINT(ZRGL(JI,JPATCH))
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZRGL(JI,JPATCH)),'.',DEC(ZRGL(JI,JPATCH)),')'
      WRITE(YSTRING6, FMT=YFMT) ZRGL(JI,JPATCH)
      YDATA_VEGPARAM(9) = YSTRING6
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NBT(ZWRMAX_CF(JI,JPATCH)),'.',DEC(ZWRMAX_CF(JI,JPATCH)),')'
      WRITE(YSTRING6, FMT=YFMT) ZWRMAX_CF(JI,JPATCH)
      YDATA_VEGPARAM(10) = YSTRING6
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZZ0_O_Z0H(JI,JPATCH)),'.',DEC(ZZ0_O_Z0H(JI,JPATCH)),')'
      WRITE(YSTRING6, FMT=YFMT) ZZ0_O_Z0H(JI,JPATCH)
      YDATA_VEGPARAM(11) = YSTRING6
      WRITE(YFMT,'(a2,i1,a1,i1,a5)') '(f',NB(10000. *ZCV(JI,JPATCH)),'.',DEC(10000. *ZCV(JI,JPATCH)),',a10)'
      WRITE(YDATA_CV, FMT=YFMT) 10000. *ZCV(JI,JPATCH),' $10^{-4}$'

      WRITE(NTEX, FMT=*) &
             JI,' & ',CNAME(JI),' & ', &
                     YDATA_VEGPARAM(4),' & ',YDATA_VEGPARAM(5),' & ', &
                     YDATA_VEGPARAM(6),' & ',YDATA_VEGPARAM(7),' & ', &
                     YDATA_VEGPARAM(8),' & ',YDATA_VEGPARAM(9),' & ', &
                     YDATA_VEGPARAM(10),' & ',YDATA_VEGPARAM(11),' & ', &
                     YDATA_CV,' '  
      WRITE(NTEX,*) '\hline'
      GLINE=.TRUE.
    END IF
    CALL HLINE(NTEX,GLINE,JI)
    IF (JIP==NLINES) EXIT
  END DO
  WRITE(NTEX,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(NTEX,*) '\clearpage'
!
ENDDO
!-------------------------------------------------------------------------------
!
!
JI=0
DO 
  IF (JI==JPCOVER) EXIT

  IF (CLANG=='en') THEN
    WRITE(NTEX,*) '{\bf other vegetation parameters} (2) (patch ',JPATCH,'/',KPATCH,') '
  ELSE
    WRITE(NTEX,*) "{\bf autres param\`etres de v\'eg\'etation} (2) (partition ",JPATCH,'/',KPATCH,') '
  END IF
  WRITE(NTEX,*) '\medskip
  WRITE(ntex,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c||}'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '&&$gm$&$B/lai$&$lai_{_{m}}$&$e_{_{fold}}$&$G_c$&$D_{max}$&$f_{2i}$&stress&Re$_{25}$\\'
  WRITE(ntex,*) '\hline'
  WRITE(ntex,*) '\hline'
  jip=0
  DO
    IF (ji==jpcover) EXIT
    ji=ji+1

    IF (zlai(ji,jpatch,1)>0. .AND. zveg(ji,jpatch,1)>0. .AND. zveg(ji,jpatch,1)/=xundef) THEN
      jip=jip+1
      WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zgmes(ji,jpatch)),'.',dec(zgmes(ji,jpatch)),')'
      WRITE(ystring6, fmt=yfmt) zgmes(ji,jpatch)
      ydata_vegparam(1) = ystring6
      WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zbslai(ji,jpatch)),'.',dec(zbslai(ji,jpatch)),')'
      WRITE(ystring6, fmt=yfmt) zbslai(ji,jpatch)
      ydata_vegparam(2) = ystring6
      WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zlaimin(ji,jpatch)),'.',dec(zlaimin(ji,jpatch)),')'
      WRITE(ystring6, fmt=yfmt) zlaimin(ji,jpatch)
      ydata_vegparam(3) = ystring6
      zsefold(ji,jpatch) = nint(zsefold(ji,jpatch)/xday)
      WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nb(zsefold(ji,jpatch)),'.',dec(zsefold(ji,jpatch)),')'
      WRITE(ystring6, fmt=yfmt) zsefold(ji,jpatch)
      ydata_vegparam(4) = ystring6
      WRITE(yfmt,'(A2,I1,A1,I1,A5)') '(F',nb(10000. *zgc(ji,jpatch)),'.',dec(10000. *zgc(ji,jpatch)),',A10)'
      WRITE(ydata_cv, fmt=yfmt) 10000. *zgc(ji,jpatch),' $10^{-4}$'
      IF (zdmax(ji,jpatch) /= xundef)THEN
         WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zdmax(ji,jpatch)),'.',dec(zdmax(ji,jpatch)),')'
         WRITE(ystring6, fmt=yfmt) zdmax(ji,jpatch)
      ELSE
         ystring6 ='  -   '
      ENDIF
      ydata_vegparam(5) = ystring6
      IF (gstress(ji,jpatch,1)) THEN
        ydata_stress='def.'
      ELSE
        ydata_stress='off.'
      END IF
      WRITE(yfmt,'(A2,I1,A1,I1,A1)') '(F',nbt(zf2i(ji,jpatch,1)),'.',dec(zf2i(ji,jpatch,1)),')'
      WRITE(ystring6, fmt=yfmt) zf2i(ji,jpatch,1)
      ydata_vegparam(6) = ystring6
      WRITE(yfmt,'(A2,I1,A1,I1,A5)') '(F',nb(10000000. *zre25(ji,jpatch)),'.',dec(10000000. *zre25(ji,jpatch)),',A10)'
      WRITE(ydata_re25, fmt=yfmt) 10000000. *zre25(ji,jpatch),' $10^{-7}$'
      WRITE(ntex, fmt=*) &
             ji,' & ',cname(ji),' & ', &
                     ydata_vegparam(1),' & ',ydata_vegparam(2),' & ', &
                     ydata_vegparam(3),' & ',ydata_vegparam(4),' & ', &
                     ydata_cv,         ' & ',ydata_vegparam(5),' & ', &
                     ydata_vegparam(6),' & ',                         &
                     ydata_stress ,' & ',ydata_re25,' \\'  
      WRITE(ntex,*) '\hline'
      gline=.true.
    END IF
    CALL hline(ntex,gline,ji)
    IF (jip==nlines) EXIT
  END DO
  WRITE(ntex,*) '\end{tabular}'
!
!-------------------------------------------------------------------------------
!
  WRITE(ntex,*) '\clearpage'
  
 END DO
!
END DO 
IF (lhook) CALL dr_hook('WRITE_COVER_TEX_ISBA_PAR',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
!
END SUBROUTINE write_cover_tex_isba_par