prep_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 prep_ocean_mercatorvergrid(HPROGRAM,OUNIF)
!   ######################################################################
!
!!****  *PREP_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_ocean_grid
USE modi_read_z1d_netcdf
USE modi_read_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
 CHARACTER(LEN=6),  INTENT(IN)  :: hprogram ! program calling
LOGICAL, INTENT(IN) :: ounif
!
!*      0.2    declarations of local variables
!
INTEGER            :: ik1,ikk1
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!       1.     Allocations
!              -----------
IF (lhook) CALL dr_hook('PREP_OCEAN_MERCATORVERGRID',0,zhook_handle)

!       1.     Principal grid
!              --------------
!CLB
nockmin=0
ik1=nockmin+1

IF (.NOT.ounif) THEN
  !
  CALL read_z1d_netcdf ! doit initialiser(NOCKMAX,XZHOC(:))
  !
ELSE
  !
  ! Initialize vertical grid for idealized case
  nockmax=40
  ikk1=nockmax-1
  ALLOCATE(xzhoc(nockmin:nockmax))
  xzhoc(0)=0.
  xzhoc(ik1)    = -1.
  xzhoc(2)      = -5.
  xzhoc(3)      = -10.
  xzhoc(4)      = -15.
  xzhoc(5)      = -20.
  xzhoc(6)      = -25.
  xzhoc(7)      = -30.
  xzhoc(8)      = -40.
  xzhoc(9)      = -50.
  xzhoc(10)     = -60.
  xzhoc(11)     = -75.
  xzhoc(12)     = -100.
  xzhoc(13)     = -125.
  xzhoc(14)     = -150.
  xzhoc(15)     = -175.
  xzhoc(16)     = -200.
  xzhoc(17)     = -225.
  xzhoc(18)     = -250.
  xzhoc(19)     = -300.
  xzhoc(20)     = -400.
  xzhoc(21)     = -500.
  xzhoc(22)     = -600.
  xzhoc(23)     = -700.
  xzhoc(24)     = -800.
  xzhoc(25)     = -900.
  xzhoc(26)     = -1000.
  xzhoc(27)     = -1100.
  xzhoc(28)     = -1200.
  xzhoc(29)     = -1300.
  xzhoc(30)     = -1400.
  xzhoc(31)     = -1500.
  xzhoc(32)     = -1750.
  xzhoc(33)     = -2000.
  xzhoc(34)     = -2250.
  xzhoc(35)     = -2500.
  xzhoc(36)     = -2750.
  xzhoc(37)     = -3000.
  xzhoc(38)     = -3250.
  xzhoc(ikk1)   = -3500.
  xzhoc(nockmax)  = -4000.
!  WRITE(0,*) 'Number of vertical levels ',NOCKMAX
!  WRITE(0,*) 'Depth of vertical level (m)',XZHOC(:)
ENDIF
!
IF (lhook) CALL dr_hook('PREP_OCEAN_MERCATORVERGRID',1,zhook_handle)
!
END SUBROUTINE prep_ocean_mercatorvergrid