interpol_ts_water_mth.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 interpol_ts_water_mth (W, &
                                        kyear,kmonth,kday,pts)
!     #######################################################
!
!!****  *INTERPOL_TS_WATER_MTH* - Interpolation of monthly TS water
!!
!!    PURPOSE
!!    -------
!!
!!
!!**  METHOD
!!    ------
!!
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      
!     B.Decharme  Meteo-France
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    28/01/10
!!      Modified    07/2015   B. Decharme : new linear interpolation
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
USE modd_watflux_n, ONLY : watflux_t
!
USE modi_interpol_quadra
USE modi_interpol_linear
!
USE modi_abor1_sfx
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declaration of arguments
!------------------------
! 
!
TYPE(watflux_t), INTENT(INOUT) :: w
!
INTEGER, INTENT(IN ) :: kyear  ! year of date
INTEGER, INTENT(IN ) :: kmonth ! month of date
INTEGER, INTENT(IN ) :: kday   ! day of date
!
REAL, DIMENSION(:), INTENT(OUT) :: pts   ! Water surface temperature at time t 
!
!*       0.2   Declaration of local variables
!              ------------------------------
!
REAL            :: zdat   ! current day in the current month
REAL            :: zndat  ! number of days in the current month
INTEGER         :: imth0  ! previous month
INTEGER         :: imth1  ! current month 
INTEGER         :: imth2  ! next month
INTEGER         :: indays ! number of days in KMONTH
!
INTEGER         :: idelta
!
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!*       1.    Number of days in a month
!              -------------------------
!
IF (lhook) CALL dr_hook('INTERPOL_TS_WATER_MTH',0,zhook_handle)
SELECT CASE (kmonth)
    CASE(4,6,9,11)
      indays=30
    CASE(1,3,5,7:8,10,12)
      indays=31
    CASE(2)
      IF( ((mod(kyear,4)==0).AND.(mod(kyear,100)/=0)) .OR. (mod(kyear,400)==0))THEN
        indays=29
      ELSE
        indays=28
      ENDIF
END SELECT
!
!
!-------------------------------------------------------------------------------
!
!*       2.    TS water Interpolation using previous, current and next month
!              -------------------------------------------------------------
!
zdat = REAL(kday)
zndat= REAL(indays)
!
! The current month correspond to the indice 2 (or 3 if next month))
!
IF (kmonth==w%TZTIME%TDATE%MONTH) THEN 
   idelta=0
ELSE
   idelta=1
END IF
!
imth0=1+idelta
imth1=2+idelta
imth2=3+idelta
!
IF(w%CINTERPOL_TS=='QUADRA')THEN
  CALL interpol_quadra(zdat,zndat,w%XTS_MTH(:,imth0),w%XTS_MTH(:,imth1),w%XTS_MTH(:,imth2),pts)
ELSEIF(w%CINTERPOL_TS=='LINEAR')THEN
  CALL interpol_linear(zdat,zndat,w%XTS_MTH(:,imth0),w%XTS_MTH(:,imth1),w%XTS_MTH(:,imth2),pts)
ELSEIF(w%CINTERPOL_TS=='UNIF')THEN
  pts(:) = w%XTS_MTH(:,imth1)
ELSE
  CALL abor1_sfx('INTERPOL_TS_WATER_MTH: interpolation method not supported')
ENDIF
!
IF (lhook) CALL dr_hook('INTERPOL_TS_WATER_MTH',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE interpol_ts_water_mth