ocean_mercatorvergrid.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 ocean_mercatorvergrid
!   ######################################################################
!
!!****  *OCEAN_MERCATORVERGRID*  
!!
!!    PURPOSE
!!    -------
!
!     Define the vertical ocean grid
!         
!     
!!**  METHOD
!!    ------
!
!
!
!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    MODD_OCEAN_CST
!!
!!      
!!    REFERENCE
!!    ---------
!!
!!      
!!    AUTHOR
!!    ------
!!
!!      C. Lebeaupin Brossier  * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2008
!       D.BARBARY 11/2014 : HPROGRAM,OUNIF in Calling OCEAN_MERCATORVERGRID
!                           Reading oceanic level and depth
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_surf_par,   ONLY : xundef
USE modd_ocean_grid
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.2    declarations of local variables
!
!
REAL                   :: zup,zdown
INTEGER            :: jloop
!
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!       1.     Allocations
!              -----------
IF (lhook) CALL dr_hook('OCEAN_MERCATORVERGRID',0,zhook_handle)

!       1.     Principal grid
!              --------------
!
!-------------------------------------------------------------------------------
!       2.     Ocean Secondary Grids 
!              ---------------------
!
!$OMP SINGLE
!
ALLOCATE(xk1(nockmin:nockmax))
ALLOCATE(xk2(nockmin:nockmax))
ALLOCATE(xk3(nockmin:nockmax))
ALLOCATE(xk4(nockmin:nockmax))
!
ALLOCATE(xz2(nockmin:nockmax))
ALLOCATE(xdz1(nockmin:nockmax))
ALLOCATE(xdz2(nockmin:nockmax))
!
ALLOCATE(xray(nockmin:nockmax))
!
DO jloop = nockmin,nockmax-1
  xz2(jloop) = (xzhoc(jloop+1) + xzhoc(jloop))/2.
  xdz1(jloop) = xzhoc(jloop) - xzhoc(jloop+1)
ENDDO
xz2(nockmax) = xzhoc(nockmax) + ((xzhoc(nockmax)-xzhoc(nockmax-1))/2.) !dernier niveau+delta/2
xdz1(nockmax) = xdz1(nockmax-1)
!
!
xdz2(nockmin) = - xz2(nockmin)
DO jloop = nockmin+1,nockmax
  xdz2(jloop) = xz2(jloop-1) - xz2(jloop)
ENDDO
!!
!!       3.     Grid Parameters
!!              ---------------
!
xk1(nockmin) = 0.
DO jloop = nockmin+1,nockmax
  xk1(jloop) = -1. / (xdz2(jloop)*xdz1(jloop-1))
ENDDO
!
zup=1.
DO jloop = nockmin,nockmax
  xk4(jloop) =  1. / ( xdz1(jloop) * xdz1(jloop) )
  zdown = rayo(xz2(jloop))
  xray(jloop) = zup - zdown
  zup = zdown
ENDDO
!
DO jloop = nockmin,nockmax-1
  xk2(jloop) = -1. / (xdz2(jloop)*xdz1(jloop))
  xk3(jloop) = -1. / (xdz1(jloop)*xdz2(jloop+1))
ENDDO
!
xk2(nockmax) = xk2(nockmax-1)
xk3(nockmax) = 0.
!
!$OMP END SINGLE
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('OCEAN_MERCATORVERGRID',1,zhook_handle)
 CONTAINS
!rayo
!-------------------------------------------------------------------------------
!
!-------------------------------------------------------------------------------
!!              #########################################
             FUNCTION rayo(Z) RESULT(RR)
!              #########################################
!
!
!!****  *RAYOFCTX* 
!!
!!    PURPOSE
!!    -------
!compute solar penetration coefficient
!    
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!     
!!    REFERENCE
!!    ---------
!!    Paulson and Simpson 1977
!!     
!!    AUTHOR
!!    ------
!!     C. Lebeaupin  *Meteo-France* (adapted from S. Belamari's code)
!!
!!    MODIFICATIONS
!!    -------------
!!      Original     27/02/2006
!!
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!
USE modd_ocean_csts,ONLY : xr,xd1,xd2
!
!*      0.1    declarations of arguments
!
REAL :: rr,z      
REAL(KIND=JPRB) :: zhook_handle

!
!-------------------------------------------------------------------------------
!
!*       1.     COMPUTE RAYO
!               ------------
!
IF (lhook) CALL dr_hook('RAYO',0,zhook_handle)
rr = xr*exp(z/xd1) + (1-xr)*exp(z/xd2)
IF (lhook) CALL dr_hook('RAYO',1,zhook_handle)
!
END FUNCTION rayo
!
!-------------------------------------------------------------------------------
END SUBROUTINE ocean_mercatorvergrid