read_file_masktopd.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_file_masktopd(KI)
!     ##########################
!
!!
!!    PURPOSE
!!    -------
!        
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!      
!!    REFERENCE
!!    ---------
!!     
!!    AUTHOR
!!    ------
!!
!!      B. Vincendon    * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      Original   11/2011
!!                 03/2014 (B. Vincendon) modification of mask_surf files format
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
USE modd_topd_par, ONLY : nunit
USE modd_topodyn
USE modd_coupling_topd, ONLY : nmaski, nnpix, nmaskt
USE modd_surf_par,        ONLY : xundef, nundef
USE modi_open_file
USE modi_close_file
USE modi_read_topd_file
!
USE mode_gridtype_conf_proj
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
INTEGER, INTENT(IN)             :: ki    ! Grid dimensions
!
!
!*      0.2    declarations of local variables
INTEGER                    :: jcat,jmesh,jpix
INTEGER,DIMENSION(KI)      :: jp_in_m
INTEGER                    :: inumpix,ieof,itmp
INTEGER :: iimax,ijmax
 CHARACTER(LEN=50) :: yname
REAL, DIMENSION(:),ALLOCATABLE    :: ztopd_read !Topgraphic variable read
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('READ_FILE_MASKTOPD',0,zhook_handle)
!
!*       0.     Initialization:
!               ---------------
!
ALLOCATE(ztopd_read(npmax))
!
DO jcat=1,nncat
  !
  yname=trim(ccat(jcat))//trim('.mask_topd')
  !
  CALL read_topd_file('OFFLIN',yname,'FORMATTED',nnpt(jcat),ztopd_read)
  !
  DO jpix=1,npmax
    !
    IF ( nline(jcat,jpix)/=0 ) THEN
      IF  (ztopd_read(jpix)==xundef ) THEN
        nmaskt(jcat,nline(jcat,jpix)) = nundef
      ELSE
        nmaskt(jcat,nline(jcat,jpix)) = floor(ztopd_read(jpix))
      ENDIF
    ENDIF
    !
  ENDDO
  !
ENDDO
!
ALLOCATE(nnpix(ki))
nnpix(:) = nundef
DO jmesh=1,ki
  nnpix(jmesh) = count(nmaskt(:,:)==jmesh)
ENDDO
inumpix=maxval(nnpix)
!
ALLOCATE(nmaski(ki,nncat,inumpix))
nmaski(:,:,:) = nundef
DO jcat=1,nncat
jp_in_m(:)=1
  !
  yname=trim(ccat(jcat))//trim('.mask_surf')
  CALL open_file('ASCII ',nunit,yname,'FORMATTED','READ')
  ! 
  ieof=0
  DO WHILE(ieof==0)
      READ(nunit,*,iostat=ieof) jmesh,itmp 
      nmaski(jmesh,jcat,jp_in_m(jmesh))=itmp
      jp_in_m(jmesh)=jp_in_m(jmesh)+1
  ENDDO
  !
  CALL close_file('ASCII ',nunit)
  !
ENDDO
!
IF (lhook) CALL dr_hook('READ_FILE_MASKTOPD',1,zhook_handle)
!
END SUBROUTINE read_file_masktopd