put_sfxcpln.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 put_sfxcpl_n (F, I, S, U, W, &
                         hprogram,ki,ksw,psw_bands,pzenith, &
                        pland_wtd,pland_fwtd,pland_fflood, &
                        pland_piflood,psea_sst,psea_ucu,   &
                        psea_vcu,pseaice_sit,pseaice_cvr,  &
                        pseaice_alb,ptsrad,                &
                        pdir_alb,psca_alb,pemis,ptsurf     )  
!     #################################################################################################
!
!!****  *PUT_SFXCPL_n* - routine to modify some variables in surfex from information coming
!                        from an ocean and/or a river routing model (but already on Surfex grid)
!
!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      B. Decharme      *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    08/2009
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
!
!
!
USE modd_flake_n, ONLY : flake_t
USE modd_isba_n, ONLY : isba_t
USE modd_seaflux_n, ONLY : seaflux_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_watflux_n, ONLY : watflux_t
!
USE modd_surf_par,   ONLY : xundef
!
USE modn_sfx_oasis,  ONLY : lwater
USE modd_sfx_oasis,  ONLY : lcpl_sea, lcpl_seaice, &
                            lcpl_land, lcpl_gw,    &
                            lcpl_flood
!                          
USE modi_get_luout
!
USE modi_abor1_sfx
USE modi_put_sfx_land
USE modi_put_sfx_sea
USE modi_update_esm_surf_atm_n
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(flake_t), INTENT(INOUT) :: f
TYPE(isba_t), INTENT(INOUT) :: i
TYPE(seaflux_t), INTENT(INOUT) :: s
TYPE(surf_atm_t), INTENT(INOUT) :: u
TYPE(watflux_t), INTENT(INOUT) :: w
!
 CHARACTER(LEN=6),        INTENT(IN)  :: hprogram
INTEGER,                 INTENT(IN)  :: ki      ! number of points
INTEGER,                 INTENT(IN)  :: ksw     ! number of bands
!
REAL, DIMENSION(KI),      INTENT(IN) :: pzenith
REAL, DIMENSION(KSW),     INTENT(IN) :: psw_bands ! mean wavelength of each shortwave band (m)
!
REAL, DIMENSION(KI),      INTENT(IN) :: pland_wtd     ! Land water table depth (m)
REAL, DIMENSION(KI),      INTENT(IN) :: pland_fwtd    ! Land grid-cell fraction of water table rise (-)
REAL, DIMENSION(KI),      INTENT(IN) :: pland_fflood  ! Land Floodplains fraction (-)
REAL, DIMENSION(KI),      INTENT(IN) :: pland_piflood ! Land Potential flood infiltration (kg/m2)
!
REAL, DIMENSION(KI),      INTENT(IN) :: psea_sst ! Sea surface temperature (K)
REAL, DIMENSION(KI),      INTENT(IN) :: psea_ucu ! Sea u-current stress (Pa)
REAL, DIMENSION(KI),      INTENT(IN) :: psea_vcu ! Sea v-current stress (Pa)
!
REAL, DIMENSION(KI),      INTENT(IN) :: pseaice_sit ! Sea-ice Temperature (K)
REAL, DIMENSION(KI),      INTENT(IN) :: pseaice_cvr ! Sea-ice cover (-)
REAL, DIMENSION(KI),      INTENT(IN) :: pseaice_alb ! Sea-ice albedo (-)
!
REAL, DIMENSION(KI),     INTENT(OUT) :: ptsrad   ! Total radiative temperature see by the atmosphere
REAL, DIMENSION(KI),     INTENT(OUT) :: ptsurf   ! Total surface temperature see by the atmosphere
REAL, DIMENSION(KI),     INTENT(OUT) :: pemis    ! Total emissivity see by the atmosphere
REAL, DIMENSION(KI,KSW), INTENT(OUT) :: pdir_alb ! Total direct albedo see by the atmosphere
REAL, DIMENSION(KI,KSW), INTENT(OUT) :: psca_alb ! Total diffus albedo see by the atmosphere
!
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
!
INTEGER :: ilu, iluout
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('PUT_SFXCL_N',0,zhook_handle)
!
 CALL get_luout(hprogram,iluout)
!-------------------------------------------------------------------------------
!
! Global argument
!
IF(ki/=u%NSIZE_FULL)THEN
  WRITE(iluout,*) 'size of field from the coupler :', ki
  WRITE(iluout,*) 'size of field in SURFEX        :', u%NSIZE_FULL
  CALL abor1_sfx('PUT_SFXCPL_N: VECTOR SIZE NOT CORRECT FOR COUPLING')
ENDIF
!
!-------------------------------------------------------------------------------
! Put variable over land tile
!-------------------------------------------------------------------------------
!
IF(lcpl_land)THEN
  CALL put_sfx_land(i, u, &
                    iluout,lcpl_gw,lcpl_flood,pland_wtd(:),       &
                    pland_fwtd(:),pland_fflood(:),pland_piflood(:))        
ENDIF
!
!-------------------------------------------------------------------------------
! Put variable over sea and/or water tile
!-------------------------------------------------------------------------------
!
IF(lcpl_sea)THEN
!
  CALL put_sfx_sea(s, u, w, &
                   iluout,lcpl_seaice,lwater,psea_sst(:),psea_ucu(:),       &
                   psea_vcu(:),pseaice_sit(:),pseaice_cvr(:),pseaice_alb(:) )
!
ENDIF
!
!-------------------------------------------------------------------------------
! Update radiative properties at time t+1 for radiative scheme
!-------------------------------------------------------------------------------
!
IF(lcpl_sea.OR.lcpl_flood)THEN
  CALL update_esm_surf_atm_n(f, i, s, u, w, &
                             hprogram, ki, ksw, pzenith, psw_bands,    &
                             ptsrad, pdir_alb, psca_alb, pemis, ptsurf )
ENDIF
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PUT_SFXCL_N',1,zhook_handle)
!
!
END SUBROUTINE put_sfxcpl_n