read_z1d_netcdf.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 read_z1d_netcdf
!     ##############################################################
!
!!**** *READ_Z1D_NETCDF* reads the vertical grid in a netcdf file
!!
!!    PURPOSE
!!    -------
!!
!!    AUTHOR
!!    ------
!!
!!    C. Lebeaupin Brossier Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    11/2014
!!      initialisation of NOCKMAX,XZHOC
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_ocean_grid
USE modd_surf_par, ONLY : nundef
USE modd_prep_seaflux, ONLY : cfile_seaflx,ctype_seaflx
USE mode_read_netcdf_mercator
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
 CHARACTER (LEN=28) :: yfilename
 CHARACTER (LEN=28)  :: yncvarname
INTEGER :: jdimension
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL(KIND=JPRB) :: zhook_handle
                                           ! and ZLON arrays
!----------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('READ_Z1D_NETCDF',0,zhook_handle)
!
nockmax=-nundef
!
!*    1.      Read the netcdf arrays dimensions
!             ---------------------------------
IF (ctype_seaflx=="NETCDF") THEN
  yfilename=trim(cfile_seaflx)
  yncvarname="depth"
  CALL read_dim_cdf(yfilename,yncvarname,jdimension)
  nockmax=jdimension
  ALLOCATE(xzhoc(0:jdimension))
  xzhoc(0)=0.
!
!*    2.      Read the array in the netcdf file
!             ---------------------------------
  CALL read_z1d_cdf(yfilename,yncvarname,xzhoc(1:jdimension))
  IF (xzhoc(2)>0) xzhoc(:)=-xzhoc(:)
!  WRITE(0,*) 'Oceanic vertical grid readed in netcdf file'
!  WRITE(0,*) 'Number of level',NOCKMAX+1
!  WRITE(0,*) 'Depth of vertical level',XZHOC(:)
!
!----------------------------------------------------------------------------
ELSE
  WRITE(*,*) 'ERROR IN READ_Z1D_NETCF: ', yfilename, ' HAS NOT A NETCDF TYPE'
  WRITE(*,*) 'CHECK CTYPE_SEAFLX IN NAM_PREP_SEAFLUX'
ENDIF
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_Z1D_NETCDF',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_z1d_netcdf