prep_hor_watflux_field.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 prep_hor_watflux_field (DTCO, U, &
                                    wg, w, &
                                   hprogram,hsurf,hatmfile,hatmfiletype,hpgdfile,hpgdfiletype)
!     #################################################################################
!
!!****  *PREP_HOR_WATFLUX_FIELD* - Reads, interpolates and prepares a water field
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     S. Malardel
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2004
!!      P. Le Moigne 10/2005, Phasage Arome
!!------------------------------------------------------------------
!
!
!
!
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modd_watflux_grid_n, ONLY : watflux_grid_t
USE modd_watflux_n, ONLY : watflux_t
!
USE modd_prep,          ONLY : cingrid_type, cinterp_type, xzs_ls, xlat_out, xlon_out, &
                               xx_out, xy_out, cmask
!
USE modi_read_prep_watflux_conf
USE modi_prep_watflux_grib
USE modi_prep_watflux_unif
USE modi_prep_watflux_buffer
USE modi_hor_interpol
USE modi_get_luout
USE modi_prep_watflux_extern
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
!
TYPE(data_cover_t), INTENT(INOUT) :: dtco
TYPE(surf_atm_t), INTENT(INOUT) :: u
!
TYPE(watflux_grid_t), INTENT(INOUT) :: wg
TYPE(watflux_t), INTENT(INOUT) :: w
!
 CHARACTER(LEN=6),   INTENT(IN)  :: hprogram  ! program calling surf. schemes
 CHARACTER(LEN=7),   INTENT(IN)  :: hsurf     ! type of field
 CHARACTER(LEN=28),  INTENT(IN)  :: hatmfile    ! name of the Atmospheric file
 CHARACTER(LEN=6),   INTENT(IN)  :: hatmfiletype! type of the Atmospheric file
 CHARACTER(LEN=28),  INTENT(IN)  :: hpgdfile    ! name of the Atmospheric file
 CHARACTER(LEN=6),   INTENT(IN)  :: hpgdfiletype! type of the Atmospheric file
!
!*      0.2    declarations of local variables
!
 CHARACTER(LEN=6)              :: yfiletype ! type of input file
 CHARACTER(LEN=28)             :: yfile     ! name of file
 CHARACTER(LEN=6)              :: yfilepgdtype ! type of input file
 CHARACTER(LEN=28)             :: yfilepgd     ! name of file
REAL, POINTER, DIMENSION(:,:) :: zfieldin  ! field to interpolate horizontally
REAL, ALLOCATABLE, DIMENSION(:,:) :: zfieldout ! field interpolated   horizontally
INTEGER                       :: iluout    ! output listing logical unit
!
LOGICAL                       :: gunif     ! flag for prescribed uniform field
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------------
!
!
!*      1.     Reading of input file name and type
!
IF (lhook) CALL dr_hook('PREP_HOR_WATFLUX_FIELD',0,zhook_handle)
 CALL get_luout(hprogram,iluout)
!
 CALL read_prep_watflux_conf(hprogram,hsurf,yfile,yfiletype,yfilepgd,yfilepgdtype,&
                            hatmfile,hatmfiletype,hpgdfile,hpgdfiletype,iluout,gunif)
!
 cmask = 'WATER'
!
!-------------------------------------------------------------------------------------
!
!*      2.     Reading of input  configuration (Grid and interpolation type)
!
IF (gunif) THEN
  CALL prep_watflux_unif(iluout,hsurf,zfieldin)
ELSE IF (yfiletype=='GRIB  ') THEN
  CALL prep_watflux_grib(hprogram,hsurf,yfile,iluout,zfieldin)
ELSE IF (yfiletype=='MESONH' .OR. yfiletype=='ASCII ' .OR. yfiletype=='LFI   '.OR. yfiletype=='FA    ') THEN
   CALL prep_watflux_extern(&
                            hprogram,hsurf,yfile,yfiletype,yfilepgd,yfilepgdtype,iluout,zfieldin)
ELSE IF (yfiletype=='BUFFER') THEN
   CALL prep_watflux_buffer(hprogram,hsurf,iluout,zfieldin)
ELSE
  CALL abor1_sfx('PREP_HOR_WATFLUX_FIELD: data file type not supported : '//yfiletype)
END IF
!
!
!*      4.     Horizontal interpolation
!
ALLOCATE(zfieldout(SIZE(wg%XLAT),SIZE(zfieldin,2)))
!
 CALL hor_interpol(dtco, u, &
                   iluout,zfieldin,zfieldout)
!
!*      5.     Return to historical variable
!
SELECT CASE (hsurf)
 CASE('ZS     ') 
  ALLOCATE(xzs_ls(SIZE(zfieldout,1)))
  xzs_ls(:) = zfieldout(:,1)
 CASE('TSWATER')
  ALLOCATE(w%XTS(SIZE(zfieldout,1)))
  w%XTS(:) = zfieldout(:,1)
END SELECT
!
!-------------------------------------------------------------------------------------
!
!*      6.     Deallocations
!
DEALLOCATE(zfieldin )
DEALLOCATE(zfieldout)
IF (lhook) CALL dr_hook('PREP_HOR_WATFLUX_FIELD',1,zhook_handle)
!
!-------------------------------------------------------------------------------------
!
END SUBROUTINE prep_hor_watflux_field