open_aux_io_surf_nc.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 open_aux_io_surf_nc (&
                                      hfile,hfiletype,hmask)
!     #######################################################
!
!!****  *OPEN_AUX_IO_SURF_ASC* - chooses the routine to OPENialize IO
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson    *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    10/2006
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
!
!
USE modd_io_surf_nc,ONLY: nmask,nfull,cmask, nluout, &
                            cfilein_nc, nid_nc, nfull_aux 
USE modi_get_luout
USE modi_read_surf
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_get_1d_mask
!
IMPLICIT NONE
!
include "netcdf.inc"
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
!
 CHARACTER(LEN=28), INTENT(IN)  :: hfile     ! file name
 CHARACTER(LEN=6),  INTENT(IN)  :: hfiletype ! main program
 CHARACTER(LEN=6),  INTENT(IN)  :: hmask
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(:),ALLOCATABLE :: zfull  ! total cover
INTEGER                        :: ilu,iret, il
INTEGER                        :: inb ! number of articles in the file
 CHARACTER(LEN=28) :: yfile
LOGICAL :: gopened
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('OPEN_AUX_IO_SURF_NC',0,zhook_handle)

yfile = trim(hfile)//".nc"
INQUIRE(file=yfile,opened=gopened)
IF (.NOT.gopened) THEN
  iret = nf_open(yfile,nf_nowrite,nid_nc)
ENDIF

 CALL get_luout('NC    ',nluout)
!
 cmask = hmask
 cfilein_nc = yfile
 CALL read_surf(&
                'NC    ','DIM_FULL',ilu,iret)
nfull_aux = ilu
!
!------------------------------------------------------------------------------
nfull = nfull_aux
!
il = nfull
ALLOCATE(zfull(il))
ALLOCATE(nmask(il))
zfull=1.
 CALL get_1d_mask(il,il,zfull,nmask)
!
!------------------------------------------------------------------------------
 cmask = hmask
IF (lhook) CALL dr_hook('OPEN_AUX_IO_SURF_NC',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE open_aux_io_surf_nc