add_forecast_to_date_surf.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 add_forecast_to_date_surf(KYEAR,KMONTH,KDAY,PSEC)
!     #######################################################
!
!!****  *ADD_FORECAST_TO_DATE* - finds the current date and hour of a forecast
!!
!!    PURPOSE
!!    -------
!!
!!                                WARNING
!!
!!      -----> Only correct for dates between 19900301 and 21000228   <-----
!!
!!  The correct test should be:
!! IF( ((MOD(KYEAR,4)==0).AND.(MOD(KYEAR,100)/=0)) .OR. (MOD(KYEAR,400)==0))THEN
!!
!!**  METHOD
!!    ------
!!
!!      A recursive method is used, adding one day ofter the other.
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    Book 2
!!
!!    AUTHOR
!!    ------
!!      
!     V.Masson  Meteo-France
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/09/95
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declaration of arguments
!              ------------------------
INTEGER, INTENT(INOUT) :: kyear  ! year of date
INTEGER, INTENT(INOUT) :: kmonth ! month of date
INTEGER, INTENT(INOUT) :: kday   ! day of date
REAL,    INTENT(INOUT) :: psec   ! number of seconds since date at 00 UTC
!
!*       0.2   Declaration of local variables
!              ------------------------------
!
INTEGER :: idays ! number of days in KMONTH
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!*       1.    Return condition: less than one day to add
!              ------------------------------------------
!
IF (lhook) CALL dr_hook('ADD_FORECAST_TO_DATE_SURF',0,zhook_handle)
DO 
  IF (86400.-psec > 1.e-6) EXIT
!
!-------------------------------------------------------------------------------
!
!*       2.    Adding one day
!              --------------
!
  psec=psec-86400.
!
!*       2.1   Number of days in a month
!              -------------------------
!
  SELECT CASE (kmonth)
    CASE(4,6,9,11)
      idays=30
    CASE(1,3,5,7:8,10,12)
      idays=31
    CASE(2)
      IF( ((mod(kyear,4)==0).AND.(mod(kyear,100)/=0)) .OR. (mod(kyear,400)==0))THEN
        idays=29
      ELSE
        idays=28
      ENDIF
  END SELECT
!
!*       2.2   Last day of the month
!              ---------------------
!
  IF (kday==idays) THEN
    IF (kmonth==12) THEN
      kday=1
      kmonth=1
      kyear=kyear+1
    ELSE
      kday=1   
      kmonth=kmonth+1
    ENDIF
! 
!*       2.3   Other days
!              ----------
  ELSE
    kday=kday+1
  ENDIF
!
!-------------------------------------------------------------------------------
!
!*       3.    Recursive call
!              --------------
!
ENDDO
IF (lhook) CALL dr_hook('ADD_FORECAST_TO_DATE_SURF',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE add_forecast_to_date_surf