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