sfx_oasis_check.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 sfx_oasis_check (IO, U, KLUOUT)
!###################################################
!
!!****  *SFX_OASIS_CHECK* - Definitions for exchange of coupling fields
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      B. Decharme   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    10/2013
!!    10/2016 B. Decharme : bug surface/groundwater coupling
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
USE modn_sfx_oasis, ONLY : ccalving, lwater
USE modd_sfx_oasis, ONLY : lcpl_lake, lcpl_calving, lcpl_gw
!
USE modi_abor1_sfx
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(surf_atm_t), INTENT(INOUT) :: U
!
INTEGER, INTENT(IN) :: KLUOUT
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('SFX_OASIS_CHECK',0,zhook_handle)
!
IF(lcpl_lake)THEN
  IF(u%CWATER/='FLAKE ')THEN
    WRITE(kluout,*)'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
    WRITE(kluout,*)'!!!    SFX - LAKE coupling is asked   !!!'
    WRITE(kluout,*)'!!!     but CWATER /= FLAKE           !!!'
    WRITE(kluout,*)'!!!                                   !!!'
    WRITE(kluout,*)'!!!    Please check your namelist     !!!'
    WRITE(kluout,*)'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
    CALL abor1_sfx('SFX_OASIS_CHECK: SFX - LAKE coupling asked but CWATER /= FLAKE')
  ENDIF
ENDIF
!
!
IF(lcpl_calving)THEN
  IF(.NOT.io%LGLACIER)THEN
    WRITE(kluout,*)'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
    WRITE(kluout,*)'Calving flux is asked by SFX - OASIS coupling      '
    WRITE(kluout,*)'CCALVING = '//trim(ccalving)//' in NAM_SFX_LAND_CPL'
    WRITE(kluout,*)'but LGLACIER in not activated in NAM_ISBAn         '
    WRITE(kluout,*)'Please check your SURFEX namelist                  '
    WRITE(kluout,*)'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
    CALL abor1_sfx('SFX_OASIS_CHECK: Calving flux is asked by SFX - OASIS coupling')            
  ENDIF  
ENDIF  
!
!
IF(lcpl_gw.AND.io%CISBA/='DIF')THEN
   WRITE(kluout,*)'SFX_OASIS_CHECK: Water table depth / surface coupling requires ISBA-DF'
   CALL abor1_sfx('SFX_OASIS_CHECK: ISBA-DF is required for SFX - Groundwater coupling')
ENDIF   
IF(.NOT.lcpl_gw.AND.io%CISBA=='DIF'.AND.io%LWTD)THEN           
      WRITE(kluout,*)'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
      WRITE(kluout,*)'!!! A groundwater map is specified and LAND coupling    !!!'
      WRITE(kluout,*)'!!!  is activated but not groundwater/surface coupling  !!!'
      WRITE(kluout,*)'!!!                                                     !!!'
      WRITE(kluout,*)'!!! ARE YOU SURE : YOU DO NOT WANT THIS COUPLING ?      !!!'
      WRITE(kluout,*)'!!!                                                     !!!'
      WRITE(kluout,*)'!!! Use NAM_SFX_LAND_CPL with CWTD and CFWTD            !!!'
      WRITE(kluout,*)'!!!  if you want this coupling via OASIS                !!!'
      WRITE(kluout,*)'!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'  
ENDIF
!
!
IF(lwater.AND.(u%CWATER=='NONE '.OR.u%CWATER=='FLAKE'))THEN
  WRITE(kluout,*)'LWATER = ',lwater,'   CWATER = ',u%CWATER
  WRITE(kluout,*)'! Inland water should not be added to sea mask in case CWATER is NONE or FLAKE !!!'     
  WRITE(kluout,*)'! Change CWATER or put LWATER=.FALSE. in NAM_SFX_SEA_CPL !!!'     
  CALL abor1_sfx('SFX_OASIS_READ_NAM: LWATER and CWATER not consistent')
ENDIF
!
IF (lhook) CALL dr_hook('SFX_OASIS_CHECK',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE sfx_oasis_check