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, 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_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, cfile_nc, &
                            nfull
!
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
!
USE netcdf
!
IMPLICIT NONE
!
!
!
TYPE(data_cover_t), INTENT(INOUT) :: DTCO

TYPE(surf_atm_t), INTENT(INOUT) :: U
!
 CHARACTER(LEN=6),  INTENT(IN)  :: HMASK    
 CHARACTER(LEN=5),  INTENT(IN)  :: HACTION 
!
INTEGER :: 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.
!
IF (haction=='READ ') THEN
  INQUIRE(file=cfilein_nc,exist=gexist)
  cfile_nc = cfilein_nc
  IF (gexist) THEN 
    iret = nf90_open(cfilein_nc,nf90_nowrite,nid_nc)
    CALL read_surf('NC    ','DIM_FULL',nfull,iret,hdir='A')
  ENDIF
ELSE 
  CALL get_dim_full_n(u%NDIM_FULL, nfull)
  IF (nrank==npio) THEN
    INQUIRE(file=cfileout_nc,exist=gexist)
    INQUIRE(file=cfileout_nc,opened=gopened)
    IF (.NOT.gopened) THEN
      iret = nf90_open(cfileout_nc,nf90_write,nid_nc)
    ENDIF
    IF (ldef) iret = nf90_redef(nid_nc)
  ENDIF
  cfile_nc = cfileout_nc
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).AND.nrank==npio) THEN
  ALLOCATE(nindex(nfull))
  nindex(:) = 0
ELSE
  CALL get_dim_full_n(u%NDIM_FULL, nfull)
ENDIF
!
! size by MPI task. NINDEX is supposed to be initialized at this step.  
 CALL get_size_full_n('OFFLIN',nfull,u%NSIZE_FULL,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
!
!------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('INIT_IO_SURF_NC_N',1,zhook_handle)
!------------------------------------------------------------------------------
!
END SUBROUTINE init_io_surf_nc_n