write_cover_tex_water.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_water
!     ##########################
!
!!**** *WRITE_COVER_TEX* writes the water data arrays into a tex file
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    08/01/98
!!    03/2011 E. Bazile (MK10) albedo from Marat Khairoutdinov
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE mode_write_cover_tex
!
USE modi_albedo_ta96
!
USE modd_csts,           ONLY : xpi
USE modd_write_cover_tex,ONLY : ntex, cname, clang, nlines
USE modd_surf_par,       ONLY : xundef
!
USE modd_data_cover,     ONLY : xdata_sea, xdata_water
USE modd_data_cover_par, ONLY : jpcover
USE modd_water_par,      ONLY : xalbwat, xalbsca_wat, xemiswat
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
!
!
INTEGER :: i,ip
!
 CHARACTER(LEN=6)  :: ystring6
!
 CHARACTER(LEN=200):: yfmt  ! fortran format
!
LOGICAL           :: gline ! flag to write an additional horizontal line
!
 CHARACTER(LEN=6), DIMENSION(6) :: ydata_water   ! water parameters
REAL,             DIMENSION(1) :: zzenith       ! zenithal angle
REAL,             DIMENSION(1) :: zalbedo       ! direct albedo
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('WRITE_COVER_TEX_WATER',0,zhook_handle)
IF (ntex==0 .AND. lhook) CALL dr_hook('WRITE_COVER_TEX_WATER',1,zhook_handle)
IF (ntex==0) RETURN
!
i=0
DO 
  IF (i==jpcover) EXIT
  IF (clang=='EN') THEN
    WRITE(ntex,*) '{\bf water parameters}\\'
  ELSE
    WRITE(ntex,*) '{\bf param\`etres aquatiques}\\'
  END IF
  WRITE(ntex,*) '\medskip\'
  WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c||}'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '&&$\alpha$&$\alpha_{dir}$&$\alpha_{sca}$&$\epsilon$'
  WRITE(NTEX,*) '\hline'
  WRITE(NTEX,*) '\hline'
  IP=0
  DO
    IF (I==JPCOVER) EXIT
    I=I+1
    IF (XDATA_SEA(I) + XDATA_WATER(I)>0.) THEN
      IP=IP+1
      WRITE(YSTRING6, FMT='(f3.2)') XALBWAT
      YDATA_WATER(1) = YSTRING6

      ZZENITH(:) = 0.
      ZALBEDO(:) = ALBEDO_TA96(ZZENITH)
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZALBEDO(1)),'.',DEC(ZALBEDO(1)),')'
      WRITE(YSTRING6, FMT=YFMT) ZALBEDO(1)
      YDATA_WATER(2) = YSTRING6

      ZZENITH(:) = XPI/2.
      ZALBEDO(:) = ALBEDO_TA96(ZZENITH)  
      WRITE(YFMT,'(a2,i1,a1,i1,a1)') '(f',NB(ZALBEDO(1)),'.',DEC(ZALBEDO(1)),')'
      WRITE(YSTRING6, FMT=YFMT) ZALBEDO(1)
      YDATA_WATER(3) = YSTRING6

      WRITE(YSTRING6, FMT='(f3.2)') XALBSCA_WAT
      YDATA_WATER(4) = YSTRING6
      WRITE(YSTRING6, FMT='(f4.2)') XEMISWAT
      YDATA_WATER(5) = YSTRING6

      WRITE(NTEX,FMT=*) &
          I,' & ',CNAME(I),' & ',YDATA_WATER(1),' & ', &
            YDATA_WATER(2),'-',YDATA_WATER(3),' & ', &
            YDATA_WATER(4),' & ', &
            YDATA_WATER(5),' '  
      WRITE(NTEX,*) '\hline'
      GLINE=.TRUE.
    END IF
    CALL HLINE(NTEX,GLINE,I)
    IF (IP==NLINES) EXIT
  END DO
  WRITE(NTEX,*) '\end{tabular}'
!-------------------------------------------------------------------------------
!
  WRITE(NTEX,*) '\clearpage'
!
ENDDO
IF (LHOOK) CALL DR_HOOK('write_cover_tex_water',1,ZHOOK_HANDLE)
!-------------------------------------------------------------------------------
!
END SUBROUTINE WRITE_COVER_TEX_WATER