modd_oceann.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.
!     #################
      MODULE modd_ocean_n
!     #################
!
!!****  *MODD_OCEAN_n - declaration of ocean varaiables 
!!                          for 1D oceanic model
!!
!!    PURPOSE
!!    -------
!     Declaration of ocean varaiables
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      C. Lebeaupin   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       04/2006
!!      Modified       07/2012, P. Le Moigne : CMO1D phasing
!
!*       0.   DECLARATIONS
!             ------------
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
TYPE ocean_t
!
!
!   Switche for interactive coupling with oceanic model
LOGICAL:: LMERCATOR   !set to .true. to initialize oceanic var. from Mercator
LOGICAL:: LCURRENT    !set to .true. to make initialize ocean state with current      
LOGICAL:: LPROGSST    !set to .true. to make SST evolve with tendance
INTEGER:: NTIME_COUPLING! coupling time frequency 
INTEGER:: NOCTCOUNT   !oceanic model counter
REAL :: XOCEAN_TSTEP       ! time step of the oceanic 1D model
!
! General surface: 
!
REAL, POINTER, DIMENSION(:,:) :: XSEAT  ! oceanic temperature profiles 
REAL, POINTER, DIMENSION(:,:) :: XSEAS  ! oceanic salinity profiles 
REAL, POINTER, DIMENSION(:,:) :: XSEAU  ! oceanic zonal current profiles 
REAL, POINTER, DIMENSION(:,:) :: XSEAV  ! oceanic meridian current profiles 
REAL, POINTER, DIMENSION(:,:) :: XSEAE  ! oceanic kinetic turbulent energy profiles (^(1/2))
REAL, POINTER, DIMENSION(:,:) :: XSEABATH !bathymetry indice
                                          !=1 for free sea water
                                          !=0 for sea-bed
REAL, POINTER, DIMENSION(:) ::   XSEAHMO! oceanic mixing lengths
!
REAL, POINTER, DIMENSION(:,:) :: XLE,XLK! oceanic mixing lengths
REAL, POINTER, DIMENSION(:,:) :: XKMEL,XKMELM  ! oceanic mixing coefficients
!
REAL, POINTER, DIMENSION(:) ::   XSEATEND! SST tendance
!
REAL, POINTER, DIMENSION(:,:) ::   XDTFSOL ! Temp tendancy due to solar flux
REAL, POINTER, DIMENSION(:)   ::   XDTFNSOL! -------------------- non solar flux
!
END TYPE ocean_t
!


 CONTAINS

!




SUBROUTINE ocean_init(YOCEAN)
TYPE(ocean_t), INTENT(INOUT) :: yocean
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_OCEAN_N:OCEAN_INIT",0,zhook_handle)
  nullify(yocean%XSEAT)
  nullify(yocean%XSEAS)
  nullify(yocean%XSEAU)
  nullify(yocean%XSEAV)
  nullify(yocean%XSEAE)
  nullify(yocean%XSEABATH)
  nullify(yocean%XSEAHMO)
  nullify(yocean%XLE)
  nullify(yocean%XLK)
  nullify(yocean%XKMEL)
  nullify(yocean%XKMELM)
  nullify(yocean%XSEATEND)
  nullify(yocean%XDTFNSOL)
  nullify(yocean%XDTFSOL)
yocean%LMERCATOR=.false.
yocean%LCURRENT=.false.
yocean%LPROGSST=.false.
yocean%NTIME_COUPLING=0
yocean%NOCTCOUNT=0
yocean%XOCEAN_TSTEP=5*60.
IF (lhook) CALL dr_hook("MODD_OCEAN_N:OCEAN_INIT",1,zhook_handle)
END SUBROUTINE ocean_init


END MODULE modd_ocean_n