sst_update.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 sst_update (DTS, S, PSST)
!   ###############################################################
!!****  *SST_UPDATE*
!!
!!    PURPOSE
!!    -------
!
!     performs the time evolution of sst
!              
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!      
!!    none
!!
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!
!!      P. Le Moigne          * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    09/2007
!!
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
USE modd_data_seaflux_n, ONLY : data_seaflux_t
USE modd_seaflux_n, ONLY : seaflux_t
!
USE modi_temporal_dists
USE modi_temporal_lts
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
!
TYPE(data_seaflux_t), INTENT(INOUT) :: DTS
TYPE(seaflux_t), INTENT(INOUT) :: S
!
REAL,   DIMENSION(:), INTENT(INOUT) :: PSST    ! sst
!
!*      0.2    declarations of local variables
!
INTEGER                                  :: JXP
REAL, DIMENSION(SIZE(PSST))              :: ZSST, ZSST0
REAL                                     :: ZSDTJX
REAL                                     :: ZDT, ZALPHA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-----------------------------------------------------------------
!
IF (lhook) CALL dr_hook('SST_UPDATE',0,zhook_handle)
!
IF (.NOT.s%LTZTIME_DONE) THEN
   s%LTZTIME_DONE = .true.
   s%JSX = 1
   s%TZTIME%TDATE%YEAR  = s%TTIME%TDATE%YEAR
   s%TZTIME%TDATE%MONTH = s%TTIME%TDATE%MONTH
   s%TZTIME%TDATE%DAY   = s%TTIME%TDATE%DAY
   s%TZTIME%TIME        = s%TTIME%TIME
ENDIF
!
zsst0(:) = s%XSST_INI(:)
!
IF ( temporal_lts( s%TTIME, dts%TDATA_SST(1) ) ) THEN
  !
  CALL temporal_dists ( dts%TDATA_SST(1)%TDATE%YEAR,dts%TDATA_SST(1)%TDATE%MONTH, &
                        dts%TDATA_SST(1)%TDATE%DAY ,dts%TDATA_SST(1)%TIME,        &
                        s%TZTIME%TDATE%YEAR   ,s%TZTIME%TDATE%MONTH,          &
                        s%TZTIME%TDATE%DAY    ,s%TZTIME%TIME,                 &
                        zsdtjx                                            )  
  !
  CALL temporal_dists ( s%TTIME%TDATE%YEAR   ,s%TTIME%TDATE%MONTH,           &
                        s%TTIME%TDATE%DAY    ,s%TTIME%TIME,                  &
                        s%TZTIME%TDATE%YEAR  ,s%TZTIME%TDATE%MONTH,          &
                        s%TZTIME%TDATE%DAY   ,s%TZTIME%TIME,                 &
                        zdt                                              )  
  !
  zalpha = zdt / zsdtjx
  !
  zsst(:)= zsst0(:)+(dts%XDATA_SST(:,1)-zsst0(:))*zalpha
  !
ELSE IF ( .NOT. temporal_lts( s%TTIME, dts%TDATA_SST(dts%NTIME) ) ) THEN
  !
  zsst(:) = dts%XDATA_SST(:,dts%NTIME)
  !
ELSE
  !
  DO
    jxp = s%JSX + 1
    IF ( temporal_lts( s%TTIME, dts%TDATA_SST(jxp)) ) EXIT
    s%JSX = s%JSX + 1
  ENDDO
  !  
  CALL temporal_dists ( dts%TDATA_SST(jxp)%TDATE%YEAR,dts%TDATA_SST(jxp)%TDATE%MONTH,   &
                        dts%TDATA_SST(jxp)%TDATE%DAY ,dts%TDATA_SST(jxp)%TIME,          &
                        dts%TDATA_SST(s%JSX)%TDATE%YEAR ,dts%TDATA_SST(s%JSX)%TDATE%MONTH,  &
                        dts%TDATA_SST(s%JSX)%TDATE%DAY  ,dts%TDATA_SST(s%JSX)%TIME,         &
                        zsdtjx                                            )  
  !
  CALL temporal_dists ( s%TTIME%TDATE%YEAR   ,s%TTIME%TDATE%MONTH,                  &
                        s%TTIME%TDATE%DAY    ,s%TTIME%TIME,                         &
                        dts%TDATA_SST(s%JSX)%TDATE%YEAR,dts%TDATA_SST(s%JSX)%TDATE%MONTH,   &
                        dts%TDATA_SST(s%JSX)%TDATE%DAY ,dts%TDATA_SST(s%JSX)%TIME,          &
                        zdt                                             )  
  !
  zalpha = zdt / zsdtjx

  !
  zsst(:)= dts%XDATA_SST(:,s%JSX)+(dts%XDATA_SST(:,jxp)-dts%XDATA_SST(:,s%JSX))*zalpha
  !                
END IF
!
psst(:) = zsst(:)
!
IF (lhook) CALL dr_hook('SST_UPDATE',1,zhook_handle)
!-----------------------------------------------------------------
!
END SUBROUTINE sst_update