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,HDIR)
!     #######################################################
!
!!****  *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
!
USE netcdf
!
IMPLICIT NONE
!
!
!*       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
 CHARACTER(LEN=1), INTENT(IN) :: HDIR  
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
REAL, DIMENSION(:),ALLOCATABLE :: ZFULL  ! total cover
INTEGER                        :: ILU,IRET, IL
 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 = nf90_open(yfile,nf90_nowrite,nid_nc)
ENDIF

 CALL get_luout('NC    ',nluout)
!
CMASK = HMASK
CFILEIN_NC = YFILE
 CALL read_surf(&
                'NC    ','DIM_FULL',ilu,iret,hdir=hdir)
nfull_aux = ilu
!
!------------------------------------------------------------------------------
nfull = nfull_aux
!
il = nfull
ALLOCATE(zfull(il))
ALLOCATE(nmask(il))
zfull=1.
 CALL get_1d_mask(il,il,zfull,nmask)
!
DEALLOCATE(zfull)
!
!------------------------------------------------------------------------------
CMASK = HMASK
IF (lhook) CALL dr_hook('OPEN_AUX_IO_SURF_NC',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE open_aux_io_surf_nc