init_io_surf_ncn.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 init_io_surf_nc_n (DTCO, U, DGU, &
                                    hmask,haction)
!     ######################
!
!!****  *INIT_IO_SURF_NC* Keep in memory the netcdf ID of the output files
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      F. Habets   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      modified 05/04 by P. LeMoigne *Meteo France*
!!      modified 06/10 by S. Faroux *Meteo France*
!!=================================================================
!
!*       0.   DECLARATIONS
!             ------------
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_n, ONLY : surf_atm_t
!
!
USE modd_surfex_mpi, ONLY : nindex, npio, nrank
!
!
USE modd_io_surf_nc, ONLY : nmask, cfilein_nc, cfileout_nc, lmask, nid_nc, &
                            cmask, lcreated, cfileout_nc_save, ldef
!
USE mode_gridtype_gauss
!
USE modi_read_surf
USE modi_get_size_full_n
USE modi_get_type_dim_n
USE modi_init_io_surf_mask_n
!
USE modi_get_dim_full_n
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
include "netcdf.inc"
!
!
TYPE(data_cover_t), INTENT(INOUT) :: dtco

TYPE(surf_atm_t), INTENT(INOUT) :: u
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
!
 CHARACTER(LEN=6),  INTENT(IN)  :: hmask    
 CHARACTER(LEN=5),  INTENT(IN)  :: haction 
!
INTEGER :: ifull, il, ilu, iluout, iret
LOGICAL          :: gexist, gopened
!
REAL(KIND=JPRB)  :: zhook_handle
!------------------------------------------------------------------------------ 
IF (lhook) CALL dr_hook('INIT_IO_SURF_NC_N',0,zhook_handle)
!
lmask = .true.
!
!$OMP BARRIER
!
IF (haction=='READ ') THEN
  INQUIRE(file=cfilein_nc,exist=gexist)
  IF (gexist) THEN 
    IF (nrank==npio) THEN
!$OMP SINGLE
      iret = nf_open(cfilein_nc,nf_nowrite,nid_nc)
!$OMP END SINGLE
    ENDIF
    CALL read_surf(&
                   'NC    ','DIM_FULL',ifull,iret,hdir='A')
  ENDIF
ELSE 
  CALL get_dim_full_n(u, &
                      ifull)
  IF (nrank==npio) THEN
!$OMP SINGLE
    INQUIRE(file=cfileout_nc,exist=gexist)
    INQUIRE(file=cfileout_nc,opened=gopened)
    IF (.NOT.gopened) THEN
      iret = nf_open(cfileout_nc,nf_write,nid_nc)
    ENDIF
    IF (ldef) iret = nf_redef(nid_nc)
!$OMP END SINGLE
  ENDIF
ENDIF
!
! nindex is needed for call to get_size_full_n. In init_index_mpi, 
! it's not initialized for first readings.  
IF (.NOT.ALLOCATED(nindex)) THEN
  ALLOCATE(nindex(ifull))
  nindex(:) = 0
ENDIF
!
! size by MPI task. NINDEX is supposed to be initialized at this step.  
 CALL get_size_full_n(u, &
                      'OFFLIN',ifull,ilu)
!
il = ilu
 CALL get_type_dim_n(dtco, u, &
                     hmask,il)
!
 CALL init_io_surf_mask_n(dtco, u, &
                          hmask, il, iluout, ilu, nmask)
!
 cmask = hmask
!
!$OMP BARRIER
!------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('INIT_IO_SURF_NC_N',1,zhook_handle)
!------------------------------------------------------------------------------
!
END SUBROUTINE init_io_surf_nc_n